[MITgcm-support] Velocity cubed sphere + pickup seaice
Dimitris Menemenlis
menemenlis at jpl.nasa.gov
Tue Jun 3 10:00:52 EDT 2008
> Is there a matlab script available to calculate the total velocity
> field from the zonal and meridional velocity components on a cubed
> sphere? I know the quiver-command, but I am not sure that this
> option gives the correct results on the cubed sphere.
velocity is output on c-grid: west face for u and south face for v, so
quiver will not be exactly correct. are you trying to compute
sqrt(u^2+v^2)? the best way to proceed is to interpolate to tracer
points before doing this computation, i.e., (u(i,j)+u(i+1,j))/2 and
(v(i,j)+v(i,j+1))/2
> Secondly, I am using the saeaice package. If I want to do a restart,
> I change nIter0 (in data) to the required value. Then, the standard
> pick-up files are read. Do I have to take an additional step to pick
> up the pickup_seaice files, too?
when you change nitero0 in data, the correct pickup_seaice files will
be used if useSeaice is .TRUE.
No additional step is required.
More information about the MITgcm-support
mailing list