multiple 2d_plots in one system (itools) [message #36319] |
Wed, 10 September 2003 06:29  |
Olaf Stetzer
Messages: 39 Registered: September 2001
|
Member |
|
|
Hello,
I am trying to plot multiple size distributions of an
aerosol in one plot. What I have tried so far: the plot
structure of the J�lihc group, but it is limited to
25 datasets, in my case there can be 100 datasets.
What I tried next is the new itools procedure iplot,
but it does not work as i expected it to. First, the
x axis should be logarithmic (values from 0.038 to 16 �m).
When i set xrange to [min(xdata),max(xdata)] the axis ranges from
~0 to 10e16, so I tried alog10() of the above resulting in
a LINEAR axis from -1.4 to 1.2 even though I have set /x_log.
I cannot see a way to get the axis I want.
Second I don't see any data in the plot although the visualisation
manager shows all datasets (show=true). It's just my first try using
itools, but it seems they are much less intuitive than I expected.
Here is the code I used:
xmin=min(Dp)
xmax=max(Dp)
ymax=max(50)
for i=0,max_nr do begin
if i eq 0 then begin
iplot, Dp[*,i], dN[*,i]/duration[i], $
/x_log,xrange=[xmin,xmax],yrange=[0,ymax]
view_nr=itgetcurrent()
endif else iplot, Dp[*,i], dN[*,i]/duration[i], $
/x_log,xrange=[xmin,xmax],yrange=[0,ymax],overplot=view_nr
endfor
Omitting xrange, yrange and x_log does not change the result
(as can be expected).
Another question:
Is there another easy way to get these plots done? I wonder
if there is a way to get them in a 3D plot as stacked xy plots
with z as the time of measurement but NOT as surface plots
(but filled to the x-axis would be OK). Like this:
y (counts, dN)
|
|
| z(time)
| / **
| / **** *** **
| / *** *** *****
| ********** ******
| / ++
| / +++ +++ ++
|++++ ++++++++++ ++++++++++++++++
-------------------------------------------x (size, Dp)
Hmmm, I hope you get my idea :-)
So far, I haven't seen a way to do this in IDL.
Thanks,
Olaf
|
|
|