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

Diff of /trunk/dyn3d/groupeun.f

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

trunk/Sources/dyn3d/groupeun.f revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC trunk/dyn3d/groupeun.f revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC
# Line 16  SUBROUTINE groupeun(jjmax, llmax, q) Line 16  SUBROUTINE groupeun(jjmax, llmax, q)
16    INTEGER ngroup    INTEGER ngroup
17    PARAMETER (ngroup=3)    PARAMETER (ngroup=3)
18    
19    REAL airen, airecn, qn    REAL airen, qn
20    REAL aires, airecs, qs    REAL aires, qs
21    
22    INTEGER i, j, l, ig, j1, j2, i0, jd    INTEGER i, j, l, ig, j1, j2, i0, jd
23    
# Line 30  SUBROUTINE groupeun(jjmax, llmax, q) Line 30  SUBROUTINE groupeun(jjmax, llmax, q)
30        DO j = j1 - jd, j2 - jd        DO j = j1 - jd, j2 - jd
31          DO i0 = 1, iim, 2**(ngroup-ig+1)          DO i0 = 1, iim, 2**(ngroup-ig+1)
32            airen = 0.            airen = 0.
           airecn = 0.  
33            qn = 0.            qn = 0.
34            aires = 0.            aires = 0.
           airecs = 0.  
35            qs = 0.            qs = 0.
36            DO i = i0, i0 + 2**(ngroup-ig+1) - 1            DO i = i0, i0 + 2**(ngroup-ig+1) - 1
37              airen = airen + aire_2d(i, j)              airen = airen + aire_2d(i, j)
# Line 41  SUBROUTINE groupeun(jjmax, llmax, q) Line 39  SUBROUTINE groupeun(jjmax, llmax, q)
39              qn = qn + q(i, j, l)              qn = qn + q(i, j, l)
40              qs = qs + q(i, jjp1-j+1-jd, l)              qs = qs + q(i, jjp1-j+1-jd, l)
41            END DO            END DO
           airecn = 0.  
           airecs = 0.  
42            DO i = i0, i0 + 2**(ngroup-ig+1) - 1            DO i = i0, i0 + 2**(ngroup-ig+1) - 1
43              q(i, j, l) = qn*aire_2d(i, j)/airen              q(i, j, l) = qn*aire_2d(i, j)/airen
44              q(i, jjp1-j+1-jd, l) = qs*aire_2d(i, jjp1-j+1)/aires              q(i, jjp1-j+1-jd, l) = qs*aire_2d(i, jjp1-j+1)/aires
Legend:
Removed from v.134  
changed lines
  Added in v.254
  ViewVC Help
Powered by ViewVC 1.1.21