Index: /trunk/LMDZ.COMMON/libf/evolution/comsoil_h_PEM.F90 =================================================================== --- /trunk/LMDZ.COMMON/libf/evolution/comsoil_h_PEM.F90 (revision 2944) +++ /trunk/LMDZ.COMMON/libf/evolution/comsoil_h_PEM.F90 (revision 2945) @@ -3,10 +3,8 @@ 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] Index: /trunk/LMDZ.COMMON/libf/evolution/constants_marspem_mod.F90 =================================================================== --- /trunk/LMDZ.COMMON/libf/evolution/constants_marspem_mod.F90 (revision 2944) +++ /trunk/LMDZ.COMMON/libf/evolution/constants_marspem_mod.F90 (revision 2945) @@ -3,10 +3,12 @@ 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 Index: /trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 =================================================================== --- /trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 (revision 2944) +++ /trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 (revision 2945) @@ -5,5 +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 Index: /trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_slope.F90 =================================================================== --- /trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_slope.F90 (revision 2944) +++ /trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_slope.F90 (revision 2945) @@ -2,7 +2,8 @@ ! $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,9 +22,10 @@ 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,12 +36,10 @@ 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 Index: /trunk/LMDZ.COMMON/libf/evolution/update_soil.F90 =================================================================== --- /trunk/LMDZ.COMMON/libf/evolution/update_soil.F90 (revision 2944) +++ /trunk/LMDZ.COMMON/libf/evolution/update_soil.F90 (revision 2945) @@ -2,7 +2,6 @@ #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