Ryan,
I was able to fix those errors. Essentially I forgot to add these
three lines of errors and the problem goes away. But another problem
arises with Tkinter errors and the matplotlib interaction.
temp = temp[:].squeeze()
u_wind = u_wind[:].squeeze() * units('m/s')
v_wind = v_wind[:].squeeze() * units('m/s')
I am not sure but maybe the code needs to include
matplotlib.use('TkAgg'). It doesn't solve the lack of display though.
I get a blank window which when I close the long trace below arises.
The dependencies are quite involved. In order for this code to run you
need to install apart from cartopy the other two Python
packages - Shapely and pyshp.
File "/usr/lib/python3/dist-packages/matplotlib/backends/backend_agg.py",
line 474, in draw
self.figure.draw(self.renderer)File
"/usr/lib/python3/dist-packages/matplotlib/artist.py", line 61, in
draw_wrapper
draw(artist, renderer, *args, **kwargs)File
"/usr/lib/python3/dist-packages/matplotlib/figure.py", line 1159, in
draw
func(*args)
File "/usr/lib/python3/dist-packages/matplotlib/artist.py", line 61,
in draw_wrapper
draw(artist, renderer, *args, **kwargs)
File "/usr/local/lib/python3.5/dist-packages/cartopy/mpl/geoaxes.py",
line 355, in draw
inframe=inframe)File
"/usr/lib/python3/dist-packages/matplotlib/artist.py", line 61, in
draw_wrapper
draw(artist, renderer, *args, **kwargs)File
"/usr/lib/python3/dist-packages/matplotlib/axes/_base.py", line 2324,
in draw
a.draw(renderer)File
"/usr/lib/python3/dist-packages/matplotlib/artist.py", line 61, in
draw_wrapper
draw(artist, renderer, *args, **kwargs)
File "/usr/local/lib/python3.5/dist-packages/cartopy/mpl/feature_artist.py",
line 103, in draw
geoms = self._feature.intersecting_geometries(extent)File
"/usr/local/lib/python3.5/dist-packages/cartopy/feature.py", line 113,
in intersecting_geometriesreturn (geom for geom in self.geometries()
if
File "/usr/local/lib/python3.5/dist-packages/cartopy/feature.py",
line 180, in geometries
name=self.name)
File "/usr/local/lib/python3.5/dist-packages/cartopy/io/shapereader.py",
line 246, in natural_earthreturn ne_downloader.path(format_dict)File
"/usr/local/lib/python3.5/dist-packages/cartopy/io/__init__.py", line
223, in path
result_path = self.acquire_resource(target_path, format_dict)File
"/usr/local/lib/python3.5/dist-packages/cartopy/io/shapereader.py",
line 301, in acquire_resource
shapefile_online = self._urlopen(url)File
"/usr/local/lib/python3.5/dist-packages/cartopy/io/__init__.py", line
262, in _urlopenreturn urlopen(url)File
"/usr/lib/python3.5/urllib/request.py", line 163, in urlopenreturn
opener.open(url, data, timeout)File
"/usr/lib/python3.5/urllib/request.py", line 472, in open
response = meth(req, response)File
"/usr/lib/python3.5/urllib/request.py", line 582, in
http_response'http', request, response, code, msg, hdrs)File
"/usr/lib/python3.5/urllib/request.py", line 504, in error
result = self._call_chain(*args)File
"/usr/lib/python3.5/urllib/request.py", line 444, in _call_chain
result = func(*args)File "/usr/lib/python3.5/urllib/request.py", line
696, in http_error_302return self.parent.open(new,
timeout=req.timeout)File "/usr/lib/python3.5/urllib/request.py", line
472, in open
response = meth(req, response)
File "/usr/lib/python3.5/urllib/request.py", line 582, in
http_response'http', request, response, code, msg, hdrs)
File "/usr/lib/python3.5/urllib/request.py", line 510, in errorreturn
self._call_chain(*args)File "/usr/lib/python3.5/urllib/request.py",
line 444, in _call_chain
result = func(*args)
File "/usr/lib/python3.5/urllib/request.py", line 590, in http_error_default
raise HTTPError(req.full_url, code, msg, hdrs, fp)
urllib.error.HTTPError: HTTP Error 404: Not FoundTraceback (most
recent call last):
File "tempAdv.py", line 72, in <module>
plt.show()
File "/usr/lib/python3/dist-packages/matplotlib/pyplot.py", line 244,
in showreturn _show(*args, **kw)
File "/usr/lib/python3/dist-packages/matplotlib/backend_bases.py",
line 192, in __call__
self.mainloop()
File "/usr/lib/python3/dist-packages/matplotlib/backends/backend_tkagg.py",
line 74, in mainloop
Tk.mainloop()
File "/usr/lib/python3.5/tkinter/__init__.py", line 408, in mainloop
_default_root.tk.mainloop(n)
AttributeError: 'NoneType' object has no attribute 'tk
On Fri, Oct 14, 2016 at 7:39 AM, ashwinD12 . <winash12@xxxxxxxxx> wrote:
> Ryan,
> Those files are just reanalysis 2 data files that are downloaded
> from NCEP NCAR. The shape of the three data sets is identical
>
> (1,1,73,144). If you wish I can send you those files offline.
> Regards,
> Ashwin.
>
> On Fri, Oct 14, 2016 at 3:47 AM, Ryan May <rmay@xxxxxxxx> wrote:
>
>> Ashwin,
>>
>> I'm responsible chiefly for MetPy so that's fine. For future reference,
>> though, I prefer public channels in general, either Unidata's python-users
>> email list or issues over at http://github.com/metpy/MetPy/issues so
>> that the rest of the community can chime in with answers or benefit from
>> the information. I'm cc'ing the python-users list so we can take the
>> conversation there.
>>
>> That error message implies that something is off with the size of one of
>> your arrays. What are the sizes of temp, u_wind, and v_wind?
>>
>> Ryan
>>
>>
>>
>> On Thu, Oct 13, 2016 at 9:53 AM, ashwinD12 . <winash12@xxxxxxxxx> wrote:
>>
>>> Hello Ryan,
>>> I am not sure if I should be contacting you
>>> regarding MetPy errors. But I ran a slightly modified example of this
>>> cookbook -http://nbviewer.jupyter.org/github/unidata/notebook-gallery
>>> /blob/master/notebooks/850mb_Temperature_Advection.ipynb and I am
>>> getting this runtime error -
>>> Any idea what that means ?
>>> Traceback (most recent call last):
>>> File "tempAdv.py", line 38, in <module>
>>> adv = advection(temp * units.kelvin, [u_wind,
>>> v_wind],(42583.44283066879 * units.meter, 55588.7366760194 * units.meter))
>>> * units('degF/sec')
>>> File
>>> "/usr/local/lib/python3.5/dist-packages/MetPy-0.3.1-py3.5.egg/metpy/calc/kinematics.py",
>>> line 166, in advection
>>> grad = _stack(_gradient(scalar, *deltas))
>>> File
>>> "/usr/local/lib/python3.5/dist-packages/MetPy-0.3.1-py3.5.egg/metpy/calc/kinematics.py",
>>> line 18, in _gradient
>>> grad = np.gradient(f, *args, **kwargs)
>>> File "/usr/lib/python3/dist-packages/numpy/lib/function_base.py",
>>> line 1440, in gradient
>>> "Shape of array too small to calculate a numerical gradient, "
>>> ValueError: Shape of array too small to calculate a numerical gradient,
>>> at least two elements are required.
>>>
>>>
>>>
>>
>>
>> --
>> Ryan May, Ph.D.
>> Software Engineer
>> UCAR/Unidata
>> Boulder, CO
>>
>
>
