/[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.19 by jmc, Sun Jul 31 22:07:48 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    
# Line 211  C--     Cosine(lat) scaling Line 214  C--     Cosine(lat) scaling
214         ENDDO ! bi         ENDDO ! bi
215        ENDDO ! bj        ENDDO ! bj
216    
217    C--   Set default (=whole domain) for where relaxation to climatology applies
218          IF ( latBandClimRelax.EQ.UNSET_RL ) THEN
219            _BEGIN_MASTER(myThid)
220            latBandClimRelax = 0.
221            DO j=1,Ny
222              latBandClimRelax = latBandClimRelax + delY(j)
223            ENDDO
224            latBandClimRelax = latBandClimRelax*3. _d 0
225            _END_MASTER(myThid)
226          ENDIF
227    
228        RETURN        RETURN
229        END        END
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22