Changeset 3538 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Dec 9, 2024, 12:35:05 PM (2 weeks ago)

Author:

jbclement

Message:

PEM:
Fixing Norbert Schorghofer's subroutine which modified soil temperatures and was unwanted in the PEM loop.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• NS_fast_subs_mars.F90 (modified) (1 diff)
• changelog.txt (modified) (1 diff)
• deftank/lib_launchPEM.sh (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/NS_fast_subs_mars.F90

@@ -282 +282 @@
 
   if (Fgeotherm>0.) then
-     Tb = Tmean1
+     !Tb = Tmean1
      typeG = 1   ! will be overwritten by depths_avmeth
      rlow = 2*rhosatav(nz)-rhosatav(nz-1)
   else
-     Tb = T(nz)
+     !Tb = T(nz)
      typeG = -9
      rlow = rhosatav(nz-1)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -510 +510 @@
 == 06/12/2024 == JBC
 Small fixes to initialize and output ice table-related variables and in the launching script.
+
+== 09/12/2024 == JBC
+Fixing Norbert Schorghofer's subroutine which modified soil temperatures and was unwanted in the PEM loop.

trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh

@@ -310 +310 @@
         # PCM relaunch during the initialization cycle
         iPEM=1
-        cleanfiles diags/diagpem .nc $iPEM
-        cleanfiles diags/diagsoilpem .nc $iPEM
-        cleanfiles "out_PEM/run" "" $iPEM
-        cleanfiles starts/restart1D_postPEM .txt $iPEM
-        cleanfiles starts/restart_postPEM .nc $iPEM
-        cleanfiles starts/restartfi_postPEM .nc $iPEM
-        cleanfiles starts/restartpem .nc $iPEM
+        cleanfiles diags/diagpem .nc $(($iPEM - 1))
+        cleanfiles diags/diagsoilpem .nc $(($iPEM - 1))
+        cleanfiles "out_PEM/run" "" $(($iPEM - 1))
+        cleanfiles starts/restart1D_postPEM .txt $(($iPEM - 1))
+        cleanfiles starts/restart_postPEM .nc $(($iPEM - 1))
+        cleanfiles starts/restartfi_postPEM .nc $(($iPEM - 1))
+        cleanfiles starts/restartpem .nc $(($iPEM - 1))
         i_myear=$irelaunch
         sed -i "1s/.*/$i_myear $n_myear $convert_years $iPCM $iPEM $nPCM $nPCM_ini/" info_PEM.txt
@@ -334 +334 @@
         iPEM=$(echo "($iPCM - $nPCM_ini)/$nPCM + 1" | bc)
         il=$(echo "($irelaunch - $nPCM_ini + 1)%$nPCM + 1" | bc)
-        cleanfiles diags/diagpem .nc $iPEM
-        cleanfiles diags/diagsoilpem .nc $iPEM
-        cleanfiles "out_PEM/run" "" $iPEM
-        cleanfiles starts/restart1D_postPEM .txt $iPEM
-        cleanfiles starts/restart_postPEM .nc $iPEM
-        cleanfiles starts/restartfi_postPEM .nc $iPEM
-        cleanfiles starts/restartpem .nc $iPEM
-        cp starts/restartpem${iPEM}.nc startpem.nc
+        cleanfiles diags/diagpem .nc $(($iPEM - 1))
+        cleanfiles diags/diagsoilpem .nc $(($iPEM - 1))
+        cleanfiles "out_PEM/run" "" $(($iPEM - 1))
+        cleanfiles starts/restart1D_postPEM .txt $(($iPEM - 1))
+        cleanfiles starts/restart_postPEM .nc $(($iPEM - 1))
+        cleanfiles starts/restartfi_postPEM .nc $(($iPEM - 1))
+        cleanfiles starts/restartpem .nc $(($iPEM - 1))
+        cp starts/restartpem$(($iPEM - 1)).nc startpem.nc
         if [ $il -eq $(($nPCM - 1)) ]; then # Second to last PCM run
             i_myear=$(($(awk "NR==$iPEM {print \$1}" "info_PEM.txt") + $il))