[MITgcm-support] obcs_sponge.F, line 330
David Hebert
david.hebert at nrlssc.navy.mil
Tue Apr 13 14:03:28 EDT 2010
Thanks Jean-Michel.
The rest of the sponge code appears to be consistent, which is a good
start. We've been trying to muddle through the relaxation settings as
well to get a better understanding of the sponge as well. I have
attached a little one page summary I wrote that multiplies lambda and
urelax to see the resulting addition to the tendency term (gU in this
case). The focus is on U but similar results would be obtained for other
boundaries. It seems that the larger the relaxation factors (tau in my
writeup), the smaller gU is affected. Comments are welcome.
David
On 04/13/10 11:15, Jean-Michel Campin wrote:
> 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
>>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
--
David A. Hebert
Naval Research Lab
Stennis Space Center, MS 39529
david.hebert at nrlssc.navy.mil
Phone: (228) 688-5846
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100413/5fff1ce0/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MITSponge.pdf
Type: application/pdf
Size: 32860 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100413/5fff1ce0/attachment-0001.pdf>
More information about the MITgcm-support
mailing list