--- MITgcm/model/src/ini_cartesian_grid.F 1998/07/02 14:16:24 1.7 +++ MITgcm/model/src/ini_cartesian_grid.F 1998/11/06 22:44:46 1.10 @@ -1,6 +1,6 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_cartesian_grid.F,v 1.7 1998/07/02 14:16:24 adcroft Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_cartesian_grid.F,v 1.10 1998/11/06 22:44:46 cnh Exp $ -#include "CPP_EEOPTIONS.h" +#include "CPP_OPTIONS.h" CStartOfInterface SUBROUTINE INI_CARTESIAN_GRID( myThid ) @@ -63,7 +63,7 @@ C I,J,K _RL xG, yG, zG _RL phi - _RL zUpper(Nz), zLower(Nz) + _RL zUpper(Nr), zLower(Nr) _RL xBase, yBase INTEGER iG, jG INTEGER bi, bj @@ -143,8 +143,8 @@ DO bi = myBxLo(myThid), myBxHi(myThid) DO J=1,sNy DO I=1,sNx - dxC(I,J,bi,bj) = (dxF(I,J,bi,bj)+dxF(I-1,J,bi,bj))*0.5 D0 - dyC(I,J,bi,bj) = (dyF(I,J,bi,bj)+dyF(I,J-1,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 + dyC(I,J,bi,bj) = (dyF(I,J,bi,bj)+dyF(I,J-1,bi,bj))*0.5 _d 0 ENDDO ENDDO ENDDO @@ -156,14 +156,14 @@ DO bi = myBxLo(myThid), myBxHi(myThid) DO J=1,sNy DO I=1,sNx - 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) tanPhiAtU(I,J,bi,bj) = 0. _d 0 tanPhiAtV(I,J,bi,bj) = 0. _d 0 ENDDO ENDDO ENDDO ENDDO - _EXCH_XY_R4 (zA , myThid ) + _EXCH_XY_R4 (rA , myThid ) _EXCH_XY_R4 (tanPhiAtU , myThid ) _EXCH_XY_R4 (tanPhiAtV , myThid )