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

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

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

revision 1.17 by cnh, Wed Sep 26 18:09:15 2001 UTC revision 1.18 by jmc, Wed Jul 13 00:36:01 2005 UTC
# Line 188  c         dYU(I,J,bi,bj) = 0.5*(dYC(I,J, Line 188  c         dYU(I,J,bi,bj) = 0.5*(dYC(I,J,
188           ENDDO           ENDDO
189          ENDDO          ENDDO
190    
191  C     Calculate vertical face area  C--     Calculate vertical face area
192          DO J=1-Oly,sNy+Oly          DO J=1-Oly,sNy+Oly
193           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
194            rA (I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj)            rA (I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj)
195            rAw(I,J,bi,bj) = dxC(I,J,bi,bj)*dyG(I,J,bi,bj)            rAw(I,J,bi,bj) = dxC(I,J,bi,bj)*dyG(I,J,bi,bj)
196            rAs(I,J,bi,bj) = dxG(I,J,bi,bj)*dyC(I,J,bi,bj)            rAs(I,J,bi,bj) = dxG(I,J,bi,bj)*dyC(I,J,bi,bj)
197            rAz(I,J,bi,bj) = dxV(I,J,bi,bj)*dyU(I,J,bi,bj)            rAz(I,J,bi,bj) = dxV(I,J,bi,bj)*dyU(I,J,bi,bj)
198    C--     Set trigonometric terms & grid orientation:
199            tanPhiAtU(I,J,bi,bj) = 0.            tanPhiAtU(I,J,bi,bj) = 0.
200            tanPhiAtV(I,J,bi,bj) = 0.            tanPhiAtV(I,J,bi,bj) = 0.
201              angleCosC(I,J,bi,bj) = 1.
202              angleSinC(I,J,bi,bj) = 0.
203           ENDDO           ENDDO
204          ENDDO          ENDDO
205    
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22