[MITgcm-support] obcs and seaice

[Tom Haine]

Hi Jinlun, 

OK, I'll try setting ice concentration and thickness at the obcs
following the salinity obc code and see what happens. I'll let you know
how it goes. Your thickness analyses look interesting. Thanks for the
link.

Cheers, Tom.

On Thu, 2005-06-09 at 11:42 -0700, Jinlun Zhang wrote:
> Hi Tom,
> 
> One way to handle sea ice open boundary is to set ice strength (P) to 
> zero at the open boundary. This has been quite successful for 
> simulations over the whole Arctic that has limited size of open 
> boundaries. I have not implemented this treatment in MITgcm yet.  In 
> your case, ice inflow is a major source, and the treatment may not be 
> appropriate. For ice dynamics, the Von Neumann obc should be used, which 
> needs some modifications to the LSR dynamics solver. But  we should 
> certainly try things like those used for ptracer obc.
> By the way, if you want to specify ice thickness at your ob, you might 
> want to use the ice thickness fields I created. The fields were created 
> by assimilating SSM/I ice concentration and have a very good match with 
> submarine observations of ice thickness. The fields from 1979 to 2004 
> can be found on http://psc.apl.washington.edu/IDAO.
> 
> Jinlun
> 
> 
> Tom Haine wrote:
> 
> >Hi Folks,
> >
> >I need to allow sea ice to be imported into my Denmark Strait model
> >across open boundaries. I see that seaice is not implemented in obcs
> >currently. Martin Losch has made some suggestions about implementing
> >ptracer in obcs (http://forge.csail.mit.edu/pipermail/mitgcm-
> >support/2004-November/002715.html). Can I simply customize this idea for
> >seaice variables?  Are there other issues I should worry about?
> >
> >Thanks, Tom.
> >  
> >
>