Update baroclinic_instability_model.jl

Author: glwagner

src/baroclinic_instability_model.jl

@@ -26,8 +26,8 @@ function baroclinic_instability_model(arch, Nx, Ny, Nz; Δt,
     buoyancy = SeawaterBuoyancy(
         equation_of_state = SeawaterPolynomials.TEOS10EquationOfState(Oceananigans.defaults.FloatType)),
 
-    closure = nothing,    
-    # closure = Oceananigans.TurbulenceClosures.CATKEVerticalDiffusivity(),
+    # closure = nothing,    
+    closure = Oceananigans.TurbulenceClosures.CATKEVerticalDiffusivity(),
     # closure = VerticalScalarDiffusivity(VerticallyImplicitTimeDiscretization(), κ=1e-5, ν=1e-4),
 
     # Coriolis forces for a rotating Earth