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

Timestamp:

Feb 8, 2024, 2:50:21 PM (11 months ago)

Author:

jbclement

Message:

Mars PCM:
Modification of r3200 for the 1D model: the 'CO2cond_ps' coefficient (now defined in "co2condens_mod.F") is applied not only to update the surface pressure but also all the tracers and the winds from now on. It is done in "physiq_mod.F" just after the call to 'co2condens'.
JBC

File:

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

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

@@ -69 +69 @@
 use compute_soiltemp_mod,       only: compute_tsoil_pem
 use writediagpem_mod,           only: writediagpem, writediagsoilpem
+use co2condens_mod,             only: CO2cond_ps
 
 #ifndef CPP_STD
@@ -136 +137 @@
 real, dimension(:,:,:), allocatable :: q             ! champs advectes
 real, dimension(ip1jmp1)            :: ps            ! pression au sol
-real                                :: CO2cond_ps    ! Coefficient to control the surface pressure change
 real, dimension(:),     allocatable :: ps_start_PCM  ! (ngrid) surface pressure
 real, dimension(:,:),   allocatable :: ps_timeseries ! (ngrid x timelen) instantaneous surface pressure
@@ -296 +296 @@
     allocate(q(ip1jmp1,llm,nqtot))
     allocate(longitude(ngrid),latitude(ngrid),cell_area(ngrid))
-    CO2cond_ps = 1.
 #else
     allocate(q(1,llm,nqtot))
@@ -316 +315 @@
     call init_testphys1d('start1D_evol.txt','startfi_evol.nc',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,CO2cond_ps)
+                         play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
     ps(2) = ps(1)
     nsplit_phys = 1