/[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 26 by guez, Tue Mar 9 15:27:15 2010 UTC trunk/libf/dyn3d/Dissipation/divgrad2.f revision 61 by guez, Fri Apr 20 14:58:43 2012 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 13  c   h ,klevel,lh et pext  sont des argum Line 13  c   h ,klevel,lh et pext  sont des argum
13  c         divgra     est  un argument  de sortie pour le s-prg  c         divgra     est  un argument  de sortie pour le s-prg
14  c  c
15        use dimens_m        use dimens_m
16          use laplacien_m, only: laplacien
17        use paramet_m        use paramet_m
18        use comgeom        use comgeom
             use inidissip_m  
19        IMPLICIT NONE        IMPLICIT NONE
20  c  c
21    
22  c    .......    variables en arguments   .......  c    .......    variables en arguments   .......
23  c  c
24        INTEGER klevel        INTEGER klevel
25        REAL h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel )        REAL, intent(in):: h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel)
26        REAL divgra( ip1jmp1,klevel)        REAL, intent(out):: divgra( ip1jmp1,klevel)
27          real, intent(in):: cdivh
28  c  c
29  c    .......    variables  locales    ..........  c    .......    variables  locales    ..........
30  c  c
# Line 35  c    ................................... Line 36  c    ...................................
36  c  c
37        signe    = (-1.)**lh        signe    = (-1.)**lh
38        nudivgrs = signe * cdivh        nudivgrs = signe * cdivh
39          divgra = h
       CALL SCOPY ( ip1jmp1 * klevel, h, 1, divgra, 1 )  
40    
41  c  c
42        CALL laplacien( klevel, divgra, divgra )        CALL laplacien( klevel, divgra )
43            
44        DO l = 1, klevel        DO l = 1, klevel
45         DO ij = 1, ip1jmp1         DO ij = 1, ip1jmp1
# Line 68  c    ................................... Line 68  c    ...................................
68          ENDDO          ENDDO
69        ENDDO        ENDDO
70  c  c
71        CALL laplacien ( klevel, divgra, divgra )        CALL laplacien ( klevel, divgra )
72  c  c
73        DO l  = 1,klevel        DO l  = 1,klevel
74        DO ij = 1,ip1jmp1        DO ij = 1,ip1jmp1
Legend:
Removed from v.26  
changed lines
  Added in v.61
  ViewVC Help
Powered by ViewVC 1.1.21