[MITgcm-support] OBCS eta bug

[Jean-Michel Campin]

Hi Gianmaria,

Thanks for your report.
I fixed it and should be OK now.

Ciao,
Jean-Michel

On Thu, Feb 09, 2006 at 04:07:31PM +0100, Gianmaria Sannino wrote:
> 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
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support