[MITgcm-devel] seaice code beyond checkoint61j blows for ECCO-GODAE

[Patrick Heimbach]

HI there,

sorry for recent silence.
I haven't had time to do more tests, runs.
Also loosing track of things
(e.g. working vs. non-working checkpoints and non-checkpoints).
F.Y.I., I've appended my data.seaice

Cheers
-p.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.gmredi
Type: application/octet-stream
Size: 1197 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20090522/32181117/attachment.obj>
-------------- next part --------------



On May 20, 2009, at 5:35 AM, Martin Losch wrote:

> Hi all,
>
> I am not sure what to do about the problems with the seaice code.  
> Currently there is potentially (most probably) broken code in the  
> repository. I have a suspicion what it is (free-slip boundary  
> conditions), but I cannot reproduce the problems that have been  
> reported (especially not the ECCO-GODAE). Also I cannot narrow down  
> the problem further, if I don't get more details.
>
> It seems that Patrick is happy with the running system with seaice  
> reverted to ckpt61j (which I can understand), but I am very  
> uncomfortable with the current code and would like to fix it.
>
> Should I go ahead and implement/check-in the fix that has worked  
> for Dimitris (setting etaMeanZ=0 on boundaries in case of free- 
> slip)? Or should I wait until more diagnostics are reported?
>
> Martin
>
> On May 14, 2009, at 8:28 AM, Martin Losch wrote:
>
>> Hi Patrick,
>>
>> do you have any further diagnosis for the seaice problem?  
>> SEAICE_no_slip = true or false (default is false)?
>>
>> I have a fix for the free slip case, that seems to work for  
>> Dimitris, and if it works for you I can check it in (I probably  
>> should), but first I'd like to hear from you about this.
>>
>> Martin
>>
>> On May 7, 2009, at 8:26 PM, Dimitris Menemenlis wrote:
>>
>>> That's correct.  In all the tests that we have done here, dozens  
>>> of years of integration on the Arctic, Weddell, CS510, and CS510- 
>>> adjoint domains, we have not experienced a single failure when:
>>>
>>> i) SEAICE_no_slip=.true.
>>>
>>> ii) SEAICE_no_slip=.false. and the etaMeanZ mask is applied.
>>>
>>> Patrick, what is your experience with the ECCO-GODAE set up?
>>>
>>> D.
>>>
>>> On May 6, 2009, at 11:16 PM, Martin Losch wrote:
>>>
>>>> I would really like to get to the bottom of the blow-ups. As far  
>>>> as I
>>>> can see they are all with free-slip boundary conditons
>>>> (SEAICE_no_slip=.false.). Correct? (Patrick?)
>>>> If that's the case, then I would change the implementation of  
>>>> the free
>>>> slip BC (include the masking of etaMeanZ).
>>>
>>> Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
>>> Jet Propulsion Lab, California Institute of Technology
>>> MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099, USA
>>> tel: 818-354-1656;  cell: 818-625-6498;  fax: 818-393-6720
>>>
>>> _______________________________________________
>>> MITgcm-devel mailing list
>>> MITgcm-devel at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>>
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel

---
Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach