300 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
301 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
302 |
C by formula |
C by formula |
303 |
lat=yC(I,J,bi,bj) |
lat =0.5 _d 0*(yGloc(I,J)+yGloc(I,J+1)) |
304 |
dlon=delX( iGl(I,bi) ) |
dlon=0.5 _d 0*( delX( iGl(I,bi) ) + delX( iGl(I-1,bi) ) ) |
305 |
dlat=0.5*( delY( jGl(J,bj) ) + delY( jGl(J-1,bj) ) ) |
dlat=0.5 _d 0*( delY( jGl(J,bj) ) + delY( jGl(J-1,bj) ) ) |
306 |
rAz(I,J,bi,bj) = rSphere*rSphere*dlon*deg2rad |
rAz(I,J,bi,bj) = rSphere*rSphere*dlon*deg2rad |
307 |
& *abs( sin(lat*deg2rad)-sin((lat-dlat)*deg2rad) ) |
& *abs( sin(lat*deg2rad)-sin((lat-dlat)*deg2rad) ) |
308 |
IF (abs(lat).GT.90..OR.abs(lat-dlat).GT.90.) rAz(I,J,bi,bj)=0. |
IF (abs(lat).GT.90..OR.abs(lat-dlat).GT.90.) rAz(I,J,bi,bj)=0. |
309 |
ENDDO |
ENDDO |
310 |
ENDDO |
ENDDO |
311 |
|
|
312 |
C-- Calculate trigonometric terms |
C-- Calculate trigonometric terms & grid orientation: |
313 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
314 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
315 |
lat=0.5*(yGloc(I,J)+yGloc(I,J+1)) |
lat=0.5*(yGloc(I,J)+yGloc(I,J+1)) |
316 |
tanPhiAtU(i,j,bi,bj)=tan(lat*deg2rad) |
tanPhiAtU(I,J,bi,bj)=tan(lat*deg2rad) |
317 |
lat=0.5*(yGloc(I,J)+yGloc(I+1,J)) |
lat=0.5*(yGloc(I,J)+yGloc(I+1,J)) |
318 |
tanPhiAtV(i,j,bi,bj)=tan(lat*deg2rad) |
tanPhiAtV(I,J,bi,bj)=tan(lat*deg2rad) |
319 |
|
angleCosC(I,J,bi,bj) = 1. |
320 |
|
angleSinC(I,J,bi,bj) = 0. |
321 |
ENDDO |
ENDDO |
322 |
ENDDO |
ENDDO |
323 |
|
|
330 |
& **cosPower |
& **cosPower |
331 |
cosFacV(J,bi,bj)=COS((yC(1,J,bi,bj)-0.5*delY(jG))*deg2rad) |
cosFacV(J,bi,bj)=COS((yC(1,J,bi,bj)-0.5*delY(jG))*deg2rad) |
332 |
& **cosPower |
& **cosPower |
333 |
|
cosFacU(J,bi,bj)=ABS(cosFacU(J,bi,bj)) |
334 |
|
cosFacV(J,bi,bj)=ABS(cosFacV(J,bi,bj)) |
335 |
sqcosFacU(J,bi,bj)=sqrt(cosFacU(J,bi,bj)) |
sqcosFacU(J,bi,bj)=sqrt(cosFacU(J,bi,bj)) |
336 |
sqcosFacV(J,bi,bj)=sqrt(cosFacV(J,bi,bj)) |
sqcosFacV(J,bi,bj)=sqrt(cosFacV(J,bi,bj)) |
337 |
ELSE |
ELSE |