Changeset 4170 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.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/phyetat0_mod.F90 (modified) (4 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.