add a lower limit for pressure in radiation

CliMA · Oct 1, 2024 · b6c1d8b · b6c1d8b
1 parent 1c94454
commit b6c1d8b
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/callbacks/callbacks.jl b/src/callbacks/callbacks.jl
@@ -89,7 +89,10 @@ NVTX.@annotate function rrtmgp_model_callback!(integrator)
 
     ᶜp = Fields.array2field(rrtmgp_model.center_pressure, axes(Y.c))
     ᶜT = Fields.array2field(rrtmgp_model.center_temperature, axes(Y.c))
-    @. ᶜp = TD.air_pressure(thermo_params, ᶜts)
+    # Here we multiply p_min by a small number to avoid exactly same
+    # pressure in the top layer and in the extra layer added in rrtmgp
+    p_min = RRTMGPI.get_p_min(rrtmgp_model)
+    @. ᶜp = max(TD.air_pressure(thermo_params, ᶜts), p_min * FT(1.01))
     # TODO: move this to RRTMGP
     @. ᶜT =
         min(max(TD.air_temperature(thermo_params, ᶜts), FT(T_min)), FT(T_max))