--- trunk/libf/dyn3d/calfis.f90	2013/07/08 18:12:18	71
+++ trunk/dyn3d/calfis.f	2014/03/11 15:09:02	88
@@ -40,6 +40,7 @@
     use nr_util, only: pi
     use physiq_m, only: physiq
     use pressure_var, only: p3d, pls
+    use pvtheta_m, only: pvtheta
 
     ! Arguments :
 
@@ -101,8 +102,8 @@
     REAL z1(iim)
     REAL pksurcp(iim + 1, jjm + 1)
 
-    ! I. Musat: diagnostic PVteta, Amip2
-    INTEGER, PARAMETER:: ntetaSTD=3
+    ! Diagnostic PVteta pour Amip2 :
+    INTEGER, PARAMETER:: ntetaSTD = 3
     REAL:: rtetaSTD(ntetaSTD) = (/350., 380., 405./)
     REAL PVteta(klon, ntetaSTD)
 
@@ -229,13 +230,13 @@
 
     forall(l= 1: llm) v(:, l) = pack(zvfi(:, :, l), dyn_phy)
 
-    !IM calcul PV a teta=350, 380, 405K
+    ! Compute potential vorticity at theta = 350, 380 and 405 K:
     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: