Index: LMDZ6/trunk/libf/phylmdiso/lscp_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/lscp_mod.F90	(revision 4003)
+++ LMDZ6/trunk/libf/phylmdiso/lscp_mod.F90	(revision 4003)
@@ -0,0 +1,1 @@
+link ../phylmd/lscp_mod.F90
Index: LMDZ6/trunk/libf/phylmdiso/lscp_tools_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/lscp_tools_mod.F90	(revision 4003)
+++ LMDZ6/trunk/libf/phylmdiso/lscp_tools_mod.F90	(revision 4003)
@@ -0,0 +1,1 @@
+link ../phylmd/lscp_tools_mod.F90
Index: LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90	(revision 4002)
+++ LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90	(revision 4003)
@@ -557,4 +557,6 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: cldemi, cldfra, cldtau, fiwc, fl, re, flwc
 !$OMP THREADPRIVATE(cldemi, cldfra, cldtau, fiwc, fl, re, flwc)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: qlth, qith
+!$OMP THREADPRIVATE(qlth, qith)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq, ref_ice, theta, zphi
 !$OMP THREADPRIVATE(ref_liq, ref_ice, theta, zphi)
@@ -898,4 +900,6 @@
       ALLOCATE(rain_lsc(klon))
       ALLOCATE(rain_num(klon))
+      ALLOCATE(qlth(klon,klev), qith(klon,klev))
+      !
 #ifdef ISO
       ALLOCATE(xtevap(ntraciso,klon))
@@ -1268,4 +1272,5 @@
       DEALLOCATE(rain_lsc)
       DEALLOCATE(rain_num)
+      DEALLOCATE(qlth, qith)
 !
       DEALLOCATE(sens_x, sens_w)
Index: LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90	(revision 4002)
+++ LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90	(revision 4003)
@@ -75,4 +75,5 @@
     USE VERTICAL_LAYERS_MOD, ONLY: aps,bps, ap, bp
     USE write_field_phy
+    USE lscp_mod, ONLY : lscp
 
     !USE cmp_seri_mod
@@ -1059,4 +1060,6 @@
     !IM cf. AM 081204 BEG
     LOGICAL ptconvth(klon,klev)
+
+    REAL picefra(klon,klev)
     !IM cf. AM 081204 END
     !
@@ -4562,4 +4565,19 @@
 !ISO
 
+    picefra(:,:)=0.
+
+    IF (ok_new_lscp) THEN
+
+    CALL lscp(phys_tstep,paprs,pplay, &
+         t_seri, q_seri,ptconv,ratqs, &
+         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, cldliq, picefra, &
+         rain_lsc, snow_lsc, &
+         pfrac_impa, pfrac_nucl, pfrac_1nucl, &
+         frac_impa, frac_nucl, beta_prec_fisrt, &
+         prfl, psfl, rhcl,  &
+         zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cld_th, &
+         iflag_ice_thermo)
+
+    ELSE
     CALL fisrtilp(phys_tstep,paprs,pplay, &
          t_seri, q_seri,ptconv,ratqs, &
@@ -4576,4 +4594,5 @@
 #endif      
      &   )
+    ENDIF
     !
     WHERE (rain_lsc < 0) rain_lsc = 0.
@@ -5295,5 +5314,5 @@
           ENDIF
           CALL newmicro (flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, &
-               paprs, pplay, t_seri, cldliq, cldfra, &
+               paprs, pplay, t_seri, cldliq, picefra, cldfra, &
                cldtau, cldemi, cldh, cldl, cldm, cldt, cldq, &
                flwp, fiwp, flwc, fiwc, &
@@ -5303,5 +5322,5 @@
        ELSE
           CALL nuage (paprs, pplay, &
-               t_seri, cldliq, cldfra, cldtau, cldemi, &
+               t_seri, cldliq, picefra, cldfra, cldtau, cldemi, &
                cldh, cldl, cldm, cldt, cldq, &
                ok_aie, &