Changeset 3129 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Nov 17, 2023, 2:38:58 PM (13 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.COMMON/libf/evolution/pem.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/pem.F90

@@ -258 +258 @@
 
     ! Dummy variables to use the subroutine 'init_testphys1d'
-    logical                             :: startfiles_1D, therestart1D, therestartfi
+    logical                             :: therestart1D, therestartfi
     integer                             :: ndt, day0
     real                                :: ptif, pks, day, gru, grv, atm_wat_profile, atm_wat_tau
@@ -308 +308 @@
     allocate(q(1,llm,nqtot))
     allocate(longitude(1),latitude(1),cell_area(1))
-    call init_testphys1d(.true.,ngrid,nlayer,610.,nq,q,time_0,ps(1),ucov,vcov,teta,startfiles_1D,therestart1D, &
-                         therestartfi,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,                   &
+
+    therestart1D = .false.
+    inquire(file = 'start1D_evol.txt',exist = therestart1D)
+    if (.not. therestart1D) then
+        write(*,*) 'There is no "start1D_evol.txt" file!'
+        error stop 'Initialization cannot be done for the 1D PEM.'
+    endif
+    therestartfi = .false.
+    inquire(file = 'startfi_evol.nc',exist = therestartfi)
+    if (.not. therestartfi) then
+        write(*,*) 'There is no "startfi_evol.nc" file!'
+        error stop 'Initialization cannot be done for the 1D PEM.'
+    endif
+
+    call init_testphys1d('start1D_evol.txt','startfi_evol.nc',.true.,therestart1D,therestartfi,ngrid,nlayer,610., &
+                         nq,q,time_0,ps(1),ucov,vcov,teta,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,  &
                          play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
     ps(2) = ps(1)