[MITgcm-devel] free drift ice model

[Dimitris Menemenlis]

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

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720