Bug in absolute vorticity? about windspharm HOT 5 CLOSED

It looks like prep_data is adding a singleton dimension to the end of the 2-d arrays. This dimension is accounted for when creating the planetary vorticity but is lost by the spherepack routine that calculates relative vorticity, thus the shape mismatch.

I'll look into this and try and work out what the best solution is. For now you can work-around this issue by squeezing the singleton dimension from the arrays returned by prep_data before inputting them into the VectorWind instance (the recovery will still work), for example:

up, uinfo = prep_data(u, 'yx')
vp, vinfo = prep_data(v, 'yx')
up = up.squeeze()
vp = vp.squeeze()
wind = VectorWind(up, vp)

from windspharm.

Thanks for the suggested work-around - it solved the problem.

from windspharm.

I think this is actually a bug in pyspharm. I'm going to look into getting it fixed there instead, I'll update on progress.

from windspharm.

The root cause of this is now fixed in pyspharm svn revision r95 and will therefore be in the next release of pyspharm. In the meantime I think it would be wise to patch windspharm to workaround the error.

from windspharm.

I've decided not to include a workaround for this in windspharm. Jeff Whitaker has already done a pyspharm v1.0.8 release that contains the fix for this issue. I've added a note to the documentation suggesting users upgrade to pyspharm v1.0.8 in order to avoid this bug.

from windspharm.