Changeset 3749 for trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

Timestamp:

Apr 30, 2025, 5:54:41 PM (8 weeks ago)

Author:

afalco

Message:

Pluto: ecritphy changed into diagfi_output_rate (followup of generic model), which defines the output rate of the diagfi in physical timesteps rather than dynamical ones.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

@@ -18 +18 @@
   use comgeomfi_h, only: totarea, totarea_planet
   use comsoil_h, only: ini_comsoil_h, nsoilmx, lay1_soil, alpha_soil
-  use time_phylmdz_mod, only: ecritphy,day_step,iphysiq, &
+  use time_phylmdz_mod, only: diagfi_output_rate, &
                               init_time, daysec, dtphys
   use comcstfi_mod, only: rad, cpp, g, r, rcp, &
@@ -83 +83 @@
   REAL SSUM
 
-  ! deprecated parameter
+  ! deprecated parameters
+  REAL :: ecritphy ! to check that this obsolete flag is no longer used...
   logical aerohaze
 
@@ -104 +105 @@
 #endif
 
+
   ! read in some parameters from "run.def" for physics,
   ! or shared between dynamics and physics.
+  ecritphy=-666 ! dummy default value
   call getin_p("ecritphy",ecritphy) ! frequency of outputs in physics,
                                     ! in dynamical steps
-  call getin_p("day_step",day_step) ! number of dynamical steps per day
-  call getin_p("iphysiq",iphysiq) ! call physics every iphysiq dyn step
+  if (ecritphy/=-666) then
+    call abort_physic(rname, &
+         "Error: parameter ecritphy is obsolete! Remove it! "//&
+         "And use diagfi_output_rate instead",1)
+  endif
 
   ! do we read a startphy.nc file? (default: .true.)
@@ -164 +170 @@
      call getin_p("tracer",tracer)
      if (is_master) write(*,*) trim(rname)//": tracer = ",tracer
+
+     if (is_master) write(*,*) trim(rname)//&
+       ": Output rate for diagfi.nc file (in physics steps) ?"
+     diagfi_output_rate=24 !default value
+     call getin_p("diagfi_output_rate",diagfi_output_rate)
+     if (is_master) write(*,*) trim(rname)//": diagfi_output_rate = ",&
+                               diagfi_output_rate
 
      if (is_master) write(*,*) trim(rname)//&