[MITgcm-devel] free drift ice model

[Jinlun Zhang]

If it is LSR, there is a third level solution that is supposed to make 
Cori. term implicit.

Jinlun

Dimitris Menemenlis wrote:

>> My model dt = 900;  Is  SEAICE_deltaTdyn = 3600 feasible.
>
>
> SEAICE_deltaTdyn = 3600 is sufficient to properly resolve Coriolis 
> term, so I
> think yes ... but of course with pkg/seaice nothing can ever be taken 
> for granted until you actually try it.
>
>> How do you find the seaice performance burden in the high-res cube 
>> results?
>
>
> about 15% of wallclock time, if you include model diagnostics
>
>> I would love to investigate how the ice moves wrt to ocean in the 
>> high-res
>> cube results.  Can you point me towards some files I could analyze?
>
>
> Matt, some results are available at 
> ftp://ecco2.jpl.nasa.gov/data2/cube/cube37.
> Descriptions of format and how to read is in 
> http://ecco2.jpl.nasa.gov/products/output/cube/output.html
> I assume that you need UICE, VICE, UVEL, VVEL, which are on the ftp 
> server. U10m and V10m are from ncep_rgau but unfortunately they are 
> not output on the rotated cubed sphere grid.  So it's some extra work 
> to diagnose what you need. Maybe easiest to go with ECCO-GODAE fields 
> after all or better yet to experiment in cutting down cost of dynamic 
> solver as per previous discussion.
>
> D.