--- trunk/libf/dyn3d/divgrad2.f 2010/03/09 15:27:15 26 +++ trunk/libf/dyn3d/Dissipation/divgrad2.f 2011/12/12 13:25:01 55 @@ -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 @@ -15,15 +15,15 @@ use dimens_m 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, intent(in):: h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel) REAL divgra( ip1jmp1,klevel) + real, intent(in):: cdivh c c ....... variables locales .......... c @@ -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 )