/[lmdze]/trunk/Sources/dyn3d/massbar.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/massbar.f

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

revision 163 by guez, Wed Apr 29 15:47:56 2015 UTC revision 164 by guez, Tue Jul 28 14:53:31 2015 UTC
# Line 12  contains Line 12  contains
12      ! Calcule les moyennes en x et y de la masse d'air dans chaque      ! Calcule les moyennes en x et y de la masse d'air dans chaque
13      ! maille. Cf. "inigeom.txt" et "massbar.txt".      ! maille. Cf. "inigeom.txt" et "massbar.txt".
14    
15        USE comgeom, ONLY: alpha1p2, alpha1p4, alpha2p3, alpha3p4
16      USE dimens_m, ONLY: iim, llm      USE dimens_m, ONLY: iim, llm
17      USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1      USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1
     USE comgeom, ONLY: alpha1p2, alpha1p4, alpha2p3, alpha3p4  
18    
19      REAL, intent(in):: masse(ip1jmp1,llm)      REAL, intent(in):: masse(ip1jmp1, llm)
20      real, intent(out):: massebx(ip1jmp1,llm), masseby(ip1jm,llm)      real, intent(out):: massebx(ip1jmp1, llm), masseby(ip1jm, llm)
21    
22      ! Local:      ! Local:
23      integer l, ij      integer l, ij
# Line 26  contains Line 26  contains
26    
27      DO l = 1, llm      DO l = 1, llm
28         DO ij = 1, ip1jmp1 - 1         DO ij = 1, ip1jmp1 - 1
29            massebx(ij,l) = masse(ij, l) * alpha1p2(ij) + &            massebx(ij, l) = masse(ij, l) * alpha1p2(ij) + &
30                 masse(ij+1, l) * alpha3p4(ij+1)                 masse(ij+1, l) * alpha3p4(ij+1)
31         ENDDO         ENDDO
32    
33         ! correction pour massebx(iip1,j)         ! correction pour massebx(iip1, j)
34         ! massebx(iip1,j)= massebx(1,j)         ! massebx(iip1, j)= massebx(1, j)
35         DO ij = iip1, ip1jmp1, iip1         DO ij = iip1, ip1jmp1, iip1
36            massebx(ij,l) = massebx(ij - iim,l)            massebx(ij, l) = massebx(ij - iim, l)
37         ENDDO         ENDDO
38    
39         DO ij = 1,ip1jm         DO ij = 1, ip1jm
40            masseby(ij,l) = masse(ij, l) * alpha2p3(ij) + &            masseby(ij, l) = masse(ij, l) * alpha2p3(ij) + &
41                 masse(ij+iip1, l) * alpha1p4(ij+iip1)                 masse(ij+iip1, l) * alpha1p4(ij+iip1)
42         ENDDO         ENDDO
43      end DO      end DO
Legend:
Removed from v.163  
changed lines
  Added in v.164
  ViewVC Help
Powered by ViewVC 1.1.21