Index: /LMDZ6/trunk/DefLists/CMIP6_ping_atmos.xml
===================================================================
--- /LMDZ6/trunk/DefLists/CMIP6_ping_atmos.xml	(revision 3120)
+++ /LMDZ6/trunk/DefLists/CMIP6_ping_atmos.xml	(revision 3121)
@@ -65,7 +65,7 @@
    <field id="CMIP6_clhcalipso"    field_ref="pclhcalipso"       /> <!-- P1 (%) cloud_area_fraction_in_atmosphere_layer : Percentage cloud cover in layer centred on 220hPa -->
    <field id="CMIP6_cli"           field_ref="iwcon"            /> <!-- P1 (kg kg-1) mass_fraction_of_cloud_ice_in_air : Includes both large-scale and convective cloud. This is calculated as the mass of cloud ice in the grid cell divided by the mass of air (including the water in all phases) in the grid cell. It includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
-   <field id="CMIP6_clic"          field_ref="dummy_XYA"        /> <!-- P2 (1.0) mass_fraction_of_convective_cloud_ice_in_air : Calculated as the mass of convective cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
+   <field id="CMIP6_clic"          field_ref="icc3dcon"        /> <!-- P2 (1.0) mass_fraction_of_convective_cloud_ice_in_air : Calculated as the mass of convective cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
    <field id="CMIP6_climodis"      field_ref="climodis"         /> <!-- P1 (%) ice_cloud_area_fraction : MODIS Ice Cloud Area Percentage -->
-   <field id="CMIP6_clis"          field_ref="dummy_XYA"        /> <!-- P2 (1.0) mass_fraction_of_stratiform_cloud_ice_in_air : Calculated as the mass of stratiform cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
+   <field id="CMIP6_clis"          field_ref="icc3dstra"        /> <!-- P2 (1.0) mass_fraction_of_stratiform_cloud_ice_in_air : Calculated as the mass of stratiform cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
    <field id="CMIP6_clisccp"       field_ref="clisccp2"         /> <!-- P1 (%) cloud_area_fraction_in_atmosphere_layer : Percentage cloud cover in optical depth categories. -->
    <field id="CMIP6_clivi"         field_ref="iwp"              /> <!-- P1 (kg m-2) atmosphere_cloud_ice_content : mass of ice water in the column divided by the area of the column (not just the area of the cloudy portion of the column). Includes precipitating frozen hydrometeors ONLY if the precipitating hydrometeor affects the calculation of radiative transfer in model. -->
Index: /LMDZ6/trunk/DefLists/field_def_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/field_def_lmdz.xml	(revision 3120)
+++ /LMDZ6/trunk/DefLists/field_def_lmdz.xml	(revision 3121)
@@ -568,4 +568,6 @@
         <field id="lcc3dcon"    long_name="Convective cloud liquid fraction"    unit="1" />
         <field id="lcc3dstra"    long_name="Stratiform cloud liquid fraction"    unit="1" />
+        <field id="icc3dcon"    long_name="Mass Fraction of Convective Cloud Ice"    unit="kg/kg" />
+        <field id="icc3dstra"    long_name="Mass Fraction of Stratiform Cloud Ice "    unit="kg/kg" />
         <field id="clwcon"    long_name="Convective Cloud Liquid water content"    unit="kg/kg" />
         <field id="Ma"    long_name="undilute adiab updraft"    unit="kg/m2/s" />
Index: /LMDZ6/trunk/DefLists/file_def_histLES_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histLES_lmdz.xml	(revision 3120)
+++ /LMDZ6/trunk/DefLists/file_def_histLES_lmdz.xml	(revision 3121)
@@ -499,4 +499,6 @@
                 <field field_ref="lcc3dcon" level="10" />
                 <field field_ref="lcc3dstra" level="10" />
+                <field field_ref="icc3dcon" level="10" />
+                <field field_ref="icc3dstra" level="10" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />
Index: /LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml	(revision 3120)
+++ /LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml	(revision 3121)
@@ -512,4 +512,6 @@
                 <field field_ref="lcc3dcon" level="6" />
                 <field field_ref="lcc3dstra" level="6" />
+                <field field_ref="icc3dcon" level="6" />
+                <field field_ref="icc3dstra" level="6" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />
Index: /LMDZ6/trunk/DefLists/file_def_histins_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histins_lmdz.xml	(revision 3120)
+++ /LMDZ6/trunk/DefLists/file_def_histins_lmdz.xml	(revision 3121)
@@ -499,4 +499,6 @@
                 <field field_ref="lcc3dcon" level="10" />
                 <field field_ref="lcc3dstra" level="10" />
+                <field field_ref="icc3dcon" level="10" />
+                <field field_ref="icc3dstra" level="10" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />
Index: /LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml	(revision 3120)
+++ /LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml	(revision 3121)
@@ -554,4 +554,6 @@
                 <field field_ref="lcc3dcon" level="2" />
                 <field field_ref="lcc3dstra" level="2" />
+                <field field_ref="icc3dcon" level="2" />
+                <field field_ref="icc3dstra" level="2" />
                 <field field_ref="clwcon" level="4" />
                 <field field_ref="Ma" level="4" />
Index: /LMDZ6/trunk/DefLists/file_def_histstn_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histstn_lmdz.xml	(revision 3120)
+++ /LMDZ6/trunk/DefLists/file_def_histstn_lmdz.xml	(revision 3121)
@@ -499,4 +499,6 @@
                 <field field_ref="lcc3dcon" level="10" />
                 <field field_ref="lcc3dstra" level="10" />
+                <field field_ref="icc3dcon" level="10" />
+                <field field_ref="icc3dstra" level="10" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />
Index: /LMDZ6/trunk/libf/phylmd/newmicro.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/newmicro.F90	(revision 3120)
+++ /LMDZ6/trunk/libf/phylmd/newmicro.F90	(revision 3121)
@@ -10,5 +10,5 @@
   USE dimphy
   USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
-    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra
+    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
@@ -586,4 +586,11 @@
         lcc3dstra(i, k) = lcc3dstra(i, k) - lcc3dcon(i, k) ! eau liquide stratiforme
         lcc3dstra(i, k) = max(lcc3dstra(i,k), 0.0)
+        !FC pour la glace (CAUSES)
+        icc3dcon(i, k) = rnebcon(i, k)*(1-phase3d(i, k))*clwcon(i, k) !  glace convective
+        icc3dstra(i, k)= pclc(i, k)*pqlwp(i, k)*(1-phase3d(i, k))
+        icc3dstra(i, k) = icc3dstra(i, k) - icc3dcon(i, k) ! glace stratiforme
+        icc3dstra(i, k) = max( icc3dstra(i, k), 0.0)
+        !FC (CAUSES)
+
         ! Compute cloud droplet radius as above in meter
         radius = 1.1*((pqlwp(i,k)*pplay(i,k)/(rd*t(i,k)))/(4./3*rpi*1000.* &
@@ -626,4 +633,8 @@
         IF (lcc3dcon(i,k)<=0.0) lcc3dcon(i, k) = 0.0
         IF (lcc3dstra(i,k)<=0.0) lcc3dstra(i, k) = 0.0
+!FC (CAUSES)
+        IF (icc3dcon(i,k)<=0.0) icc3dcon(i, k) = 0.0
+        IF (icc3dstra(i,k)<=0.0) icc3dstra(i, k) = 0.0
+!FC (CAUSES)
       END DO
       IF (reffclwtop(i)<=0.0) reffclwtop(i) = 0.0
Index: /LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90	(revision 3120)
+++ /LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90	(revision 3121)
@@ -169,4 +169,8 @@
       REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:)
       !$OMP THREADPRIVATE(lcc3dstra)
+      REAL, SAVE, ALLOCATABLE :: icc3dcon(:,:)
+      !$OMP THREADPRIVATE(icc3dcon)
+      REAL, SAVE, ALLOCATABLE :: icc3dstra(:,:)
+      !$OMP THREADPRIVATE(icc3dstra)
       REAL, SAVE, ALLOCATABLE :: od443aer(:) 
       !$OMP THREADPRIVATE(od443aer) 
@@ -593,4 +597,6 @@
       ALLOCATE(lcc3dcon(klon, klev))
       ALLOCATE(lcc3dstra(klon, klev))
+      ALLOCATE(icc3dcon(klon, klev))
+      ALLOCATE(icc3dstra(klon, klev))
       ALLOCATE(od443aer(klon))
       ALLOCATE(od550aer(klon))
@@ -868,4 +874,6 @@
       DEALLOCATE(lcc3dcon)
       DEALLOCATE(lcc3dstra)
+      DEALLOCATE(icc3dcon)
+      DEALLOCATE(icc3dstra)
       DEALLOCATE(od443aer)
       DEALLOCATE(od550aer)
Index: /LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3120)
+++ /LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3121)
@@ -809,5 +809,5 @@
     's2', 'Surface moyenne des panaches de type 2', 'm2', (/ ('', i=1, 10) /))
               
-!-------Déclenchement stochastiqu
+!-------Déclenchement stochastique
   TYPE(ctrl_out), SAVE :: o_proba_notrig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), &
     'proba_notrig', 'Probabilite de non-declenchement', ' ', (/ ('', i=1, 10) /))
@@ -1447,4 +1447,8 @@
   TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_icc3dcon = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'icc3dcon', 'Convective cloud ice fraction', '1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_icc3dstra = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'icc3dstra', 'Stratiform cloud ice fraction', '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 3120)
+++ /LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 3121)
@@ -14,5 +14,5 @@
 CONTAINS 
 
-  ! ug Routine pour définir (los du premier passageà) ET sortir les variables
+  ! ug Routine pour définir (lors du premier passage) ET sortir les variables
   SUBROUTINE phys_output_write(itap, pdtphys, paprs, pphis, &
        pplay, lmax_th, aerosol_couple,         &
@@ -122,5 +122,6 @@
          o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
          o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
-         o_lcc3dstra, o_reffclwtop, o_ec550aer, &
+         o_lcc3dstra, o_icc3dcon, o_icc3dstra, &
+         o_reffclwtop, o_ec550aer, &
          o_lwcon, o_iwcon, o_temp, o_theta, &
          o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop, &
@@ -278,5 +279,6 @@
          sollwad0_aero, toplwai_aero, sollwai_aero, &
          scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&
-         lcc, lcc3d, lcc3dcon, lcc3dstra, reffclwtop, &
+         lcc, lcc3d, lcc3dcon, lcc3dstra, &
+         icc3dcon, icc3dstra, reffclwtop, &
          ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
          ql_seri, qs_seri, tr_seri, &
@@ -1406,4 +1408,6 @@
           CALL histwrite_phy(o_lcc3dcon, lcc3dcon)
           CALL histwrite_phy(o_lcc3dstra, lcc3dstra)
+          CALL histwrite_phy(o_icc3dcon, icc3dcon)
+          CALL histwrite_phy(o_icc3dstra, icc3dstra)
           CALL histwrite_phy(o_reffclwtop, reffclwtop)
        ENDIF