--- trunk/libf/dyn3d/integrd.f90	2011/07/01 15:00:48	47
+++ trunk/libf/dyn3d/integrd.f90	2012/10/02 15:50:56	67
@@ -11,10 +11,11 @@
     ! Author: P. Le Van 
     ! Objet: incrémentation des tendances dynamiques
 
-    USE comvert, ONLY : ap, bp
     USE comgeom, ONLY : aire, apoln, apols
     USE dimens_m, ONLY : iim, jjm, llm
+    USE disvert_m, ONLY : ap, bp
     USE filtreg_m, ONLY : filtreg
+    use massdair_m, only: massdair
     use nr_util, only: assert
     USE paramet_m, ONLY : iip1, iip2, ip1jm, ip1jmp1, jjp1, llmp1
 
@@ -27,7 +28,8 @@
     REAL masse((iim + 1) * (jjm + 1), llm)
 
     REAL vcovm1(ip1jm, llm), ucovm1((iim + 1) * (jjm + 1), llm)
-    REAL tetam1((iim + 1) * (jjm + 1), llm), psm1((iim + 1) * (jjm + 1))
+    REAL, intent(inout):: tetam1((iim + 1) * (jjm + 1), llm)
+    REAL, intent(inout):: psm1((iim + 1) * (jjm + 1))
     real massem1((iim + 1) * (jjm + 1), llm)
 
     REAL dv(ip1jm, llm), dudyn((iim + 1) * (jjm + 1), llm)
@@ -99,7 +101,7 @@
     CALL massdair(p, masse)
 
     finvmasse = masse
-    CALL filtreg(finvmasse, jjp1, llm, -2, 2, .TRUE., 1)
+    CALL filtreg(finvmasse, jjp1, llm, -2, 2, .TRUE.)
 
     ! integration de ucov, vcov, h
 
@@ -114,11 +116,9 @@
           vcov(ij, l) = vcovm1(ij, l) + dt*dv(ij, l)
        END DO
 
-       DO ij = 1, (iim + 1) * (jjm + 1)
-          hscr(ij) = teta(ij, l)
-          teta(ij, l) = tetam1(ij, l) * massem1(ij, l) / masse(ij, l) &
-               + dt * dteta(ij, l) / masse(ij, l)
-       END DO
+       hscr = teta(:, l)
+       teta(:, l) = tetam1(:, l) * massem1(:, l) / masse(:, l) &
+            + dt * dteta(:, l) / masse(:, l)
 
        ! Calcul de la valeur moyenne, unique aux poles pour teta
 
@@ -135,9 +135,9 @@
        END DO
 
        IF (leapf) THEN
-          CALL scopy((iim + 1) * (jjm + 1), uscr(1), 1, ucovm1(1, l), 1)
-          CALL scopy(ip1jm, vscr(1), 1, vcovm1(1, l), 1)
-          CALL scopy((iim + 1) * (jjm + 1), hscr(1), 1, tetam1(1, l), 1)
+          ucovm1(:, l)  =uscr
+          vcovm1(:, l) = vscr
+          tetam1(:, l) = hscr
        END IF
     END DO
 
@@ -167,13 +167,13 @@
        END DO
     END DO
 
-    CALL scopy((iim + 1) * (jjm + 1) * llm, finvmasse, 1, finvmaold, 1)
+    finvmaold = finvmasse
 
     ! Fin de l'integration de q
 
     IF (leapf) THEN
-       CALL scopy((iim + 1) * (jjm + 1), pscr, 1, psm1, 1)
-       CALL scopy((iim + 1) * (jjm + 1)*llm, massescr, 1, massem1, 1)
+       psm1 = pscr
+       massem1 = massescr
     END IF
 
   END SUBROUTINE integrd