/[lmdze]/trunk/dyn3d/divergf.f
ViewVC logotype

Diff of /trunk/dyn3d/divergf.f

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

revision 57 by guez, Mon Jan 30 12:54:02 2012 UTC revision 65 by guez, Thu Sep 20 09:57:03 2012 UTC
# Line 6  contains Line 6  contains
6    
7    SUBROUTINE divergf(klevel, x, y, div)    SUBROUTINE divergf(klevel, x, y, div)
8    
9      ! From libf/dyn3d/divergf.F, v 1.1.1.1 2004/05/19 12:53:05      ! From libf/dyn3d/divergf.F, version 1.1.1.1 2004/05/19 12:53:05
10    
11      ! P. Le Van      ! P. Le Van
12      ! Calcule la divergence à tous les niveaux d'un vecteur de      ! Calcule la divergence à tous les niveaux d'un vecteur de
13      ! composantes x et y. x et y sont des composantes covariantes.      ! composantes x et y. x et y sont des composantes covariantes.
14    
     USE dimens_m, ONLY: iim  
     USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmi1, ip1jmp1, jjp1  
15      USE comgeom, ONLY: apoln, apols, cuvsurcv, cvusurcu, unsaire      USE comgeom, ONLY: apoln, apols, cuvsurcv, cvusurcu, unsaire
16        USE dimens_m, ONLY: iim
17      USE filtreg_m, ONLY: filtreg      USE filtreg_m, ONLY: filtreg
18        USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmi1, ip1jmp1, jjp1
     ! div est un argument de sortie pour le s-prog  
   
     ! variables en arguments  
19    
20      INTEGER, intent(in):: klevel      INTEGER, intent(in):: klevel
21      REAL, intent(in):: x(ip1jmp1, klevel), y(ip1jm, klevel)      REAL, intent(in):: x(ip1jmp1, klevel), y(ip1jm, klevel)
22      real div(ip1jmp1, klevel)      real, intent(out):: div(ip1jmp1, klevel) ! in (unit of x, y) m-2
23    
24      ! variables locales      ! Variables locales :
25    
26      INTEGER l, ij      INTEGER l, ij
27      REAL aiy1(iip1) , aiy2(iip1)      REAL aiy1(iim) , aiy2(iim)
28      REAL sumypn, sumyps      REAL sumypn, sumyps
29    
     REAL SSUM  
   
30      !------------------------------------------------------------      !------------------------------------------------------------
31    
32      DO l = 1, klevel      DO l = 1, klevel
# Line 52  contains Line 46  contains
46            aiy1(ij) = cuvsurcv(ij) * y(ij , l)            aiy1(ij) = cuvsurcv(ij) * y(ij , l)
47            aiy2(ij) = cuvsurcv(ij+ ip1jmi1) * y(ij+ ip1jmi1, l)            aiy2(ij) = cuvsurcv(ij+ ip1jmi1) * y(ij+ ip1jmi1, l)
48         ENDDO         ENDDO
49         sumypn = SSUM (iim, aiy1, 1) / apoln         sumypn = SUM(aiy1) / apoln
50         sumyps = SSUM (iim, aiy2, 1) / apols         sumyps = SUM(aiy2) / apols
51    
52         DO ij = 1, iip1         DO ij = 1, iip1
53            div(ij , l) = - sumypn            div(ij , l) = - sumypn
# Line 61  contains Line 55  contains
55         ENDDO         ENDDO
56      end DO      end DO
57    
58      CALL filtreg(div, jjp1, klevel, 2, 2, .TRUE., 1)      CALL filtreg(div, jjp1, klevel, 2, 2, .TRUE.)
59    
60      DO l = 1, klevel      DO l = 1, klevel
61         DO ij = iip2, ip1jm         DO ij = iip2, ip1jm
Legend:
Removed from v.57  
changed lines
  Added in v.65
  ViewVC Help
Powered by ViewVC 1.1.21