Changeset 4160 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Mar 31, 2026, 3:38:07 PM (10 days ago)

Author:

jbclement

Message:

Mars PCM:
Giving standard explicit names to variables related to "paleoclimate" module + adding yearly average of flux exchanged with subsurface ice in XIOS output file intended for the PEM.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90 (modified) (3 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, h2o_ice_depth, lag_co2_ice, d_coef, zdqsdif_ssi_tot
+use paleoclimate_mod,    only: paleoclimate, h2oice_depth, co2ice_depth, coef_ssdif, zdqsdif_ssi_tot
 use comcstfi_h,          only: pi
 use geometry_mod,        only: latitude
@@ -908 +908 @@
   if (startphy_file) then
 ! Depth of H2O lag
-   call get_field("h2o_ice_depth",h2o_ice_depth,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: Failed loading <h2o_ice_depth> : ", &
-                          "<h2o_ice_depth> is set as -1 (no subsurface ice)"
-     h2o_ice_depth(:,:) = -1.
+   call get_field("h2oice_depth",h2oice_depth,found,indextime)
+   if (.not.found) then
+     write(*,*) "phyetat0: Failed loading <h2oice_depth> : ", &
+                          "<h2oice_depth> is set as -1 (no subsurface ice)"
+     h2oice_depth(:,:) = -1.
    endif
    
 ! Total flux with SSI
-   call get_field("zdqsdif_ssi_tot",zdqsdif_ssi_tot,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: Failed loading <zdqsdif_ssi_tot> : ", &
-                          "<zdqsdif_ssi_tot> is set as -1 (no subsurface ice)"
+   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 -1 (no subsurface ice)"
      zdqsdif_ssi_tot(:,:) = 0.
    endif
 
 ! Diffusion coeficent
-   call get_field("d_coef",d_coef,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: Failed loading <d_coef> : ", &
-                          "<d_coef> is set as 4e-4 (defualt value)"
-     d_coef(:,:) = 4e-4
+   call get_field("coef_ssdif",coef_ssdif,found,indextime)
+   if (.not.found) then
+     write(*,*) "phyetat0: Failed loading <coef_ssdif> : ", &
+                          "<coef_ssdif> is set as 4e-4 (default value)"
+     coef_ssdif(:,:) = 4.e-4
    endif
 
 ! Depth of CO2 lag
-   call get_field("lag_co2_ice",lag_co2_ice,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: Failed loading <lag_co2_ice> : ", &
-                          "<lag_co2_ice> is set as -1 (no subsurface ice)"
-     lag_co2_ice(:,:) = -1.
+   call get_field("co2ice_depth",co2ice_depth,found,indextime)
+   if (.not.found) then
+     write(*,*) "phyetat0: Failed loading <co2ice_depth> : ", &
+                          "<co2ice_depth> is set as -1 (no subsurface ice)"
+     co2ice_depth(:,:) = -1.
    endif
 
@@ -948 +948 @@
     endif ! not found
   else ! no startphy_file
-    h2o_ice_depth = -1.
-    lag_co2_ice = -1.
-    zdqsdif_ssi_tot = 0.
-    d_coef = 4.e-4
-    perennial_co2ice = 0.
+    h2oice_depth(:,:) = -1.
+    co2ice_depth(:,:) = -1.
+    zdqsdif_ssi_tot(:,:) = 0.
+    coef_ssdif(:,:) = 4.e-4
+    perennial_co2ice(:,:) = 0.
     if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) perennial_co2ice(ngrid,:) = 10*1.6e3 ! 10m which is convert to kg/m^2
   endif !startphy_file
 else
-   h2o_ice_depth = -1.
-   lag_co2_ice = -1.
-   zdqsdif_ssi_tot = 0.
-   d_coef = 4.e-4
-   perennial_co2ice = 0.
+   h2oice_depth(:,:) = -1.
+   co2ice_depth(:,:) = -1.
+   zdqsdif_ssi_tot(:,:) = 0.
+   coef_ssdif(:,:) = 4.e-4
+   perennial_co2ice(:,:) = 0.
 endif !paleoclimate