Changeset 3728

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

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• libf/phystd/tabfi_mod.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/changelog.txt

@@ -2072 +2072 @@
 This should makes things (hopefully) clearer for users and also better
 enforces a cleaner and clearer separation between dynamics and physics.
+
+== 18/04/2025 == EM
+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).

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)