--- trunk/Sources/dyn3d/diverg_gam.f	2016/08/30 12:52:46	206
+++ trunk/Sources/dyn3d/diverg_gam.f	2016/09/01 10:30:53	207
@@ -1,80 +1,84 @@
+module diverg_gam_m
 
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/diverg_gam.F,v 1.1.1.1 2004/05/19
-! 12:53:05 lmdzadmin Exp $
+  IMPLICIT NONE
 
-SUBROUTINE diverg_gam(klevel, cuvscvgam, cvuscugam, unsairegam, unsapolnga, &
-    unsapolsga, x, y, div)
+contains
 
-  ! P. Le Van
+  SUBROUTINE diverg_gam(klevel, cuvscvgam, cvuscugam, unsairegam, unsapolnga, &
+       unsapolsga, x, y, div)
 
-  ! *********************************************************************
-  ! ... calcule la divergence a tous les niveaux d'1 vecteur de compos.
-  ! x et y...
-  ! x et y  etant des composantes covariantes   ...
-  ! *********************************************************************
-  USE dimens_m
-  USE paramet_m
-  USE comgeom
-  IMPLICIT NONE
+    ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/diverg_gam.F,v 1.1.1.1 2004/05/19
+    ! 12:53:05 lmdzadmin Exp $
 
-  ! x  et  y  sont des arguments  d'entree pour le s-prog
-  ! div      est  un argument  de sortie pour le s-prog
+    ! P. Le Van
 
+    ! *********************************************************************
+    ! ... calcule la divergence a tous les niveaux d'1 vecteur de compos.
+    ! x et y...
+    ! x et y  etant des composantes covariantes   ...
+    ! *********************************************************************
+    USE dimens_m
+    USE paramet_m
+    USE comgeom
 
-  ! ---------------------------------------------------------------------
+    ! x  et  y  sont des arguments  d'entree pour le s-prog
+    ! div      est  un argument  de sortie pour le s-prog
 
-  ! ATTENTION : pendant ce s-pg , ne pas toucher au COMMON/scratch/  .
 
-  ! ---------------------------------------------------------------------
+    ! ---------------------------------------------------------------------
 
-  ! ..........          variables en arguments    ...................
+    ! ATTENTION : pendant ce s-pg , ne pas toucher au COMMON/scratch/  .
 
-  INTEGER, INTENT (IN) :: klevel
-  REAL x(ip1jmp1, klevel), y(ip1jm, klevel), div(ip1jmp1, klevel)
-  REAL cuvscvgam(ip1jm), cvuscugam(ip1jmp1), unsairegam(ip1jmp1)
-  REAL unsapolnga, unsapolsga
+    ! ---------------------------------------------------------------------
 
-  ! ...............     variables  locales   .........................
+    ! ..........          variables en arguments    ...................
 
-  REAL aiy1(iip1), aiy2(iip1)
-  REAL sumypn, sumyps
-  INTEGER l, ij
-  ! ...................................................................
+    INTEGER, INTENT (IN) :: klevel
+    REAL x(ip1jmp1, klevel), y(ip1jm, klevel), div(ip1jmp1, klevel)
+    REAL cuvscvgam(ip1jm), cvuscugam(ip1jmp1), unsairegam(ip1jmp1)
+    REAL unsapolnga, unsapolsga
 
-  REAL ssum
+    ! ...............     variables  locales   .........................
 
+    REAL aiy1(iip1), aiy2(iip1)
+    REAL sumypn, sumyps
+    INTEGER l, ij
+    ! ...................................................................
 
-  DO l = 1, klevel
+    REAL ssum
 
-    DO ij = iip2, ip1jm - 1
-      div(ij+1, l) = (cvuscugam(ij+1)*x(ij+1,l)-cvuscugam(ij)*x(ij,l)+ &
-        cuvscvgam(ij-iim)*y(ij-iim,l)-cuvscvgam(ij+1)*y(ij+1,l))* &
-        unsairegam(ij+1)
-    END DO
 
-    ! ....  correction pour  div( 1,j,l)  ......
-    ! ....   div(1,j,l)= div(iip1,j,l) ....
+    DO l = 1, klevel
 
-    ! DIR$ IVDEP
-    DO ij = iip2, ip1jm, iip1
-      div(ij, l) = div(ij+iim, l)
-    END DO
+       DO ij = iip2, ip1jm - 1
+          div(ij+1, l) = (cvuscugam(ij+1)*x(ij+1,l)-cvuscugam(ij)*x(ij,l)+ &
+               cuvscvgam(ij-iim)*y(ij-iim,l)-cuvscvgam(ij+1)*y(ij+1,l))* &
+               unsairegam(ij+1)
+       END DO
 
-    ! ....  calcul  aux poles  .....
+       ! ....  correction pour  div( 1,j,l)  ......
+       ! ....   div(1,j,l)= div(iip1,j,l) ....
 
-    DO ij = 1, iim
-      aiy1(ij) = cuvscvgam(ij)*y(ij, l)
-      aiy2(ij) = cuvscvgam(ij+ip1jmi1)*y(ij+ip1jmi1, l)
-    END DO
-    sumypn = ssum(iim, aiy1, 1)*unsapolnga
-    sumyps = ssum(iim, aiy2, 1)*unsapolsga
+       ! DIR$ IVDEP
+       DO ij = iip2, ip1jm, iip1
+          div(ij, l) = div(ij+iim, l)
+       END DO
+
+       ! ....  calcul  aux poles  .....
+
+       DO ij = 1, iim
+          aiy1(ij) = cuvscvgam(ij)*y(ij, l)
+          aiy2(ij) = cuvscvgam(ij+ip1jmi1)*y(ij+ip1jmi1, l)
+       END DO
+       sumypn = ssum(iim, aiy1, 1)*unsapolnga
+       sumyps = ssum(iim, aiy2, 1)*unsapolsga
 
-    DO ij = 1, iip1
-      div(ij, l) = -sumypn
-      div(ij+ip1jm, l) = sumyps
+       DO ij = 1, iip1
+          div(ij, l) = -sumypn
+          div(ij+ip1jm, l) = sumyps
+       END DO
     END DO
-  END DO
 
+  END SUBROUTINE diverg_gam
 
-  RETURN
-END SUBROUTINE diverg_gam
+end module diverg_gam_m