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

[Menemenlis, Dimitris (3248)]

Martin I am not familiar with embm.  Is this another energy balance atmosphere model?

Yes switching to thsice is also on my nodo (er, I mean todo) list.

Initially, best to check in ifdefed and off by default until it's been vetted a little bit, no?

Dimitris

On Jun 23, 2010, at 6:50 AM, Martin Losch wrote:

> 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