/[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 82 by guez, Wed Mar 5 14:57:53 2014 UTC 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.82  
changed lines
  Added in v.254
  ViewVC Help
Powered by ViewVC 1.1.21