[MITgcm-support] obcs_sponge.F, line 330

[Jean-Michel Campin]

Hi David,

Thanks for the report. You are right, and this is clearly a bug. 
Will check-in this fix into the CVS repository.
However, I have to admit that I don't fully understand 
this relaxobcsinner business in obcs_sponge.F, so it's difficult 
for me to tell if it will be bug free once this one is fixed.

Cheers,
Jean-Michel

On Tue, Apr 13, 2010 at 10:38:43AM -0500, David Hebert wrote:
> In combing through obcs_sponge.F, in obcs_sponge_v, for the eastern  
> boundary near line 330, has the following
>
> 327     IF ((i.ge.imin).and.(i.le.imax)) THEN
> 328          vrelax=(
> 329 &      float(spongeThickness-isl)*OBEv(j,kLev,bi,bj)
> 330 &    + float(isl-1)*vVel(i,j,kLev,bi,bj) )
> 331 &    / float(spongeThickness)
>
> Should that (isl-1) be (isl)? This is the only place in obcs_sponge.F  
> where isl-1 (or jsl-1) appears.
>
> Thanks,
>
> David
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support