22 |
|
|
23 |
C === Routine arguments === |
C === Routine arguments === |
24 |
C nSepBnd :: Number of latitude bands |
C nSepBnd :: Number of latitude bands |
25 |
C ySepBnd :: latitude band egdes (from 1 to nSepBnd-1) |
C ySepBnd :: Southern latitude egde (from 2 to nSepBnd, 1 is not used) |
26 |
INTEGER myNr, mskNr, kLoc |
INTEGER myNr, mskNr, kLoc |
27 |
INTEGER nSepBnd |
INTEGER nSepBnd |
28 |
_RS ySepBnd(nSepBnd) |
_RS ySepBnd(nSepBnd) |
83 |
noPnts(n)=.FALSE. |
noPnts(n)=.FALSE. |
84 |
ENDIF |
ENDIF |
85 |
IF (arrMask(i,j,km,bi,bj).NE.0.) THEN |
IF (arrMask(i,j,km,bi,bj).NE.0.) THEN |
86 |
theMin(n)=min(theMin(n),tmpVal) |
theMin(n)=MIN(theMin(n),tmpVal) |
87 |
theMax(n)=max(theMax(n),tmpVal) |
theMax(n)=MAX(theMax(n),tmpVal) |
88 |
tmpVol = arrArea(i,j,bi,bj)*arrhFac(i,j,km,bi,bj)*arrDr(k) |
tmpVol = arrArea(i,j,bi,bj)*arrhFac(i,j,km,bi,bj)*arrDr(k) |
89 |
& *arrMask(i,j,km,bi,bj) |
& *arrMask(i,j,km,bi,bj) |
90 |
theVol(n) = theVol(n) + tmpVol |
theVol(n) = theVol(n) + tmpVol |
131 |
C *==========================================================* |
C *==========================================================* |
132 |
C |
C |
133 |
C !INPUT PARAMETERS: |
C !INPUT PARAMETERS: |
134 |
C nBnd :: Nb of latitude band |
C nBnd :: Number of latitude bands |
135 |
C yBnd :: latitude of band boundaries |
C yBnd :: latitude of southern boundary (for each lat. band) |
136 |
C yLoc :: current latitude |
C yLoc :: current latitude |
137 |
INTEGER nBnd |
INTEGER nBnd |
138 |
_RS yBnd(nBnd) |
_RS yBnd(nBnd) |
141 |
C !LOCAL VARIABLES: |
C !LOCAL VARIABLES: |
142 |
INTEGER n |
INTEGER n |
143 |
|
|
144 |
NLATBND = nBnd |
NLATBND = 1 |
145 |
DO n=1,nBnd-1 |
DO n=2,nBnd |
146 |
IF (yLoc .LE. yBnd(n)) NLATBND = n |
IF (yLoc .GT. yBnd(n)) NLATBND = n |
147 |
ENDDO |
ENDDO |
148 |
|
|
149 |
RETURN |
RETURN |