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

Timestamp:

Jan 29, 2024, 4:18:28 PM (11 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.COMMON/libf/evolution/pem.F90 (modified) (4 diffs)

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

@@ -137 +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
@@ -300 +301 @@
     allocate(q(ip1jmp1,llm,nqtot))
     allocate(longitude(ngrid),latitude(ngrid),cell_area(ngrid))
+    CO2cond_ps = 1.
 #else
     allocate(q(1,llm,nqtot))
@@ -319 +321 @@
     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)
+                         play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau,CO2cond_ps)
     ps(2) = ps(1)
     nsplit_phys = 1
@@ -643 +645 @@
     do i = 1,ngrid
         do islope = 1,nslope
-            global_avg_press_new = global_avg_press_new - g*cell_area(i)*tend_co2_ice(i,islope)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180.)/Total_surface
+            global_avg_press_new = global_avg_press_new - CO2cond_ps*g*cell_area(i)*tend_co2_ice(i,islope)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180.)/Total_surface
         enddo
     enddo