Changeset 3207

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

Location:

trunk

Files:

• LMDZ.COMMON/libf/evolution/pem.F90 (modified) (4 diffs)
• LMDZ.MARS/changelog.txt (modified) (1 diff)
• LMDZ.MARS/libf/phymars/co2condens_mod.F (modified) (1 diff)
• LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (4 diffs)
• LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90 (modified) (3 diffs)
• LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (2 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

trunk/LMDZ.MARS/changelog.txt

@@ -4473 +4473 @@
 == 08/02/2024 == JBC
 Fix the situation where the 1D model read the "startfi.nc" or "start1D.txt" if they were present even though 'startfiles=.false', instead of initializing by default.
+
+== 08/02/2024 == JBC
+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'.

trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F

@@ -5 +5 @@
       logical, save :: scavco2cond = .false. ! flag for using scavenging_by_co2
 !$OMP THREADPRIVATE(scavco2cond)
+      real,    save :: CO2cond_ps = 1.       ! Coefficient to control the surface pressure change
 
       CONTAINS

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

@@ -7 +7 @@
 SUBROUTINE init_testphys1d(start1Dname,startfiname,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,CO2cond_ps)
+                           play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
 
 use ioipsl_getincom,          only: getin ! To use 'getin'
@@ -43 +43 @@
 use paleoclimate_mod,         only: paleoclimate, ini_paleoclimate_h, end_paleoclimate_h
 use comslope_mod,             only: nslope, subslope_dist, ini_comslope_h, end_comslope_h
+use co2condens_mod,           only: CO2cond_ps
 ! Mostly for XIOS outputs:
 use mod_const_mpi,            only: COMM_LMDZ
@@ -79 +80 @@
 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                   ! Coefficient to control the surface pressure change
 
 !=======================================================================
@@ -308 +308 @@
 ! Coefficient to control the surface pressure change
 ! To be defined in "callphys.def" so that the PEM can read it asa well
+! If CO2cond_ps = 0, surface pressure is kept constant; If CO2cond_ps = 1, surface pressure varies normally
+! CO2cond_ps = 300*400/144.37e6 = 8.31e-4 (ratio of polar cap surface over planetary surface) is a typical value for tests
 ! --------------------------------------------------------------------
 CO2cond_ps = 1. ! Default value

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

@@ -72 +72 @@
 real, dimension(nlayer + 1)         :: plev          ! intermediate pressure levels (Pa)
 real                                :: psurf         ! Surface pressure
-real                                :: CO2cond_ps    ! Coefficient to control the surface pressure change
 real, dimension(nlayer)             :: u, v          ! zonal, meridional wind
 real                                :: gru, grv      ! prescribed "geostrophic" background wind
@@ -147 +146 @@
 call init_testphys1d('start1D.txt','startfi.nc',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,CO2cond_ps)
+                     play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
 
 ! Write a "startfi" file
@@ -259 +258 @@
     ! Compute pressure for next time step
     ! -----------------------------------
-    ! CO2cond_ps is a coefficient to control the surface pressure change
-    ! If CO2cond_ps = 0, psurf is kept constant; If CO2cond_ps = 1, psurf varies normally
-    ! CO2cond_ps = 300*400/144.37e6 = 8.31e-4 (ratio of polar cap surface over planetary surface) is a typical value for tests
-    psurf = psurf + CO2cond_ps*dttestphys*dpsurf(1) ! Surface pressure change
+    psurf = psurf + dttestphys*dpsurf(1) ! Surface pressure change
     plev = ap + psurf*bp
     play = aps + psurf*bps
 
     ! Increment tracers
-    if (abs(CO2cond_ps) < 1.e-10) dq(1,:,igcm_co2) = 0.
     q = q + dttestphys*dq
 enddo ! End of time stepping loop (idt=1,ndt)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -17 +17 @@
       use calchim_mod, only: calchim, ichemistry, zdqchim, zdqschim
       use watersat_mod, only: watersat
-      use co2condens_mod, only: co2condens
+      use co2condens_mod, only: co2condens, CO2cond_ps
       use co2cloud_mod, only: co2cloud
       use callradite_mod, only: callradite
@@ -2231 +2231 @@
      $              zdqssed_co2,zcondicea_co2microp,
      &              zdqsc)
+
+          if (ngrid == 1) then ! For the 1D model
+          ! CO2cond_ps is a coefficient to control the surface pressure change
+              pdpsrf = CO2cond_ps*pdpsrf
+              zduc   = CO2cond_ps*zduc
+              zdvc   = CO2cond_ps*zdvc
+              zdqc   = CO2cond_ps*zdqc
+          endif
 
          DO iq=1, nq