--- trunk/dyn3d/calfis.f90	2013/11/15 18:45:49	76
+++ trunk/dyn3d/calfis.f	2014/03/11 16:03:19	89
@@ -5,7 +5,7 @@
 contains
 
   SUBROUTINE calfis(rdayvrai, time, ucov, vcov, teta, q, ps, pk, phis, phi, &
-       dudyn, dv, w, dufi, dvfi, dtetafi, dqfi, dpfi, lafin)
+       dudyn, w, dufi, dvfi, dtetafi, dqfi, dpfi, lafin)
 
     ! From dyn3d/calfis.F, version 1.3 2005/05/25 13:10:09
     ! Authors: P. Le Van, F. Hourdin
@@ -69,7 +69,6 @@
     REAL, intent(in):: phis(iim + 1, jjm + 1)
     REAL, intent(in):: phi(iim + 1, jjm + 1, llm)
     REAL dudyn(iim + 1, jjm + 1, llm)
-    REAL dv(iim + 1, jjm, llm)
     REAL, intent(in):: w(iim + 1, jjm + 1, llm)
 
     REAL, intent(out):: dufi(iim + 1, jjm + 1, llm)
@@ -83,7 +82,7 @@
 
     ! Local variables :
 
-    INTEGER i, j, l, ig0, ig, iq, iiq
+    INTEGER i, j, l, ig0, iq, iiq
     REAL zpsrf(klon)
     REAL paprs(klon, llm+1), play(klon, llm)
     REAL pphi(klon, llm), pphis(klon)
@@ -101,11 +100,6 @@
     REAL z1(iim)
     REAL pksurcp(iim + 1, jjm + 1)
 
-    ! I. Musat: diagnostic PVteta, Amip2
-    INTEGER, PARAMETER:: ntetaSTD=3
-    REAL:: rtetaSTD(ntetaSTD) = (/350., 380., 405./)
-    REAL PVteta(klon, ntetaSTD)
-
     !-----------------------------------------------------------------------
 
     !!print *, "Call sequence information: calfis"
@@ -229,13 +223,9 @@
 
     forall(l= 1: llm) v(:, l) = pack(zvfi(:, :, l), dyn_phy)
 
-    !IM calcul PV a teta=350, 380, 405K
-    CALL PVtheta(klon, llm, ucov, vcov, teta, t, play, paprs, ntetaSTD, &
-         rtetaSTD, PVteta)
-
     ! Appel de la physique :
     CALL physiq(lafin, rdayvrai, time, dtphys, paprs, play, pphi, pphis, u, &
-         v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps, dudyn, PVteta)
+         v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps, dudyn)
 
     ! transformation des tendances physiques en tendances dynamiques: