--- trunk/libf/dyn3d/divgrad2.f	2011/12/06 15:07:04	54
+++ trunk/libf/dyn3d/Dissipation/divgrad2.f	2012/01/10 19:02:02	56
@@ -21,8 +21,8 @@
 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    ..........
@@ -35,8 +35,7 @@
 c
       signe    = (-1.)**lh
       nudivgrs = signe * cdivh
-
-      CALL SCOPY ( ip1jmp1 * klevel, h, 1, divgra, 1 )
+      divgra = h
 
 c
       CALL laplacien( klevel, divgra, divgra )