--- trunk/Sources/phylmd/Thermcell/thermcell.f	2015/07/20 16:01:49	157
+++ trunk/Sources/phylmd/Thermcell/thermcell.f	2016/03/11 18:47:26	178
@@ -53,13 +53,10 @@
 
     real zmax(klon), zw, zw2(klon, klev+1), ztva(klon, klev)
 
-    real zlev(klon, klev+1), zlay(klon, klev)
+    real zlev(klon, klev+1)
     REAL zh(klon, klev), zdhadj(klon, klev)
     REAL ztv(klon, klev)
     real zu(klon, klev), zv(klon, klev), zo(klon, klev)
-    real zla(klon, klev+1)
-    real zwa(klon, klev+1)
-    real zld(klon, klev+1)
     real zva(klon, klev)
     real zua(klon, klev)
     real zoa(klon, klev)
@@ -72,18 +69,13 @@
     real thetath2(klon, klev), wth2(klon, klev)
     common/comtherm/thetath2, wth2
 
-    integer isplit, nsplit
+    integer nsplit
     parameter (nsplit=10)
-    data isplit/0/
-    save isplit
 
-    logical sorties
     real rho(klon, klev), rhobarz(klon, klev+1), masse(klon, klev)
     real zpspsk(klon, klev)
 
     real wmax(klon), wmaxa(klon)
-    real wa(klon, klev, klev+1)
-    real wd(klon, klev+1)
     real fracd(klon, klev+1)
     real xxx(klon, klev+1)
     real larg_cons(klon, klev+1)
@@ -104,7 +96,6 @@
 
     ! initialisation:
 
-    sorties=.true.
     IF(ngrid.NE.klon) THEN
        PRINT *
        PRINT *, 'STOP dans convadj'
@@ -141,11 +132,6 @@
        zlev(ig, 1)=0.
        zlev(ig, nlay+1)=(2.*pphi(ig, klev)-pphi(ig, klev-1))/RG
     enddo
-    do l=1, nlay
-       do ig=1, ngrid
-          zlay(ig, l)=pphi(ig, l)/RG
-       enddo
-    enddo
 
     ! Calcul des densites
 
@@ -161,14 +147,6 @@
        enddo
     enddo
 
-    do k=1, nlay
-       do l=1, nlay+1
-          do ig=1, ngrid
-             wa(ig, k, l)=0.
-          enddo
-       enddo
-    enddo
-
     ! Calcul de w2, quarre de w a partir de la cape
     ! a partir de w2, on calcule wa, vitesse de l'ascendance
 
@@ -545,9 +523,6 @@
        do ig=1, ngrid
           if(fracd(ig, l).lt.0.1) then
              stop'fracd trop petit'
-          else
-             ! vitesse descendante "diagnostique"
-             wd(ig, l)=fm(ig, l)/(fracd(ig, l)*rhobarz(ig, l))
           endif
        enddo
     enddo
@@ -619,25 +594,6 @@
        enddo
     enddo
 
-    print *, '14 OK convect8'
-
-    ! Calculs pour les sorties
-
-    if(sorties) then
-       do l=1, nlay
-          do ig=1, ngrid
-             zla(ig, l)=(1.-fracd(ig, l))*zmax(ig)
-             zld(ig, l)=fracd(ig, l)*zmax(ig)
-             if(1.-fracd(ig, l).gt.1.e-10) &
-                  zwa(ig, l)=wd(ig, l)*fracd(ig, l)/(1.-fracd(ig, l))
-          enddo
-       enddo
-
-       isplit=isplit+1
-    endif
-
-    print *, '19 OK convect8'
-
   end SUBROUTINE thermcell
 
 end module thermcell_m