[Thread Prev][Thread Next][Index]

shaded area between two lines



Hi,

I would like to shade the area between two lines.
Is there a straightforward way to do this?

I started using the 3 argument shade command, but I do not
quite get what I want.
I want to represent the shaded area +/- one standard 
deviation around something along my time axis. So I calculated 
the 2d-variables:
TT=(t-dt, t+dt) XX=(x-std, x+std) and made up a 2d-dummy FF
of the same size containing only 1's.
Then
shade FF,TT,XX 
does look right, except that I do not have the time axis I had 
before any more. Somehow, I can't get the right time axis 
into that plot.
This would be nice, because one could them overlay some other
data one wants to compare to the indicated range.
However, I fail getting any overlay to work...


Here is what I do:

I have 15 min data:

define axis/t="15-AUG-1998 00:07:30":"31-DEC-1999
23:52:30":0.010416667/unit="days since 01-Jan-1998 00:00:00" timtax

define axis/x=1:2:1 timxax2
define axis/x=1:4:1 timxax4

define grid/t=timtax/x=timxax2 twgrd2
define grid/t=timtax/x=timxax4 twgrd4

! time series
set data/ez/var=w/grid=twgrd4/col=4 wrccrf_swc_199899_tim_d.dat

! data file with the data for the 3 argument shade
set data/ez/var=ff,xx,tt/grid=twgrd2/col=3 tt

shade ff,tt,xx

plot/over w[d=1,i=@ave]

This plot over now does not understand the axes from the previous shade.

Vice versa,

plot w[d=1,i=@ave]

then

shade/over/noaxis  ff,tt,xx

draws a new plot.



Any ideas would be most appreciated.

Cheers,
Joerg



--
Jörg Kaduk                            Tel.: 1 650 325 1521 x 416
Carnegie Institution of Washington    FAX: 1 650 325 6857
Dept. of Plant Biology
260 Panama Street                     joerg@jasper.stanford.edu
Stanford, CA 94305-1297               http://Jasper.Stanford.EDU/joerg/


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement