Re: [Pflexible-users] plot accumulative with log10 for multiple levels
A python tool for data analysis of FLEXPART transport modeling.
Status: Beta
Brought to you by:
jfburkhart
|
From: John B. <jfb...@gm...> - 2014-11-14 12:47:11
|
Hello Albert, It would be helpful if you could describe the differences a bit more. However, one thing I see immediately is that level=0 should give you the total column. Why are you breaking the plot ranges into tuples? It looks like you are making individual plots of each levels. Also, you are titling your plots with minlvl and not i... So I need to know if you are trying to plot each individual model level or an average or summation of the minlvl to maxlvl, and then I could provide some code. --john On Thu, Nov 13, 2014 at 6:11 PM, Albert Jornet Puig <alb...@ic...> wrote: > Hi, > > I am having difficulties to plot flexpart outputs with pflexible. We have > several script in Matlab and we would like to move to python. > For this purpose, we are make things work with pflexible if possible. > > As an example, I want to plot levels from 0 to 4, 5 to 15 and the last one > accumulated and then multiply values to log10 to be able to improve > visualitzation. I have seen a couple of example and I managed to get the > code found below but it is not producing the same outputs as in Matlab. > > lvl_values = [(0,4),(5,15),(16,17)] > > data = H.C[(0,0)] > data_range=[-1, 4] > > for minlvl, maxlvl in lvl_values: > TC = None > for i in range(minlvl, (maxlvl)): > print i > TC = pf.plot_at_level(H, H.C[(0,0)], level=i, map_region='Europe', > overlay=True, FIGURE=TC) > filename = '%s_tc_%s_%d_%d_overlay.png' % (data.species, > data.timestamp, minlvl, maxlvl) > TC.fig.savefig(filename) > > > I also checked source code but it is not clear to me a the moment. > > Does somebody have any idea or any way to improve the code? > > Thanks > > > ------------------------------------------------------------------------------ > Comprehensive Server Monitoring with Site24x7. > Monitor 10 servers for $9/Month. > Get alerted through email, SMS, voice calls or mobile push notifications. > Take corrective actions from your mobile device. > > http://pubads.g.doubleclick.net/gampad/clk?id=154624111&iu=/4140/ostg.clktrk > _______________________________________________ > Pflexible-users mailing list > Pfl...@li... > https://lists.sourceforge.net/lists/listinfo/pflexible-users > > |