[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