/[MITgcm]/MITgcm/model/src/ini_spherical_polar_grid.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_spherical_polar_grid.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.20 by adcroft, Thu Sep 27 18:14:52 2001 UTC revision 1.25 by mlosch, Fri Feb 8 13:01:25 2008 UTC
# Line 300  C--     Calculate vertical face area (vo Line 300  C--     Calculate vertical face area (vo
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    
# Line 328  C--     Cosine(lat) scaling Line 330  C--     Cosine(lat) scaling
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
# Line 341  C--     Cosine(lat) scaling Line 345  C--     Cosine(lat) scaling
345         ENDDO ! bi         ENDDO ! bi
346        ENDDO ! bj        ENDDO ! bj
347    
348          IF ( rotateGrid ) THEN
349           CALL ROTATE_SPHERICAL_POLAR_GRID( xC, yC, myThid )
350           CALL ROTATE_SPHERICAL_POLAR_GRID( xG, yG, myThid )
351           CALL CALC_ANGLES( myThid )
352          ENDIF
353    
354        RETURN        RETURN
355        END        END
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22