/[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.11 by adcroft, Mon Nov 30 23:45:25 1998 UTC revision 1.14 by adcroft, Mon Mar 27 22:25:44 2000 UTC
# Line 33  C     | Under the spherical polar grid m Line 33  C     | Under the spherical polar grid m
33  C     | in X and Y are in degrees. Distance in Z are in m or Pa  |  C     | in X and Y are in degrees. Distance in Z are in m or Pa  |
34  C     | depending on the vertical gridding mode.                 |  C     | depending on the vertical gridding mode.                 |
35  C     \==========================================================/  C     \==========================================================/
36          IMPLICIT NONE
37    
38  C     === Global variables ===  C     === Global variables ===
39  #include "SIZE.h"  #include "SIZE.h"
# Line 47  CEndOfInterface Line 48  CEndOfInterface
48    
49  C     == Local variables ==  C     == Local variables ==
50  C     xG, yG - Global coordinate location.  C     xG, yG - Global coordinate location.
 C     zG  
51  C     xBase  - South-west corner location for process.  C     xBase  - South-west corner location for process.
52  C     yBase  C     yBase
 C     zUpper - Work arrays for upper and lower  
 C     zLower   cell-face heights.  
 C     phi    - Temporary scalar  
53  C     iG, jG - Global coordinate index. Usually used to hold  C     iG, jG - Global coordinate index. Usually used to hold
54  C              the south-west global coordinate of a tile.  C              the south-west global coordinate of a tile.
55  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
# Line 63  C     yBase Line 60  C     yBase
60  C     lat, latN, - Temporary variables used to hold latitude  C     lat, latN, - Temporary variables used to hold latitude
61  C     latS         values.  C     latS         values.
62  C     I,J,K  C     I,J,K
63        _RL    xG, yG, zG        _RL    xG, yG
       _RL    phi  
       _RL    zUpper(Nr), zLower(Nr)  
64        _RL    xBase, yBase        _RL    xBase, yBase
65        INTEGER iG, jG        INTEGER iG, jG
66        INTEGER bi, bj        INTEGER bi, bj
67        INTEGER  I,  J, K        INTEGER  I,  J
68        _RL lat, latS, latN        _RL lat, latS, latN
69    
70  C--   Example of inialisation for spherical polar grid  C--   Example of inialisation for spherical polar grid
# Line 171  C     Calculate vertical face area and t Line 166  C     Calculate vertical face area and t
166          DO J=1,sNy          DO J=1,sNy
167           DO I=1,sNx           DO I=1,sNx
168            jG = myYGlobalLo + (bj-1)*sNy + J-1            jG = myYGlobalLo + (bj-1)*sNy + J-1
169              iG = myXGlobalLo + (bi-1)*sNx + I-1
170            latS = yc(i,j,bi,bj)-delY(jG)*0.5 _d 0            latS = yc(i,j,bi,bj)-delY(jG)*0.5 _d 0
171            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0
172  #ifdef OLD_UV_GEOMETRY  #ifdef OLD_UV_GEOMETRY
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22