Index: LMDZ6/trunk/libf/phylmd/newmicro.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/newmicro.F90	(revision 3112)
+++ LMDZ6/trunk/libf/phylmd/newmicro.F90	(revision 3113)
@@ -10,5 +10,5 @@
   USE dimphy
   USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
-    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra
+    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra, zfice
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
@@ -152,5 +152,4 @@
   REAL zrho(klon, klev) !--rho pour la couche
   REAL dh(klon, klev) !--dz pour la couche
-  REAL zfice(klon, klev)
   REAL rad_chaud(klon, klev) !--rayon pour les nuages chauds
   REAL rad_chaud_pi(klon, klev) !--rayon pour les nuages chauds pre-industriels
Index: LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90	(revision 3112)
+++ LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90	(revision 3113)
@@ -175,4 +175,6 @@
       REAL, SAVE, ALLOCATABLE :: icc3dcon(:,:)
       !$OMP THREADPRIVATE(icc3dcon)
+      REAL, SAVE, ALLOCATABLE :: zfice(:,:)
+      !$OMP THREADPRIVATE(zfice)
       REAL, SAVE, ALLOCATABLE :: od443aer(:) 
       !$OMP THREADPRIVATE(od443aer) 
@@ -603,4 +605,5 @@
       ALLOCATE(icc3dcon(klon, klev))
       ALLOCATE(icc3dstra(klon, klev))
+      ALLOCATE(zfice(klon, klev))
       ALLOCATE(od443aer(klon))
       ALLOCATE(od550aer(klon))
@@ -883,4 +886,5 @@
       DEALLOCATE(icc3dcon)
       DEALLOCATE(icc3dstra)
+      DEALLOCATE(zfice)
       DEALLOCATE(od443aer)
       DEALLOCATE(od550aer)
Index: LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3112)
+++ LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3113)
@@ -1453,4 +1453,10 @@
     'icc3dstra', 'Mass Fraction of Stratiform Cloud Ice', &
     'kg kg-1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_cldicemxrat = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'cldicemxrat', 'Cloud Ice Mixing Ratio', &
+    '1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_cldwatmxrat = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'cldwatmxrat', 'Cloud Water Mixing Ratio', &
+    '1', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_solbnd = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'solbnd', 'Top-of-Atmosphere Solar Insolation for each band', 'W m-2', (/ ('', i=1, 10) /))
Index: LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 3112)
+++ LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 3113)
@@ -123,5 +123,6 @@
          o_rsdsdiff, o_rsdscsdiff, o_stratomask,&
          o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
-         o_lcc3dstra, o_icc3dcon, o_icc3dstra, o_reffclwtop, o_ec550aer, &
+         o_lcc3dstra, o_icc3dcon, o_icc3dstra, &
+         o_cldicemxrat, o_cldwatmxrat, o_reffclwtop, o_ec550aer, &
          o_lwcon, o_iwcon, o_temp, o_theta, &
          o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop, &
@@ -280,5 +281,6 @@
          scdnc, cldncl, reffclws, reffclwc, cldnvi, &
          PSFSWDIR, PSFSWDIF, stratomask,&
-         lcc, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra, reffclwtop, &
+         lcc, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra, &
+         zfice, reffclwtop, &
          ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
          ql_seri, qs_seri, tr_seri, &
@@ -1410,4 +1412,7 @@
           CALL histwrite_phy(o_icc3dcon, icc3dcon)
           CALL histwrite_phy(o_icc3dstra, icc3dstra)
+          CALL histwrite_phy(o_cldicemxrat, zfice)
+          zx_tmp_fi3d(:,:)=1-zfice(:,:)
+          CALL histwrite_phy(o_cldwatmxrat, zx_tmp_fi3d)
           CALL histwrite_phy(o_reffclwtop, reffclwtop)
        ENDIF
