[MITgcm-support] Calculating curl on a cube-sphere (calc_vort_cs.m)

[Shatwell, Peter A]

Hi all,

I am using a coupled configuration of the MITgcm with a cube-sphere grid (C24), and I am trying to calculate the curl of the surface wind-stress in order to understand Ekman layer dynamics in the model ocean.

I found the script calc_vort_cs.m in the MATLAB MITgcm utils and provide it with surface wind stress fields, however I’m unsure how to interpret the outputs, which come in ‘compressed form’ (nc*6*nc+2, 1) or ‘face-splitted form’ (nc+1, nc+1, 1, 6). I’d like to be able to get them into a familiar form of (6*nc, nc, 1), which I can then diagnose with other utils.

Is anyone familiar with converting between these forms for cube-sphere grid data?

Thanks and all the best,

Peter

--
Peter Shatwell
PhD student at MPECDT
Space and Atmospheric Physics
Imperial College London
Huxley 714
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200519/b8bfbfc2/attachment.html>