/[lmdze]/trunk/dyn3d/Dissipation/divgrad2.f
ViewVC logotype

Diff of /trunk/dyn3d/Dissipation/divgrad2.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/dyn3d/divgrad2.f revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC trunk/libf/dyn3d/Dissipation/divgrad2.f revision 55 by guez, Mon Dec 12 13:25:01 2011 UTC
# Line 1  Line 1 
1  !  !
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad2.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad2.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3  !  !
4        SUBROUTINE divgrad2 ( klevel, h, deltapres, lh, divgra )        SUBROUTINE divgrad2 ( klevel, h, deltapres, lh, divgra, cdivh )
5  c  c
6  c     P. Le Van  c     P. Le Van
7  c  c
# Line 17  c Line 17  c
17        use comgeom        use comgeom
18        IMPLICIT NONE        IMPLICIT NONE
19  c  c
       include "comdissipn.h"  
20    
21  c    .......    variables en arguments   .......  c    .......    variables en arguments   .......
22  c  c
23        INTEGER klevel        INTEGER klevel
24        REAL h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel )        REAL, intent(in):: h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel)
25        REAL divgra( ip1jmp1,klevel)        REAL divgra( ip1jmp1,klevel)
26          real, intent(in):: cdivh
27  c  c
28  c    .......    variables  locales    ..........  c    .......    variables  locales    ..........
29  c  c
30        REAL     signe, nudivgrs, sqrtps( ip1jmp1,llm )        REAL     signe, nudivgrs, sqrtps( ip1jmp1,llm )
31        INTEGER  l,ij,iter,lh        INTEGER  l,ij,iter
32          integer, intent(in):: lh
33  c    ...................................................................  c    ...................................................................
34    
35  c  c
36        signe    = (-1.)**lh        signe    = (-1.)**lh
37        nudivgrs = signe * cdivh        nudivgrs = signe * cdivh
38          divgra = h
       CALL SCOPY ( ip1jmp1 * klevel, h, 1, divgra, 1 )  
39    
40  c  c
41        CALL laplacien( klevel, divgra, divgra )        CALL laplacien( klevel, divgra, divgra )
Legend:
Removed from v.3  
changed lines
  Added in v.55
  ViewVC Help
Powered by ViewVC 1.1.21