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

Timestamp:

Apr 18, 2023, 3:16:47 PM (20 months ago)

Author:

llange

Message:

PEM
Modification in recomp_tend_co2_slope. The emissivity of ice was assumed to be 0.95, but it can be lower. It is now set to the values given by the PCM
Remove unused variables (n_1km)

LL

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• comsoil_h_PEM.F90 (modified) (1 diff)
• constants_marspem_mod.F90 (modified) (1 diff)
• pemetat0.F90 (modified) (1 diff)
• recomp_tend_co2_slope.F90 (modified) (3 diffs)
• update_soil.F90 (modified) (1 diff)

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

@@ -3 +3 @@
 implicit none
   integer, parameter :: nsoilmx_PEM = 27               ! number of layers in the PEM
-  integer, parameter :: n_1km = 23                     ! index at which we have overcome 1km
-
   real,save,allocatable,dimension(:) :: layer_PEM      ! soil layer depths   [m]
   real,save,allocatable,dimension(:) :: mlayer_PEM     ! soil mid-layer depths [m]
   real,save,allocatable,dimension(:,:,:) :: TI_PEM     ! soil thermal inertia [SI]
-    real,save,allocatable,dimension(:,:) :: inertiedat_PEM ! soil thermal inertia saved as reference for current climate [SI]
+  real,save,allocatable,dimension(:,:) :: inertiedat_PEM ! soil thermal inertia saved as reference for current climate [SI]
   ! variables (FC: built in firstcall in soil.F)
   REAL,SAVE,ALLOCATABLE :: tsoil_PEM(:,:,:)       ! sub-surface temperatures [K]

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

@@ -3 +3 @@
 
     IMPLICIT NONE
-
+! Duration of a year and day 
+      INTEGER,PARAMETER :: sols_per_my =669 ! Number of Sols per year
+      REAL, PARAMETER :: sec_per_sol=88775.  ! Duration of a sol, in seconds
 
 ! Molecular masses for CO2,H2O and non condensible gaz, following Franz et al. 2017    
-       REAL,PARAMETER :: m_co2 = 44.01E-3  ! CO2 molecular mass (kg/mol)   
-       REAL,PARAMETER :: m_noco2 = 33.37E-3  ! Non condensible mol mass (kg/mol)   
-       REAL,PARAMETER :: m_h2o = 18.01528E-3      ! Molecular weight of h2o (kg/mol)
+      REAL,PARAMETER :: m_co2 = 44.01E-3  ! CO2 molecular mass (kg/mol)   
+      REAL,PARAMETER :: m_noco2 = 33.37E-3  ! Non condensible mol mass (kg/mol)   
+      REAL,PARAMETER :: m_h2o = 18.01528E-3      ! Molecular weight of h2o (kg/mol)
 
 !     Coefficient for Clapeyron law for CO2 condensation temperature (Tco2 = beta/(alpha-log(vmr)),following James et al. 1992

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

@@ -5 +5 @@
    
    use iostart_PEM, only:  open_startphy, close_startphy, get_field, get_var
-   use comsoil_h_PEM, only: soil_pem,layer_PEM, mlayer_PEM,n_1km,fluxgeo,inertiedat_PEM,water_reservoir_nom,depth_breccia,depth_bedrock,index_breccia,index_bedrock
+   use comsoil_h_PEM, only: soil_pem,layer_PEM, mlayer_PEM,fluxgeo,inertiedat_PEM,water_reservoir_nom,depth_breccia,depth_bedrock,index_breccia,index_bedrock
    use comsoil_h, only:  volcapa,inertiedat 
    use adsorption_mod, only : regolith_adsorption,adsorption_pem

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

@@ -2 +2 @@
 ! $Id $
 !
-SUBROUTINE recomp_tend_co2_slope(tendencies_co2_ice_phys,tendencies_co2_ice_phys_ini,co2ice_slope,vmr_co2_gcm,vmr_co2_pem,ps_GCM_2,global_ave_press_GCM,global_ave_press_new,timelen,ngrid,nslope)
+SUBROUTINE recomp_tend_co2_slope(ngrid,nslope,timelen,tendencies_co2_ice_phys,tendencies_co2_ice_phys_ini,co2ice_slope, emissivity_slope, &
+                                 vmr_co2_gcm,vmr_co2_pem,ps_GCM_2,global_ave_press_GCM,global_ave_press_new)
 
-      use constants_marspem_mod, only : alpha_clap_co2, beta_clap_co2, sigmaB,Lco2
+      use constants_marspem_mod, only : alpha_clap_co2, beta_clap_co2, sigmaB,Lco2, sols_per_my, sec_per_sol
 
       IMPLICIT NONE
@@ -21 +22 @@
   REAL, INTENT(in) ::  vmr_co2_gcm(ngrid,timelen)                ! physical point field : Volume mixing ratio of co2 in the first layer
   REAL, INTENT(in) ::  vmr_co2_pem(ngrid,timelen)                ! physical point field : Volume mixing ratio of co2 in the first layer
-  REAL, intent(in) :: ps_GCM_2(ngrid,timelen)                 ! physical point field : Surface pressure in the GCM
-  REAL, intent(in) :: global_ave_press_GCM
-  REAL, intent(in) :: global_ave_press_new
+  REAL, intent(in) :: ps_GCM_2(ngrid,timelen)                    ! physical point field : Surface pressure in the GCM
+  REAL, intent(in) :: global_ave_press_GCM                       ! global averaged pressure at previous timestep
+  REAL, intent(in) :: global_ave_press_new                       ! global averaged pressure at current timestep
   REAL, intent(in) ::  tendencies_co2_ice_phys_ini(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year
-  REAL, intent(in) :: co2ice_slope(ngrid,nslope)
+  REAL, intent(in) :: co2ice_slope(ngrid,nslope)                 ! CO2 ice per mesh and sub-grid slope(kg/m^2)
+  REAL, intent(in) :: emissivity_slope(ngrid,nslope)            ! Emissivity per mesh and sub-grid slope(1)
 
 !   OUTPUT
@@ -34 +36 @@
 
   INTEGER :: i,t,islope
-  REAL :: eps, coef, ave
-
-  eps=0.95
-  coef=669*88875*eps*sigmaB/Lco2
+  REAL ::  coef, ave
 
 ! Evolution of the water ice for each physical point
   do i=1,ngrid
     do islope=1,nslope
+      coef=sols_per_my*sec_per_sol*emissivity_slope(i,islope)*sigmaB/Lco2
       ave=0.
 !      if(abs(tendencies_co2_ice_phys(i,islope)).gt.1e-4) then

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

@@ -2 +2 @@
 #ifndef CPP_STD
  USE comsoil_h, only:  inertiedat, volcapa
- USE comsoil_h_PEM, only: layer_PEM,n_1km,inertiedat_PEM,depth_breccia,depth_bedrock,index_breccia,index_bedrock,reg_thprop_dependp
+ USE comsoil_h_PEM, only: layer_PEM,inertiedat_PEM,depth_breccia,depth_bedrock,index_breccia,index_bedrock,reg_thprop_dependp
  USE vertical_layers_mod, ONLY: ap,bp
- USE comsoil_h_PEM, only: n_1km
  USE constants_marspem_mod,only: TI_breccia,TI_bedrock, TI_regolith_avg
  implicit none