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

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

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

revision 56 by guez, Tue Jan 10 19:02:02 2012 UTC revision 65 by guez, Thu Sep 20 09:57:03 2012 UTC
# Line 4  module gradiv2_m Line 4  module gradiv2_m
4    
5  contains  contains
6    
7    SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy, cdivu)    SUBROUTINE gradiv2(xcov, ycov, ld, gdx, gdy, cdivu)
8    
9      ! From LMDZ4/libf/dyn3d/gradiv2.F, version 1.1.1.1 2004/05/19 12:53:07      ! From LMDZ4/libf/dyn3d/gradiv2.F, version 1.1.1.1 2004/05/19 12:53:07
10      ! P. Le Van      ! P. Le Van
11      ! Calcul de grad div du vecteur v.      ! Calcul du gradient de la divergence du vecteur v.
12    
13      USE dimens_m, ONLY : llm      USE comgeom, ONLY: cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, &
     USE paramet_m, ONLY : ip1jm, ip1jmp1, jjp1  
     USE comgeom, ONLY : cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, &  
14           unsapolsga1           unsapolsga1
15      USE filtreg_m, ONLY : filtreg      USE dimens_m, ONLY: iim, jjm, llm
16        use divergf_m, only: divergf
17      INTEGER, intent(in):: klevel      USE filtreg_m, ONLY: filtreg
18        use grad_m, only: grad
19      ! composantes covariantes de v:      use laplacien_m, only: laplacien
20      REAL, intent(in):: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel)      use nr_util, only: assert_eq, assert
21    
22        ! Composantes covariantes de v :
23        REAL, intent(in):: xcov(:, :, :) ! (iim + 1, jjm + 1, klevel)
24        REAL, intent(in):: ycov(:, :, :) ! (iim + 1, jjm, klevel)
25    
26      integer, intent(in):: ld      integer, intent(in):: ld
27      REAL, intent(out):: gdx(ip1jmp1, klevel), gdy(ip1jm, klevel)      REAL, intent(out):: gdx(:, :, :) ! (iim + 1, jjm + 1, klevel)
28        REAL, intent(out):: gdy(:, :, :) ! (iim + 1, jjm, klevel)
29      real, intent(in):: cdivu      real, intent(in):: cdivu
30    
31      ! Variables locales :      ! Variables locales :
32      REAL div(ip1jmp1, llm)      REAL nugrads, div(iim + 1, jjm + 1, llm)
33      REAL nugrads      INTEGER iter, klevel
     INTEGER iter  
34    
35      !--------------------------------------------------------------      !--------------------------------------------------------------
36    
37      gdx = xcov      call assert((/size(xcov, 1), size(ycov, 1), size(gdx, 1), size(gdy, 1)/) &
38      gdy = ycov           == iim + 1, "gradiv2 iim")
39        call assert((/size(xcov, 2) - 1, size(ycov, 2), size(gdx, 2) - 1, &
40             size(gdy, 2)/) == jjm, "gradiv2 iim")
41        klevel = assert_eq(size(xcov, 3), size(ycov, 3), size(gdx, 3), &
42             size(gdy, 3), "gradiv2 klevel")
43    
44      CALL divergf(klevel, gdx, gdy, div)      CALL divergf(klevel, xcov, ycov, div)
45    
46      IF (ld > 1) THEN      IF (ld > 1) THEN
47         CALL laplacien(klevel, div, div)         CALL laplacien(klevel, div)
48    
49         ! Itération de l'opérateur laplacien_gam         ! Itération de l'opérateur laplacien_gam
50         DO iter = 1, ld -2         DO iter = 1, ld -2
# Line 47  contains Line 53  contains
53         END DO         END DO
54      ENDIF      ENDIF
55    
56      CALL filtreg(div, jjp1, klevel, 2, 1, .TRUE., 1)      CALL filtreg(div, jjm + 1, klevel, 2, 1, .TRUE.)
57      CALL grad(klevel, div, gdx, gdy)      CALL grad(klevel, div, gdx, gdy)
58      nugrads = (-1.)**ld * cdivu      nugrads = (-1.)**ld * cdivu
59    
Legend:
Removed from v.56  
changed lines
  Added in v.65
  ViewVC Help
Powered by ViewVC 1.1.21