Changeset 3129 for trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.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/init_testphys1d_mod.F90 (modified) (10 diffs)

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

@@ -5 +5 @@
 contains
 
-SUBROUTINE init_testphys1d(pem1d,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,        &
+SUBROUTINE init_testphys1d(start1Dname,startfiname,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)
 
@@ -38 +38 @@
 use phys_state_var_init_mod,  only: phys_state_var_init
 use turb_mod,                 only: q2
-use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd
-use conf_phys_mod, only: conf_phys
+use nonoro_gwd_ran_mod,       only: du_nonoro_gwd, dv_nonoro_gwd
+use conf_phys_mod,            only: conf_phys
 ! Mostly for XIOS outputs:
 use mod_const_mpi,            only: COMM_LMDZ
@@ -51 +51 @@
 ! Arguments
 !=======================================================================
-integer, intent(in) :: ngrid, nlayer
-real,    intent(in) :: odpref        ! DOD reference pressure (Pa)
-logical, intent(in) :: pem1d         ! If initialization for the 1D PEM
+integer,      intent(in) :: ngrid, nlayer
+real,         intent(in) :: odpref                                    ! DOD reference pressure (Pa)
+character(*), intent(in) :: start1Dname, startfiname                  ! Name of starting files for 1D
+logical,      intent(in) :: startfiles_1D, therestart1D, therestartfi ! Use of starting files for 1D
 
 integer, intent(inout) :: nq
@@ -62 +63 @@
 real, dimension(nlayer),             intent(out) :: u, v  ! zonal, meridional wind
 real, dimension(nlayer),             intent(out) :: temp  ! temperature at the middle of the layers
-logical,                             intent(out) :: startfiles_1D, therestart1D, therestartfi ! Use of starting files for 1D
 integer,                             intent(out) :: ndt
 real,                                intent(out) :: ptif, pks
@@ -91 +91 @@
 ! RV & JBC: Use of starting files for 1D
 logical                :: found
-character(len = 30)    :: header
+character(30)          :: header
 real, dimension(100)   :: tab_cntrl
 real, dimension(1,2,1) :: albedo_read ! surface albedo
@@ -100 +100 @@
 
 ! MVals: isotopes as in the dynamics (CRisi)
-integer                                        :: ifils, ipere, generation, ierr0
-character(len = 30), dimension(:), allocatable :: tnom_transp ! transporting fluid short name
-character(len = 80)                            :: line        ! to store a line of text
-logical                                        :: continu, there
+integer                                  :: ifils, ipere, generation, ierr0
+character(30), dimension(:), allocatable :: tnom_transp ! transporting fluid short name
+character(80)                            :: line        ! to store a line of text
+logical                                  :: continu
 
 ! LL: Possibility to add subsurface ice
@@ -113 +113 @@
 real :: flux_geo_tmp
 
-! JBC: To initialize the 1D PEM
-character(:), allocatable :: start1Dname, startfiname ! Name of starting files for 1D
-
 !=======================================================================
 ! Code
 !=======================================================================
-if (.not. pem1d) then
-    start1Dname = 'start1D.txt'
-    startfiname = 'startfi.nc'
-    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
-
-    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
-else
-    start1dname = 'start1D_evol.txt'
-    startfiname = 'startfi_evol.nc'
-    startfiles_1D = .true.
-endif
-
-therestart1D = .false.
-therestartfi = .false.
-inquire(file = start1Dname,exist = therestart1D)
-if (startfiles_1D .and. .not. therestart1D) then
-    write(*,*) 'There is no "'//start1Dname//'" file!'
-    if (.not. pem1d) then
-        write(*,*) 'Initialization is done with default values.'
-    else
-        error stop 'Initialization cannot be done for the 1D PEM.'
-    endif
-endif
-inquire(file = startfiname,exist = therestartfi)
-if (.not. therestartfi) then
-    write(*,*) 'There is no "'//startfiname//'" file!'
-    if (.not. pem1d) then
-        write(*,*) 'Initialization is done with default values.'
-    else
-        error stop 'Initialization cannot be done for the 1D PEM.'
-    endif
-endif
-
 !------------------------------------------------------
 ! Prescribed constants to be set here
@@ -309 +259 @@
     time = time/24. ! convert time (hours) to fraction of sol
 else
-    call open_startphy(startfiname)
+    call open_startphy(trim(startfiname))
     call get_var("controle",tab_cntrl,found)
     if (.not. found) then
@@ -351 +301 @@
     call getin("psurf",psurf)
 else
-    open(3,file = start1Dname,status = "old",action = "read")
+    open(3,file = trim(start1Dname),status = "old",action = "read")
     read(3,*) header, psurf
 endif
@@ -683 +633 @@
     ! check if "h2o_vap" has already been initialized
     ! (it has been if there is a "profile_h2o_vap" file around)
-    inquire(file = "profile_h2o_vap",exist = there)
-    if (there) then
+    inquire(file = "profile_h2o_vap",exist = found)
+    if (found) then
         flagh2o = 0 ! 0: do not initialize h2o_vap
     else