Due to the current gap in continued funding from the U.S. National Science Foundation (NSF), the NSF Unidata Program Center has temporarily paused most operations. See NSF Unidata Pause in Most Operations for details.

[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: