--- trunk/libf/dyn3d/Dissipation/gradiv2.f90	2012/01/30 14:37:26	60
+++ trunk/dyn3d/Dissipation/gradiv2.f	2014/03/05 14:57:53	82
@@ -10,12 +10,13 @@
     ! P. Le Van
     ! Calcul du gradient de la divergence du vecteur v.
 
-    USE dimens_m, ONLY : iim, jjm, llm
-    use divergf_m, only: divergf
-    USE comgeom, ONLY : cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, &
+    USE comgeom, ONLY: cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, &
          unsapolsga1
-    USE filtreg_m, ONLY : filtreg
+    USE dimens_m, ONLY: iim, jjm, llm
+    use divergf_m, only: divergf
+    USE filtreg_m, ONLY: filtreg
     use grad_m, only: grad
+    use laplacien_m, only: laplacien
     use nr_util, only: assert_eq, assert
 
     ! Composantes covariantes de v :
@@ -43,7 +44,7 @@
     CALL divergf(klevel, xcov, ycov, div)
 
     IF (ld > 1) THEN
-       CALL laplacien(klevel, div, div)
+       CALL laplacien(klevel, div)
 
        ! Itération de l'opérateur laplacien_gam
        DO iter = 1, ld -2
@@ -52,7 +53,7 @@
        END DO
     ENDIF
 
-    CALL filtreg(div, jjm + 1, klevel, 2, 1, .TRUE., 1)
+    CALL filtreg(div, jjm + 1, klevel, 2, 1, .TRUE.)
     CALL grad(klevel, div, gdx, gdy)
     nugrads = (-1.)**ld * cdivu