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

Timestamp:

May 29, 2023, 12:41:48 PM (20 months ago)

Author:

emillour

Message:

Generic PCM:
Make number of scatterers fully dynamic (i.e. set in callphys.def
and no longer by compilation option "-s #").
One should now specify
naerkind = #
in callphys.def (default is 0).
EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (4 diffs)

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

@@ -29 +29 @@
       USE vertical_layers_mod, ONLY: init_vertical_layers
       USE logic_mod, ONLY: hybrid
+      use radinc_h, only: naerkind
       use regular_lonlat_mod, only: init_regular_lonlat
       use planete_mod, only: ini_planete_mod
@@ -154 +155 @@
 c INITIALISATION
 c=======================================================================
+! check if 'rcm1d.def' file is around 
+      open(90,file='rcm1d.def',status='old',form='formatted',
+     &     iostat=ierr)
+      if (ierr.ne.0) then
+        write(*,*) 'Cannot find required file "rcm1d.def"'
+        write(*,*) 'which should contain some input parameters'
+        write(*,*) ' ... might as well stop here ...'
+        stop
+      else
+        close(90)
+      endif
+
+! now, run.def is needed anyway. so we create a dummy temporary one
+! for ioipsl to work. if a run.def is already here, stop the
+! program and ask the user to do a bit of cleaning
+      open(90,file='run.def',status='old',form='formatted',
+     &     iostat=ierr)
+      if (ierr.eq.0) then
+        close(90)
+        write(*,*) 'There is already a run.def file.'
+        write(*,*) 'This is not compatible with 1D runs.'
+        write(*,*) 'Please remove the file and restart the run.'
+        write(*,*) 'Runtime parameters are supposed to be in rcm1d.def'
+        stop
+      else
+        call system('touch run.def')
+        call system("echo 'INCLUDEDEF=callphys.def' >> run.def")
+        call system("echo 'INCLUDEDEF=rcm1d.def' >> run.def")
+      endif
 
       ! read nq from traceur.def
@@ -187 +217 @@
       ! Initialize dimphy module
       call init_dimphy(1,llm) 
+      
       ! now initialize arrays using phys_state_var_init
+      ! but first initialise naerkind (from callphys.def)
+      naerkind=0 !default
+      call getin("naerkind",naerkind)
+      
       call phys_state_var_init(nq)
       
@@ -221 +256 @@
 c ------------------------------------------------------
 
-! check if 'rcm1d.def' file is around 
-      open(90,file='rcm1d.def',status='old',form='formatted',
-     &     iostat=ierr)
-      if (ierr.ne.0) then
-        write(*,*) 'Cannot find required file "rcm1d.def"'
-        write(*,*) 'which should contain some input parameters'
-        write(*,*) ' ... might as well stop here ...'
-        stop
-      else
-        close(90)
-      endif
-
-! now, run.def is needed anyway. so we create a dummy temporary one
-! for ioipsl to work. if a run.def is already here, stop the
-! program and ask the user to do a bit of cleaning
-      open(90,file='run.def',status='old',form='formatted',
-     &     iostat=ierr)
-      if (ierr.eq.0) then
-        close(90)
-        write(*,*) 'There is already a run.def file.'
-        write(*,*) 'This is not compatible with 1D runs.'
-        write(*,*) 'Please remove the file and restart the run.'
-        write(*,*) 'Runtime parameters are supposed to be in rcm1d.def'
-        stop
-      else
-        call system('touch run.def')
-        call system("echo 'INCLUDEDEF=callphys.def' >> run.def")
-        call system("echo 'INCLUDEDEF=rcm1d.def' >> run.def")
-      endif
 
 ! check if we are going to run with or without tracers