/[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.7 by adcroft, Thu Jul 2 14:16:24 1998 UTC revision 1.8 by cnh, Sat Aug 22 17:51:08 1998 UTC
# Line 65  C     latS         values. Line 65  C     latS         values.
65  C     I,J,K  C     I,J,K
66        _RL    xG, yG, zG        _RL    xG, yG, zG
67        _RL    phi        _RL    phi
68        _RL    zUpper(Nz), zLower(Nz)        _RL    zUpper(Nr), zLower(Nr)
69        _RL    xBase, yBase        _RL    xBase, yBase
70        INTEGER iG, jG        INTEGER iG, jG
71        INTEGER bi, bj        INTEGER bi, bj
# Line 168  C     Calculate vertical face area and t Line 168  C     Calculate vertical face area and t
168            jG = myYGlobalLo + (bj-1)*sNy + J-1            jG = myYGlobalLo + (bj-1)*sNy + J-1
169            latS = yc(i,j,bi,bj)-delY(jG)*0.5 _d 0            latS = yc(i,j,bi,bj)-delY(jG)*0.5 _d 0
170            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0
171            zA(I,J,bi,bj) = dyF(I,J,bi,bj)            rA(I,J,bi,bj) = dyF(I,J,bi,bj)
172       &    *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad))       &    *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad))
173            tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad)            tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad)
174            tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad)            tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad)
# Line 176  C     Calculate vertical face area and t Line 176  C     Calculate vertical face area and t
176          ENDDO          ENDDO
177         ENDDO         ENDDO
178        ENDDO        ENDDO
179        _EXCH_XY_R4 (zA       , myThid )        _EXCH_XY_R4 (rA       , myThid )
180        _EXCH_XY_R4 (tanPhiAtU , myThid )        _EXCH_XY_R4 (tanPhiAtU , myThid )
181        _EXCH_XY_R4 (tanPhiAtV , myThid )        _EXCH_XY_R4 (tanPhiAtV , myThid )
182    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22