--- trunk/Sources/dyn3d/laplacien_gam.f	2016/08/30 12:52:46	206
+++ trunk/Sources/dyn3d/laplacien_gam.f	2016/09/01 10:30:53	207
@@ -1,48 +1,47 @@
+module laplacien_gam_m
 
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_gam.F,v 1.1.1.1
-! 2004/05/19 12:53:06 lmdzadmin Exp $
-
-SUBROUTINE laplacien_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, &
-    unsapolsga, teta, divgra)
-
-  ! P. Le Van
-
-  ! ************************************************************
+  IMPLICIT NONE
 
-  ! ....   calcul de  (div( grad ))   de   teta  .....
-  ! ************************************************************
-  ! klevel et teta  sont des arguments  d'entree pour le s-prog
-  ! divgra     est  un argument  de sortie pour le s-prog
+contains
 
-  USE grad_m, ONLY: grad
-  USE dimens_m
-  USE paramet_m
-  USE comgeom
-  IMPLICIT NONE
+  SUBROUTINE laplacien_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, &
+       unsapolsga, teta, divgra)
 
+    ! From LMDZ4/libf/dyn3d/laplacien_gam.F,v 1.1.1.1, 2004/05/19 12:53:06
 
+    ! P. Le Van
 
-  ! ............     variables  en arguments    ..........
+    ! ************************************************************
 
-  INTEGER, INTENT (IN) :: klevel
-  REAL teta(ip1jmp1, klevel), divgra(ip1jmp1, klevel)
-  REAL cuvsga(ip1jm), cvusga(ip1jmp1), unsaigam(ip1jmp1), unsapolnga, &
-    unsapolsga
+    ! ....   calcul de  (div( grad ))   de   teta  .....
+    ! ************************************************************
+    ! klevel et teta  sont des arguments  d'entree pour le s-prog
+    ! divgra     est  un argument  de sortie pour le s-prog
 
-  ! ...........    variables  locales    .................
+    USE comgeom
+    USE dimens_m
+    use diverg_gam_m, only: diverg_gam
+    USE grad_m, ONLY: grad
+    USE paramet_m
 
-  REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
-  ! ......................................................
 
-  CALL scopy(ip1jmp1*klevel, teta, 1, divgra, 1)
+    ! ............     variables  en arguments    ..........
 
-  CALL grad(klevel, divgra, ghx, ghy)
+    INTEGER, INTENT (IN) :: klevel
+    REAL teta(ip1jmp1, klevel), divgra(ip1jmp1, klevel)
+    REAL cuvsga(ip1jm), cvusga(ip1jmp1), unsaigam(ip1jmp1), unsapolnga, &
+         unsapolsga
 
-  CALL diverg_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, unsapolsga, &
-    ghx, ghy, divgra)
+    ! ...........    variables  locales    .................
 
+    REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
+    ! ......................................................
 
+    CALL scopy(ip1jmp1*klevel, teta, 1, divgra, 1)
+    CALL grad(klevel, divgra, ghx, ghy)
+    CALL diverg_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, unsapolsga, &
+         ghx, ghy, divgra)
 
+  END SUBROUTINE laplacien_gam
 
-  RETURN
-END SUBROUTINE laplacien_gam
+end module laplacien_gam_m