# Re: [ferret_users] More specific question on time series overlay wrt Grid Transformation and Definition

 Hi Ag, You can overlay time series plots even if the time axes are not the same. For instance, try this quick script.  monthly_navy_winds is one of the datasets that comes with the Ferret distribution. The time axis of variable uwnd from monthly_navy_winds has units of days, and a time origin of 14-Jan-1980; the time axis I define has units of months and a time origin of 15-jan-1901. yes? use monthly_navy_winds yes? let var1 = uwnd[x=300,y=1] yes? define axis/t=1-jan-1982:1-jan-1991:1/units=months tmonth  *** NOTE: /UNIT=MONTHS is ambiguous ... using 1/12 of 365.2425 days yes? let var2 = 2*cos(t[gt=tmonth]) yes? plot var1,var2 For your script, I would suggest you plot each of your time series separately. Define each variable, don't regrid to the same time axis, but just create each on the time axis that makes sense for it.  Then look at the individual plots and see if they look correct.  Are the commands to read from your ascii file really skipping the first block and reading the second set of data? You could do this: SET VIEW UPPER PLOT new_var_11 SET VIEW LOWER new_var_21 ag wrote: Dear Ferret User Group: OK ...after a more thorough search....I understand that in order to place (overlay) a second time series, I need to re-grid the second variable to the time grid that is specified by the first time axis...here are the commands I am using to do so and I am still getting only one plot in the figure. Any pointers are greatly appreciated. !For first data block  in ASCII file...(\$3 is the variable I want to plot specified by Column name) DEFINE axis/T/T0=1-jan-2001/UNIT=DAYS ttaxis_days_1=tstep_1[L=1:`nntime_1`]    !nntime_1 is read from file LET dummy_1   = t[GT=ttaxis_days_1] LET new_var_1 = RESHAPE(\$3,dummy_1) list new_var_1  ! Prints out the correct values DEFINE GRID/T=ttaxis_days_1 gttaxis_1 LET dummy_11  = t[GT=gttaxis_1] LET new_var_11 = RESHAPE(\$3,dummy_11) ! For Second Data Block in ASCII file ! Read File for Station = S2 and DataProduct/Dataset =D1 ! Move Back to beginning of ASCII File CANCEL DATA "\$1" !...SKIP the first data block and read second one etc! DEFINE axis/T/T0=1-jan-2001/UNIT=DAYS ttaxis_days_2=tstep_2[L=1:`nntime_2`]  !nntime_2 is also read from file and is NE nntime_1 LET dummy_2   = t[GT=ttaxis_days_2] LET new_var_2 = RESHAPE(\$3,dummy_2) ! Now re-gridding second Time Series Variable new_var_2 to time axis ttaxis_days_1 (defined above) DEFINE GRID/T=ttaxis_days_1 gttaxis_2 LET dummy_21  = t[G=gttaxis_2] LET new_var_21 = RESHAPE(\$3,dummy_21) list new_var_11, new_var_21 ! Unfortunately, the above diagnostic shows identical values for new_var_11 and new_var_21...this is should not be the case ! Make the Time Series Plots as unconnected symbols with different colors PLOT /NOLAB/THICKNESS=3/SIZE=0.2/SET new_var_11, new_var_21[g=gttaxis_1, gt=@asn] !PPL TXLINT,1,0 PPL pen 1,2                   ! set pen 1 color to 2 (red) PPL line 1,19,3               ! set line 1 to symbol #19, unconnected PPL pen 2,3                   ! set pen 2 color to 3 (green) PPL line 2,17,3               ! set line 2 to symbol #17, unconnected PPL PLOT ! creat gif file ! Create a GIF file in local directory FRAME/FILE="\$2" SPAWN chmod ugo+rw "\$2" quit Thanks to any insights. agdc2005

• References: