/[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 60 by guez, Mon Jan 30 14:37:26 2012 UTC
# Line 17  contains Line 17  contains
17      USE comgeom, ONLY: apoln, apols, cuvsurcv, cvusurcu, unsaire      USE comgeom, ONLY: apoln, apols, cuvsurcv, cvusurcu, unsaire
18      USE filtreg_m, ONLY: filtreg      USE filtreg_m, ONLY: filtreg
19    
     ! div est un argument de sortie pour le s-prog  
   
     ! variables en arguments  
   
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
Legend:
Removed from v.57  
changed lines
  Added in v.60
  ViewVC Help
Powered by ViewVC 1.1.21