Changeset 3161 for trunk/LMDZ.COMMON/libf/evolution
- Timestamp:
- Dec 19, 2023, 3:14:48 PM (11 months ago)
- Location:
- trunk/LMDZ.COMMON/libf/evolution
- Files:
-
- 8 edited
-
changelog.txt (modified) (1 diff)
-
comsoil_h_PEM.F90 (modified) (1 diff)
-
conf_pem.F90 (modified) (3 diffs)
-
constants_marspem_mod.F90 (modified) (1 diff)
-
deftank/run_PEM.def (modified) (3 diffs)
-
glaciers_mod.F90 (modified) (1 diff)
-
pem.F90 (modified) (15 diffs)
-
pemetat0.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/changelog.txt
r3159 r3161 179 179 - Addition of (CO2 and H2O) ice metamorphism: if the frost is above a given value, then the excess frost is transferred to the perennial ice. 180 180 - Thresholds value for ice management can now be set in the "run_PEM.def". 181 182 == 19/12/2023 == JBC 183 PEM: 184 - Addition of flags defined in the "run_PEM.def" to decide to do or not CO2 & H2O ice metamorphism: 'metam_co2ice' and 'metam_h2oice' (default is false). 185 - The variations of infinite reservoirs ('watercap') during the PCM years are now taken into account to update H2O ice at the PEM initialization. 186 - 'ini_h2o_bigreservoir' is renamed into 'ini_huge_h2oice'. 187 - Some cleanings, in particular for the main program "pem.F90". -
trunk/LMDZ.COMMON/libf/evolution/comsoil_h_PEM.F90
r3149 r3161 3 3 implicit none 4 4 5 integer, parameter :: nsoilmx_PEM = 68 6 real, allocatable, dimension(:), save :: layer_PEM 7 real, allocatable, dimension(:), save :: mlayer_PEM 8 real, allocatable, dimension(:,:,:), save :: TI_PEM 9 real, allocatable, dimension(:,:), save :: inertiedat_PEM 5 integer, parameter :: nsoilmx_PEM = 68 ! number of layers in the PEM 6 real, allocatable, dimension(:), save :: layer_PEM ! soil layer depths [m] 7 real, allocatable, dimension(:), save :: mlayer_PEM ! soil mid-layer depths [m] 8 real, allocatable, dimension(:,:,:), save :: TI_PEM ! soil thermal inertia [SI] 9 real, allocatable, dimension(:,:), save :: inertiedat_PEM ! soil thermal inertia saved as reference for current climate [SI] 10 10 ! Variables (FC: built in firstcall in soil.F) 11 real, allocatable, dimension(:,:,:), save :: tsoil_PEM 12 real, allocatable, dimension(:,:), save :: mthermdiff_PEM 13 real, allocatable, dimension(:,:), save :: thermdiff_PEM 14 real, allocatable, dimension(:), save :: coefq_PEM 15 real, allocatable, dimension(:,:), save :: coefd_PEM 16 real, allocatable, dimension(:,:), save :: alph_PEM 17 real, allocatable, dimension(:,:), save :: beta_PEM 18 real, save :: mu_PEM 19 real, save :: fluxgeo 20 real, save :: depth_breccia 21 real, save :: depth_bedrock 22 integer, save :: index_breccia 23 integer, save :: index_bedrock 24 logical :: soil_pem 25 real, save :: ini_h 2o_bigreservoir ! Initial value for the large reservoir of water[kg/m^2]26 logical, save :: reg_thprop_dependp 11 real, allocatable, dimension(:,:,:), save :: tsoil_PEM ! sub-surface temperatures [K] 12 real, allocatable, dimension(:,:), save :: mthermdiff_PEM ! (FC) mid-layer thermal diffusivity [SI] 13 real, allocatable, dimension(:,:), save :: thermdiff_PEM ! (FC) inter-layer thermal diffusivity [SI] 14 real, allocatable, dimension(:), save :: coefq_PEM ! (FC) q_{k+1/2} coefficients [SI] 15 real, allocatable, dimension(:,:), save :: coefd_PEM ! (FC) d_k coefficients [SI] 16 real, allocatable, dimension(:,:), save :: alph_PEM ! (FC) alpha_k coefficients [SI] 17 real, allocatable, dimension(:,:), save :: beta_PEM ! beta_k coefficients [SI] 18 real, save :: mu_PEM ! mu coefficient [SI] 19 real, save :: fluxgeo ! Geothermal flux [W/m^2] 20 real, save :: depth_breccia ! Depth at which we have breccia [m] 21 real, save :: depth_bedrock ! Depth at which we have bedrock [m] 22 integer, save :: index_breccia ! last index of the depth grid before having breccia 23 integer, save :: index_bedrock ! last index of the depth grid before having bedrock 24 logical :: soil_pem ! True by default, to run with the subsurface physic. Read in pem.def 25 real, save :: ini_huge_h2oice ! Initial value for huge reservoirs of H2O ice [kg/m^2] 26 logical, save :: reg_thprop_dependp ! thermal properites of the regolith vary with the surface pressure 27 27 28 28 !======================================================================= -
trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90
r3159 r3161 25 25 use time_evol_mod, only: year_bp_ini, dt_pem, h2o_ice_crit, co2_ice_crit, ps_criterion, Max_iter_pem, & 26 26 evol_orbit_pem, var_obl, var_ecc, var_lsp, convert_years 27 use comsoil_h_pem, only: soil_pem, fluxgeo, ini_h 2o_bigreservoir, depth_breccia, depth_bedrock, reg_thprop_dependp27 use comsoil_h_pem, only: soil_pem, fluxgeo, ini_huge_h2oice, depth_breccia, depth_bedrock, reg_thprop_dependp 28 28 use adsorption_mod, only: adsorption_pem 29 use glaciers_mod, only: co2_ice_flow, h2o_ice_flow29 use glaciers_mod, only: h2oice_flow, co2ice_flow, inf_h2oice_threshold, metam_co2ice_threshold, metam_h2oice_threshold, metam_h2oice, metam_co2ice 30 30 use ice_table_mod, only: icetable_equilibrium, icetable_dynamic 31 31 use time_phylmdz_mod, only: ecritphy 32 use constants_marspem_mod, only: inf_h2oice_threshold, metam_co2ice_threshold, metam_h2oice_threshold33 32 34 33 implicit none … … 101 100 call getin('adsorption_pem',adsorption_pem) 102 101 103 co2_ice_flow = .true.104 call getin('co2_ice_flow',co2_ice_flow)105 106 h2o_ice_flow = .true.107 call getin('h2o_ice_flow',h2o_ice_flow)108 109 102 reg_thprop_dependp = .false. 110 103 call getin('reg_thprop_dependp',reg_thprop_dependp) … … 158 151 159 152 !#---------- Ice management parameters ----------# 160 ini_h 2o_bigreservoir= 1.e4161 call getin('ini_h 2o_bigreservoir',ini_h2o_bigreservoir)153 ini_huge_h2oice = 1.e4 154 call getin('ini_huge_h2oice',ini_huge_h2oice) 162 155 163 156 inf_h2oice_threshold = 2.e3 164 157 call getin('inf_h2oice_threshold',inf_h2oice_threshold) 165 158 159 metam_h2oice = .false. 160 call getin('metam_h2oice',metam_h2oice) 161 166 162 metam_h2oice_threshold = 5.e-2 167 163 call getin('metam_h2oice_threshold',metam_h2oice_threshold) 164 165 h2oice_flow = .true. 166 call getin('h2oice_flow',h2oice_flow) 167 168 metam_co2ice = .false. 169 call getin('metam_co2ice',metam_co2ice) 168 170 169 171 metam_co2ice_threshold = 16.e3 170 172 call getin('metam_co2ice_threshold',metam_co2ice_threshold) 171 173 174 co2ice_flow = .true. 175 call getin('co2ice_flow',co2ice_flow) 176 172 177 END SUBROUTINE conf_pem 173 178 -
trunk/LMDZ.COMMON/libf/evolution/constants_marspem_mod.F90
r3159 r3161 37 37 real, parameter :: Lco2 = 5.71e5 ! Pilorget and Forget 2016 38 38 39 ! Thresholds for ice management40 real, save :: inf_h2oice_threshold ! To consider the amount of H2O ice as an infinite reservoir41 real, save :: metam_h2oice_threshold ! To consider frost is becoming perennial H2O ice42 real, save :: metam_co2ice_threshold ! To consider frost is becoming perennial CO2 ice43 44 39 END MODULE constants_marspem_mod -
trunk/LMDZ.COMMON/libf/evolution/deftank/run_PEM.def
r3159 r3161 28 28 # Max_iter_pem=100000000 29 29 30 # Acceptance rate of sublimating H2 oice surface change? Default = 0.230 # Acceptance rate of sublimating H2O ice surface change? Default = 0.2 31 31 # h2o_ice_crit=0.2 32 32 … … 47 47 # adsorption_pem=.true. 48 48 49 # Do you want the glaciers to flow along subslope inside a mesh? Default = .true.50 # co2glaciersflow=.true.51 52 49 # Do you want to modify the soil thermal properties with the pressure? Default = .false. 53 50 # reg_thprop_dependp=.false. 54 51 55 52 #---------- Layering parameters ----------# 56 # Value of the geothermal flux? Default = 0.53 # Value of the geothermal flux? Default = 0. 57 54 # fluxgeo=0. 58 55 … … 70 67 71 68 #---------- Ice management parameters ----------# 72 # Amount of h2o ice to initialize the bigreservoir if the variable is not present in "startfi_PEM.nc"? Default = 1.e473 # ini_h 2o_bigreservoir=1.e469 # Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startfi_PEM.nc"? Default = 1.e4 70 # ini_huge_h2oice=1.e4 74 71 75 72 # Threshold to consider the amount of H2O ice as an infinite reservoir? Default = 2.e3 76 73 # inf_h2oice_threshold=2.e3 77 74 75 # Do you want H2O frost to transform into perennial H2O ice? Default = .false. 76 # metam_h2oice=.false. 77 78 78 # Threshold to consider frost is becoming perennial H2O ice? Default = 5.e-2 79 79 # metam_h2oice_threshold=5.e-2 80 81 # Do you want the H2O ice to flow along subslope inside a cell? Default = .true. 82 # h2oice_flow=.true. 83 84 # Do you want CO2 frost to transform into perennial CO2 ice? Default = .false. 85 # metam_co2ice=.false. 80 86 81 87 # Threshold to consider frost is becoming perennial CO2 ice? Default = 16.e3 82 88 # metam_co2ice_threshold=16.e3 83 89 90 # Do you want the CO2 ice to flow along subslope inside a cell? Default = .true. 91 # co2ice_flow=.true. 92 84 93 # Some definitions for the physics, in file 'callphys.def' 85 94 INCLUDEDEF=callphys.def -
trunk/LMDZ.COMMON/libf/evolution/glaciers_mod.F90
r3149 r3161 3 3 implicit none 4 4 5 logical :: co2_ice_flow ! True by default, to compute co2 ice flow. Read in run_PEM.def 6 logical :: h2o_ice_flow ! True by default, to compute h2o ice flow. Read in run_PEM.def 5 6 ! Flags for ice management 7 logical :: h2oice_flow ! True by default, to compute H2O ice flow. Read in "run_PEM.def" 8 logical :: co2ice_flow ! True by default, to compute CO2 ice flow. Read in "run_PEM.def" 9 logical :: metam_h2oice ! False by default, to compute H2O ice metamorphism. Read in "run_PEM.def" 10 logical :: metam_co2ice ! False by default, to compute CO2 ice metamorphism. Read in "run_PEM.def" 11 12 ! Thresholds for ice management 13 real, save :: inf_h2oice_threshold ! To consider the amount of H2O ice as an infinite reservoir 14 real, save :: metam_h2oice_threshold ! To consider frost is becoming perennial H2O ice 15 real, save :: metam_co2ice_threshold ! To consider frost is becoming perennial CO2 ice 7 16 8 17 !======================================================================= -
trunk/LMDZ.COMMON/libf/evolution/pem.F90
r3159 r3161 1 1 !------------------------ 2 2 ! I Initialization 3 ! I_a Read run.def4 ! I_b Read of start_evol.nc and starfi_evol.nc3 ! I_a Read the "run.def" 4 ! I_b Read the "start_evol.nc" and "startfi_evol.nc" 5 5 ! I_c Subslope parametrisation 6 ! I_d Read PCM data and convertto the physical grid6 ! I_d Read the PCM data and convert them to the physical grid 7 7 ! I_e Initialization of the PEM variable and soil 8 8 ! I_f Compute tendencies 9 9 ! I_g Save initial PCM situation 10 ! I_h Read the startpem.nc10 ! I_h Read the "startpem.nc" 11 11 ! I_i Compute orbit criterion 12 12 … … 22 22 ! III Output 23 23 ! III_a Update surface value for the PCM start files 24 ! III_b Write restart_evol.nc and restartfi_evol.nc25 ! III_c Write restartpem.nc24 ! III_b Write the "restart_evol.nc" and "restartfi_evol.nc" 25 ! III_c Write the "restartpem.nc" 26 26 !------------------------ 27 27 28 28 PROGRAM pem 29 29 30 use phyetat0_mod, 31 use phyredem, 32 use netcdf, 33 use turb_mod, 34 use comslope_mod, 35 use logic_mod, 36 use mod_const_mpi, 30 use phyetat0_mod, only: phyetat0 31 use phyredem, only: physdem0, physdem1 32 use netcdf, only: nf90_open, NF90_NOWRITE, nf90_get_var, nf90_inq_varid, nf90_close 33 use turb_mod, only: q2, wstar 34 use comslope_mod, only: nslope, def_slope, def_slope_mean, subslope_dist, iflat, major_slope, ini_comslope_h 35 use logic_mod, only: iflag_phys 36 use mod_const_mpi, only: COMM_LMDZ 37 37 use infotrac 38 use geometry_mod, 39 use conf_pem_mod, 40 use pemredem, 41 use glaciers_mod, only: flow_co2glaciers, flow_h2oglaciers, co2_ice_flow, h2o_ice_flow42 use stopping_crit_mod, only: stopping_crit_h2o_ice, stopping_crit_co2 43 use constants_marspem_mod, only: alpha_clap_co2, beta_clap_co2, alpha_clap_h2o, beta_clap_h2o, m_co2, &44 m_noco2, inf_h2oice_threshold, metam_co2ice_threshold, metam_h2oice_threshold 45 use evol_ice_mod, 46 use comsoil_h_PEM, 47 48 49 50 use adsorption_mod, 51 52 53 use time_evol_mod, 54 use orbit_param_criterion_mod, 55 use recomp_orb_param_mod, 56 use ice_table_mod, 57 58 use soil_thermalproperties_mod, 59 use time_phylmdz_mod, 60 use abort_pem_mod, 61 use soil_settings_PEM_mod, 62 use compute_tend_mod, 63 use info_PEM_mod, 64 use interpol_TI_PEM2PCM_mod, only: interpol_TI_PEM2PCM65 use nb_time_step_PCM_mod, 66 use pemetat0_mod, 67 use read_data_PCM_mod, 68 use recomp_tend_co2_slope_mod, 69 use soil_pem_compute_mod, 70 use writediagpem_mod, 38 use geometry_mod, only: latitude_deg 39 use conf_pem_mod, only: conf_pem 40 use pemredem, only: pemdem0, pemdem1 41 use glaciers_mod, only: flow_co2glaciers, flow_h2oglaciers, co2ice_flow, h2oice_flow, inf_h2oice_threshold, & 42 metam_h2oice_threshold, metam_co2ice_threshold, metam_h2oice, metam_co2ice 43 use stopping_crit_mod, only: stopping_crit_h2o_ice, stopping_crit_co2 44 use constants_marspem_mod, only: alpha_clap_co2, beta_clap_co2, alpha_clap_h2o, beta_clap_h2o, m_co2, m_noco2 45 use evol_ice_mod, only: evol_co2_ice, evol_h2o_ice 46 use comsoil_h_PEM, only: soil_pem, ini_comsoil_h_PEM, end_comsoil_h_PEM, nsoilmx_PEM, & 47 TI_PEM, inertiedat_PEM, & ! soil thermal inertia 48 tsoil_PEM, mlayer_PEM, layer_PEM, & ! Soil temp, number of subsurface layers, soil mid layer depths 49 fluxgeo ! Geothermal flux for the PEM and PCM 50 use adsorption_mod, only: regolith_adsorption, adsorption_pem, & ! Bool to check if adsorption, main subroutine 51 ini_adsorption_h_PEM, end_adsorption_h_PEM, & ! Allocate arrays 52 co2_adsorbded_phys, h2o_adsorbded_phys ! Mass of co2 and h2O adsorbded 53 use time_evol_mod, only: dt_pem, evol_orbit_pem, Max_iter_pem, convert_years, year_bp_ini 54 use orbit_param_criterion_mod, only: orbit_param_criterion 55 use recomp_orb_param_mod, only: recomp_orb_param 56 use ice_table_mod, only: porefillingice_depth, porefillingice_thickness, end_ice_table_porefilling, & 57 ini_ice_table_porefilling, computeice_table_equilibrium,compute_massh2o_exchange_ssi 58 use soil_thermalproperties_mod, only: update_soil_thermalproperties 59 use time_phylmdz_mod, only: daysec, dtphys, day_end 60 use abort_pem_mod, only: abort_pem 61 use soil_settings_PEM_mod, only: soil_settings_PEM 62 use compute_tend_mod, only: compute_tend 63 use info_PEM_mod, only: info_PEM 64 use interpol_TI_PEM2PCM_mod, only: interpol_TI_PEM2PCM 65 use nb_time_step_PCM_mod, only: nb_time_step_PCM 66 use pemetat0_mod, only: pemetat0 67 use read_data_PCM_mod, only: read_data_PCM 68 use recomp_tend_co2_slope_mod, only: recomp_tend_co2_slope 69 use soil_pem_compute_mod, only: soil_pem_compute 70 use writediagpem_mod, only: writediagpem 71 71 72 72 #ifndef CPP_STD … … 291 291 !------------------------ 292 292 ! I Initialization 293 ! I_a Read run.def293 ! I_a Read the "run.def" 294 294 !------------------------ 295 295 #ifndef CPP_1D … … 328 328 !------------------------ 329 329 ! I Initialization 330 ! I_b Read of start_evol.ncand starfi_evol.nc331 !------------------------ 332 ! I_b.1 Read start_evol.nc330 ! I_b Read of the "start_evol.nc" and starfi_evol.nc 331 !------------------------ 332 ! I_b.1 Read "start_evol.nc" 333 333 allocate(ps_start_PCM(ngrid)) 334 334 #ifndef CPP_1D … … 355 355 356 356 ! In the PCM, these values are given to the physic by the dynamic. 357 ! Here we simply read them in the startfi_evol.ncfile357 ! Here we simply read them in the "startfi_evol.nc" file 358 358 status = nf90_open(startfi_name, NF90_NOWRITE, ncid) 359 359 … … 372 372 status = nf90_close(ncid) 373 373 374 ! I_b.2 Read startfi_evol.nc374 ! I_b.2 Read the "startfi_evol.nc" 375 375 ! First we read the initial state (starfi.nc) 376 376 #ifndef CPP_STD … … 468 468 !------------------------ 469 469 ! I Initialization 470 ! I_d Read PCM data and convertto the physical grid470 ! I_d Read the PCM data and convert them to the physical grid 471 471 !------------------------ 472 472 ! First we read the evolution of water and co2 ice (and the mass mixing ratio) over the first year of the PCM run, saving only the minimum value … … 579 579 !------------------------ 580 580 ! I Initialization 581 ! I_h Read the startpem.nc581 ! I_h Read the "startpem.nc" 582 582 !------------------------ 583 583 allocate(co2_ice(ngrid,nslope),h2o_ice(ngrid,nslope)) … … 655 655 write(*,*) 'Global average pressure new time step',global_avg_press_new 656 656 657 ! II.a.2. Old pressure levels for the timeseries, this value is deleted when unused and recreated each time (big memory consu ption)657 ! II.a.2. Old pressure levels for the timeseries, this value is deleted when unused and recreated each time (big memory consumption) 658 658 allocate(zplev_old_timeseries(ngrid,nlayer + 1,timelen)) 659 659 write(*,*) "Recomputing the old pressure levels timeserie adapted to the old pressure..." … … 725 725 ! II_c Flow of glaciers 726 726 !------------------------ 727 if (co2 _ice_flow) call flow_co2glaciers(timelen,ngrid,nslope,iflat,subslope_dist,def_slope_mean,vmr_co2_PEM_phys,ps_timeseries, &727 if (co2ice_flow) call flow_co2glaciers(timelen,ngrid,nslope,iflat,subslope_dist,def_slope_mean,vmr_co2_PEM_phys,ps_timeseries, & 728 728 global_avg_press_PCM,global_avg_press_new,co2_ice,flag_co2flow,flag_co2flow_mesh) 729 if (h2o _ice_flow) call flow_h2oglaciers(timelen,ngrid,nslope,iflat,subslope_dist,def_slope_mean,tsurf_ave,h2o_ice,flag_h2oflow,flag_h2oflow_mesh)729 if (h2oice_flow) call flow_h2oglaciers(timelen,ngrid,nslope,iflat,subslope_dist,def_slope_mean,tsurf_ave,h2o_ice,flag_h2oflow,flag_h2oflow_mesh) 730 730 731 731 !------------------------ … … 944 944 do ig = 1,ngrid 945 945 ! H2O ice metamorphism 946 if ( sum(qsurf(ig,igcm_h2o_ice,:)*subslope_dist(ig,:)/cos(pi*def_slope_mean(:)/180.)) > metam_h2oice_threshold) then946 if (metam_h2oice .and. sum(qsurf(ig,igcm_h2o_ice,:)*subslope_dist(ig,:)/cos(pi*def_slope_mean(:)/180.)) > metam_h2oice_threshold) then 947 947 h2o_ice(ig,:) = h2o_ice(ig,:) + qsurf(ig,igcm_h2o_ice,:) - metam_h2oice_threshold/subslope_dist(ig,:)*cos(pi*def_slope_mean(:)*180.) 948 948 qsurf(ig,igcm_h2o_ice,:) = metam_h2oice_threshold/subslope_dist(ig,:)*cos(pi*def_slope_mean(:)*180.) … … 961 961 962 962 ! CO2 ice metamorphism 963 if ( sum(qsurf(ig,igcm_co2,:)*subslope_dist(ig,:)/cos(pi*def_slope_mean(:)/180.)) > metam_co2ice_threshold) then963 if (metam_co2ice .and. sum(qsurf(ig,igcm_co2,:)*subslope_dist(ig,:)/cos(pi*def_slope_mean(:)/180.)) > metam_co2ice_threshold) then 964 964 perennial_co2ice(ig,:) = perennial_co2ice(ig,:) + qsurf(ig,igcm_co2,:) - metam_co2ice_threshold/subslope_dist(ig,:)*cos(pi*def_slope_mean(:)*180.) 965 965 qsurf(ig,igcm_co2,:) = metam_co2ice_threshold/subslope_dist(ig,:)*cos(pi*def_slope_mean(:)*180.) … … 1022 1022 !------------------------ 1023 1023 ! III Output 1024 ! III_b Write restart_evol.nc and restartfi_evol.nc1025 !------------------------ 1026 ! III_b.1 Write restart_evol.nc1024 ! III_b Write "restart_evol.nc" and "restartfi_evol.nc" 1025 !------------------------ 1026 ! III_b.1 Write "restart_evol.nc" 1027 1027 ptimestep = iphysiq*daysec/real(day_step)/nsplit_phys ! dtphys/nsplit_phys 1028 1028 pday = day_ini … … 1041 1041 #endif 1042 1042 1043 ! III_b.2 Write restartfi_evol.nc1043 ! III_b.2 Write the "restartfi_evol.nc" 1044 1044 #ifndef CPP_STD 1045 1045 call physdem0("restartfi_evol.nc",longitude,latitude,nsoilmx,ngrid, & … … 1063 1063 !------------------------ 1064 1064 ! III Output 1065 ! III_c Write restartpem.nc1065 ! III_c Write the "restartpem.nc" 1066 1066 !------------------------ 1067 1067 call pemdem0("restartpem.nc",longitude,latitude,cell_area,nsoilmx_PEM,ngrid, & -
trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90
r3159 r3161 13 13 14 14 use iostart_PEM, only: open_startphy, close_startphy, get_field, get_var 15 use comsoil_h_PEM, only: soil_pem, layer_PEM, mlayer_PEM, fluxgeo, inertiedat_PEM, ini_h 2o_bigreservoir, depth_breccia, depth_bedrock, index_breccia, index_bedrock15 use comsoil_h_PEM, only: soil_pem, layer_PEM, mlayer_PEM, fluxgeo, inertiedat_PEM, ini_huge_h2oice, depth_breccia, depth_bedrock, index_breccia, index_bedrock 16 16 use comsoil_h, only: volcapa, inertiedat 17 17 use adsorption_mod, only: regolith_adsorption, adsorption_pem … … 23 23 use soil_pem_compute_mod, only: soil_pem_compute 24 24 use soil_pem_ini_mod, only: soil_pem_ini 25 use comslope_mod, only: def_slope_mean, subslope_dist 25 26 26 27 #ifndef CPP_STD 27 use comcstfi_h, only: r, mugaz 28 use surfdat_h, only: watercaptag, perennial_co2ice28 use comcstfi_h, only: r, mugaz, pi 29 use surfdat_h, only: watercaptag, watercap, perennial_co2ice 29 30 #else 30 use comcstfi_mod, only: r, mugaz 31 use comcstfi_mod, only: r, mugaz, pi 31 32 #endif 32 33 … … 120 121 write(*,*)'Pemetat0: failed loading <h2o_ice>' 121 122 write(*,*)'will reconstruct the values from watercaptag' 122 write(*,*)'with default value ''ini_h 2o_bigreservoir'''123 write(*,*)'with default value ''ini_huge_h2oice''' 123 124 do ig = 1,ngrid 124 if (watercaptag(ig)) h2o_ice(ig,:) = ini_h2o_bigreservoir 125 enddo 125 if (watercaptag(ig)) h2o_ice(ig,:) = ini_huge_h2oice/subslope_dist(ig,:)*cos(pi*def_slope_mean(:)*180.) 126 enddo 127 else 128 ! The variations of infinite reservoirs during the PCM years are taken into account 129 h2o_ice = h2o_ice + watercap 126 130 endif 127 131 … … 319 323 ! h2o ice 320 324 h2o_ice = 0. 321 write(*,*)'There is no "startpem.nc" so ''h2o_ice'' is reconstructed from watercaptag with default value ''ini_h 2o_bigreservoir''.'325 write(*,*)'There is no "startpem.nc" so ''h2o_ice'' is reconstructed from watercaptag with default value ''ini_huge_h2oice''.' 322 326 do ig = 1,ngrid 323 if (watercaptag(ig)) h2o_ice(ig,:) = ini_h 2o_bigreservoir327 if (watercaptag(ig)) h2o_ice(ig,:) = ini_huge_h2oice/subslope_dist(ig,:)*cos(pi*def_slope_mean(:)*180.) 324 328 enddo 325 329
Note: See TracChangeset
for help on using the changeset viewer.
