NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.

To learn about what's going on, see About the Archive Site.

[python-users] Need clarification on shiftgrid and addcylic

Hello,
          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 -

<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
</code>

Why do I get an IndexError: index out of range on shiftgrid call ?

BR,
Ashwin.
  • 2016 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the python-users archives: