Index: LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 4069)
+++ LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 4070)
@@ -1826,20 +1826,27 @@
 
 !--aviation & supersaturation
-       CALL histwrite_phy(o_oclr, qclr)
-       CALL histwrite_phy(o_ocld, qcld)
-       CALL histwrite_phy(o_oss, qss)
-       CALL histwrite_phy(o_ovc, qvc)
-       CALL histwrite_phy(o_rnebclr, rnebclr)
-       CALL histwrite_phy(o_rnebss, rnebss)
-       CALL histwrite_phy(o_rnebseri, rneb_seri)
-       CALL histwrite_phy(o_gammass, gamma_ss)
-       CALL histwrite_phy(o_N1_ss, N1_ss)
-       CALL histwrite_phy(o_N2_ss, N2_ss)
-       CALL histwrite_phy(o_drnebsub, drneb_sub)
-       CALL histwrite_phy(o_drnebcon, drneb_con)
-       CALL histwrite_phy(o_drnebtur, drneb_tur)
-       CALL histwrite_phy(o_drnebavi, drneb_avi)
-       CALL histwrite_phy(o_qsatl, zqsatl)
-       CALL histwrite_phy(o_qsats, zqsats)
+       IF (ok_ice_sursat) THEN
+         CALL histwrite_phy(o_oclr, qclr)
+         CALL histwrite_phy(o_ocld, qcld)
+         CALL histwrite_phy(o_oss, qss)
+         CALL histwrite_phy(o_ovc, qvc)
+         CALL histwrite_phy(o_rnebclr, rnebclr)
+         CALL histwrite_phy(o_rnebss, rnebss)
+         CALL histwrite_phy(o_rnebseri, rneb_seri)
+         CALL histwrite_phy(o_gammass, gamma_ss)
+         CALL histwrite_phy(o_N1_ss, N1_ss)
+         CALL histwrite_phy(o_N2_ss, N2_ss)
+         CALL histwrite_phy(o_drnebsub, drneb_sub)
+         CALL histwrite_phy(o_drnebcon, drneb_con)
+         CALL histwrite_phy(o_drnebtur, drneb_tur)
+         CALL histwrite_phy(o_drnebavi, drneb_avi)
+         CALL histwrite_phy(o_qsatl, zqsatl)
+         CALL histwrite_phy(o_qsats, zqsats)
+         CALL histwrite_phy(o_Tcontr, Tcontr)
+         CALL histwrite_phy(o_qcontr, qcontr)
+         CALL histwrite_phy(o_qcontr2, qcontr2)
+         CALL histwrite_phy(o_fcontrN, fcontrN)
+         CALL histwrite_phy(o_fcontrP, fcontrP)
+       ENDIF
        IF (ok_plane_contrail) THEN
          CALL histwrite_phy(o_flight_m, flight_m)
@@ -1848,9 +1855,4 @@
          CALL histwrite_phy(o_flight_h2o, flight_h2o)
        ENDIF
-       CALL histwrite_phy(o_Tcontr, Tcontr)
-       CALL histwrite_phy(o_qcontr, qcontr)
-       CALL histwrite_phy(o_qcontr2, qcontr2)
-       CALL histwrite_phy(o_fcontrN, fcontrN)
-       CALL histwrite_phy(o_fcontrP, fcontrP)
        
        IF (vars_defined) zx_tmp_fi3d = wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd