Changeset 3129 for trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

Timestamp:

Nov 17, 2023, 2:38:58 PM (14 months ago)

Author:

jbclement

Message:

Mars PCM/PEM:
Cleaning of the 1D initialization: any reference of the PEM has been removed from "init_testphys1D_mod.F90". This way is much cleaner even though it needs more code.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90 (modified) (4 diffs)

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

@@ -15 +15 @@
 use turb_mod,            only: q2
 use write_output_mod,    only: write_output
+use ioipsl_getincom,     only: getin ! To use 'getin'
 use init_testphys1d_mod, only: init_testphys1d
 use writerestart1D_mod,  only: writerestart1D
@@ -89 +90 @@
 integer                 :: nq = 1 ! number of tracers
 real, dimension(1)      :: latitude, longitude, cell_area
-logical                 :: there
 character(2)            :: str2
 character(7)            :: str7
@@ -95 +95 @@
 
 ! RV & JBC: Use of starting files for 1D
-logical :: startfiles_1D, therestart1D, therestartfi
+logical :: startfiles_1D, therestart1D, therestartfi, there
 
 ! JN & JBC: Force atmospheric water profiles
@@ -115 +115 @@
 !call initcomgeomphy
 
-call init_testphys1d(.false.,ngrid,nlayer,odpref,nq,q,time,psurf,u,v,temp,startfiles_1D,therestart1D, &
-                     therestartfi,ndt,ptif,pks,dttestphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,          &
+startfiles_1D = .false.
+!------------------------------------------------------
+! Loading run parameters from "run.def" file
+!------------------------------------------------------
+! check if 'run.def' file is around. Otherwise reading parameters
+! from callphys.def via getin() routine won't work.
+inquire(file = 'run.def',exist = there)
+if (.not. there) then
+    write(*,*) 'Cannot find required file "run.def"'
+    write(*,*) '  (which should contain some input parameters along with the following line: INCLUDEDEF=callphys.def)'
+    write(*,*) ' ... might as well stop here ...'
+    error stop
+endif
+
+startfiles_1D = .false. ! Default value
+write(*,*) 'Do you want to use starting files and/or to write restarting files?'
+call getin('startfiles_1D',startfiles_1D)
+write(*,*) 'startfiles_1D =', startfiles_1D
+
+therestart1D = .false. ! Default value
+inquire(file = 'start1D.txt',exist = therestart1D)
+if (startfiles_1D .and. .not. therestart1D) then
+    write(*,*) 'There is no "start1D.txt" file!'
+    write(*,*) 'Initialization is done with default values.'
+endif
+therestartfi = .false. ! Default value
+inquire(file = 'startfi.nc',exist = therestartfi)
+if (.not. therestartfi) then
+    write(*,*) 'There is no "startfi.nc" file!'
+    write(*,*) 'Initialization is done with default values.'
+endif
+
+call init_testphys1d('start1D.txt','startfi.nc',startfiles_1D,therestart1D,therestartfi,ngrid,nlayer,odpref, &
+                     nq,q,time,psurf,u,v,temp,ndt,ptif,pks,dttestphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,     &
                      play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)