/[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.22 by jmc, Sun Jun 19 21:40:08 2005 UTC revision 1.23 by jmc, Wed Jul 13 00:36:01 2005 UTC
# Line 309  C         by formula Line 309  C         by formula
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    
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22