Changeset 5171 for LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Aug 31, 2024, 8:39:33 PM (2 months ago)

Author:

evignon

Message:

adding a condition to prescribe surface temperature in 1D cases forced in flux
this is needed for the new MOSAI 1D case to avoid spurious crashes

File:

• LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -2715 +2715 @@
             yt1_new=(1./RCPD)*(AcoefH(j)+BcoefH(j)*y_flux_t1(j)*dtime)
             ytsurf_new(j)=yt1_new-y_flux_t1(j)/(Kech_h(j)*RCPD)
+            ! for cases forced in flux and for which forcing in Ts is needed
+            ! to prevent the latter to reach unrealistic value (even if not used,
+            ! Ts is calculated and hgardfou can appear during the calculation
+            ! of surface saturation humidity for example
+            if (ok_forc_tsurf) ytsurf_new(j)=tg
           ENDDO