--- trunk/libf/dyn3d/divgrad.f	2010/03/05 16:43:45	25
+++ trunk/libf/dyn3d/Dissipation/divgrad.f	2012/01/10 19:02:02	56
@@ -1,12 +1,12 @@
 !
 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
 !
-      SUBROUTINE divgrad (klevel,h, lh, divgra )
+      SUBROUTINE divgrad (klevel,h, lh, divgra, cdivh )
       use dimens_m
       use paramet_m
       use logic
       use comgeom
-            use comdissipn
+      use filtreg_m, only: filtreg
       IMPLICIT NONE
 c
 c=======================================================================
@@ -26,15 +26,18 @@
 c
 c
       INTEGER klevel
-      REAL h( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
+      REAL, intent(in):: h( ip1jmp1,klevel )
+      real, intent(out):: divgra( ip1jmp1,klevel )
+      real, intent(in):: cdivh
 c
       REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
 
-      INTEGER  l,ij,iter,lh
+      INTEGER  l,ij,iter
+      integer, intent(in):: lh
 c
 c
 c
-      CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )
+      divgra = h
 c
       DO 10 iter = 1,lh