65 |
C I,J,K |
C I,J,K |
66 |
_RL xG, yG, zG |
_RL xG, yG, zG |
67 |
_RL phi |
_RL phi |
68 |
_RL zUpper(Nz), zLower(Nz) |
_RL zUpper(Nr), zLower(Nr) |
69 |
_RL xBase, yBase |
_RL xBase, yBase |
70 |
INTEGER iG, jG |
INTEGER iG, jG |
71 |
INTEGER bi, bj |
INTEGER bi, bj |
168 |
jG = myYGlobalLo + (bj-1)*sNy + J-1 |
jG = myYGlobalLo + (bj-1)*sNy + J-1 |
169 |
latS = yc(i,j,bi,bj)-delY(jG)*0.5 _d 0 |
latS = yc(i,j,bi,bj)-delY(jG)*0.5 _d 0 |
170 |
latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0 |
latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0 |
171 |
zA(I,J,bi,bj) = dyF(I,J,bi,bj) |
rA(I,J,bi,bj) = dyF(I,J,bi,bj) |
172 |
& *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad)) |
& *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad)) |
173 |
tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad) |
tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad) |
174 |
tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad) |
tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad) |
176 |
ENDDO |
ENDDO |
177 |
ENDDO |
ENDDO |
178 |
ENDDO |
ENDDO |
179 |
_EXCH_XY_R4 (zA , myThid ) |
_EXCH_XY_R4 (rA , myThid ) |
180 |
_EXCH_XY_R4 (tanPhiAtU , myThid ) |
_EXCH_XY_R4 (tanPhiAtU , myThid ) |
181 |
_EXCH_XY_R4 (tanPhiAtV , myThid ) |
_EXCH_XY_R4 (tanPhiAtV , myThid ) |
182 |
|
|