Index: LMDZ6/trunk/libf/phylmd/lmdz_call_lscp.f90
===================================================================
--- LMDZ6/trunk/libf/phylmd/lmdz_call_lscp.f90	(revision 6106)
+++ LMDZ6/trunk/libf/phylmd/lmdz_call_lscp.f90	(revision 6107)
@@ -305,8 +305,4 @@
       ELSE
 
-         ! old (pre-cmip7) code
-         !$gpum nocall
-         CALL fisrtilp_first(klon, klev, dtime, pfrac_impa, pfrac_nucl, pfrac_1nucl)
-
          CALL fisrtilp(klon, klev, dtime, paprs, pplay, &
                        temp, qt, ptconv, ratqs, sigma_qtherm, &
Index: LMDZ6/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 6106)
+++ LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 6107)
@@ -56,4 +56,6 @@
     USE assert_m, only: assert
     USE change_srf_frac_mod
+    USE lmdz_lscp_old, ONLY: fisrtilp_first
+
 
     USE conf_phys_m, ONLY : conf_phys
@@ -1179,5 +1181,5 @@
     !======================================================================!
     if (debut) then                                                        !
-       print*,'DEBUT PHYSIQ_MOD TRAVAIL DEBUT 2'
+       print*,'DEBUT PHYSIQ_MOD TRAVAIL DEBUT 5'
        iflag_physiq=0
        call getin_p('iflag_physiq', iflag_physiq)                          !
@@ -1722,7 +1724,17 @@
        CALL thermcell_ini(iflag_thermals,prt_level,tau_thermals,lunout, &
             &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)
+       IF (iflag_thermals>=1) THEN
+          CALL alpale_th_first()
+       ENDIF
+
        CALL lscp_ini(pdtphys,klon,klev,iflag_thermals,iflag_cld_th,iflag_ice_thermo, & 
                   lunout,prt_level,ok_ice_supersat,fl_cor_ebil,                       &
                   RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RV,RG,RPI,EPS_W)
+       IF (ok_new_lscp) THEN
+           ! old (pre-cmip7) code
+           !$gpum nocall
+           CALL fisrtilp_first(klon, klev, pdtphys, pfrac_impa, pfrac_nucl, pfrac_1nucl)
+       ENDIF
+
        CALL blosno_ini(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG, RV, RPI)
        CALL gwd_ini(klon, klev, pplay, paprs, lunout, prt_level, RCPD, RD, RG, RPI, & 
@@ -1895,4 +1907,5 @@
              CALL ini_wake(0.,0.,it_wape_prescr,wape_prescr,fip_prescr &
                   ,alp_bl_prescr, ale_bl_prescr)
+             CALL calWAKE_first(phys_tstep)
              ! 11/09/06 rajout initialisation ALE et ALP du wake et PBL(YU)
              !        print*,'apres ini_wake iflag_cld_th=', iflag_cld_th
@@ -1914,4 +1927,9 @@
              d_dens_wk(:) = 0.
              d_dens_a_wk(:) = 0.
+
+             IF (iflag_alp_wk_cond .GT. 0.) THEN
+                CALL alpale_wk_first(cell_area)
+             ENDIF
+
           ENDIF  !  (iflag_wake>=1)
 
@@ -3500,5 +3518,4 @@
        IF (iflag_alp_wk_cond .GT. 0.) THEN
 
-          CALL alpale_wk_first(cell_area)
           CALL alpale_wk(phys_tstep, cell_area, wake_k, wake_s, wake_dens, wake_fip_0, &
                wake_fip)
@@ -3564,4 +3581,5 @@
        ENDDO
        !cc fin nrlmd le 10/04/2012
+
 
        IF (iflag_thermals>=1) THEN
@@ -3684,5 +3702,4 @@
           !
           !
-          CALL alpale_th_first()
           CALL alpale_th( phys_tstep, lmax_th, t_seri, cell_area,  &
                cin, s2, n2, strig, &
Index: LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90	(revision 6106)
+++ LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90	(revision 6107)
@@ -1258,4 +1258,5 @@
     !======================================================================!
     IF (debut) THEN                                                        !
+       print*,'DEBUT PHYSIQ_MOD ISO TRAVAIL DEBUT 5'
        iflag_physiq=0                                                      !
        CALL getin_p('iflag_physiq', iflag_physiq)                          !
@@ -1838,7 +1839,17 @@
        CALL thermcell_ini(iflag_thermals,prt_level,tau_thermals,lunout, &
             &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)
+       IF (iflag_thermals>=1) THEN
+          CALL alpale_th_first()
+       ENDIF
+
        CALL lscp_ini(pdtphys,klon,klev,iflag_thermals,iflag_cld_th,iflag_ice_thermo, &
                   lunout,prt_level,ok_ice_supersat,fl_cor_ebil,                       &
                   RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RV,RG,RPI,EPS_W)
+       IF (ok_new_lscp) THEN
+           ! old (pre-cmip7) code
+           !$gpum nocall
+           CALL fisrtilp_first(klon, klev, pdtphys, pfrac_impa, pfrac_nucl, pfrac_1nucl)
+       ENDIF
+
        CALL blosno_ini(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG, RV, RPI)
        CALL gwd_ini(klon, klev, pplay, paprs, lunout, prt_level, RCPD, RD, RG, RPI, ROMEGA, RA, RKAPPA, RLVTT, ok_strato)
@@ -2028,4 +2039,6 @@
              CALL ini_wake(0.,0.,it_wape_prescr,wape_prescr,fip_prescr &
                   ,alp_bl_prescr, ale_bl_prescr)
+             CALL calWAKE_first(phys_tstep)
+
              ! 11/09/06 rajout initialisation ALE et ALP du wake et PBL(YU)
              !        print*,'apres ini_wake iflag_cld_th=', iflag_cld_th
@@ -2053,4 +2066,9 @@
              d_deltaxt_ajs_cv(:,:,:) = 0.
 #endif
+             IF (iflag_alp_wk_cond .GT. 0.) THEN
+                CALL alpale_wk_first(cell_area)
+             ENDIF
+
+
           ENDIF  !  (iflag_wake>=1)
 
@@ -4352,5 +4370,4 @@
         !
         !calcul caracteristiques de la poche froide
-        CALL calWAKE_first(phys_tstep)
         CALL calWAKE (iflag_wake_tend, paprs, pplay, phys_tstep, &
                t_seri, q_seri, omega,  &
@@ -4425,5 +4442,4 @@
       IF (iflag_alp_wk_cond .GT. 0.) THEN
 
-        CALL alpale_wk_first(cell_area)
         CALL alpale_wk(phys_tstep, cell_area, wake_k, wake_s, wake_dens, wake_fip_0, &
                         wake_fip)
@@ -4752,5 +4768,4 @@
           !
 !
-          CALL alpale_th_first()
           CALL alpale_th( phys_tstep, lmax_th, t_seri, cell_area,  &
                           cin, s2, n2, strig, &