/[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.13 by adcroft, Mon Mar 27 22:25:44 2000 UTC revision 1.14 by adcroft, Fri Feb 2 21:04:48 2001 UTC
# Line 50  C     == Local variables == Line 50  C     == Local variables ==
50  C     xG, yG - Global coordinate location.  C     xG, yG - Global coordinate location.
51  C     xBase  - South-west corner location for process.  C     xBase  - South-west corner location for process.
52  C     yBase  C     yBase
53    C     zUpper - Work arrays for upper and lower
54    C     zLower   cell-face heights.
55    C     phi    - Temporary scalar
56  C     xBase  - Temporaries for lower corner coordinate  C     xBase  - Temporaries for lower corner coordinate
57  C     yBase  C     yBase
58  C     iG, jG - Global coordinate index. Usually used to hold  C     iG, jG - Global coordinate index. Usually used to hold
# Line 58  C     bi,bj  - Loop counters Line 61  C     bi,bj  - Loop counters
61  C     zUpper - Temporary arrays holding z coordinates of  C     zUpper - Temporary arrays holding z coordinates of
62  C     zLower   upper and lower faces.  C     zLower   upper and lower faces.
63  C     I,J,K  C     I,J,K
64        _RL    xG, yG        _RL    xGloc, yGloc
65        _RL    xBase, yBase        _RL    xBase, yBase
66        INTEGER iG, jG        INTEGER iG, jG
67        INTEGER bi, bj        INTEGER bi, bj
# Line 84  C     Set up my local grid first Line 87  C     Set up my local grid first
87          DO j=1,jG-1          DO j=1,jG-1
88           yBase = yBase + delY(j)           yBase = yBase + delY(j)
89          ENDDO          ENDDO
90          yG = yBase          yGloc = yBase
91          DO J=1,sNy          DO J=1,sNy
92           xG = xBase           xGloc = xBase
93           DO I=1,sNx           DO I=1,sNx
94            xc(I,J,bi,bj)  = xG + delX(iG+i-1)*0.5 _d 0            xG(I,J,bi,bj)  = xGloc
95            yc(I,J,bi,bj)  = yG + delY(jG+j-1)*0.5 _d 0            yG(I,J,bi,bj)  = yGloc
96            xG = xG + delX(iG+I-1)            xc(I,J,bi,bj)  = xGloc + delX(iG+i-1)*0.5 _d 0
97              yc(I,J,bi,bj)  = yGloc + delY(jG+j-1)*0.5 _d 0
98              xGloc = xGloc + delX(iG+I-1)
99            dxF(I,J,bi,bj) = delX(iG+i-1)            dxF(I,J,bi,bj) = delX(iG+i-1)
100            dyF(I,J,bi,bj) = delY(jG+j-1)            dyF(I,J,bi,bj) = delY(jG+j-1)
101           ENDDO           ENDDO
102           yG = yG + delY(jG+J-1)           yGloc = yGloc + delY(jG+J-1)
103          ENDDO          ENDDO
104         ENDDO         ENDDO
105        ENDDO        ENDDO
# Line 154  C     Calculate vertical face area Line 159  C     Calculate vertical face area
159            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)
160            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)
161            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)
162              rAz(I,J,bi,bj) = dxV(I,J,bi,bj)*dyU(I,J,bi,bj)
163            tanPhiAtU(I,J,bi,bj) = 0. _d 0            tanPhiAtU(I,J,bi,bj) = 0. _d 0
164            tanPhiAtV(I,J,bi,bj) = 0. _d 0            tanPhiAtV(I,J,bi,bj) = 0. _d 0
165           ENDDO           ENDDO
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22