--- trunk/libf/dyn3d/gradiv2.f	2010/03/05 16:43:45	25
+++ trunk/libf/dyn3d/Dissipation/gradiv2.f90	2011/12/12 13:25:01	55
@@ -1,79 +1,65 @@
-!
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/gradiv2.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
-!
-      SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy )
-c
-c     P. Le Van
-c
-c   **********************************************************
-c                                ld
-c       calcul  de  (grad (div) )   du vect. v ....
-c
-c     xcov et ycov etant les composant.covariantes de v
-c   **********************************************************
-c     xcont , ycont et ld  sont des arguments  d'entree pour le s-prog
-c      gdx   et  gdy       sont des arguments de sortie pour le s-prog
-c
-c
-      use dimens_m
-      use paramet_m
-      use comgeom
-            use comdissipn
-      IMPLICIT NONE
-c
-c
-c     ........    variables en arguments      ........
-
-      INTEGER klevel
-      REAL  xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel )
-      REAL   gdx( ip1jmp1,klevel ),  gdy( ip1jm,klevel )
-c
-c     ........       variables locales       .........
-c
-      REAL div(ip1jmp1,llm)
-      REAL signe, nugrads
-      INTEGER l,ij,iter,ld
-      
-c    ........................................................
-c
-c
-      CALL SCOPY( ip1jmp1 * klevel, xcov, 1, gdx, 1 )
-      CALL SCOPY(   ip1jm * klevel, ycov, 1, gdy, 1 )
-c
-c
-      signe   = (-1.)**ld
-      nugrads = signe * cdivu
-c
-
-
-      CALL    divergf( klevel, gdx,   gdy , div )
-
-      IF( ld.GT.1 )   THEN
-
-        CALL laplacien ( klevel, div,  div     )
-
-c    ......  Iteration de l'operateur laplacien_gam   .......
-
-        DO iter = 1, ld -2
-         CALL laplacien_gam ( klevel,cuvscvgam1,cvuscugam1,unsair_gam1,
-     *                       unsapolnga1, unsapolsga1,  div, div       )
-        ENDDO
-
-      ENDIF
-
-
-       CALL filtreg( div   , jjp1, klevel, 2, 1, .TRUE., 1 )
-       CALL  grad  ( klevel,  div,   gdx,  gdy             )
-
-c
-       DO   l = 1, klevel
-         DO  ij = 1, ip1jmp1
-          gdx( ij,l ) = gdx( ij,l ) * nugrads
-         ENDDO
-         DO  ij = 1, ip1jm
-          gdy( ij,l ) = gdy( ij,l ) * nugrads
-         ENDDO
+module gradiv2_m
+
+  IMPLICIT NONE
+
+contains
+
+  SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy, cdivu)
+
+    ! From LMDZ4/libf/dyn3d/gradiv2.F, version 1.1.1.1 2004/05/19 12:53:07
+    ! P. Le Van
+    ! Calcul de grad div du vecteur v.
+
+    USE dimens_m, ONLY : llm
+    USE paramet_m, ONLY : ip1jm, ip1jmp1, jjp1
+    USE comgeom, ONLY : cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, &
+         unsapolsga1
+    USE filtreg_m, ONLY : filtreg
+
+    INTEGER, intent(in):: klevel
+
+    ! composantes covariantes de v:
+    REAL, intent(in):: xcov(ip1jmp1,klevel), ycov(ip1jm,klevel)
+
+    integer, intent(in):: ld
+    REAL, intent(out):: gdx(ip1jmp1,klevel), gdy(ip1jm,klevel)
+    real, intent(in):: cdivu
+
+    ! Variables locales :
+    REAL div(ip1jmp1,llm)
+    REAL nugrads
+    INTEGER l,ij,iter
+
+    !--------------------------------------------------------------
+
+    gdx = xcov
+    gdy = ycov
+
+    CALL divergf(klevel, gdx, gdy, div)
+
+    IF(ld.GT.1) THEN
+       CALL laplacien (klevel, div, div)
+
+       ! Iteration de l'operateur laplacien_gam
+       DO iter = 1, ld -2
+          CALL laplacien_gam (klevel,cuvscvgam1,cvuscugam1,unsair_gam1, &
+               unsapolnga1, unsapolsga1, div, div)
        ENDDO
-c
-       RETURN
-       END
+    ENDIF
+
+    CALL filtreg(div, jjp1, klevel, 2, 1, .TRUE., 1)
+    CALL grad (klevel, div, gdx, gdy)
+    nugrads = (-1.)**ld * cdivu
+
+    DO l = 1, klevel
+       DO ij = 1, ip1jmp1
+          gdx(ij,l) = gdx(ij,l) * nugrads
+       ENDDO
+       DO ij = 1, ip1jm
+          gdy(ij,l) = gdy(ij,l) * nugrads
+       ENDDO
+    ENDDO
+
+  END SUBROUTINE gradiv2
+
+end module gradiv2_m