[MITgcm-devel] seaice growth

[Matthew Mazloff]

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