Changeset 3042 for trunk

Timestamp:

Sep 13, 2023, 2:33:49 PM (22 months ago)

Author:

jbclement

Message:

Mars PEM:
The date and time are redefined at the end of PEM to write the file "restarfi_evol.nc". Like so, the next GCM runs can restart at the date and time given at the beginning of the simulation since the PEM is only progressing by step of one Martian year.
JBC

Location:

trunk

Files:

• LMDZ.COMMON/libf/evolution/pem.F90 (modified) (2 diffs)
• LMDZ.MARS/changelog.txt (modified) (2 diffs)
• LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (1 diff)
• LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (1 diff)

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

@@ -116 +116 @@
 integer :: day_ini   ! First day of the simulation
 real    :: pday      ! Physical day
-real    :: ptime     ! Physical time
 real    :: time_phys ! Same as GCM
 real    :: ptimestep ! Same as GCM
@@ -1136 +1135 @@
 !------------------------
 ! III_b.1 Write restart_evol.nc
-ptimestep = iphysiq*daysec/real(day_step)/nsplit_phys
-print*, ptimestep,dtphys,iphysiq*daysec/real(day_step),nsplit_phys
+ptimestep = iphysiq*daysec/real(day_step)/nsplit_phys ! dtphys/nsplit_phys
 pday = day_ini
-ptime = 0.
-ztime_fin = pday - day_end + ptime + ptimestep/(float(iphysiq)*daysec)
+ztime_fin = time_phys
 
 allocate(p(ip1jmp1,nlayer + 1))

trunk/LMDZ.MARS/changelog.txt

@@ -4184 +4184 @@
 * Some improvements for the script modify_startfi_orbit.sh;
 * Some changes in run_PEM.def and xml files.
-JBC
 
 == 11/09/2023 == JBC
@@ -4195 +4194 @@
 == 12/09/2023 == JBC
 The variable 'timeperi' (defined in "planete_h.F90" and computed in "iniorbit.F") is renamed into 'lsperi' and thus slightly changed to be coherent to the solar longitude of perihelion in radian. It can now be used out of the box by other subroutines/programs like the PEM.
+
+== 13/09/2023 == JBC
+The date and time are redefined at the end of PEM to write the file "restarfi_evol.nc". Like so, the next GCM runs can restart at the date and time given at the beginning of the simulation since the PEM is only progressing by step of one Martian year.

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -97 +97 @@
       REAL day            ! date during the run
       REAL time           ! time (0<time<1 ; time=0.5 a midi)
-      REAL dttestphys         ! testphys1d timestep
+      REAL dttestphys     ! testphys1d timestep
       REAL play(nlayer)   ! Pressure at the middle of the layers (Pa)
       REAL plev(nlayer+1) ! intermediate pressure levels (pa)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -258 +258 @@
 c ----------------------
       INTEGER,SAVE :: day_ini ! Initial date of the run (sol since Ls=0)
-      INTEGER,SAVE :: icount     ! counter of calls to physiq during the run.
+      INTEGER,SAVE :: icount  ! Counter of calls to physiq during the run.
       REAL,SAVE :: time_phys