Has anyone written code to calculate the z-component of the vorticity using gcmfaces objects? I know this has been done for the cubed sphere in /MITgcm/utils/matlab/cs_grid/calc_vort_cs.m. Specifically, I'm using the lat-lon-cap grid, but using gcmfaces, one can write a generic program. thanks, Martha