[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