/[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.7 by adcroft, Thu Jul 2 14:16:24 1998 UTC revision 1.10 by cnh, Fri Nov 6 22:44:46 1998 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  CStartOfInterface  CStartOfInterface
6        SUBROUTINE INI_CARTESIAN_GRID( myThid )        SUBROUTINE INI_CARTESIAN_GRID( myThid )
# Line 63  C     zLower   upper and lower faces. Line 63  C     zLower   upper and lower faces.
63  C     I,J,K  C     I,J,K
64        _RL    xG, yG, zG        _RL    xG, yG, zG
65        _RL    phi        _RL    phi
66        _RL    zUpper(Nz), zLower(Nz)        _RL    zUpper(Nr), zLower(Nr)
67        _RL    xBase, yBase        _RL    xBase, yBase
68        INTEGER iG, jG        INTEGER iG, jG
69        INTEGER bi, bj        INTEGER bi, bj
# Line 143  C     dxC, dyC is separation between cel Line 143  C     dxC, dyC is separation between cel
143         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
144          DO J=1,sNy          DO J=1,sNy
145           DO I=1,sNx           DO I=1,sNx
146            dxC(I,J,bi,bj)    = (dxF(I,J,bi,bj)+dxF(I-1,J,bi,bj))*0.5 D0            dxC(I,J,bi,bj)    = (dxF(I,J,bi,bj)+dxF(I-1,J,bi,bj))*0.5 _d 0
147            dyC(I,J,bi,bj)    = (dyF(I,J,bi,bj)+dyF(I,J-1,bi,bj))*0.5 D0            dyC(I,J,bi,bj)    = (dyF(I,J,bi,bj)+dyF(I,J-1,bi,bj))*0.5 _d 0
148           ENDDO           ENDDO
149          ENDDO          ENDDO
150         ENDDO         ENDDO
# Line 156  C     Calculate vertical face area Line 156  C     Calculate vertical face area
156         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
157          DO J=1,sNy          DO J=1,sNy
158           DO I=1,sNx           DO I=1,sNx
159            zA(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            tanPhiAtU(I,J,bi,bj) = 0. _d 0            tanPhiAtU(I,J,bi,bj) = 0. _d 0
161            tanPhiAtV(I,J,bi,bj) = 0. _d 0            tanPhiAtV(I,J,bi,bj) = 0. _d 0
162           ENDDO           ENDDO
163          ENDDO          ENDDO
164         ENDDO         ENDDO
165        ENDDO        ENDDO
166        _EXCH_XY_R4 (zA       , myThid )        _EXCH_XY_R4 (rA       , myThid )
167        _EXCH_XY_R4 (tanPhiAtU , myThid )        _EXCH_XY_R4 (tanPhiAtU , myThid )
168        _EXCH_XY_R4 (tanPhiAtV , myThid )        _EXCH_XY_R4 (tanPhiAtV , myThid )
169    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22