[MITgcm-support] OBCS eta bug

[Gianmaria Sannino]

Hello,

There is a bug in obcs_apply_eta.F

On the western boundary condition is applied OBEeta instead of OBWeta

C     Set model variables to OB values on East/West Boundaries
        DO j=1-Oly,sNy+Oly
C Eastern boundary
          i_obc = OB_Ie(j,bi,bj)
          IF (i_obc.NE.0) THEN
            IF ( ksurfW(i_obc,j,bi,bj).LE.Nr )
     &        etaFld(i_obc,j,bi,bj) = OBEeta(j,bi,bj)
          ENDIF
C Western boundary
          i_obc=OB_Iw(j,bi,bj)
          IF (i_obc.NE.0) THEN
            IF ( ksurfW(i_obc+1,j,bi,bj).LE.Nr )
***************************************************************
     &        etaFld(i_obc,j,bi,bj) = OBEeta(j,bi,bj)
***************************************************************
          ENDIF
        ENDDO



ciao
gianmaria
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20060209/c733b695/attachment.htm>