Index: /LMDZ6/trunk/libf/phylmd/lmdz_wake.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/lmdz_wake.F90	(revision 4813)
+++ /LMDZ6/trunk/libf/phylmd/lmdz_wake.F90	(revision 4814)
@@ -708,4 +708,5 @@
       ! IM v3JYG; ptop_provis(i).LT. ph(i,1)
 
+! LJYF : mettre cette ligne à la place de la suivante      IF (dth(i,k)>=-delta_t_min .AND. dth(i,k-1)<-delta_t_min .AND. &
       IF (dth(i,k)>-delta_t_min .AND. dth(i,k-1)<-delta_t_min .AND. &
           ptop_provis(i)==ph(i,1)) THEN
@@ -766,6 +767,4 @@
   ! -5/ Determination de ktop et kupper
 
-  CALL pkupper (klon, klev, ptop, ph, pupper, kupper)
-  
   DO k = klev, 1, -1
     DO i = 1, klon
@@ -785,4 +784,5 @@
     DO i = 1, klon
       IF (k<=ktop(i) .AND. ptop_new(i)==ptop(i) .AND. &
+! LJYF changer par cette ligne          dth(i,k)>=-delta_t_min .AND. dth(i,k-1)<-delta_t_min) THEN
           dth(i,k)>-delta_t_min .AND. dth(i,k-1)<-delta_t_min) THEN
         ptop_new(i) = ((dth(i,k)+delta_t_min)*p(i,k-1)-(dth(i, &
@@ -802,4 +802,6 @@
   END DO
 
+  CALL pkupper (klon, klev, ptop, ph, pupper, kupper)
+  
   IF (prt_level>=10) THEN
     PRINT *, 'wake-3, ktop(igout), kupper(igout) ', ktop(igout), kupper(igout)
@@ -989,9 +991,4 @@
   !
   ! ------------------------------------------------------------------------
-  !
-    CALL pkupper (klon, klev, ptop, ph, pupper, kupper)
- 
-    !print*, 'ptop, pupper, ktop, kupper', ptop, pupper, ktop, kupper
-
     ! wk_adv is the logical flag enabling wake evolution in the time advance
     ! loop
@@ -1843,4 +1840,5 @@
       DO i = 1, klon
         IF (wk_adv(i) .AND. ptop_provis(i)==ph(i,1) .AND. &
+! LJYF changer :           dth(i,k)>=-delta_t_min .AND. dth(i,k-1)<-delta_t_min) THEN
             dth(i,k)>-delta_t_min .AND. dth(i,k-1)<-delta_t_min) THEN
           ptop_provis(i) = ((dth(i,k)+delta_t_min)*p(i,k-1) - &
@@ -1916,4 +1914,5 @@
         ! IM v3JYG; IF (k .GE. ktop(i)
         IF (wk_adv(i) .AND. k<=ktop(i) .AND. ptop_new(i)==ptop(i) .AND. &
+! LJYF changer :           dth(i,k)>=-delta_t_min .AND. dth(i,k-1)<-delta_t_min) THEN
             dth(i,k)>-delta_t_min .AND. dth(i,k-1)<-delta_t_min) THEN
           ptop_new(i) = ((dth(i,k)+delta_t_min)*p(i,k-1) - &
@@ -1937,4 +1936,7 @@
       END DO
     END DO
+
+    CALL pkupper (klon, klev, ptop, ph, pupper, kupper)
+ 
 
     ! 5/ Set deltatw & deltaqw to 0 above kupper
@@ -2087,5 +2089,6 @@
   ! ------------------------------------------------------------------------
   !
-  END DO                  ! end sub-timestep loop
+  END DO   ! isubstep end sub-timestep loop
+
   !
   ! ------------------------------------------------------------------------