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

Timestamp:

Apr 3, 2026, 4:34:51 PM (5 weeks 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

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• phyetat0_mod.F90 (modified) (4 diffs)
• phyredem.F90 (modified) (2 diffs)

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

@@ -37 +37 @@
 use comsoil_h,           only: flux_geo
 use comslope_mod,        only: nslope, major_slope
-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
 use comcstfi_h,          only: pi
 use geometry_mod,        only: latitude
@@ -914 +914 @@
      h2oice_depth(:,:) = -999.
    endif
-   
-! Total flux with SSI
-   call get_field("flux_ssice",zdqsdif_ssi_tot,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: Failed loading <flux_ssice> : ", &
-                          "<flux_ssice> is set as 0 (no subsurface ice)"
-     zdqsdif_ssi_tot(:,:) = 0.
-   endif
 
 ! Diffusion coeficent
@@ -950 +942 @@
     h2oice_depth(:,:) = -999.
     co2ice_depth(:,:) = -999.
-    zdqsdif_ssi_tot(:,:) = 0.
     coef_ssdif(:,:) = 4.e-4
     perennial_co2ice(:,:) = 0.
@@ -958 +949 @@
    h2oice_depth(:,:) = -999.
    co2ice_depth(:,:) = -999.
-   zdqsdif_ssi_tot(:,:) = 0.
    coef_ssdif(:,:) = 4.e-4
    perennial_co2ice(:,:) = 0.

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