Changeset 3728 for trunk/LMDZ.GENERIC/libf

Timestamp:

Apr 18, 2025, 2:53:06 PM (2 months ago)

Author:

emillour

Message:

Generic PCM:
Minor fix in tabfi: the physics time step stored in the startfi file should
not overwrite the real value of dtphys (computed in and passed by the dynamics).
EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/tabfi_mod.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/tabfi_mod.F90 ¶

@@ -87 +87 @@
       LOGICAL :: found
       CHARACTER(len=5) :: modname="tabfi"
+      real :: cntrl_daysec, cntrl_dtphys
       
       write(*,*)"tabfi: nid=",nid," tab0=",tab0," Lmodif=",Lmodif
@@ -181 +182 @@
       rcp = tab_cntrl(tab0+9)
       cpp=(8.314463/(mugaz/1000.0))/rcp
-      daysec = tab_cntrl(tab0+10)
-      dtphys = tab_cntrl(tab0+11)
+! daysec and dtphys are already set by inifis
+      cntrl_daysec = tab_cntrl(tab0+10)
+      cntrl_dtphys = tab_cntrl(tab0+11)
 ! Informations about planet for the physics only
       year_day = tab_cntrl(tab0+14)
@@ -201 +203 @@
 ! soil properties
       volcapa = tab_cntrl(tab0+35) ! volumetric heat capacity
+! warn user if there has been a change in daysec or dtphys
+      if (daysec.ne.cntrl_daysec) then
+        write(*,*) modname//" Warning: lenght of day daysec changed from ", &
+                   cntrl_daysec," to ",daysec
+      endif
+      if (dtphys.ne.cntrl_dtphys) then
+        write(*,*) modname//" Warning: time step dtphys changed from ", &
+                   cntrl_dtphys," to ",dtphys
+      endif
 !-----------------------------------------------------------------------
 !       Save some constants for later use (as routine arguments)