Index: /trunk/LMDZ.MARS/changelog.txt
===================================================================
--- /trunk/LMDZ.MARS/changelog.txt	(revision 3963)
+++ /trunk/LMDZ.MARS/changelog.txt	(revision 3964)
@@ -4995,2 +4995,6 @@
 == 30/10/2025 == JM
 Update run.def.1d with slow_diagfi option.
+
+== 19/11/2025 == EM
+Minor cleaning: don't systematically read/write "totcloudfrac" (water ice
+cloud fraction) in startfi/restartfi files. Do it only if CLFvarying is .true.
Index: /trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 3963)
+++ /trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 3964)
@@ -38,4 +38,5 @@
 use tabfi_mod, only: tabfi
 use callkeys_mod, only: startphy_file, rdstorm, hdo
+use callkeys_mod, only: CLFvarying, CLFfixval
 
 implicit none
@@ -483,13 +484,14 @@
             minval(dtau), maxval(dtau)
 
-! Sub-grid cloud fraction
-if (startphy_file) then
+! Sub-grid water ice cloud fraction
+totcloudfrac(:)=1.0 ! default value; no cloud fraction
+if (startphy_file) then
+  if (CLFvarying) then
    call get_field("totcloudfrac",totcloudfrac,found,indextime)
    if (.not.found) then
-     write(*,*) "phyetat0: <totcloudfrac> not in file WARNING put to 1"
-     totcloudfrac(:) = 1.0 !valeur par defaut (CLFfixval par defaut)
-   endif
-else
-   totcloudfrac(:)=1.0
+     write(*,*) "phyetat0: <totcloudfrac> not in file; set to ",CLFfixval
+   endif
+   totcloudfrac(:)=CLFfixval
+  endif ! of if (CLFvarying)
 endif ! if (startphy_file)
 write(*,*) "phyetat0: total cloud fraction <totcloudfrac> range:", &
Index: /trunk/LMDZ.MARS/libf/phymars/phyredem.F90
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/phyredem.F90	(revision 3963)
+++ /trunk/LMDZ.MARS/libf/phymars/phyredem.F90	(revision 3964)
@@ -207,4 +207,5 @@
   use paleoclimate_mod,    only: paleoclimate
   use callkeys_mod,        only: calltherm, dustinjection, calllott_nonoro
+  use callkeys_mod, only: CLFvarying
   use paleoclimate_mod,    only: paleoclimate, h2o_ice_depth, lag_co2_ice, d_coef, zdqsdif_ssi_tot
 
@@ -274,5 +275,7 @@
 
   ! Sub-grid cloud fraction
-  call put_field("totcloudfrac","Total cloud fraction",totcloudfrac,time) 
+  if (CLFvarying) then
+    call put_field("totcloudfrac","Total cloud fraction",totcloudfrac,time)
+  endif 
  
   ! Dust conversion factor