Index: LMDZ6/branches/LMDZISOdevSN/libf/phylmd/phys_output_write_mod.F90
===================================================================
--- LMDZ6/branches/LMDZISOdevSN/libf/phylmd/phys_output_write_mod.F90	(revision 6108)
+++ LMDZ6/branches/LMDZISOdevSN/libf/phylmd/phys_output_write_mod.F90	(revision 6109)
@@ -3224,5 +3224,7 @@
         CALL histwrite_phy(o_xtpluc(ixt), zx_tmp_fi2d)
         CALL histwrite_phy(o_xtevap(ixt),   xtevap(ixt,:))
-        CALL histwrite_phy(o_Rland_ice(ixt), Rland_ice(ixt,:))
+        IF (ok_veget) THEN
+          CALL histwrite_phy(o_Rland_ice(ixt), Rland_ice(ixt,:))
+        ENDIF
         CALL histwrite_phy(o_xtovap(ixt),  xt_seri(ixt,:,:))
         CALL histwrite_phy(o_xtoliq(ixt), xtl_seri(ixt,:,:))
@@ -3230,7 +3232,9 @@
         !! runoff land bucket + soil isotopes ratio + river and costal flows PRSN 2025
         CALL histwrite_phy(o_xtrunoff_diag(ixt), xtrunoff_diag(ixt,:))
-        CALL histwrite_phy(o_xtriverflow(ixt),   xtriverflow(ixt,:))
-        CALL histwrite_phy(o_xtcoastalflow(ixt),   xtcoastalflow(ixt,:))
-        CALL histwrite_phy(o_Rsol(ixt),   Rsol(ixt,:))
+        IF (ok_veget) THEN
+          CALL histwrite_phy(o_xtriverflow(ixt),   xtriverflow(ixt,:))
+          CALL histwrite_phy(o_xtcoastalflow(ixt),   xtcoastalflow(ixt,:))
+          CALL histwrite_phy(o_Rsol(ixt),   Rsol(ixt,:))
+        ENDIF
 
         DO nsrf = 1, nbsrf ! ajout Camille 8 mai 2023