A "CubeSphere object" is simply an xarray DataSet with a specific format:
area (tile, y, x)
Below is an example from the GFDL-FV3 output.
With coordinate information available, we can plot the data with one function, without regridding to lat-lon.
plotCS_quick( ) simply calls plt.colormesh( ) 6 times. It only works with the PlateCarree projection and does not show the cells near map boundaries. But it is very fast and can be used for rough checking.
Note: For the rest of this notebook we use the wrapper function plotCS( ) for convenience. Will need to decide which method should be default based on both appearance and performance.
All details including colormap, colorbar and gridline can be tuned as usual. There's nothing special with the cubesphere plot.
xarray's multi-dimensional groupby allows us to take zonal/meridonal average without regridding to lat-lon. This is much faster and more accurate than first regridding to lat-lon and then taking average.
If the "area" coordinate variable is present, the binning will be weighted by area; otherwise uniform weights will be used.
We can set up subplots as usual -- nothing special with cubed-sphere. However, as Cartopy documentation is not very comprehensive (especially about subplots), we make some examples here.