Update twolayerdns.jl

jbisits · jbisits · commit 0154aa9bd463 · 2023-10-30T20:03:45.000+11:00
diff --git a/src/twolayerdns.jl b/src/twolayerdns.jl
@@ -187,7 +187,7 @@ function TLDNS_simulation_setup(tldns::TwoLayerDNS, Δt::Number,
 
     # Custom saved output
     # Potential density
-    σ = seawater_density(model, geopotential_height = 0)
+    σ = seawater_density(model, geopotential_height = density_reference_gp_height)
 
     # Inferred vertical temperature diffusivity
     # T_mean = Average(T)
