--- trunk/libf/dyn3d/divgrad2.f	2010/03/09 15:27:15	26
+++ trunk/libf/dyn3d/Dissipation/divgrad2.f	2012/04/20 14:58:43	61
@@ -1,7 +1,7 @@
 !
 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad2.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
 !
-      SUBROUTINE divgrad2 ( klevel, h, deltapres, lh, divgra )
+      SUBROUTINE divgrad2 ( klevel, h, deltapres, lh, divgra, cdivh )
 c
 c     P. Le Van
 c
@@ -13,17 +13,18 @@
 c         divgra     est  un argument  de sortie pour le s-prg
 c
       use dimens_m
+      use laplacien_m, only: laplacien
       use paramet_m
       use comgeom
-            use inidissip_m
       IMPLICIT NONE
 c
 
 c    .......    variables en arguments   .......
 c
       INTEGER klevel
-      REAL h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel )
-      REAL divgra( ip1jmp1,klevel)
+      REAL, intent(in):: h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel)
+      REAL, intent(out):: divgra( ip1jmp1,klevel)
+      real, intent(in):: cdivh
 c
 c    .......    variables  locales    ..........
 c
@@ -35,11 +36,10 @@
 c
       signe    = (-1.)**lh
       nudivgrs = signe * cdivh
-
-      CALL SCOPY ( ip1jmp1 * klevel, h, 1, divgra, 1 )
+      divgra = h
 
 c
-      CALL laplacien( klevel, divgra, divgra )
+      CALL laplacien( klevel, divgra )
      
       DO l = 1, klevel
        DO ij = 1, ip1jmp1
@@ -68,7 +68,7 @@
         ENDDO
       ENDDO
 c
-      CALL laplacien ( klevel, divgra, divgra )
+      CALL laplacien ( klevel, divgra )
 c
       DO l  = 1,klevel
       DO ij = 1,ip1jmp1