[MITgcm-support] obcs and seaice

[Jinlun Zhang]

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

-- 

Jinlun Zhang
Polar Science Center, Applied Physics Laboratory
University of Washington, 1013 NE 40th St, Seattle, WA 98105-6698

Phone: (206)-543-5569    Fax:    (206)-616-3142

zhang at apl.washington.edu
http://psc.apl.washington.edu/IDAO