/[lmdze]/trunk/libf/dyn3d/Dissipation/divgrad.f
ViewVC logotype

Diff of /trunk/libf/dyn3d/Dissipation/divgrad.f

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

trunk/libf/dyn3d/divgrad.f revision 26 by guez, Tue Mar 9 15:27:15 2010 UTC trunk/libf/dyn3d/Dissipation/divgrad.f revision 64 by guez, Wed Aug 29 14:47:17 2012 UTC
# Line 1  Line 1 
1  !  !
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3  !  !
4        SUBROUTINE divgrad (klevel,h, lh, divgra )        SUBROUTINE divgrad (klevel,h, lh, divgra, cdivh )
5          use grad_m, only: grad
6        use dimens_m        use dimens_m
7        use paramet_m        use paramet_m
8        use logic        use conf_gcm_m
9        use comgeom        use comgeom
10              use inidissip_m        use filtreg_m, only: filtreg
11        IMPLICIT NONE        IMPLICIT NONE
12  c  c
13  c=======================================================================  c=======================================================================
# Line 26  c   ------------- Line 27  c   -------------
27  c  c
28  c  c
29        INTEGER klevel        INTEGER klevel
30        REAL h( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )        REAL, intent(in):: h( ip1jmp1,klevel )
31          real, intent(out):: divgra( ip1jmp1,klevel )
32          real, intent(in):: cdivh
33  c  c
34        REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)        REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
35    
# Line 35  c Line 38  c
38  c  c
39  c  c
40  c  c
41        CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )        divgra = h
42  c  c
43        DO 10 iter = 1,lh        DO 10 iter = 1,lh
44    
45        CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1  )        CALL filtreg ( divgra,jjp1,klevel,2,1,.true.)
46    
47        CALL    grad (klevel,divgra, ghx  , ghy          )        CALL    grad (klevel,divgra, ghx  , ghy          )
48        CALL  diverg (klevel,  ghx , ghy  , divgra       )        CALL  diverg (klevel,  ghx , ghy  , divgra       )
49    
50        CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1)        CALL filtreg ( divgra,jjp1,klevel,2,1,.true.)
51    
52        DO 5 l = 1,klevel        DO 5 l = 1,klevel
53        DO 4  ij = 1, ip1jmp1        DO 4  ij = 1, ip1jmp1
Legend:
Removed from v.26  
changed lines
  Added in v.64
  ViewVC Help
Powered by ViewVC 1.1.21