matplotlib-users Mailing List for matplotlib (Page 79)

I apologize if this has been fixed already, I can only check different
versions at home. However, the documentation of mpl
1.3.1<http://matplotlib.org/1.3.1/api/pyplot_api.html#matplotlib.pyplot.errorbar>.
has the same information. So unless the code changed to reflect the
documentation, this is still present.

When using errorbar, the documentation says the color of the errorbar lines
will match with the color of the markers if ecolor=None. That’s not what I
found. Apparently it takes over the color of the Line2D instance which
interconnects the markers.

Short, Self Contained, Correct Example:

from pylab import *
plt.ion()  # saves typing show

x = np.arange(10)
y = np.random.rand(10)
xerr, yerr = y/4., y/4.

# Markers in red, but errorlines assume the color of the "trendline"
(default rcparams: blue).
errorbar(x, y, yerr=yerr, mfc='r', marker='o', ecolor=None)

# Errorlines get color green now - documentation not in line with results

figure(); errorbar(x, y, yerr=xerr, mfc='r', marker='o', ecolor=None, color='g')
# Errorlines get color blue now, because it can be specified -
expected behaviour
figure(); errorbar(x, y, yerr=xerr, mfc='r', marker='o', ecolor='b', color='g')

Is this an oversight mistake?

Hello everyone,

Just a quick reminder that the EuroScipy call for abstracts closes on the
14th: don't forget to submit your talk proposal! It is in four days only!

In short, EuroScipy is a cross-disciplinary gathering focused on the use
and development of the Python language in scientific research. This event
strives to bring together both users and developers of scientific tools, as
well as academic research and state of the art industry.

EuroSciPy 2014, the Seventh Annual Conference on Python in Science, takes
place in *Cambridge, UK on 27 - 30 August 2014*. The conference features
two days of tutorials followed by two days of scientific talks. The day
after the main conference, developer sprints will be organized on projects
of interest to attendees.

The topics presented at EuroSciPy are very diverse, with a focus on
advanced software engineering and original uses of Python and its
scientific libraries, either in theoretical or experimental research, from
both academia and the industry. The program includes keynotes, contributed
talks and posters.

Submissions for talks and posters are welcome on our website (
http://www.euroscipy.org/2014/). In your abstract, please provide details
on what Python tools are being employed, and how. The deadline for
submission is 14 April 2014.

Also until 14 April 2014, you can apply for a sprint session on 31 August
2014. See https://www.euroscipy.org/2014/calls/sprints/ for details.

Thanks,

N

 I want to plot audio data (samples versus time), and then animate a marker
moving across the plot while the audio is played back. The marker should
stay in sync with the audio data that is being played at that moment. Is
there a way to do this using the animation capabilities of matplotlib? If I
could sync the movement of the market with an accurate timer (like a QTimer
in Qt) it should be good enough.

Mike

Thanks Paul, I will try it out.


On Wed, Apr 9, 2014 at 12:21 PM, Paul Hobson <pmh...@gm...> wrote:

>
>
>
> On Wed, Apr 9, 2014 at 9:00 AM, Adam Hughes <hug...@gm...>wrote:
>
>> Thanks.  That's probably the way I'll go.  At first, I thought creating
>> separate legend markers and removing them from the plot seemed hacky, but I
>> guess there's no way that matplotlib could know which legend size I want.
>>  I wonder if there'd be any interest in a PR to add a keyword to legend to
>> handle this situation?
>>
>
> Why not just work the other way around with proxy artists. IOW, make the
> artists but never add them to the plot.
>
>
> http://matplotlib.org/users/legend_guide.html?highlight=proxy%20artists#using-proxy-artist
> (works with Line2D artists)
>
> -p
>
>
>
>>
>> On Wed, Apr 9, 2014 at 1:44 AM, Sterling Smith <sm...@fu...>wrote:
>>
>>> Adam,
>>>
>>> I haven't investigated, but does the discussion of the legend marker at
>>> [1] help?
>>>
>>> -Sterling
>>>
>>> [1]
>>> https://www.mail-archive.com/mat...@li.../msg25200.html
>>>
>>> On Apr 8, 2014, at 3:44PM, Adam Hughes wrote:
>>>
>>> > Hello,
>>> >
>>> > I've been searching but can't seem to find this topic addressed
>>> (perhaps wrong search terms)
>>> >
>>> > Simply put, I have a scatter plot with variable size markers, and I'd
>>> like to have the markers all be a single size in the legend.  Is there a
>>> standard way to do this?
>>> >
>>> > Thanks.
>>> >
>>> ------------------------------------------------------------------------------
>>> > Put Bad Developers to Shame
>>> > Dominate Development with Jenkins Continuous Integration
>>> > Continuously Automate Build, Test & Deployment
>>> > Start a new project now. Try Jenkins in the cloud.
>>> >
>>> http://p.sf.net/sfu/13600_Cloudbees_______________________________________________
>>> > Matplotlib-users mailing list
>>> > Mat...@li...
>>> > https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>>
>>>
>>
>>
>> ------------------------------------------------------------------------------
>> Put Bad Developers to Shame
>> Dominate Development with Jenkins Continuous Integration
>> Continuously Automate Build, Test & Deployment
>> Start a new project now. Try Jenkins in the cloud.
>> http://p.sf.net/sfu/13600_Cloudbees
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>
>>
>

On Wed, Apr 9, 2014 at 9:00 AM, Adam Hughes <hug...@gm...> wrote:

> Thanks.  That's probably the way I'll go.  At first, I thought creating
> separate legend markers and removing them from the plot seemed hacky, but I
> guess there's no way that matplotlib could know which legend size I want.
>  I wonder if there'd be any interest in a PR to add a keyword to legend to
> handle this situation?
>

Why not just work the other way around with proxy artists. IOW, make the
artists but never add them to the plot.

http://matplotlib.org/users/legend_guide.html?highlight=proxy%20artists#using-proxy-artist
(works with Line2D artists)

-p



>
> On Wed, Apr 9, 2014 at 1:44 AM, Sterling Smith <sm...@fu...>wrote:
>
>> Adam,
>>
>> I haven't investigated, but does the discussion of the legend marker at
>> [1] help?
>>
>> -Sterling
>>
>> [1]
>> https://www.mail-archive.com/mat...@li.../msg25200.html
>>
>> On Apr 8, 2014, at 3:44PM, Adam Hughes wrote:
>>
>> > Hello,
>> >
>> > I've been searching but can't seem to find this topic addressed
>> (perhaps wrong search terms)
>> >
>> > Simply put, I have a scatter plot with variable size markers, and I'd
>> like to have the markers all be a single size in the legend.  Is there a
>> standard way to do this?
>> >
>> > Thanks.
>> >
>> ------------------------------------------------------------------------------
>> > Put Bad Developers to Shame
>> > Dominate Development with Jenkins Continuous Integration
>> > Continuously Automate Build, Test & Deployment
>> > Start a new project now. Try Jenkins in the cloud.
>> >
>> http://p.sf.net/sfu/13600_Cloudbees_______________________________________________
>> > Matplotlib-users mailing list
>> > Mat...@li...
>> > https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>
>>
>
>
> ------------------------------------------------------------------------------
> Put Bad Developers to Shame
> Dominate Development with Jenkins Continuous Integration
> Continuously Automate Build, Test & Deployment
> Start a new project now. Try Jenkins in the cloud.
> http://p.sf.net/sfu/13600_Cloudbees
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Thanks.  That's probably the way I'll go.  At first, I thought creating
separate legend markers and removing them from the plot seemed hacky, but I
guess there's no way that matplotlib could know which legend size I want.
 I wonder if there'd be any interest in a PR to add a keyword to legend to
handle this situation?


On Wed, Apr 9, 2014 at 1:44 AM, Sterling Smith <sm...@fu...>wrote:

> Adam,
>
> I haven't investigated, but does the discussion of the legend marker at
> [1] help?
>
> -Sterling
>
> [1]
> https://www.mail-archive.com/mat...@li.../msg25200.html
>
> On Apr 8, 2014, at 3:44PM, Adam Hughes wrote:
>
> > Hello,
> >
> > I've been searching but can't seem to find this topic addressed (perhaps
> wrong search terms)
> >
> > Simply put, I have a scatter plot with variable size markers, and I'd
> like to have the markers all be a single size in the legend.  Is there a
> standard way to do this?
> >
> > Thanks.
> >
> ------------------------------------------------------------------------------
> > Put Bad Developers to Shame
> > Dominate Development with Jenkins Continuous Integration
> > Continuously Automate Build, Test & Deployment
> > Start a new project now. Try Jenkins in the cloud.
> >
> http://p.sf.net/sfu/13600_Cloudbees_______________________________________________
> > Matplotlib-users mailing list
> > Mat...@li...
> > https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Adam,

I haven't investigated, but does the discussion of the legend marker at [1] help?

-Sterling

[1] https://www.mail-archive.com/mat...@li.../msg25200.html 

On Apr 8, 2014, at 3:44PM, Adam Hughes wrote:

> Hello,
> 
> I've been searching but can't seem to find this topic addressed (perhaps wrong search terms)
> 
> Simply put, I have a scatter plot with variable size markers, and I'd like to have the markers all be a single size in the legend.  Is there a standard way to do this?
> 
> Thanks.
> ------------------------------------------------------------------------------
> Put Bad Developers to Shame
> Dominate Development with Jenkins Continuous Integration
> Continuously Automate Build, Test & Deployment 
> Start a new project now. Try Jenkins in the cloud.
> http://p.sf.net/sfu/13600_Cloudbees_______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Hello,

I've been searching but can't seem to find this topic addressed (perhaps
wrong search terms)

Simply put, I have a scatter plot with variable size markers, and I'd like
to have the markers all be a single size in the legend.  Is there a
standard way to do this?

Thanks.

Which version of matplotlib are you using? There used to be some issues
with this, but I could have sworn I resolved this a while back.


Ben Root


On Mon, Apr 7, 2014 at 2:40 PM, Emilia Petrisor
<emi...@gm...>wrote:

> Dear all,
>
>
> I want to visualize the image through the stereographic projection (that
> maps
> the extended complex plane to the unit sphere), of a circle and a line
> from the plane.
>
> I imported:
>
> import numpy as np
> import matplotlib.pyplot as plt
> from mpl_toolkits.mplot3d import Axes3D
>
> fig = plt.figure()
> ax = fig.add_subplot(111,  projection='3d')
>
> If I plot only the  the sphere and  a circle in the plane xoy and set:
>
> ax.set_aspect('equal')
>
> all is OK, but inserting one more plot for the image of the circle onto
> the sphere,
> the aspect is changed and the sphere is displayed as an ellipsoid.
>
> My question is:  how could I lock the aspect in order to keep the geometry
> of the sphere?
>
> Thanks,
>
> Em
>
>
>
>
>
>
> ------------------------------------------------------------------------------
> Put Bad Developers to Shame
> Dominate Development with Jenkins Continuous Integration
> Continuously Automate Build, Test & Deployment
> Start a new project now. Try Jenkins in the cloud.
> http://p.sf.net/sfu/13600_Cloudbees
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Hi list!

Has anybody ported the example at,

http://matplotlib.org/examples/user_interfaces/mpl_with_glade.html

to the new versions of GTK (3.10) and Glade (3.16) yet?

It would be really helpful to me to see a basic working example.

Thank you all!
-Tobias

Dear all,


I want to visualize the image through the stereographic projection (that
maps
the extended complex plane to the unit sphere), of a circle and a line from
the plane.

I imported:

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

fig = plt.figure()
ax = fig.add_subplot(111,  projection='3d')

If I plot only the  the sphere and  a circle in the plane xoy and set:

ax.set_aspect('equal')

all is OK, but inserting one more plot for the image of the circle onto the
sphere,
the aspect is changed and the sphere is displayed as an ellipsoid.

My question is:  how could I lock the aspect in order to keep the geometry
of the sphere?

Thanks,

Em

dear all, 
I get it.

I have only to set:

cm = matplotlib.cm.get_cmap()
and then I have again the default colormap 




--
View this message in context: http://matplotlib.1069221.n5.nabble.com/colormap-default-tp43196p43197.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Dear all,

I am trying to plot different subplot with different colormap colors.

this is the first one, with colormap in red (YlOrRd):

cm = matplotlib.cm.get_cmap('*YlOrRd*')
sc=ax1.scatter(xp[::N],yp[::N], s=20, c=QA[::N],
marker='o',edgecolor='none',alpha=0.8,
                                                cmap=cm,
                                                vmin=vmin,vmax=vmax) 

How can set up again the default matplotlib colormap for the other. I do not
know the name of it.

thanks you all





--
View this message in context: http://matplotlib.1069221.n5.nabble.com/colormap-default-tp43196.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

I tried to get the resizing working, but to no avail.

Without modifying your code, when you click on multiple squares, and then
resize, all the previously displayed squares are dawn, as if there was a
merge of all the previous blits...

Any idea ?



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Question-about-copy-from-bbox-tp43182p43194.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

This worked brilliantly! Thank you so much!

Secondary question now, is this:
I have a contourf and a polygon made with

ax.add_collection3d(a3.art3d.Poly3DCollection(verts, facecolors=[.5,.5,.5],
linewidths=1,edgecolor='k',zorder=0))

How do I change it so the contourf is on top of the polygon? By default it
seems to plot the polygon on top. I tried zorder but that did nothing.

Ta

Hi

thanks for the help, this is exactly what I needed. I updated the code in my
app and it worked.

Some minor observations though :
- This does not play well with tight_layout. I think tight_layout tries to
readjust the plots even after the singleShot call. So I needed to disable
it.
- I have plots in multiple tabs (in a QtGui.QTabWidget()). When moving from
one tab to another, I update the plots. I had to set the delay to 50 ms (on
faster computers the delay could be made smaller). 0 ms is too short in this
specific case and the copy_from_bbox will copy the background of the
previous tab.
- I have to work on the resizing of the window. I have to catch the
resizeEvent in the plot and update it to get a new fresh background. Don't
know what's the best solution for this.

Thanks very much for the singleShot trick, this was very helpful.

Michka



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Question-about-copy-from-bbox-tp43182p43192.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

So, you would like to plot a surface, and then have some slices through it?
Perhaps you are looking for something like this?
http://matplotlib.org/examples/mplot3d/contourf3d_demo2.html

In your case, I would use the contourf()'s "zdir" argument to set the
direction of the slice's normal vector, and the "offset" argument to set
the location of the slice. In this particular example, the three
contourf()'s were purposely placed on the edges of the axes, but they can
be placed anywhere you want them to be. You will need to supply contourf()
with the appropriate surface data to contour.

A word of caution. mplot3d's rendering is very crude. When you start to
intersect planes together, you will very likely get visual artifacts that I
tend to refer to as "Escher-like" in that they are physically impossible.
If you want true 3d volume rendering, then you may want to take a look at
glumpy or vis3d or mayavi2.

I hope that helps!
Ben Root



On Thu, Apr 3, 2014 at 12:50 AM, james <ja...@gm...> wrote:

> Hi All,
>
> I wish to make a 3d volume with a series of slices through it. I have X-Z
> data at 15 different Y planes. In MATLAB, I would make a pcolor plot, then
> set the Z data to the original ydata, and the Y data to a constant. This
> works fine.
>
> I would like to achieve a similar result in matplotlib, but I cannot figure
> out anyway to do it.
>
> Any help would be appreciated.
>
> James
>
>
>
> --
> View this message in context:
> http://matplotlib.1069221.n5.nabble.com/Set-X-Y-and-Z-data-for-a-pcolormesh-tp43187.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

If I understand your intent well, then you want to make slices through a 3D
volume and show such slice planes.

I’ve always done this with the various functions of an Axes3D instance that
allow you to specify the slice position. On the matplotlib examples page,
there is a good example that showed me how to do what I wanted (even though
I only had a series of 1D data):
polys3d_demo<http://matplotlib.org/1.3.1/examples/mplot3d/polys3d_demo.html>
.

You could use the same technique with contourf, which is probably closer to
what you want. Have a look at this stackoverflow
question<https://stackoverflow.com/questions/15582105/python-plot-stacked-image-slices>.


I hope this helps.


2014-04-03 6:50 GMT+02:00 james <ja...@gm...>:

> Hi All,
>
> I wish to make a 3d volume with a series of slices through it. I have X-Z
> data at 15 different Y planes. In MATLAB, I would make a pcolor plot, then
> set the Z data to the original ydata, and the Y data to a constant. This
> works fine.
>
> I would like to achieve a similar result in matplotlib, but I cannot figure
> out anyway to do it.
>
> Any help would be appreciated.
>
> James
>
>
>
> --
> View this message in context:
> http://matplotlib.1069221.n5.nabble.com/Set-X-Y-and-Z-data-for-a-pcolormesh-tp43187.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

Hi All,

I wish to make a 3d volume with a series of slices through it. I have X-Z
data at 15 different Y planes. In MATLAB, I would make a pcolor plot, then
set the Z data to the original ydata, and the Y data to a constant. This
works fine.

I would like to achieve a similar result in matplotlib, but I cannot figure
out anyway to do it.

Any help would be appreciated.

James

Hi All,

I wish to make a 3d volume with a series of slices through it. I have X-Z
data at 15 different Y planes. In MATLAB, I would make a pcolor plot, then
set the Z data to the original ydata, and the Y data to a constant. This
works fine.

I would like to achieve a similar result in matplotlib, but I cannot figure
out anyway to do it.

Any help would be appreciated.

James



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Set-X-Y-and-Z-data-for-a-pcolormesh-tp43187.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Hello,

I noticed that in the matplotlib source code in lib/matplotlib/rcsetup.py
there is a parameter called axes.formatter.useoffset that defaults to True.
However, in the documentation the example rc file doesn't include this
parameter, and if I try to use it I get the exception "Bad key
axes.formatter.useoffset".

I am constantly using ax.ticklabel_format(useoffset=False) to get at this
parameter and I would love to just have my rc file take care of it instead.
How can I do this? If I can't, then can this be added? Why is the parameter
listed in rcsetup.py but not in the example matplotlibrc?

Thanks,
Julian

I should have posted this question on the numpy list of course, even though
floor is imported with pylab.

Off-list, Felix pointed me to the following discussion where some of the
reasons are explained why a float is returned rather than an integer:

http://stackoverflow.com/questions/8582741/why-do-pythons-math-ceil-and-math-floor-operations-return-floats-instead-of



On Mon, Mar 31, 2014 at 5:05 PM, Mark Bakker <ma...@gm...> wrote:

> I expected that floor(x) would return an integer especially since the docs
> state:
>
> Return the floor of the input, element-wise.
>
>
> The floor of the scalar `x` is the largest integer `i`, such that
>
> `i <= x`. It is often denoted as :math:`\lfloor x \rfloor`.
>
>
> Any reason why it returns a float? Bug/feature?
>
>
> Thanks, Mark
>

Hi Michka,

I haven't practiced PyQt for some time, but I think I remember there is
a common practice of using a "0 ms" timer to launch a function after the
Gui setup.

I've modified your gist here :
https://gist.github.com/pierre-haessig/9909708
(for some reason the Github fork button printed "The change you wanted
was rejected.")

I don't know if it does what you want, but at least it doesn't move anymore.

best,
Pierre

Hi

I am using blitting and am not able to get it to work as I want.

I am using copy_from_bbox during the setup of my PyQt GUI to save a meshgrid as a starting empty background.
Then, once the GUI is loaded, the user can click on the meshgrid and a square is displayed on the selected pixel with blitting, so that no re-plotting needs to be done.

The problem is that between the copy_from_bbox call and the final GUI setup, there is some movement of the widgets (don’t really know why, surely PyQt layout adjustments).
This leads to a shift of the restored region, because the saved one (after the draw() call) is not the one which is finally displayed in the GUI. This can lead to strange display errors …

I tried to trim down my code to the essential, I hope it’s short enough.
I am using Matplotlib 1.3.1 on OS X, but I am having this problem one or two years and had never time to tackle it (I have it also on Windows and Fedora)

Here is my code, just run it and click on the meshgrid, you will see it move down, so that the top of the red square is not drawn correctly.
https://gist.github.com/iMichka/9901318

In the case of my example the shift is quite small (I could live with a shift of one pixel …), but in my real app there is more stuff and the shift is bigger … the plot even overlaps with PyQt buttons …

Is there any way to call copy_from_bbox when the PyQt GUI has stabilized ? This means after the PyQt event loop has finished refreshing, and not before the refreshing ?

Any help would be much appreciated :)


Michka

Hello Marry,

I am highly interested in getting "wrf_cape_3d" wrapped to be accessible in
Python. So far, that's how I calculate CAPE and CIN for my WRF outputs.
wrf_cape_3d is more robust comparing to the function in the SkewT script.
For some reason, I have no luck getting wrf_cape_2d working properly as it
throws a NaN error.


On Mon, Mar 31, 2014 at 2:25 PM, Mary Haley <ha...@uc...> wrote:

> Hi all,
>
> I’ve been following this thread somewhat peripherally.
>
> I’ve slowly started creating Python wrappers of some of the WRF Fortran
> calculation functions (not the graphical ones) that are used in NCL.
>
> You can see the list of the NCL ones at:
>
> http://www.ncl.ucar.edu/Document/Functions/wrf.shtml
>
> So far I’ve only wrapped these six: wrf_avo, wrf_tk, wrf_td, wrf_slp,
> wrf_rh, wrf_dbz.
>
> Would the wrf_cape_2d and wrf_cape_3d routines be of interest? These are
> specific to WRF data.  I believe these are the same ones that Wanli is
> referring to.
>
> We also have the ones that we use for the basic Skew-T code in NCL, that
> Gökhan has been corresponding with Dennis on.
> I could wrap these as well. These routines are not advertised in NCL, but
> they are used by the Skew-T examples you see at:
>
> http://www.ncl.ucar.edu/Applications/skewt.shtml
>
> --Mary
>
> On Mar 31, 2014, at 11:41 AM, Wanli Wu <wu...@gm...> wrote:
>
> All,
>
> Another good example of Skew-T with all Parcel stability info including
> CIN, CAPE is produced through RIP4 (see this example:
> http://www.mmm.ucar.edu/mm5/WRF_post/RIP4/pages/rip_sample_cgm030_gif.htm).
> If this one can be duplicated with python, it'd great for the community.
>
> Wanli Wu
>
>
>
> On Mon, Mar 31, 2014 at 7:49 AM, Gökhan Sever <gok...@gm...>wrote:
>
>> Hi James,
>>
>> I have managed to run CLIMT's thermodyn.py . Most of the functions I
>> tested from within the Driver.f90 works fine except the CAPE and CIN
>> routines. I sent an e-mail to the author regarding this but nothing back
>> from him so far. Would you give a test if I send you a simple sounding
>> data? My system is Window 7 (x64) and using Python(XY). f2py uses the
>> gfortran provided in MinGW32 folder.
>>
>> Could you provide an example (with some test data) for Kerry Emanuel's
>> code? That code has definitely more functions than I need but it might be a
>> valuable source.
>>
>> As for the NCL, it is easy to interface to a WRF output, it also includes
>>  a SkewT/LogP [https://www.ncl.ucar.edu/Applications/skewt.shtml], but
>> the CAPE estimation in this script is very sensitive to the number of
>> data-points, which I have bitten a couple of times. Dennis Shea has
>> provided some CAPE calculation routines coded in fortran. [Check under
>> http://www.cgd.ucar.edu/~shea/ for the files starting with cape*]. Yet I
>> have no luck wrapping them via f2py.
>>
>>
>> On Sat, Mar 29, 2014 at 7:11 PM, James Boyle <jsb...@gm...>wrote:
>>
>>> I have used the CAPE and CIN from ClMT - a bit of overkill but many
>>> useful functions:
>>> http://people.su.se/~rcaba/climt/
>>>
>>> I have also wrapped using f2py the the Fortran CAPE and CIN of Kerry
>>> Emanuel ( a prestigious source) in his convect code:
>>> http://eaps4.mit.edu/faculty/Emanuel/products
>>>
>>> If you prowl about in the NCL source distribution, you will find the
>>> fortran that the NCL skew - T uses.
>>> If you ask, Dennis Shea of NCAR might break the code out for you. It is
>>> trivial to wrap using f2py ( f77).
>>>
>>>
>>> On Mar 29, 2014, at 3:32 PM, Gökhan Sever <gok...@gm...> wrote:
>>>
>>> Hello,
>>>
>>> Lately, I am working on plotting sounding profiles on a SkewT/LogP
>>> diagram. The SkewT package which is located at
>>> https://github.com/tchubb/SkewT has a nice feature to lift a parcel on
>>> dry/moist adiabats. This is very useful to demonstrate the regions of CIN
>>> and CAPE overlaid with the full sounding.
>>>
>>> However, the package misses these diagnostic calculations. This is the
>>> only step holding me back to use Python only (migrating from NCL) for my
>>> plotting tasks.  I am aware that these calculations are usually performed
>>> in fortran. Are there any routines wrapped in Python to calculate CAPE and
>>> CIN parameters?  Any suggestions or comments would be really appreciated.
>>>
>>> --
>>> Gökhan
>>>
>>> ------------------------------------------------------------------------------
>>> _______________________________________________
>>> Matplotlib-users mailing list
>>> Mat...@li...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>>
>>>
>>>
>>
>>
>> --
>> Gökhan
>>
>> _______________________________________________
>> Pyaos mailing list
>> Py...@li...
>> http://lists.johnny-lin.com/listinfo.cgi/pyaos-johnny-lin.com
>>
>>
> _______________________________________________
> Pyaos mailing list
> Py...@li...
> http://lists.johnny-lin.com/listinfo.cgi/pyaos-johnny-lin.com
>
>
>


-- 
Gökhan