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

Timestamp:

Oct 2, 2023, 2:30:47 PM (15 months ago)

Author:

jbclement

Message:

PEM:
Initialization of the PEM in 1D through the subroutine "init_testphys1d_mod.F90" + Some adaptations of the Mars PCM in 1D + Update of "launch_pem.sh" in deftank.
JBC

File:

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

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

@@ -1 +1 @@
 PROGRAM testphys1d
 
-use ioipsl_getincom,         only: getin ! To use 'getin'
 use comsoil_h,               only: inertiedat, inertiesoil, nsoilmx
 use surfdat_h,               only: albedodat, perenial_co2ice, watercap
@@ -100 +99 @@
 character(len = 44)     :: txt
 
-! RV & JBC: Use of "start1D.txt" and "startfi.nc" files
+! RV & JBC: Use of starting files for 1D
 logical :: startfiles_1D, therestart1D, therestartfi
 
@@ -121 +120 @@
 !call initcomgeomphy
 
-!------------------------------------------------------
-! 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 ...'
-    stop
-endif
-
-write(*,*)'Do you want to use "start1D.txt" and "startfi.nc" files?'
-startfiles_1D = .false.
-therestart1D = .false.
-therestartfi = .false.
-call getin("startfiles_1D",startfiles_1D)
-write(*,*) " startfiles_1D = ", startfiles_1D
-
-if (startfiles_1D) then
-    inquire(file = 'start1D.txt',exist = therestart1D)
-    if (.not. therestart1D) then
-        write(*,*) 'There is no "start1D.txt" file!'
-        write(*,*) 'Initialization is done with default values.'
-    endif
-    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
-endif
-
-call init_testphys1d(startfiles_1D,therestart1D,therestartfi,ngrid,nlayer,odpref,nq,ndt,ptif,pks,dttestphys, &
-                     q,zqsat,qsurf,dq,dqdyn,day0,day,time,psurf,tsurf,gru,grv,u,v,w,q2,play,plev,tsoil,temp, &
+call init_testphys1d(.false.,ngrid,nlayer,odpref,nq,q,time,psurf,u,v,temp,startfiles_1D,therestart1D,therestartfi, &
+                     ndt,ptif,pks,dttestphys,zqsat,qsurf,dq,dqdyn,day0,day,tsurf,gru,grv,w,q2,play,plev,tsoil,     &
                      albedo,emis,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
 
@@ -265 +231 @@
 
 ! Writing the "restart1D.txt" file for the next run
-if (startfiles_1D) call writerestart1D(psurf,tsurf,nlayer,temp,u,v,nq,noms,qsurf,q)
+if (startfiles_1D) call writerestart1D('restart1D.txt',psurf,tsurf,nlayer,temp,u,v,nq,noms,qsurf,q)
 
 write(*,*) "testphys1d: everything is cool."