Changeset 3725 for trunk/LMDZ.GENERIC/libf/phystd/dyn1d

Timestamp:

Apr 15, 2025, 9:34:31 AM (9 months ago)

Author:

emillour

Message:

Generic 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 "diagfi_output_rate" to specify output rate (in physics steps) 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.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -21 +21 @@
       use comcstfi_mod, only: pi, cpp, rad, g, r,
      &                        mugaz, rcp, omeg
-      use time_phylmdz_mod, only: daysec, dtphys, day_step, ecritphy,
-     &                            nday, iphysiq
+      use time_phylmdz_mod, only: daysec, dtphys, nday,
+     &                            diagfi_output_rate
       use callkeys_mod, only: tracer,rings_shadow,
      &                        specOLR,water,pceil,ok_slab_ocean,photochem
@@ -71 +71 @@
 !==================================================================
 
-#include "dimensions.h"
-#include "paramet.h"
-!include "dimphys.h"
-#include "netcdf.inc"
-#include "comgeom.h"
+      include "dimensions.h"
+      include "paramet.h"
+      include "comgeom.h"
 
 c --------------------------------------------------------------
@@ -91 +89 @@
 c
       INTEGER day0          ! date initial (sol ; =0 a Ls=0)
-      REAL day              ! date durant le run
+      REAL day              ! date during the run
+      INTEGER day_step ! number of time steps per day
       REAL time             ! time (0<time<1 ; time=0.5 a midi)
       REAL play(llm)        ! Pressure at the middle of the layers (Pa)
@@ -596 +595 @@
       write(*,*) " day_step = ",day_step
 
-      iphysiq=1 ! in 1D model physics are called evry time step
-      ecritphy=day_step ! default value for ecritphy
+      diagfi_output_rate=day_step ! default value for diagfi_output_rate
       PRINT *,'Nunber of steps between writediagfi ?'
-      call getin("ecritphy",ecritphy)
-      write(*,*) " ecritphy = ",ecritphy
+      call getin("diagfi_output_rate",diagfi_output_rate)
+      write(*,*) " diagfi_output_rate = ",diagfi_output_rate
 
       ndt=10 ! default value for ndt