Changeset 4170 for trunk/LMDZ.MARS/libf/phymars/phyredem.F90

Timestamp:

Apr 3, 2026, 4:34:51 PM (7 days ago)

Author:

jbclement

Message:

PEM:

• Deletion of 'flux_ssice' ('zqdsdif_tot') from the "startfi.nc" as it is not needed.
• Using the yearly average flux for the sublimating subsurface ice instead of the value at last timestep.
• Keeping a clear separation between the subsurface ice flux and the surface ice tendency.
• Making sure that subsurface ice depth is well given to the PCM at the end of the PEM (ice table VS layering).

JBC

File:

• trunk/LMDZ.MARS/libf/phymars/phyredem.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/phyredem.F90

@@ -210 +210 @@
   use callkeys_mod,        only: calltherm, dustinjection, calllott_nonoro
   use callkeys_mod, only: CLFvarying
-  use paleoclimate_mod,    only: paleoclimate, h2oice_depth, co2ice_depth, coef_ssdif, zdqsdif_ssi_tot
+  use paleoclimate_mod,    only: paleoclimate, h2oice_depth, co2ice_depth, coef_ssdif
 
   implicit none
@@ -390 +390 @@
     call put_field("h2oice_depth","Depth of the shallowest H2O ice layer",h2oice_depth)
     call put_field("co2ice_depth","Depth of the shallowest CO2 ice layer",co2ice_depth)
-    call put_field("flux_ssice","Total flux exchanged with subsurface water ice",zdqsdif_ssi_tot)
     call put_field("coef_ssdif","Diffusion coefficent for subsurface water",coef_ssdif)
   endif