185 |
_EXCH_XY_R4(rDyF, myThid ) |
_EXCH_XY_R4(rDyF, myThid ) |
186 |
_EXCH_XY_R4(rDxV, myThid ) |
_EXCH_XY_R4(rDxV, myThid ) |
187 |
_EXCH_XY_R4(rDyU, myThid ) |
_EXCH_XY_R4(rDyU, myThid ) |
188 |
C Calculate vertical face area |
C Calculate vertical face area and trigonometric terms |
189 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
190 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
191 |
DO J=1,sNy |
DO J=1,sNy |
195 |
latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0 |
latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0 |
196 |
zA(I,J,bi,bj) = dyF(I,J,bi,bj) |
zA(I,J,bi,bj) = dyF(I,J,bi,bj) |
197 |
& *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad)) |
& *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad)) |
198 |
|
tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad) |
199 |
|
tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad) |
200 |
ENDDO |
ENDDO |
201 |
ENDDO |
ENDDO |
202 |
ENDDO |
ENDDO |
239 |
_EXCH_XYZ_R4(maskW , myThid ) |
_EXCH_XYZ_R4(maskW , myThid ) |
240 |
_EXCH_XYZ_R4(maskS , myThid ) |
_EXCH_XYZ_R4(maskS , myThid ) |
241 |
_EXCH_XY_R4 (zA , myThid ) |
_EXCH_XY_R4 (zA , myThid ) |
|
|
|
|
CcnhDebugStarts |
|
|
tanPhiAtU = 0. _d 0 |
|
|
tanPhiAtV = 0. _d 0 |
|
242 |
_EXCH_XY_R4 (tanPhiAtU , myThid ) |
_EXCH_XY_R4 (tanPhiAtU , myThid ) |
243 |
_EXCH_XY_R4 (tanPhiAtV , myThid ) |
_EXCH_XY_R4 (tanPhiAtV , myThid ) |
|
CcnhDebugEnds |
|
|
|
|
244 |
|
|
245 |
C |
C |
246 |
RETURN |
RETURN |