Changeset 4166 for trunk

Timestamp:

Apr 2, 2026, 2:54:31 PM (8 days ago)

Author:

jbclement

Message:

Mars PCM:
Small correction regarding old fixed Fortran formatting following r4160 + Add missing 'watercaptag' field in "field_def_physics.xml" for XIOS following r4161.
JBC

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• deftank/xios/field_def_physics_mars.xml (modified) (1 diff)
• libf/dynphy_lonlat/phymars/newstart.F (modified) (1 diff)
• libf/phymars/phyetat0_mod.F90 (modified) (5 diffs)
• libf/phymars/soil_settings.F (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -5155 +5155 @@
 Initialize mass_predyn to -999 so that newstart writes a defined value (-999) to mass_predyn_<tracer> in restartfi.nc, avoiding uninitialized values (e.g., 0 or random memory).
 If mass_predyn_<tracer> <= 0 (such as -999 written by newstart) in startfi.nc, the mass fixer for dynamics is skipped at the first physiq step.
+
+== 02/04/2026 == JBC
+Small correction regarding old fixed Fortran formatting following r4160 + Add missing 'watercaptag' field in "field_def_physics.xml" for XIOS following r4161.

trunk/LMDZ.MARS/deftank/xios/field_def_physics_mars.xml

@@ -427 +427 @@
 
             <!-- water cycle -->
+            <field id="watercaptag"
+                   long_name="Perennial water ice location"
+                   unit="T/F" />
             <field id="watercap"
                    long_name="Perennial water ice thickness"

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -60 +60 @@
       use comslope_mod, ONLY: nslope,def_slope,def_slope_mean,
      &             subslope_dist,end_comslope_h,ini_comslope_h
-      use paleoclimate_mod, only: h2oice_depth, co2ice_depth, coef_ssdif,
-     &             ini_paleoclimate, end_paleoclimate
+      use paleoclimate_mod, only: h2oice_depth, co2ice_depth,
+     &             coef_ssdif, ini_paleoclimate, end_paleoclimate
       use subslope_mola_mod, ONLY: subslope_mola
       use parse_args_mod, only: parse_args

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

@@ -911 +911 @@
    if (.not.found) then
      write(*,*) "phyetat0: Failed loading <h2oice_depth> : ", &
-                          "<h2oice_depth> is set as -1 (no subsurface ice)"
-     h2oice_depth(:,:) = -1.
+                          "<h2oice_depth> is set as -999 (no subsurface ice)"
+     h2oice_depth(:,:) = -999.
    endif
    
@@ -919 +919 @@
    if (.not.found) then
      write(*,*) "phyetat0: Failed loading <flux_ssice> : ", &
-                          "<flux_ssice> is set as -1 (no subsurface ice)"
+                          "<flux_ssice> is set as 0 (no subsurface ice)"
      zdqsdif_ssi_tot(:,:) = 0.
    endif
@@ -935 +935 @@
    if (.not.found) then
      write(*,*) "phyetat0: Failed loading <co2ice_depth> : ", &
-                          "<co2ice_depth> is set as -1 (no subsurface ice)"
-     co2ice_depth(:,:) = -1.
+                          "<co2ice_depth> is set as -999 (no subsurface ice)"
+     co2ice_depth(:,:) = -999.
    endif
 
@@ -948 +948 @@
     endif ! not found
   else ! no startphy_file
-    h2oice_depth(:,:) = -1.
-    co2ice_depth(:,:) = -1.
+    h2oice_depth(:,:) = -999.
+    co2ice_depth(:,:) = -999.
     zdqsdif_ssi_tot(:,:) = 0.
     coef_ssdif(:,:) = 4.e-4
@@ -956 +956 @@
   endif !startphy_file
 else
-   h2oice_depth(:,:) = -1.
-   co2ice_depth(:,:) = -1.
+   h2oice_depth(:,:) = -999.
+   co2ice_depth(:,:) = -999.
    zdqsdif_ssi_tot(:,:) = 0.
    coef_ssdif(:,:) = 4.e-4

trunk/LMDZ.MARS/libf/phymars/soil_settings.F

@@ -477 +477 @@
                qsoil(:,:,igcm_h2o_ice_soil,:)= 0. 
              else
-               write(*,*) "phyetat0: suburface tracer <",trim(txt),
+               write(*,*) "phyetat0: subsurface tracer <",trim(txt),
      &                     "> range:", 
      &                     minval(qsoil(:,:,igcm_h2o_ice_soil,:)),
@@ -493 +493 @@
                qsoil(:,:,igcm_h2o_vap_ads,:)= 0. 
              else
-               write(*,*) "phyetat0: suburface tracer <",trim(txt),">
+               write(*,*) "phyetat0: subsurface tracer <",trim(txt),">
      &                     range:", 
      &                     minval(qsoil(:,:,igcm_h2o_vap_ads,:)),