Changeset 4650

Timestamp:

Aug 28, 2023, 8:17:17 AM (8 months ago)

Author:

evignon

Message:

petite correction de l'initialisation de la Ts forcee dans les cas au format standard

Location:

LMDZ6/trunk/libf/phylmd/dyn1d

Files:

• 1DUTILS.h (modified) (1 diff)
• 1D_read_forc_cases.h (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -266 +266 @@
 !Config  Def  = 0.
 !Config  Help = 
-       tsurf = 0.
+       time_ini = 0.
        CALL getin('time_ini',time_ini)
 

LMDZ6/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h

@@ -73 +73 @@
 
 ! initial and boundary conditions :
+      ! initial pressure
       psurf = ps_prof_cas
-      if (tskin_prof_cas .NE. 0.) THEN
+
+      !initial surface temperature
+      if (tskin_prof_cas .NE. 0.) then
+      ! we take the first value of the prescribed ts
           tsurf=tskin_prof_cas
+      else if (ts_prof_cas .NE. 0) then
+      ! if an initial ts value is present, we take it
+          tsurf=ts_prof_cas
       endif
 
+      ! ts forcing during the run (if any)
       tg = ts_prof_cas
       if ((tg .eq. 0.) .and. (tskin_prof_cas .NE. 0.)) THEN