Did you get this worked out? My first inclination is that the problem might
be that the hgt array doesn't have lons as the "right-most dimension"
(which is needed according to the docs
In other words, hgt should have a shape of (M, N) where N is the size of
On Tue, Jun 21, 2016 at 5:06 AM, ashwinD12 . <winash12@xxxxxxxxx> wrote:
> I am just starting to learn how to plot NCEP reanalysis data
> using Python3 and matplotlib and basemap. I have some clarifications on how
> to go about doing it.
> My dataset is between latitude 5 N to 40 N and 65 E and 100 E.
> I am trying to plot geopotential height from NCEP reanalysis data. Here is
> my code -
> level = 5
> time = 0
> meters_per_grid = 277830
> d = Dataset("hgt_500_2014_12_5_00Z.nc")
> hgt = d.get_variables_by_attributes(axis="hgt")
> lon = d.get_variables_by_attributes(axis="lon")
> lat =d.get_variables_by_attributes(axis="lat")
> hgt = hgt[::-1]
> hgt, lon = shiftgrid(180, hgt, lon, start = False
> Why do I get an IndexError: index out of range on shiftgrid call ?
> NOTE: All exchanges posted to Unidata maintained email lists are
> recorded in the Unidata inquiry tracking system and made publicly
> available through the web. Users who post to any of the lists we
> maintain are reminded to remove any personal information that they
> do not want to be made public.
> python-users mailing list
> For list information, to unsubscribe, or change your membership options,
> visit: http://www.unidata.ucar.edu/mailing_lists/
Ryan May, Ph.D.