[MITgcm-support] Processing and visualising cubed sphere output

[Flynn Ames]

Dear MITgcm community,

Hello - I've been attempting to use the cubed sphere topology in the MITgcm. I've managed to successfully run the model and produce output. However, I am having trouble processing the output for visualisation, namely with the velocity outputs.

To cut to the chase, I was wondering if anyone could help with the following questions:

  *
What are the steps to be taken to process outputs for ocean velocity (e.g., zonal, meridional components) so that they can be visualised?
  *
Do translations for the existing MATLAB MITgcmutils cubed sphere library exist for Python?

As I understand it, velocity outputs need to first be re-orientated, because as outputted, they are orientated relative to the cube face that they are on. I have tried using the function 'rotate_uv2uvEN' (from MITgcmutils) for this. However, I find the function requires the velocity arrays to be a different shape than that outputted by the model. I'm therefore unsure if there are some preprocessing steps I need to take before this? Or if I have the wrong idea entirely?

Secondly, I noticed that the python library for cubed sphere (in MITgcmutils) has one function which is to plot a field (which in contrast to the MATLAB function I tried, did accept the velocity arrays in the shape outputted by the model). I'm curious to ask if any further python translations exist for the rest of the cubed sphere library?

If any more information is required about my setup/intentions, I can definitely provide it!

Thanks in advance!
Flynn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20241009/0b25209d/attachment.html>