Index: trunk/LMDZ.MARS/changelog.txt
===================================================================
--- trunk/LMDZ.MARS/changelog.txt	(revision 4165)
+++ trunk/LMDZ.MARS/changelog.txt	(revision 4166)
@@ -5155,2 +5155,5 @@
 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.
Index: trunk/LMDZ.MARS/deftank/xios/field_def_physics_mars.xml
===================================================================
--- trunk/LMDZ.MARS/deftank/xios/field_def_physics_mars.xml	(revision 4165)
+++ trunk/LMDZ.MARS/deftank/xios/field_def_physics_mars.xml	(revision 4166)
@@ -427,4 +427,7 @@
 
             <!-- water cycle -->
+            <field id="watercaptag"
+                   long_name="Perennial water ice location"
+                   unit="T/F" />
             <field id="watercap"
                    long_name="Perennial water ice thickness"
Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F	(revision 4165)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F	(revision 4166)
@@ -60,6 +60,6 @@
       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
Index: trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 4165)
+++ trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 4166)
@@ -911,6 +911,6 @@
    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,5 +919,5 @@
    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,6 +935,6 @@
    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,6 +948,6 @@
     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,6 +956,6 @@
   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
Index: trunk/LMDZ.MARS/libf/phymars/soil_settings.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/soil_settings.F	(revision 4165)
+++ trunk/LMDZ.MARS/libf/phymars/soil_settings.F	(revision 4166)
@@ -477,5 +477,5 @@
                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,5 +493,5 @@
                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,:)),