[MITgcm-devel] [MITgcm-support] maintaining sea ice coverage over multi-year run on cs32

[Martin Losch]

Hi Dimitris,

I do have this code around, but it is in a version which is taylored to the embm that I added. I had to introduce it because the ice model kept creating water (ice) volume, because the sublimation part is accounted for in the latent heat flux (seaice_budget_ice) but not in the volume part. In this context, I realized again, how poorly written this seaice_growth routine is. It's such a drag to try and change anything in there.

I really recommend switching to the thsice package (which is even harder to understand, but probably has a much lower problem-density). We should work on streamlining thsice for the adjoint (and vectorize it) and then use it as the default!

If I get to check the sublimation code for seaice in (I could try this next week), it will change the results, but clearly it's a bug fix. I guess we should have CPP flags around it, have them defined by default and then turned off in verification experiments? Or rather vice versa (you'd have to turn them on yourself, ie. you have to know about this)?

Martin

On Jul 16, 2010, at 2:32 PM, Menemenlis, Dimitris (3248) wrote:

> Martin, you have sublimation code available for pkg/seaice?
> That would be awesome.
> It's been on my todo list for quite some time.
> (I should really call it my "nodo" list because I never seem to get around to it ).
> 
> D.
> 
> On Jun 23, 2010, at 4:00 AM, Martin Losch wrote:
> 
>> I need to check in a contribution of ice sublimation to the mass balance
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel