MODULE planet
!-----------------------------------------------------------------------
! NAME
!     planet
!
! DESCRIPTION
!     Mars-specific constants kept hard-coded for the PEM.
!
! AUTHORS & DATE
!     R. Vandemeulebrouck
!     JB Clement, 2023-2025
!
! NOTES
!     Only valid for the Mars planet.
!-----------------------------------------------------------------------

! DEPENDENCIES
! ------------
use numerics, only: dp

! DECLARATION
! -----------
implicit none

! PARAMETERS
! ----------
! Molecular masses for CO2, H2O and non condensible gaz, following Franz et al. 2017
real(dp), parameter :: m_co2 = 44.01e-3_dp               ! CO2 molecular mass [kg/mol]
real(dp), parameter :: m_noco2 = 33.37e-3_dp             ! Non condensible mol mass [kg/mol]
real(dp), parameter :: m_h2o = 18.01528e-3_dp            ! Molecular weight of h2o [kg/mol]
real(dp), parameter :: A = (1._dp/m_co2 - 1._dp/m_noco2) ! Intermediate parameter for computation
real(dp), parameter :: B = 1._dp/m_noco2                 ! Intermediate parameter for computation

! Thermal inertia for breccia and bedrock following Mellon et al. 2000; Wood et al. 2008
real(dp), parameter :: TI_breccia = 750._dp  ! Thermal inertia of Breccia following Wood 2009 [SI]
real(dp), parameter :: TI_bedrock = 2300._dp ! Thermal inertia of Bedrock following Wood 2009 [SI]

! Coefficient for Clapeyron law for psat (psat = exp(beta/Th2o + alpha)), following Murphy and Koop 2005
real(dp), parameter :: alpha_clap_h2o = 28.9074_dp ! Uniteless, Murphy and Koop 2005
real(dp), parameter :: beta_clap_h2o = -6143.7_dp  ! Kelvin, Murphy and Koop 2005

! Coefficient for Clapeyron law for CO2 condensation temperature (Tco2 = beta/(alpha-log(vmr)),following James et al. 1992
real(dp), parameter :: alpha_clap_co2 = 23.3494_dp ! Unitless, James et al. 1992
real(dp), parameter :: beta_clap_co2 = 3182.48_dp  ! Kelvin, James et al. 1992

! Density of the regolith (Zent et al., 1995; Buhler and Piqueux 2021)
real(dp), parameter :: rho_regolith = 2000._dp ! kg/m^3

! Average thermal inertia of the surface, breccia, bedrock, following Mellon et al., 2000., Wood et al., 2008
real(dp), parameter :: TI_regolith_avg = 250._dp ! Averaged of the observed thermal inertia for regolith following Mellon et al., 2000[SI]

! Porosity of the soil
real(dp), parameter :: porosity = 0.4_dp ! porosity of the martian soil, correspond to the value for a random loose packing of monodiperse sphere (Scott, 1960)

! Stefan Boltzmann constant
real(dp), parameter :: sigmaB = 5.678e-8_dp

! Latent heat of CO2
real(dp), parameter :: Lco2 =  5.71e5_dp ! Pilorget and Forget 2016

! Current average pressure [Pa]
real(dp), parameter :: P610 = 610._dp

contains
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

END MODULE planet