188 |
ENDDO |
ENDDO |
189 |
ENDDO |
ENDDO |
190 |
|
|
191 |
C Calculate vertical face area |
C-- Calculate vertical face area |
192 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
193 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
194 |
rA (I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj) |
rA (I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj) |
195 |
rAw(I,J,bi,bj) = dxC(I,J,bi,bj)*dyG(I,J,bi,bj) |
rAw(I,J,bi,bj) = dxC(I,J,bi,bj)*dyG(I,J,bi,bj) |
196 |
rAs(I,J,bi,bj) = dxG(I,J,bi,bj)*dyC(I,J,bi,bj) |
rAs(I,J,bi,bj) = dxG(I,J,bi,bj)*dyC(I,J,bi,bj) |
197 |
rAz(I,J,bi,bj) = dxV(I,J,bi,bj)*dyU(I,J,bi,bj) |
rAz(I,J,bi,bj) = dxV(I,J,bi,bj)*dyU(I,J,bi,bj) |
198 |
|
C-- Set trigonometric terms & grid orientation: |
199 |
tanPhiAtU(I,J,bi,bj) = 0. |
tanPhiAtU(I,J,bi,bj) = 0. |
200 |
tanPhiAtV(I,J,bi,bj) = 0. |
tanPhiAtV(I,J,bi,bj) = 0. |
201 |
|
angleCosC(I,J,bi,bj) = 1. |
202 |
|
angleSinC(I,J,bi,bj) = 0. |
203 |
ENDDO |
ENDDO |
204 |
ENDDO |
ENDDO |
205 |
|
|
214 |
ENDDO ! bi |
ENDDO ! bi |
215 |
ENDDO ! bj |
ENDDO ! bj |
216 |
|
|
217 |
|
C-- Set default (=whole domain) for where relaxation to climatology applies |
218 |
|
IF ( latBandClimRelax.EQ.UNSET_RL ) THEN |
219 |
|
_BEGIN_MASTER(myThid) |
220 |
|
latBandClimRelax = 0. |
221 |
|
DO j=1,Ny |
222 |
|
latBandClimRelax = latBandClimRelax + delY(j) |
223 |
|
ENDDO |
224 |
|
latBandClimRelax = latBandClimRelax*3. _d 0 |
225 |
|
_END_MASTER(myThid) |
226 |
|
ENDIF |
227 |
|
|
228 |
RETURN |
RETURN |
229 |
END |
END |