Changeset 3369 for trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

Timestamp:

Jun 12, 2024, 11:00:16 AM (6 months ago)

Author:

emillour

Message:

Mars PCM:
Change the way the rate of outputs for diagfi.nc files is specified:
IMPORTANT: Specifying "ecritphy" no longer possible and will trigger an error.
Use "outputs_per_sol" to specify output rate instead.
This should makes things (hopefully) clearer for users and also better
enforces a cleaner and clearer separation between dynamics and physics.
EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -11 +11 @@
 use ioipsl_getincom,          only: getin ! To use 'getin'
 use comcstfi_h,               only: pi, rad, omeg, g, mugaz, rcp, r, cpp
-use time_phylmdz_mod,         only: daysec, day_step, ecritphy, iphysiq
+use time_phylmdz_mod, only: daysec, steps_per_sol, outputs_per_sol, iphysiq
 use planete_h,                only: year_day, periheli, aphelie, peri_day, obliquit, emin_turb, lmixmin
 use surfdat_h,                only: albedodat, z0_default, z0, emissiv, emisice, albedice, iceradius, dtemisice,      &
@@ -287 +287 @@
 nsoil = nsoilmx
 
-day_step = 48 ! Default value for day_step
+steps_per_sol = 48 ! Default value for day_step (old name for steps_per_sol)
 write(*,*)'Number of time steps per sol?'
-call getin("day_step",day_step)
-write(*,*) " day_step = ",day_step
-
-ecritphy = day_step ! Default value for ecritphy, output every time step
+call getin("day_step",steps_per_sol)
+write(*,*) " steps_per_sol (aka day_step) = ",steps_per_sol
+
+outputs_per_sol = steps_per_sol ! Default value for outputs_per_sol
 
 ndt = 10 ! Default value for ndt
@@ -300 +300 @@
 
 dayn = day0 + ndt
-ndt = ndt*day_step
-dttestphys = daysec/day_step
+ndt = ndt*steps_per_sol
+dttestphys = daysec/steps_per_sol
 
 ! Imposed surface pressure