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*rSphere*COS(yc(I,J,bi,bj)*deg2rad) |
dxF(I,J,bi,bj) = delX(iG+i-1)*deg2rad* |
106 |
|
& 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) |