[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