Index: LMDZ6/trunk/libf/phylmd/coef_diff_turb_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/coef_diff_turb_mod.F90	(revision 3575)
+++ LMDZ6/trunk/libf/phylmd/coef_diff_turb_mod.F90	(revision 3577)
@@ -182,4 +182,7 @@
        ycoefh(1:knon,2:klev)=ykmn(1:knon,2:klev)
                 
+    ELSE
+       ! No TKE for Standard Physics
+       yq2=0.
     ENDIF !(iflag_pbl.ge.3)
 
Index: LMDZ6/trunk/libf/phylmd/cv30_routines.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cv30_routines.F90	(revision 3575)
+++ LMDZ6/trunk/libf/phylmd/cv30_routines.F90	(revision 3577)
@@ -3050,7 +3050,8 @@
   ! variables pour tracer dans precip de l'AA et des mel
   ! local variables:
-  INTEGER i, j, k
+  INTEGER i, j, k, nam1
   REAL epm(nloc, na, na)
 
+  nam1=na-1 ! Introduced because ep is not defined for j=na
   ! variables d'Emanuel : du second indice au troisieme
   ! --->    tab(i,k,j) -> de l origine k a l arrivee j
@@ -3082,5 +3083,5 @@
   ! fraction deau condensee dans les melanges convertie en precip : epm
   ! et eau condensée précipitée dans masse d'air saturé : l_m*dM_m/dzdz.dzdz
-  DO j = 1, na
+  DO j = 1, nam1
     DO k = 1, j - 1
       DO i = 1, ncum
@@ -3095,6 +3096,6 @@
   END DO
 
-  DO j = 1, na
-    DO k = 1, na
+  DO j = 1, nam1
+    DO k = 1, nam1
       DO i = 1, ncum
         IF (k>=icb(i) .AND. k<=inb(i)) THEN
@@ -3106,5 +3107,5 @@
   END DO
 
-  DO j = 1, na
+  DO j = 1, nam1
     DO k = 1, j - 1
       DO i = 1, ncum
@@ -3117,6 +3118,6 @@
 
   ! matrices pour calculer la tendance des concentrations dans cvltr.F90
-  DO j = 1, na
-    DO k = 1, na
+  DO j = 1, nam1
+    DO k = 1, nam1
       DO i = 1, ncum
         da(i, j) = da(i, j) + (1.-sij(i,k,j))*ment(i, k, j)
@@ -3127,5 +3128,5 @@
   END DO
 
-  DO j = 1, na
+  DO j = 1, nam1
     DO k = 1, j - 1
       DO i = 1, ncum
Index: LMDZ6/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 3575)
+++ LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 3577)
@@ -1701,5 +1701,9 @@
 
 !jyg<
-       IF (klon_glo==1) THEN
+       IF (iflag_pbl<=1) THEN
+          ! No TKE for Standard Physics
+          pbl_tke(:,:,:)=0.
+
+       ELSE IF (klon_glo==1) THEN
           pbl_tke(:,:,is_ave) = 0.
           DO nsrf=1,nbsrf