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

Diff of /trunk/dyn3d/groupeun.f

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

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 178 by guez, Fri Mar 11 18:47:26 2016 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.178
  ViewVC Help
Powered by ViewVC 1.1.21