Index: /LMDZ5/trunk/libf/phylmd/climb_hq_mod.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/climb_hq_mod.F90	(revision 2851)
+++ /LMDZ5/trunk/libf/phylmd/climb_hq_mod.F90	(revision 2852)
@@ -221,5 +221,6 @@
 !****************************************************************************************
 !!! jyg le 07/02/2012
-       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+!!jyg       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+       IF (mod(iflag_pbl_split,10) .ge.1) THEN
 !!! nrlmd le 02/05/2011
     DO k= 1, klev
@@ -240,5 +241,5 @@
     ENDDO
 !!!      
-       ENDIF  ! (mod(iflag_pbl_split,2) .eq.1)
+       ENDIF  ! (mod(iflag_pbl_split,2) .ge.1)
 !!!
 
@@ -377,5 +378,6 @@
 
 !!! jyg le 07/02/2012
-       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+!!jyg       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+       IF (mod(iflag_pbl_split,10) .ge.1) THEN
 !!! nrlmd le 02/05/2011
     DO i = 1, knon
@@ -399,5 +401,5 @@
     ENDDO
 !!!      
-       ENDIF  ! (mod(iflag_pbl_split,2) .eq.1)
+       ENDIF  ! (mod(iflag_pbl_split,2) .ge.1)
 !!!
 
Index: /LMDZ5/trunk/libf/phylmd/climb_wind_mod.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/climb_wind_mod.F90	(revision 2851)
+++ /LMDZ5/trunk/libf/phylmd/climb_wind_mod.F90	(revision 2852)
@@ -172,5 +172,6 @@
 !****************************************************************************************
 !!! jyg le 07/02/2012
-       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+!!jyg       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+       IF (mod(iflag_pbl_split,10) .ge.1) THEN
 !!! nrlmd le 02/05/2011
     DO k= 1, klev
@@ -188,5 +189,5 @@
     ENDDO
 !!!      
-       ENDIF  ! (mod(iflag_pbl_split,2) .eq.1)
+       ENDIF  ! (mod(iflag_pbl_split,2) .ge.1)
 !!!
 
@@ -308,5 +309,6 @@
 
 !!! jyg le 07/02/2012
-       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+!!jyg       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+       IF (mod(iflag_pbl_split,10) .ge.1) THEN
 !!! nrlmd le 02/05/2011
     DO i = 1, knon
@@ -326,5 +328,5 @@
     ENDDO
 !!!
-       ENDIF  ! (mod(iflag_pbl_split,2) .eq.1)
+       ENDIF  ! (mod(iflag_pbl_split,2) .ge.1)
 !!!
 
Index: /LMDZ5/trunk/libf/phylmd/conf_phys_m.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/conf_phys_m.F90	(revision 2851)
+++ /LMDZ5/trunk/libf/phylmd/conf_phys_m.F90	(revision 2852)
@@ -1369,10 +1369,10 @@
     !
     !Config Key  = iflag_pbl_split
-    !Config Desc = binary flag: least signif bit = split vdf; next bit = split thermals
+    !Config Desc = decimal flag: least signif digit = split vdf; next digit = split thermals
     !Config Def  = 0
-    !Config Help = 0-> no splitting; 1-> vdf splitting; 2-> thermals splitting; 3-> full splitting
+    !Config Help = 0-> no splitting; 1-> vdf splitting; 10-> thermals splitting; 11-> full splitting
     !
     iflag_pbl_split_omp = 0
-    CALL getin('iflag_pbl_split',iflag_pbl_split_omp)
+    call getin('iflag_pbl_split',iflag_pbl_split_omp)
     !
     !Config Key  = iflag_order2_sollw
Index: /LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 2851)
+++ /LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 2852)
@@ -6,5 +6,5 @@
 ! Planetary Boundary Layer and Surface module
 !
-! This module manage the calculation of turbulent diffusion in the boundary layer 
+! This module manages the calculation of turbulent diffusion in the boundary layer 
 ! and all interactions towards the differents sub-surfaces.
 !
@@ -778,5 +778,6 @@
       IF (prt_level >=10) print *,' -> pbl_surface, itap ',itap
 !
-      iflag_split = mod(iflag_pbl_split,2)
+!!jyg      iflag_split = mod(iflag_pbl_split,2)
+      iflag_split = mod(iflag_pbl_split,10)
 
 !****************************************************************************************
@@ -1239,5 +1240,5 @@
         ENDDO
 !
-       IF (iflag_split .eq.1) THEN
+       IF (iflag_split .ge.1) THEN
 !!! nrlmd le 02/05/2011
         DO k = 1, klev
@@ -1287,5 +1288,5 @@
         ENDDO
 !!!
-       ENDIF  ! (iflag_split .eq.1)
+       ENDIF  ! (iflag_split .ge.1)
 !!!
        DO k = 1, nsoilmx
@@ -2047,5 +2048,5 @@
          print*,'effusivity,facteur,cstar', effusivity, facteur,wake_cstar(j)
         ENDDO
-       ENDIF
+       ENDIF  ! (prt_level >=10)
 
 !!! jyg le 07/02/2012
@@ -2282,5 +2283,5 @@
 
 !!! jyg le 07/02/2012
-       IF (iflag_split .eq.1) THEN
+       IF (iflag_split .ge.1) THEN
 !!!
 !!! nrlmd le 02/05/2011
@@ -2302,5 +2303,5 @@
         END DO
 !!!      
-       ENDIF  ! (iflag_split .eq.1)
+       ENDIF  ! (iflag_split .ge.1)
 !!!
 !!! nrlmd le 02/05/2011
@@ -2337,5 +2338,5 @@
         END DO
 
-       ELSE
+       ELSE  ! (iflag_split .eq.0)
         DO k = 1, klev
           DO j = 1, knon
@@ -2384,5 +2385,5 @@
        
 !!! jyg le 07/02/2012
-       IF (iflag_split .eq.1) THEN
+       IF (iflag_split .ge.1) THEN
 !!!
 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
@@ -2407,5 +2408,5 @@
         END DO
 !!!
-       ENDIF  ! (iflag_split .eq.1)
+       ENDIF  ! (iflag_split .ge.1)
 !!!
        
@@ -2752,5 +2753,5 @@
 
 !!! jyg le 07/02/2012
-       IF (iflag_split .eq.1) THEN
+       IF (iflag_split .ge.1) THEN
 !!!
 !!! nrlmd & jyg les 02/05/2011, 05/02/2012
@@ -2777,5 +2778,5 @@
     END DO
 !!!
-       ENDIF  ! (iflag_split .eq.1)
+       ENDIF  ! (iflag_split .ge.1)
 !!!
 
Index: /LMDZ5/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2851)
+++ /LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2852)
@@ -2137,5 +2137,6 @@
 
        !jyg+nrlmd<
-       IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,2) .eq. 1) THEN
+!!jyg       IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,2) .eq. 1) THEN
+       IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,10) .ge. 1) THEN
           print *,'debut du splitting de la PBL'
        ENDIF
@@ -2217,5 +2218,6 @@
        !
        !  Add turbulent diffusion tendency to the wake difference variables
-       IF (mod(iflag_pbl_split,2) .NE. 0) THEN
+!!jyg       IF (mod(iflag_pbl_split,2) .NE. 0) THEN
+       IF (mod(iflag_pbl_split,10) .NE. 0) THEN
 !jyg<
           d_deltat_vdf(:,:) = d_t_vdf_w(:,:)-d_t_vdf_x(:,:)
@@ -2829,5 +2831,6 @@
        IF (iflag_thermals>=1) THEN
           !jyg<
-          IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
+!!       IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
+       IF (mod(iflag_pbl_split/10,10) .GE. 1) THEN
              !  Appel des thermiques avec les profils exterieurs aux poches
              DO k=1,klev
@@ -2874,5 +2877,6 @@
           !
           !jyg<
-          IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
+!!jyg          IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN
+          IF (mod(iflag_pbl_split/10,10) .GE. 1) THEN
              !  Si les thermiques ne sont presents que hors des
              !  poches, la tendance moyenne associ\'ee doit etre
@@ -2896,5 +2900,5 @@
              call prt_enerbil('the',itap)
           !
-          ENDIF  ! (mod(iflag_pbl_split/2,2) .EQ. 1)
+          ENDIF  ! (mod(iflag_pbl_split/10,10) .GE. 1)
           !
           CALL add_phys_tend(d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs,  &