[MITgcm-devel] seaice growth

[Ian Fenty]

Matt,

If the LSR or EVP solvers can't handle ice which has "piled-up" during 
dynamical convergence and is too thick then I suggest that you hack the 
dynamics routines and leave seaice_growth alone!  Limiting thermodynamic 
growth when HEFF > HEFF_MAX probably isn't going to solve your problem 
and arbitrarily capping thickness (  HEFF = MIN(HEFF_MAX,HEFF)  ) in 
seaice_growth will lead to violations of energy and water conservation.

"Hacking" the dynamical solvers could be as simple as making sure that 
the ice strength function returns a sufficiently large value for "thick" 
ice.

Roughly, if: ice_strength(HEFF) \approx \infinity when HEFF >= HEFF_MAX, 
then your "piling-up" problems will be over.

I'm sure Martin can guide you in the details of how to implement it.  I 
believe you would be modifying the 'P' function.

-Ian


On 7/27/2011 4:50 PM, Matthew Mazloff wrote:
> Hi Jean-Michel,
>
> Ooops, sorry, missed it at the top there!
>
> OK, so I see now that #define SEAICE_CAP_HEFF in SEAICE_OPTIONS did 
> nothing because capping HEFF is only allowed in LEGACY code.
>
> Are there any "hacks" that I can use in the new sea ice code to keep 
> the model well conditioned and running.  I would like to diagnose how 
> the ice is piling up, and not having it crash as soon as it starts to 
> pile up would be nice
>
> Thanks,
> Matt
>
>
>
>
> On Jul 27, 2011, at 9:23 AM, Jean-Michel Campin wrote:
>
>> Hi Matt,
>>
>> I don't know which version of seaice_growth.F you are looking at,
>> since:
>>> #include "SEAICE_OPTIONS.h" ?
>> has been there from the beginning (seaice_growth.F version 1.1)
>>
>> Cheers,
>> Jean-Michel
>>
>> On Wed, Jul 27, 2011 at 09:16:27AM -0700, Matthew Mazloff wrote:
>>> Hello,
>>>
>>> Shouldn't seaice_growth.F include the header file:
>>>
>>> #include "SEAICE_OPTIONS.h" ?
>>>
>>> Thanks
>>> Matt
>>>
>>>
>>> ps> This comes up because I am having issues of too much ice
>>> accumulation in the Southern Ocean with the latest seaice package.
>>> It may very well be due to forcing or ocean state, not a seaice
>>> model problem.  When the thickness gets too extreme, however, the
>>> LSR solver crashes.  So I can make a little progress diagnosing
>>> this, I want to set
>>> #define SEAICE_CAP_HEFF
>>> and felt it should go in
>>> SEAICE_OPTIONS.h
>>> but realized that it cannot go there...thus the question
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> 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