source: trunk/LMDZ.COMMON/libf/evolution/constants_marspem_mod.F90 @ 3152

Last change on this file since 3152 was 3149, checked in by jbclement, 19 months ago

PEM:

  • Simplification of the algorithm managing the stopping criteria;
  • Complete rework of the ice management in the PEM (H2O & CO2);

    Subroutines to evolve the H2O and CO2 ice are now in the same module "evol_ice_mod.F90".
    Tendencies are computed from the variation of "ice + frost" between the 2 PCM runs.
    Evolving ice in the PEM is now called 'h2o_ice' or 'co2_ice' (not anymore in 'qsurf' and free of 'water_reservoir').
    Default value 'ini_h2o_bigreservoir' (= 10 m) initializes the H2O ice of the first PEM run where there is 'watercap'. For the next PEM runs, initialization is done with the value kept in "startpem.nc". CO2 ice is taken from 'perennial_co2ice' of the PCM (paleoclimate flag must be true).
    Simplification of the condition to compute the surface ice cover needed for the stopping criteria.
    Frost ('qsurf') is not evolved by the PEM and given back to the PCM.
    New default threshold value 'inf_h2oice_threshold' (= 2 m) to decide at the end of the PEM run if the H2O ice should be 'watercap' or not for the next PCM runs. If H2O ice cannot be 'watercap', then the remaining H2O ice is transferred to the frost ('qsurf').

  • Renaming of variables/subroutines for clarity;
  • Some cleanings throughout the code;
  • Small updates in files of the deftank.

JBC
File size: 2.1 KB
Line 
1MODULE constants_marspem_mod
2
3implicit none
4
5! Duration of a year and day
6integer, parameter :: sols_per_my = 669    ! Number of Sols per year
7real,    parameter :: sec_per_sol = 88775. ! Duration of a sol, in seconds
8
9! Molecular masses for CO2,H2O and non condensible gaz, following Franz et al. 2017   
10real, parameter :: m_co2 = 44.01e-3    ! CO2 molecular mass (kg/mol)   
11real, parameter :: m_noco2 = 33.37e-3  ! Non condensible mol mass (kg/mol)   
12real, parameter :: m_h2o = 18.01528e-3 ! Molecular weight of h2o (kg/mol)
13
14! Coefficient for Clapeyron law for CO2 condensation temperature (Tco2 = beta/(alpha-log(vmr)),following James et al. 1992
15real, parameter :: alpha_clap_co2 = 23.3494 ! Uniteless, James et al. 1992
16real, parameter :: beta_clap_co2 = 3182.48  ! Kelvin, James et al. 1992
17
18! Coefficient for Clapeyron law for psat (psat = exp(beta/Th2o+alpha)),following Murphy and Koop 2005
19real, parameter :: alpha_clap_h2o = 28.9074 ! Uniteless, Murphy and Koop 2005
20real, parameter :: beta_clap_h2o = -6143.7  ! Kelvin, Murphy and Koop 2005
21
22! Density of the regolith (Zent et al., 1995, Buhler and Piqueux 2021)     
23real, parameter :: rho_regolith = 2000. ! kg/m^3
24
25! Average thermal inertia of the surface, breccia, bedrock, following Mellon et al., 2000., Wood et al., 2008
26real, parameter :: TI_regolith_avg = 250. ! Averaged of the observed thermal inertia for regolith following Mellon et al., 2000[SI]
27real, parameter :: TI_breccia = 750.      ! Thermal inertia of Breccia following Wood 2009 [SI]
28real, parameter :: TI_bedrock = 2300.     ! Thermal inertia of Bedrock following Wood 2009 [SI]
29
30! Porosity of the soil
31real, parameter :: porosity = 0.4 ! porosity of the martian soil, correspond to the value for a random loose packing of monodiperse sphere (Scott, 1960)
32
33! Stefan Boltzmann constant
34real, parameter :: sigmaB = 5.678e-8
35
36! Latent heat of CO2
37real, parameter :: Lco2 =  5.71e5 ! Pilorget and Forget 2016
38
39! Threshold to consider the amount of H2O ice as an infinite reservoir
40real, parameter :: inf_h2oice_threshold = 2000.  !~ 1 m
41
42END MODULE constants_marspem_mod
Note: See TracBrowser for help on using the repository browser.