Changeset 3076 for trunk/LMDZ.COMMON/libf/evolution/constants_marspem_mod.F90

Timestamp:

Oct 6, 2023, 5:32:11 PM (14 months ago)

Author:

jbclement

Message:

PEM:
Big cleaning/improvements of the PEM:

• Conversion of "abort_pem.F" and "soil_settings_PEM.F" into Fortran 90;
• Transformation of every PEM subroutines into module;
• Rewriting of many subroutines with modern Fortran syntax;
• Correction of a bug in "pem.F90" when calling 'recomp_tend_co2_slope'. The arguments were given in disorder and emissivity was missing;
• Update of "launch_pem.sh" in deftank.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/constants_marspem_mod.F90 (modified) (1 diff)

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

@@ -1 @@
-    MODULE constants_marspem_mod
+MODULE constants_marspem_mod
 
+implicit none
 
-    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
+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
-      REAL,PARAMETER :: alpha_clap_co2 = 23.3494  !Uniteless,James et al. 1992
-      REAL,PARAMETER :: beta_clap_co2 = 3182.48   !Kelvin, James et al. 1992
+real, parameter :: alpha_clap_co2 = 23.3494 ! Uniteless, James et al. 1992
+real, parameter :: beta_clap_co2 = 3182.48  ! Kelvin, James et al. 1992
 
 !     Coefficient for Clapeyron law for psat (psat = exp(beta/Th2o+alpha)),following Murphy and Koop 2005
-      REAL,PARAMETER :: alpha_clap_h2o = 28.9074  ! Uniteless, Murphy and Koop 2005
-      REAL,PARAMETER :: beta_clap_h2o = -6143.7   ! Kelvin, Murphy and Koop 2005
+real, parameter :: alpha_clap_h2o = 28.9074 ! Uniteless, Murphy and Koop 2005
+real, parameter :: beta_clap_h2o = -6143.7  ! Kelvin, Murphy and Koop 2005
 
 !     Density of the regolith (Zent et al., 1995, Buhler and Piqueux 2021)     
-      REAL,PARAMETER ::  rho_regolith = 2000.     ! kg/m^3
+real, parameter :: rho_regolith = 2000. ! kg/m^3
 
 !     Average  Thermal inertia of the surface, breccia, bedrock, following Mellon et al., 2000., Wood et al., 2008
-      REAL,PARAMETER :: TI_regolith_avg = 250.        ! Averaged of the observed thermal inertia for regolith following Mellon et al., 2000[SI]
-      REAL,PARAMETER :: TI_breccia = 750.             ! Thermal inertia of Breccia following Wood 2009 [SI]
-      REAL,PARAMETER :: TI_bedrock = 2300.            ! Thermal inertia of Bedrock following Wood 2009 [SI]
+real, parameter :: TI_regolith_avg = 250. ! Averaged of the observed thermal inertia for regolith following Mellon et al., 2000[SI]
+real, parameter :: TI_breccia = 750.      ! Thermal inertia of Breccia following Wood 2009 [SI]
+real, parameter :: TI_bedrock = 2300.     ! Thermal inertia of Bedrock following Wood 2009 [SI]
 
 !     Porosity of the soil
-      REAL,PARAMETER :: porosity = 0.4                ! porosity of the martian soil, correspond to the value for a random loose packing of monodiperse sphere (Scott, 1960)
+real, parameter :: porosity = 0.4 ! porosity of the martian soil, correspond to the value for a random loose packing of monodiperse sphere (Scott, 1960)
 
 !     Stefan Boltzmann constant
-      REAL,PARAMETER :: sigmaB=5.678E-8 
+real, parameter :: sigmaB = 5.678e-8 
 
 !     Latent heat of CO2
-      REAL,PARAMETER ::  Lco2 =  5.71E5                ! Pilorget and Forget 2016
+real, parameter :: Lco2 =  5.71e5 ! Pilorget and Forget 2016
 
 ! Conversion H2O/CO2 frost to perenial frost and vice versa
-      REAL,PARAMETER :: threshold_water_frost2perenial = 1000. !~ 1m
-      REAL,PARAMETER :: threshold_co2_frost2perenial = 3200.   !~ 2m
+real, parameter :: threshold_water_frost2perenial = 1000. !~ 1m
+real, parameter :: threshold_co2_frost2perenial = 3200.   !~ 2m
 
+END MODULE constants_marspem_mod
 
-      
-    END MODULE constants_marspem_mod
-