The output of shtns._analys is a 1 d array

I am using the SHTNS python interface - getvrtdivspec. I am passing in 2d arrays (u and v). They are both functions of latitude and longitude. I should expect vrt and divergence to be also 2d arrays. But when I print the shape of vrtspec and divspec I get a 1 d array.

‌

My nlats = 73 and nlons = 144

‌

The shape of vrtspec is 1225. It should be 73 and 144. It maybe a usage question. Where am I going wrong ?

Comments (2)