[Mitgcm-support] EXChanges in ini_depths.F
mitgcm-support at dev.mitgcm.org
mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:15:16 EDT 2003
The reason the following code was in this order was so that
the overlap regions contained land on the lat-lon grid.
With the switch to exchanging Rlo/Rhigh instead of H the
overlaps contain land again...
C _EXCH_XY_R4( H, myThid )
IF (usingSphericalPolarGrid) THEN
DO bj = myByLo(myThid), myByHi(myThid)
DO bi = myBxLo(myThid), myBxHi(myThid)
DO j=1-Oly,sNy+Oly
DO i=1-Olx,sNx+Olx
IF (abs(yC(I,J,bi,bj)).GE.90.) H(I,J,bi,bj)=0.
ENDDO
ENDDO
ENDDO
ENDDO
ENDIF
<sigh>
A.
More information about the MITgcm-support
mailing list