DamienIrving
commented
9 years ago The problem was the windspharm requires 3 things:
- A global field (because of the spherical harmonics it uses)
- Latitude and longitude to be the leading axis
- Latitude from north to south
The cdms2 implementation of windspharm seemed to be able to deal if lat and lon weren't the leading axes, but not if latitude went from south to north. I switched to using the standard implementation (where you have to use their convenience functions to switch axes etc) and it all worked fine.
When I calculate the streamfunction zonal anomaly the data seems erroneous.
In the meantime I'm using the geopotential height zonal anomaly instead, but in the long run I might want to get to the bottom of this issue.