[MITgcm-support] obcs_sponge.F, line 330
David Hebert
david.hebert at nrlssc.navy.mil
Tue Apr 13 11:38:43 EDT 2010
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
More information about the MITgcm-support
mailing list