[MITgcm-support] seaice setup for idealized ice leads experiment

[Mehmet Ilicak]

Hi Martin,

Thanks for your suggestions.

First of all, hacking the code didn't work. When I set the 
stefan-boltzmann constant to 0,
the ice melts and I have negative salinity bias at the top of the model 
domain.

So I decided to prescribe downward long-wave radiation.
I remove the lead and put the ice everywhere, so that I can check the 
results easily.
I assume Tice=-2C and compute stefanBoltzmann*Tice^4 and run the simulation.
It worked so far.

Now I have two questions since I don't have any experience in ice 
modelling.
First, is it proper to assume Tice=-2?
Second, if I run the simulation with lead, should I prescribe the LW 
radiation over the lead also?
Because if I do it, that will reduce the forcing and less water will sink.

Best,
Mehmet


On 09/17/2012 06:00 PM, mitgcm-support-request at mitgcm.org wrote:
> Message: 1
> Date: Mon, 17 Sep 2012 11:57:51 +0200
> From: Martin Losch<Martin.Losch at awi.de>
> To:mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] seaice setup for idealized ice leads
> 	experiment
> Message-ID:<239DE773-EAC5-4308-BFF7-B51663307E7C at awi.de>
> Content-Type: text/plain; CHARSET=US-ASCII
>
> Mehmet,
>
> I think you need to specify some radiation fields, at least non-zero downward longwave radiation to balance the upward radiation by the ice (stefanBoltzmann*Tice**4). I think (without checking your set up very closely) you see ice formation because of this strong imbalance.
>
> Martin