Index: trunk/LMDZ.COMMON/libf/evolution/backup.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/backup.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/backup.F90 (revision 4135) @@ -114,37 +114,30 @@ ! LOCAL VARIABLES ! --------------- -real(dp), dimension(:,:), allocatable :: h2o_ice4PCM, co2_ice4PCM, tsurf4PCM, flux_geo4PCM, albedo4PCM, emissivity4PCM -real(dp), dimension(:,:), allocatable :: teta4PCM, air_mass4PCM -real(dp), dimension(:,:,:), allocatable :: tsoil4PCM, inertiesoil4PCM, q4PCM -real(dp), dimension(:), allocatable :: ps4PCM -real(dp) :: pa4PCM, preff4PCM -logical(k4), dimension(:), allocatable :: is_h2o_perice +real(dp), dimension(ngrid,nslope) :: h2o_ice4PCM, co2_ice4PCM, tsurf4PCM, flux_geo4PCM, albedo4PCM, emissivity4PCM +real(dp), dimension(ngrid,nlayer) :: teta4PCM, air_mass4PCM +real(dp), dimension(ngrid,nsoil_PCM,nslope) :: tsoil4PCM, inertiesoil4PCM +real(dp), dimension(ngrid,nlayer,nq) :: q4PCM +real(dp), dimension(ngrid) :: ps4PCM +logical(k4), dimension(ngrid) :: is_h2o_perice +real(dp) :: pa4PCM, preff4PCM ! CODE ! ---- ! Build ice for the PCM -allocate(h2o_ice4PCM(ngrid,nslope),co2_ice4PCM(ngrid,nslope),is_h2o_perice(ngrid)) call build4PCM_perice(h2o_ice,co2_ice,is_h2o_perice,h2o_ice4PCM,co2_ice4PCM) ! Build surface temperature for the PCM -allocate(tsurf4PCM(ngrid,nslope)) call build4PCM_tsurf(tsurf_avg,tsurf_dev,tsurf4PCM) ! Build soil for the PCM -if (do_soil) then - allocate(tsoil4PCM(ngrid,nsoil_PCM,nslope),inertiesoil4PCM(ngrid,nsoil_PCM,nslope),flux_geo4PCM(ngrid,nslope)) - call build4PCM_soil(tsoil_avg,tsoil_dev,inertiesoil4PCM,tsoil4PCM,flux_geo4PCM) -end if +if (do_soil) call build4PCM_soil(tsoil_avg,tsoil_dev,inertiesoil4PCM,tsoil4PCM,flux_geo4PCM) ! Build atmosphere for the PCM -allocate(ps4PCM(ngrid),teta4PCM(ngrid,nlayer),air_mass4PCM(ngrid,nlayer)) call build4PCM_atmosphere(ps_avg,ps_dev,ps_avg_glob,ps_avg_glob_ini,ps4PCM,pa4PCM,preff4PCM,teta4PCM,air_mass4PCM) ! Build tracers for the PCM -allocate(q4PCM(ngrid,nlayer,nq)) call build4PCM_tracers(ps4PCM,q4PCM) ! Build surface radiative properties state for the PCM -allocate(albedo4PCM(ngrid,nslope),emissivity4PCM(ngrid,nslope)) call build4PCM_surf_rad_prop(h2o_ice,co2_ice,albedo4PCM,emissivity4PCM) @@ -158,24 +151,9 @@ end if -! Deallocation -if (allocated(emissivity4PCM)) deallocate(emissivity4PCM) -if (allocated(albedo4PCM)) deallocate(albedo4PCM) -if (allocated(q4PCM)) deallocate(q4PCM) -if (allocated(ps4PCM)) deallocate(ps4PCM) -if (allocated(teta4PCM)) deallocate(teta4PCM) -if (allocated(air_mass4PCM)) deallocate(air_mass4PCM) -if (allocated(tsoil4PCM)) deallocate(tsoil4PCM) -if (allocated(inertiesoil4PCM)) deallocate(inertiesoil4PCM) -if (allocated(flux_geo4PCM)) deallocate(flux_geo4PCM) -if (allocated(tsurf4PCM)) deallocate(tsurf4PCM) -if (allocated(co2_ice4PCM)) deallocate(co2_ice4PCM) -if (allocated(h2o_ice4PCM)) deallocate(h2o_ice4PCM) -if (allocated(is_h2o_perice)) deallocate(is_h2o_perice) - END SUBROUTINE save_clim_state !======================================================================= !======================================================================= -SUBROUTINE backup_restarts(idt) +SUBROUTINE backup_restarts(backup_idt) !----------------------------------------------------------------------- ! NAME @@ -205,5 +183,5 @@ ! ARGUMENTS ! --------- -integer(di), intent(in) :: idt +integer(di), intent(in) :: backup_idt ! LOCAL VARIABLES @@ -213,7 +191,7 @@ ! CODE ! ---- -suffix = '_ts'//int2str(idt) - -call print_msg('> Backup at dt = '//int2str(idt),LVL_NFO) +suffix = '_ts'//int2str(backup_idt) + +call print_msg('> Backup of "restart" files at dt = '//int2str(backup_idt),LVL_NFO) call copy_restart_if_present('re'//startevo_name,suffix2filename(startevo_name,suffix)) call copy_restart_if_present('re'//startfi_name,suffix2filename(startfi_name,suffix)) Index: trunk/LMDZ.COMMON/libf/evolution/changelog.txt =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/changelog.txt (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/changelog.txt (revision 4135) @@ -911,2 +911,5 @@ - Add a periodic backup based on a single wrapper which builds PCM-compatible climate state and writes "restart" files. - Remove condition "if (.not. allocated(*))" for more strict safeguard. + +== 16/03/2026 == JBC +Relocate Mars-specific parameters out of "planet" module and into the modules that own their physics domain. Index: trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90 (revision 4135) @@ -319,11 +319,10 @@ use geometry, only: ngrid, nslope, nsoil, nsoil_PCM use evolution, only: dt -use physics, only: mugaz, r -use planet, only: TI_breccia, TI_bedrock, alpha_clap_h2o, beta_clap_h2o +use physics, only: mugaz, r, alpha_clap_h2o, beta_clap_h2o use frost, only: h2ofrost_PCM, h2o_frost4PCM, co2frost_PCM, co2_frost4PCM use surf_ice, only: is_h2o_perice_PCM, h2oice_huge_ini, co2_perice_PCM, threshold_h2oice_cap use soil, only: do_soil, TI, depth_breccia, depth_bedrock, index_breccia, index_bedrock, inertiedat, layer, inertiedat_PCM use soil_temp, only: ini_tsoil_pem, compute_tsoil -use soil_therm_inertia, only: update_soil_TI +use soil_therm_inertia, only: update_soil_TI, TI_breccia, TI_bedrock use ice_table, only: icetable_equilibrium, icetable_dynamic, computeice_table_equilibrium use sorption, only: do_sorption, evolve_regolith_adsorption Index: trunk/LMDZ.COMMON/libf/evolution/glaciers.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/glaciers.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/glaciers.F90 (revision 4135) @@ -339,6 +339,6 @@ ! DEPENDENCIES ! ------------ -use geometry, only: ngrid, nday -use planet, only: alpha_clap_co2, beta_clap_co2 +use geometry, only: ngrid, nday +use physics, only: alpha_clap_co2, beta_clap_co2 ! DECLARATION Index: trunk/LMDZ.COMMON/libf/evolution/ice_table.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/ice_table.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/ice_table.F90 (revision 4135) @@ -276,7 +276,6 @@ ! ------------ use geometry, only: ngrid, nslope, nsoil -use soil, only: mlayer +use soil, only: mlayer, regolith_porosity use slopes, only: subslope_dist, def_slope_mean -use planet, only: porosity use maths, only: pi @@ -306,5 +305,5 @@ do islope = 1,nslope call compute_Tice(tsoil(ig,:,islope),tsurf(ig,islope),icetable_depth(ig,islope),Tice) - delta_icetable(ig) = delta_icetable(ig) + porosity*rho_ice(Tice,'H2O')*(icetable_thickness(ig,islope) - icetable_thickness_old(ig,islope)) & ! convention > 0. <=> it condenses + delta_icetable(ig) = delta_icetable(ig) + regolith_porosity*rho_ice(Tice,'H2O')*(icetable_thickness(ig,islope) - icetable_thickness_old(ig,islope)) & ! convention > 0. <=> it condenses *subslope_dist(ig,islope)/cos(def_slope_mean(islope)*pi/180._dp) end do Index: trunk/LMDZ.COMMON/libf/evolution/layered_deposits.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/layered_deposits.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/layered_deposits.F90 (revision 4135) @@ -36,10 +36,7 @@ real(dp), parameter, private :: dry_lag_porosity = 0.2_dp ! Porosity of dust lag real(dp), parameter, private :: wet_lag_porosity = 0.15_dp ! Porosity of water ice lag -real(dp), parameter, private :: regolith_porosity = 0.4_dp ! Porosity of regolith -real(dp), parameter, private :: breccia_porosity = 0.1_dp ! Porosity of breccia real(dp), parameter, private :: co2ice_porosity = 0.05_dp ! Porosity of CO2 ice real(dp), parameter, private :: h2oice_porosity = 0.1_dp ! Porosity of H2O ice real(dp), parameter, private :: rho_dust = 2500._dp ! Density of dust [kg.m-3] -real(dp), parameter, private :: rho_rock = 3200._dp ! Density of rock [kg.m-3] real(dp), parameter, private :: h_patchy_dust = 5.e-4_dp ! Height under which a dust lag is considered patchy [m] real(dp), parameter, private :: h_patchy_ice = 5.e-4_dp ! Height under which a H2O ice lag is considered patchy [m] @@ -426,5 +423,5 @@ ! DEPENDENCIES ! ------------ -use soil, only: do_soil, layer, index_breccia, index_bedrock +use soil, only: do_soil, layer, index_breccia, index_bedrock, regolith_porosity, breccia_porosity use geometry, only: nsoil @@ -1250,4 +1247,5 @@ use utility, only: real2str use display, only: print_msg, LVL_ERR +use soil, only: regolith_porosity ! DECLARATION Index: trunk/LMDZ.COMMON/libf/evolution/physics.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/physics.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/physics.F90 (revision 4135) @@ -24,4 +24,20 @@ ! PARAMETERS ! ---------- +real(dp), parameter :: sigmaB = 5.678e-8_dp ! Stefan-Boltzmann constant [W/m^2/K^4] +real(dp), parameter :: m_co2 = 44.01e-3_dp ! CO2 molecular mass [kg/mol], Franz et al. 2017 +real(dp), parameter :: m_noco2 = 33.37e-3_dp ! Non-condensible molecular mass [kg/mol], Franz et al. 2017 +real(dp), parameter :: m_h2o = 18.01528e-3_dp ! H2O molecular mass [kg/mol], Franz et al. 2017 +real(dp), parameter :: A = 1._dp/m_co2 - 1._dp/m_noco2 ! MMR/VMR conversion factor +real(dp), parameter :: B = 1._dp/m_noco2 ! MMR/VMR conversion factor + +! H2O coefficients for Clausius-Clapeyron law for psat (psat = exp(beta/Th2o + alpha)), following Murphy and Koop 2005 +real(dp), parameter :: alpha_clap_h2o = 28.9074_dp +real(dp), parameter :: beta_clap_h2o = -6143.7_dp ! [K] +! CO2 coefficients for Clausius-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 +real(dp), parameter :: beta_clap_co2 = 3182.48_dp ! [K] + +real(dp), parameter :: Lco2 = 5.71e5_dp ! CO2 latent heat [J/kg], Pilorget and Forget 2016 + real(dp), protected :: g ! Gravity [m/s2] real(dp), protected :: r ! Reduced gas constant,r = 8.314511/[mugaz/1000.0] Index: trunk/LMDZ.COMMON/libf/evolution/planet.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/planet.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/planet.F90 (revision 4135) @@ -5,5 +5,5 @@ ! ! DESCRIPTION -! Mars-specific constants kept hard-coded for the PEM. +! Legacy compatibility module for planet-specific definitions. ! ! AUTHORS & DATE @@ -15,50 +15,10 @@ !----------------------------------------------------------------------- -! 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 +! NOTE +! ---- contains Index: trunk/LMDZ.COMMON/libf/evolution/soil.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/soil.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/soil.F90 (revision 4135) @@ -25,4 +25,9 @@ ! PARAMETERS ! ---------- +! Density and porosity of regolith and breccia +real(dp), parameter :: rho_regolith = 2000._dp ! [kg/m^3], Zent et al. 1995; Buhler and Piqueux 2021 +real(dp), parameter :: regolith_porosity = 0.4_dp ! Porosity of the martian regolith for a random loose packing of monodisperse sphere, Scott 1960 +real(dp), parameter :: breccia_porosity = 0.1_dp ! Porosity of breccia + logical(k4), protected :: do_soil ! To run with subsurface physics logical(k4), protected :: reg_thprop_dependp ! Regolith thermal properties depend on surface pressure Index: trunk/LMDZ.COMMON/libf/evolution/soil_temp.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/soil_temp.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/soil_temp.F90 (revision 4135) @@ -569,12 +569,11 @@ ! DEPENDENCIES ! ------------ -use geometry, only: ngrid, nsoil, nslope, nday -use planet, only: alpha_clap_h2o, beta_clap_h2o -use tracers, only: mmol, iPCM_qh2o -use physics, only: mugaz, r -use evolution, only: dt -use soil, only: TI -use stoppage, only: stop_clean -use display, only: print_msg, LVL_NFO +use geometry, only: ngrid, nsoil, nslope, nday +use tracers, only: mmol, iPCM_qh2o +use physics, only: mugaz, r, alpha_clap_h2o, beta_clap_h2o +use evolution, only: dt +use soil, only: TI +use stoppage, only: stop_clean +use display, only: print_msg, LVL_NFO ! DECLARATION Index: trunk/LMDZ.COMMON/libf/evolution/soil_therm_inertia.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/soil_therm_inertia.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/soil_therm_inertia.F90 (revision 4135) @@ -26,4 +26,8 @@ ! PARAMETERS ! ---------- +real(dp), parameter :: TI_regolith_avg = 250._dp ! Average of observed thermal inertia for regolith [SI], Mellon et al. 2000 +real(dp), parameter :: TI_breccia = 750._dp ! Thermal inertia for breccia [SI], Mellon et al. 2000; Wood et al. 2008 +real(dp), parameter :: TI_bedrock = 2300._dp ! Thermal inertia for bedrock [SI], Mellon et al. 2000; Wood et al. 2008 +real(dp), parameter, private :: P610 = 610._dp ! Reference mean surface pressure [Pa] real(dp), parameter, private :: reg_inertie_thresold = 370._dp ! Above this thermal inertia, the regolith has too much cementation to be dependant on the pressure [J/m^2/K/s^1/2] real(dp), parameter, private :: reg_inertie_minvalue = 50._dp ! Minimum value of the Thermal Inertia at low pressure (Piqueux & Christensen 2009) [J/m^2/K/s^1/2] @@ -56,5 +60,5 @@ ! DEPENDENCIES ! ------------ -use planet, only: porosity +use soil, only: regolith_porosity ! DECLARATION @@ -78,5 +82,5 @@ ice_thermalinertia = inertie_purewaterice else - ice_thermalinertia = sqrt(surf_thermalinertia**2 + porosity*pore_filling*inertie_purewaterice**2) ! Siegler et al., 2012 + ice_thermalinertia = sqrt(surf_thermalinertia**2 + regolith_porosity*pore_filling*inertie_purewaterice**2) ! Siegler et al., 2012 end if @@ -106,5 +110,4 @@ use geometry, only: ngrid, nslope, nsoil use soil, only: volcapa, layer, inertiedat, depth_breccia, depth_bedrock, index_breccia, index_bedrock, reg_thprop_dependp -use planet, only: TI_breccia, TI_bedrock, TI_regolith_avg, P610 use display, only: print_msg, LVL_NFO Index: trunk/LMDZ.COMMON/libf/evolution/sorption.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/sorption.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/sorption.F90 (revision 4135) @@ -139,8 +139,8 @@ ! ------------ use geometry, only: ngrid, nslope, nsoil, nday -use soil, only: layer, index_breccia +use soil, only: layer, index_breccia, rho_regolith use slopes, only: subslope_dist, def_slope_mean use atmosphere, only: ap, bp -use planet, only: alpha_clap_h2o, beta_clap_h2o, m_h2o, m_co2,m_noco2, rho_regolith +use physics, only: alpha_clap_h2o, beta_clap_h2o, m_h2o, m_co2, A, B use maths, only: pi @@ -177,5 +177,4 @@ real(dp), dimension(ngrid,nslope) :: deltam_reg_slope ! Difference density of H2O adsorbed per slope [kg/m^3] real(dp), dimension(ngrid,nsoil,nslope) :: dm_h2o_regolith_slope ! Elementary H2O mass adsorded per mesh per slope -real(dp) :: A, B ! Used to compute the mean mass above the surface real(dp) :: p_sat ! Saturated vapor pressure of ice real(dp), dimension(:,:), allocatable :: mass_mean ! Mean mass above the surface @@ -188,6 +187,4 @@ ! 0. Some initializations allocate(mass_mean(ngrid,nday),zplev_mean(ngrid,nday),pvapor(ngrid,nday),pvapor_avg(ngrid)) -A = 1._dp/m_co2 - 1._dp/m_noco2 -B = 1._dp/m_noco2 theta_h2o_ads = 0._dp dm_h2o_regolith_slope = 0._dp @@ -273,8 +270,8 @@ ! ------------ use geometry, only: ngrid, nslope, nsoil, nday -use soil, only: layer, index_breccia +use soil, only: layer, index_breccia, rho_regolith use slopes, only: subslope_dist, def_slope_mean use atmosphere, only: ap, bp -use planet, only: m_co2, m_noco2, rho_regolith +use physics, only: m_co2, A, B use maths, only: pi @@ -303,5 +300,4 @@ ! real(dp) :: as = 18.9e3_dp ! Specific area, Buhler & Piqueux 2021 real(dp) :: as = 9.48e4_dp ! Same as previous but from zent -real(dp) :: A, B ! Used to compute the mean mass above the surface integer(di) :: ig, islope, iloop, it ! For loops real(dp), dimension(ngrid,nsoil,nslope) :: dm_co2_regolith_slope ! Elementary mass adsorded per mesh per slope @@ -322,6 +318,4 @@ allocate(mass_mean(ngrid,nday),zplev_mean(ngrid,nday),pco2(ngrid,nday),pco2_avg(ngrid)) m_h2o_adsorbed = 0._dp -A = 1._dp/m_co2 - 1._dp/m_noco2 -B = 1._dp/m_noco2 dm_co2_regolith_slope = 0._dp delta_co2_ads = 0._dp Index: trunk/LMDZ.COMMON/libf/evolution/tendencies.F90 =================================================================== --- trunk/LMDZ.COMMON/libf/evolution/tendencies.F90 (revision 4134) +++ trunk/LMDZ.COMMON/libf/evolution/tendencies.F90 (revision 4135) @@ -88,9 +88,9 @@ ! DEPENDENCIES ! ------------ -use geometry, only: ngrid, nslope, nday -use planet, only: alpha_clap_co2, beta_clap_co2, sigmaB, Lco2, m_co2, A, B -use orbit, only: yr_len_sol, sol_len_s -use display, only: print_msg, LVL_NFO -use utility, only: real2str +use geometry, only: ngrid, nslope, nday +use physics, only: sigmaB, alpha_clap_co2, beta_clap_co2, Lco2, m_co2, A, B +use orbit, only: yr_len_sol, sol_len_s +use display, only: print_msg, LVL_NFO +use utility, only: real2str ! DECLARATION