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.
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.
Ryan,
I was able to fix and I am able to see a wonderful plot using
matplotlib and cartopy ! I had a older version of cartopy and when I
upgraded to the latest version all errors were resolved and a beautiful
plot was generated.
Best,
Ashwin.
On Fri, Oct 14, 2016 at 1:16 PM, ashwinD12 . <winash12@xxxxxxxxx> wrote:
> 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
>>>
>>
>>
>
python-users archives: