[MITgcm-support] seaice_salinity

[Patrick Heimbach]

On Nov 1, 2008, at 7:13 PM, Matthew Mazloff wrote:

> Hi Patrick,
>
> Yes you have it correct.  The MITgcm_contrib/high_res_cube/ 
> cs510_adjoint/ settings work fine without dynamics, and model blows  
> up when I turn on dynamics.
>
> So I can turn on ATMOSPHERIC_LOADING (any preference on what header  
> file it goes in, and is there a runtime flag I need to turn on too)  
> and try without adjustments.  Any other advice?  Are there  
> parameters I can play with to stabilize the solution?

Should go into CPP_OPTIONS.h
(but it is already there if you're using the one from high_res_cube/ 
cs510_adjoint/ ).

You could try increase solver accuracy initially
to make sure you have converged solution:
LSR_ERROR          = 1.e-8,

Other parameters to play with re. dynamics are
SEAICE_strength
SEAICE_waterDrag
SEAICE_drag

but not sure, what stabiling values ought to be.

I would really try a pure NCEP run first, and with high solver accuracy.

-p.



> And by the way, the blow up is gradual, the model runs for about 15  
> days before crashing....I'll look more into where the instability  
> is occurring
>
> -Matt
>
>
>
> On Nov 1, 2008, at 3:47 PM, Patrick Heimbach wrote:
>
>> If easy, can you re-try with un-adjusted atmos. state?
>> THere's a possibility that the adjusted fields you have
>> reflect some compensation for lack of seaiceDynamics in prior runs.
>>
>> Another thing is use of ATMOSPHERIC_LOADING
>> to include sea ice loading on the surface.
>> But if I understand you correctly the
>> MITgcm_contrib/high_res_cube/cs510_adjoint/
>> settings work fine without dynamics, and model only blows with  
>> dynamics(?)
>>
>> -p.
>>
>>
>>
>> On Nov 1, 2008, at 6:34 PM, Matthew Mazloff wrote:
>>
>>> Hi Patrick,
>>>
>>> I am using the adjusted atmos. state and the same  
>>> SEAICE_OPTIONS.h and data.seaice as in the cube sphere adjoint.
>>> http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm_contrib/ 
>>> high_res_cube/cs510_adjoint/
>>>
>>> -matt
>>>
>>>
>>>
>>> On Nov 1, 2008, at 2:52 PM, Patrick Heimbach wrote:
>>>
>>>>
>>>> Matt,
>>>> is this using adjusted atmos. state, or pure NCEP forcings?
>>>> Also, are you using B-grid or C-grid,
>>>> and which of the two solvers, LSR or EVP?
>>>> -p.
>>>>
>>>> On Nov 1, 2008, at 4:07 PM, Matthew Mazloff wrote:
>>>>
>>>>> sorry, it wasn't the seaice_salinity....i had seaice dynamics  
>>>>> on...meant to turn them off at runtime.
>>>>>
>>>>> So, new question:  using seaice dynamics makes my model blow  
>>>>> up.  Any advice on how to increase stability when using seaice  
>>>>> dynamics?
>>>>>
>>>>> -matt
>>>>>
>>>>>
>>>>>
>>>>> On Oct 31, 2008, at 12:56 PM, Matthew Mazloff wrote:
>>>>>
>>>>>> Hello,
>>>>>>
>>>>>> My model crashes when I  #define SEAICE_SALINITY
>>>>>>
>>>>>> I just turned it on to keep up with ECCO2...but is this a flag  
>>>>>> that is even worth using.  And if so, does anyone have any  
>>>>>> advice (besides to turn it back off)?
>>>>>>
>>>>>> Thanks
>>>>>> Matt
>>>>>>
>>>>>> ps> and could someone also let me know an expected value of  
>>>>>> seaice_hsalt_max so I know if I have initialized it properly
>>>>>>
>>>>>> _______________________________________________
>>>>>> 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
>>>>
>>>> ---
>>>> 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
>>>>
>>>>
>>>> _______________________________________________
>>>> 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
>>
>> ---
>> 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
>>
>>
>> _______________________________________________
>> 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

---
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