Changeset 3912 for trunk/LMDZ.COMMON/libf

Timestamp:

Sep 5, 2025, 2:46:19 PM (5 months ago)

Author:

jbclement

Message:

Mars PCM 1D:
Prescription of the atmospheric water ice profile in 1D. The boolean 'ctrl_h2oice' activates the option, the file "profile_ref_h2ovap" defines the prescribed profile and the value 'relaxtime_h2oice', if positive, activates the relaxation and gives the constant of time.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (2 diffs)

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

@@ -276 +276 @@
 
     ! Dummy variables to use the subroutine 'init_testphys1d'
-    logical                             :: therestart1D, therestartfi, prescribed_h2ovap
+    logical                             :: therestart1D, therestartfi, ctrl_h2ovap, ctrl_h2oice
     integer                             :: ndt, day0
-    real                                :: ptif, pks, day, gru, grv, h2ovap_relax_time
+    real                                :: ptif, pks, day, gru, grv, relaxtime_h2ovap, relaxtime_h2oice
     real, dimension(:),     allocatable :: zqsat
     real, dimension(:,:,:), allocatable :: dq, dqdyn
-    real, dimension(nlayer)             :: play, w, q_prescribed_h2o_vap
+    real, dimension(nlayer)             :: play, w, qref_h2ovap, qref_h2oice
     real, dimension(nlayer + 1)         :: plev
 #else
@@ -380 +380 @@
     call init_testphys1d('start1D.txt','startfi.nc',therestart1D,therestartfi,ngrid,nlayer,610.,nq,q,         &
                          time_0,pdyn(1),ucov,vcov,teta,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w, &
-                         play,plev,latitude,longitude,cell_area,prescribed_h2ovap,h2ovap_relax_time,q_prescribed_h2o_vap)
+                         play,plev,latitude,longitude,cell_area,                                              &
+                         ctrl_h2ovap,relaxtime_h2ovap,qref_h2ovap,ctrl_h2oice,relaxtime_h2oice,qref_h2oice)
     nsplit_phys = 1
     day_step = steps_per_sol