Changeset 3188 for trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

Timestamp:

Jan 29, 2024, 4:18:28 PM (13 months ago)

Author:

jbclement

Message:

Mars PCM:
Addition of the parameter 'CO2cond_ps' (= 1 by default) for 1D. This coefficient controls the surface pressure change. If 'relaxcoef_ps = 1', then surface pressure varies normally. If 'relaxcoef_ps = 0', then surface pressure is kept constant. The ratio of polar cap surface over planetary surface is a typical value (8.3e-4) for tests. To be defined in "callphys.def" so that both PCM and PEM can read it.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (3 diffs)

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

@@ -7 +7 @@
 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)
+                           play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau,CO2cond_ps)
 
 use ioipsl_getincom,          only: getin ! To use 'getin'
@@ -78 +78 @@
 real, dimension(1),                  intent(out) :: latitude, longitude, cell_area
 real,                                intent(out) :: atm_wat_profile, atm_wat_tau ! Force atmospheric water profiles
+real,                                intent(out) :: CO2cond_ps                 ! Relaxation coefficient for psurf
 
 !=======================================================================
@@ -308 +309 @@
 write(*,*) " psurf = ",psurf
 
+! Coefficient to control the surface pressure change
+! To be defined in "callphys.def" so that the PEM can read it asa well
+CO2cond_ps = 1. ! Default value
+write(*,*) 'Coefficient to control the surface pressure change?'
+call getin("CO2cond_ps",CO2cond_ps)
+if (CO2cond_ps < 0. .or. CO2cond_ps > 1.) then
+    write(*,*) 'Value for ''CO2cond_ps'' is not between 0 and 1.'
+    error stop 'Please, specify a correct value!'
+endif
+
 ! Reference pressures
 pa = 20.     ! transition pressure (for hybrid coord.)