I am having a tough time looking (did search User Forum) for hints on how to overlay two time series (that are represented as unconnected
symbols on the same plot). I have tried various combinations for overlay (including using various PPL commands) and am obviously not doing
something properly. Any help is much appreciated.
For Instance I am trying to plot Column 4 (AMOD_mean_AOT0550img) of the attached ASCII file and ferret output (see attached GIF file) contains only data for the second
time series (for the second station). Please note that the ASCII text file is tab-delimited (not space so that there are only 4 Columns in the file separated
by tabs and the second column "2007-03-01 14:35:00" is actually one string and cannot be broken up. I am manipulating it as below)
I am using the following code snippet to read this file and set up the time axis for the station data for each case
!>>
LET var_year_1 = STRFLOAT( SUBSTRING(begindatetime, 1, 4) )
LET var_month_1 = STRFLOAT( SUBSTRING(begindatetime, 6, 2) )
LET var_day_1 = STRFLOAT( SUBSTRING(begindatetime, 9, 2) )
LET var_hour_1 = STRFLOAT( SUBSTRING(begindatetime, 12, 2) )
LET var_minute_1 = STRFLOAT( SUBSTRING(begindatetime, 15, 2) )
!<<
! Construct the Time Axis using Numeric Variables created above
LET days_2001_1 = DAYS1900(2001,1,1)
LET day_frac_1 = (var_hour_1 + var_minute_1/60 ) / 24
LET tstep_1 = DAYS1900(var_year_1,var_month_1,var_day_1) - days_2001_1 + day_frac_1
The plotting code snippet I am using is below between the !>> !<< symbols. new_var_1 and new_var_2 are the transformed variables from original time series
They are constructed using e.g commands like (note the UNIT=hours below does not give me the proper time-axis output but the UNIT=days at least gives the proper
output even if only for 1 station):
DEFINE axis/T/T0=1-jan-2001/UNIT=DAYS ttaxis_days_1=tstep_1[L=1:`nntime_1`] ! nntime_1 is read from record 3 of attached ascii file
LET dummy_1 = t[GT=ttaxis_days_1]
LET new_var_1 = RESHAPE($3,dummy_1) ! $3 is an argument containing column name e.g AMOD_mean_AOT0550img
DEFINE axis/T/T0=1-jan-2001/UNIT=DAYS ttaxis_days_2=tstep_2[L=1:`nntime_2`] !nntime_2 is also read from record 3 of corresponding block of data fileLET new_var_1 = RESHAPE($3,dummy_1) ! $3 is an argument containing column name e.g AMOD_mean_AOT0550img
LET dummy_2 = t[GT=ttaxis_days_2]
LET new_var_2 = RESHAPE($3,dummy_2)
!>>>>>
! ppl commands with the first plot
! ppl commands with the first plot
plot/symbol=19/size=.2/color=red/nolab/vlim=0:1/set_up new_var_1
ppl plot
plot/symbol=17/size=.2/color=green/nolab/over new_var_2
! ppl commands with the second overlay plot
plot/symbol=19/size=.2/color=red/nolab/vlim=0:1 new_var_1
plot/symbol=17/size=.2/color=green/nolab/over/set_up new_var_2
! insert ppl commands here
ppl plot/over
!<<<<<
plot/symbol=17/size=.2/color=green/nolab/over new_var_2
! ppl commands with the second overlay plot
plot/symbol=17/size=.2/color=green/nolab/over/set_up new_var_2
! insert ppl commands here
ppl plot/over
!<<<<<
Thanks in advance for any pointers and for taking the time.
I did not have any luck using the overlapping viewports suggestions that I read somewhere. BTW is it possible to obtain horizontal and vertical gridlines
on the plot to make it more readable.
Sincerely
agdc2005
On Mon, Sep 28, 2009 at 1:38 PM, Ansley Manke <Ansley.B.Manke@xxxxxxxx> wrote:
Hi -
That depends on how you made the symbols. Are they done with PLOT/SYMBOL ? If so, you'll get a basic plot key at the bottom of the plot, if they're done by first making a blank plot, and then overlaying each line on that. Try this example:
! Define some example variables.
yes? let xx = x[i=1:5]
yes? let variable_A = xx-1
yes? let variable_B = 2*xx
! Plot the variables just to see what their x and y axis extent will be
yes? plot variable_a, variable_b
! Now use /HLIM and /VLIM and plot any variable with /NOLAB and /COLOR=white
! This makes a blank plot that the variables will fit into.
yes? plot/nolab/color=white/hlim=0:6/vlim=-1:11 variable_a
! PLOT/OVERLAY will make a basic key below the plot
yes? plot/over/sym=20/color=green/siz=0.12 variable_a
yes? plot/over/sym=18/color=red/siz=0.12 variable_b
If you want more ideas, maybe try sending this example to the Users List and ask if anyone has a scrip or ideas for making a nice key with symbols. (I see that this could be made better in the internal Ferret code by just drawing one symbol for the key at the bottom, not three symbols!)
Ansley
ag wrote:Thanks again for your help. I did a search of the User Archives to look for legend with symbols. But I could not find how to put a legend of symbols. Please see attached gif file. I was trying to put the respective symbols next to the station names on top right.
On Mon, Sep 28, 2009 at 12:26 PM, Ansley Manke <Ansley.B.Manke@xxxxxxxx> wrote:
Hi Ag,
Use one of the ASCII fonts, @AC or @AS
LABEL/NOUSER `my_xlen`, `my_ylen2`, 0, 0, 0.15 @AC(($dataprodname1)
