/[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.9 by cnh, Mon Nov 2 03:34:11 1998 UTC revision 1.15 by cnh, Sun Feb 4 14:38:47 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CStartOfInterface
7        SUBROUTINE INI_CARTESIAN_GRID( myThid )        SUBROUTINE INI_CARTESIAN_GRID( myThid )
# Line 33  C     | Under the cartesian grid mode pr Line 34  C     | Under the cartesian grid mode pr
34  C     | and Y are in metres. Disktance in Z are in m or Pa       |  C     | and Y are in metres. Disktance in Z are in m or Pa       |
35  C     | depending on the vertical gridding mode.                 |  C     | depending on the vertical gridding mode.                 |
36  C     \==========================================================/  C     \==========================================================/
37          IMPLICIT NONE
38    
39  C     === Global variables ===  C     === Global variables ===
40  #include "SIZE.h"  #include "SIZE.h"
# Line 47  CEndOfInterface Line 49  CEndOfInterface
49    
50  C     == Local variables ==  C     == Local variables ==
51  C     xG, yG - Global coordinate location.  C     xG, yG - Global coordinate location.
 C     zG  
52  C     xBase  - South-west corner location for process.  C     xBase  - South-west corner location for process.
53  C     yBase  C     yBase
54  C     zUpper - Work arrays for upper and lower  C     zUpper - Work arrays for upper and lower
# Line 61  C     bi,bj  - Loop counters Line 62  C     bi,bj  - Loop counters
62  C     zUpper - Temporary arrays holding z coordinates of  C     zUpper - Temporary arrays holding z coordinates of
63  C     zLower   upper and lower faces.  C     zLower   upper and lower faces.
64  C     I,J,K  C     I,J,K
65        _RL    xG, yG, zG        _RL    xGloc, yGloc
       _RL    phi  
       _RL    zUpper(Nr), zLower(Nr)  
66        _RL    xBase, yBase        _RL    xBase, yBase
67        INTEGER iG, jG        INTEGER iG, jG
68        INTEGER bi, bj        INTEGER bi, bj
69        INTEGER  I,  J, K        INTEGER  I,  J
70    
71  C--   Simple example of inialisation on cartesian grid  C--   Simple example of inialisation on cartesian grid
72  C--   First set coordinates of cell centers  C--   First set coordinates of cell centers
# Line 89  C     Set up my local grid first Line 88  C     Set up my local grid first
88          DO j=1,jG-1          DO j=1,jG-1
89           yBase = yBase + delY(j)           yBase = yBase + delY(j)
90          ENDDO          ENDDO
91          yG = yBase          yGloc = yBase
92          DO J=1,sNy          DO J=1,sNy
93           xG = xBase           xGloc = xBase
94           DO I=1,sNx           DO I=1,sNx
95            xc(I,J,bi,bj)  = xG + delX(iG+i-1)*0.5 _d 0            xG(I,J,bi,bj)  = xGloc
96            yc(I,J,bi,bj)  = yG + delY(jG+j-1)*0.5 _d 0            yG(I,J,bi,bj)  = yGloc
97            xG = xG + delX(iG+I-1)            xc(I,J,bi,bj)  = xGloc + delX(iG+i-1)*0.5 _d 0
98              yc(I,J,bi,bj)  = yGloc + delY(jG+j-1)*0.5 _d 0
99              xGloc = xGloc + delX(iG+I-1)
100            dxF(I,J,bi,bj) = delX(iG+i-1)            dxF(I,J,bi,bj) = delX(iG+i-1)
101            dyF(I,J,bi,bj) = delY(jG+j-1)            dyF(I,J,bi,bj) = delY(jG+j-1)
102           ENDDO           ENDDO
103           yG = yG + delY(jG+J-1)           yGloc = yGloc + delY(jG+J-1)
104          ENDDO          ENDDO
105         ENDDO         ENDDO
106        ENDDO        ENDDO
# Line 156  C     Calculate vertical face area Line 157  C     Calculate vertical face area
157         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
158          DO J=1,sNy          DO J=1,sNy
159           DO I=1,sNx           DO I=1,sNx
160            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)
161              rAw(I,J,bi,bj) = dxC(I,J,bi,bj)*dyG(I,J,bi,bj)
162              rAs(I,J,bi,bj) = dxG(I,J,bi,bj)*dyC(I,J,bi,bj)
163              rAz(I,J,bi,bj) = dxV(I,J,bi,bj)*dyU(I,J,bi,bj)
164            tanPhiAtU(I,J,bi,bj) = 0. _d 0            tanPhiAtU(I,J,bi,bj) = 0. _d 0
165            tanPhiAtV(I,J,bi,bj) = 0. _d 0            tanPhiAtV(I,J,bi,bj) = 0. _d 0
166           ENDDO           ENDDO
# Line 164  C     Calculate vertical face area Line 168  C     Calculate vertical face area
168         ENDDO         ENDDO
169        ENDDO        ENDDO
170        _EXCH_XY_R4 (rA       , myThid )        _EXCH_XY_R4 (rA       , myThid )
171          _EXCH_XY_R4 (rAw      , myThid )
172          _EXCH_XY_R4 (rAs      , myThid )
173        _EXCH_XY_R4 (tanPhiAtU , myThid )        _EXCH_XY_R4 (tanPhiAtU , myThid )
174        _EXCH_XY_R4 (tanPhiAtV , myThid )        _EXCH_XY_R4 (tanPhiAtV , myThid )
175    
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22