1 |
C $Header$ |
C $Header$ |
2 |
|
|
3 |
#include "CPP_EEOPTIONS.h" |
#include "CPP_OPTIONS.h" |
4 |
|
|
5 |
CStartOfInterface |
CStartOfInterface |
6 |
SUBROUTINE INI_SPHERICAL_POLAR_GRID( myThid ) |
SUBROUTINE INI_SPHERICAL_POLAR_GRID( myThid ) |
102 |
xc(I,J,bi,bj) = xG + delX(iG+i-1)*0.5 _d 0 |
xc(I,J,bi,bj) = xG + delX(iG+i-1)*0.5 _d 0 |
103 |
yc(I,J,bi,bj) = yG + delY(jG+j-1)*0.5 _d 0 |
yc(I,J,bi,bj) = yG + delY(jG+j-1)*0.5 _d 0 |
104 |
xG = xG + delX(iG+I-1) |
xG = xG + delX(iG+I-1) |
105 |
dxF(I,J,bi,bj) = delX(iG+i-1)*deg2rad* |
dxF(I,J,bi,bj) = delX(iG+i-1)*deg2rad |
106 |
& rSphere*COS(yc(I,J,bi,bj)*deg2rad) |
& *rSphere*COS(yc(I,J,bi,bj)*deg2rad) |
107 |
dyF(I,J,bi,bj) = delY(jG+j-1)*deg2rad*rSphere |
dyF(I,J,bi,bj) = delY(jG+j-1)*deg2rad*rSphere |
108 |
ENDDO |
ENDDO |
109 |
yG = yG + delY(jG+J-1) |
yG = yG + delY(jG+J-1) |
171 |
latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0 |
latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0 |
172 |
rA(I,J,bi,bj) = dyF(I,J,bi,bj) |
rA(I,J,bi,bj) = dyF(I,J,bi,bj) |
173 |
& *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad)) |
& *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad)) |
174 |
|
C Area cannot be zero but sin can be if lat if < -90. |
175 |
|
IF ( rA(I,J,bi,bj) .LT. 0. ) rA(I,J,bi,bj) = -rA(I,J,bi,bj) |
176 |
tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad) |
tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad) |
177 |
tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad) |
tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad) |
178 |
ENDDO |
ENDDO |