--- trunk/Sources/dyn3d/caldyn0.f	2015/07/24 14:27:59	161
+++ trunk/Sources/dyn3d/caldyn0.f	2018/01/22 15:02:56	252
@@ -4,19 +4,19 @@
 
 contains
 
-  SUBROUTINE caldyn0(ucov, vcov, teta, ps, masse, pk, phis, phi, w, &
-       pbaru, pbarv)
+  SUBROUTINE caldyn0(ucov, vcov, teta, ps, pk, phis, phi)
 
     ! From dyn3d/caldyn0.F, version 1.1.1.1, 2004/05/19 12:53:07
     ! Authors:  P. Le Van, F. Forget
     ! Objet : calcul des tendances dynamiques
 
     use bernoui_m, only: bernoui
-    USE comgeom, ONLY: airesurg
+    USE comgeom, ONLY: airesurg_2d
     use convmas_m, only: convmas
     use covcont_m, only: covcont
     USE dimens_m, ONLY: iim, jjm, llm
     USE disvert_m, ONLY: ap, bp
+    use enercin_m, only: enercin
     use flumass_m, only: flumass
     use massbar_m, only: massbar
     use massbarxy_m, only: massbarxy
@@ -30,20 +30,20 @@
     REAL, INTENT(IN):: vcov(:, :, :) ! (iim + 1, jjm, llm) ! vent covariant
     REAL, INTENT(IN):: teta(ip1jmp1, llm)
     REAL, INTENT (IN):: ps(ip1jmp1)
-    real, intent(out):: masse(ip1jmp1, llm)
     REAL, INTENT (IN):: pk(iip1, jjp1, llm)
     REAL, INTENT (IN):: phis(ip1jmp1)
     REAL, INTENT (IN):: phi(iim + 1, jjm + 1, llm)
-    REAL, INTENT(out):: w(:, :, :) ! (iim + 1, jjm + 1, llm)
-    REAL, intent(out):: pbaru(ip1jmp1, llm), pbarv((iim + 1) * jjm, llm)
 
     ! Local:
+    real masse(ip1jmp1, llm)
+    REAL w(iim + 1, jjm + 1, llm)
+    REAL pbaru(ip1jmp1, llm), pbarv((iim + 1) * jjm, llm)
     REAL vcont((iim + 1) * jjm, llm), ucont(ip1jmp1, llm)
     REAL p(ip1jmp1, llmp1)
     REAL massebx(ip1jmp1, llm), masseby((iim + 1) * jjm, llm)
     REAL vorpot(iim + 1, jjm, llm)
-    real ecin(iim + 1, jjm + 1, llm), convm(ip1jmp1, llm)
-    REAL massebxy(iim + 1, jjm, llm), dp(ip1jmp1)
+    real ecin(iim + 1, jjm + 1, llm), convm(iim + 1, jjm + 1, llm)
+    REAL massebxy(iim + 1, jjm, llm), dp(iim + 1, jjm + 1)
     INTEGER l
 
     !-----------------------------------------------------------------------
@@ -57,8 +57,8 @@
     CALL massbarxy(masse, massebxy)
     CALL flumass(massebx, masseby, vcont, ucont, pbaru, pbarv)
     CALL convmas(pbaru, pbarv, convm)
-    dp = convm(:, 1) / airesurg
-    CALL vitvert(convm, w)
+    dp = convm(:, :, 1) / airesurg_2d
+    w = vitvert(convm)
     CALL tourpot(vcov, ucov, massebxy, vorpot)
     CALL enercin(vcov, ucov, vcont, ucont, ecin)
     CALL sortvarc(ucov, teta, ps, masse, pk, phis, vorpot, phi, &