Changeset 3983 for trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90

Timestamp:

Dec 8, 2025, 11:27:43 AM (7 days ago)

Author:

jbclement

Message:

PEM:

• Removing completely the ice metamorphism computed by a threshold at the end of the PCM (which was commented).
• Addition of a module "metamorphism" to compute the PCM frost at the PEM beginning and give it back to the PCM at the PEM end. The frost is considered as the ice given by the PCM "startfi.nc" which is above the yearly minimum. Thereby, metamorphism is performed through this operation.
• Ice reservoirs representation in the PEM is modernized.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90 (modified) (2 diffs)

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

@@ -27 +27 @@
 use comsoil_h_pem,         only: soil_pem, fluxgeo, ini_huge_h2oice, depth_breccia, depth_bedrock, reg_thprop_dependp
 use adsorption_mod,        only: adsorption_pem
-use glaciers_mod,          only: h2oice_flow, co2ice_flow, inf_h2oice_threshold, metam_co2ice_threshold, metam_h2oice_threshold, metam_h2oice, metam_co2ice
+use glaciers_mod,          only: h2oice_flow, co2ice_flow, inf_h2oice_threshold
 use ice_table_mod,         only: icetable_equilibrium, icetable_dynamic
 use layering_mod,          only: layering_algo, d_dust, impose_dust_ratio, dust2ice_ratio
@@ -166 +166 @@
 call getin('inf_h2oice_threshold',inf_h2oice_threshold)
 
-metam_h2oice = .false.
-call getin('metam_h2oice',metam_h2oice)
-
-metam_h2oice_threshold = 460. ! kg.m-2 (= 0.5 m)
-call getin('metam_h2oice_threshold',metam_h2oice_threshold)
-
 h2oice_flow = .true.
 call getin('h2oice_flow',h2oice_flow)
-
-metam_co2ice = .false.
-call getin('metam_co2ice',metam_co2ice)
-
-metam_co2ice_threshold = 16500. ! kg.m-2 (= 10 m)
-call getin('metam_co2ice_threshold',metam_co2ice_threshold)
 
 co2ice_flow = .true.