[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