Re: overlaying years without regrid ?

[Mick Spillane <Mick.Spillane@xxxxxxxx>]

Hi Jaison,
	You are quite right, and thanks for pointing out what I should
have thought to mention in my response to Domingos.  By hiding the
axes of overplots - time or otherwise - you lose one reality check that
what you do plot makes sense.

	Actually there is a typo in my demo compared to what I had
intended:
  def axis/t="1-jan-1981:12":"31-dec-1984:12":1/units=days tax
i.e. a one day increment, is what I had intended.  The problem I had
thought to mention, but forgot, is that the fourth year in my example,
being a leap year, would have 366 points crammed into a space designed
for 365.  Also the first time I ran through it I forgot to include the
"ppl yaxis" command for years 2,3,4 and their amplitude was fouled up
as a result.

There is another "almost" non-regrid way to solve the problem using @SHF

! overlaying years using the @SHF operation
def axis/t="1-jan-1981:12":"31-dec-1984:12":1/units=days tax
def grid/t=tax grd
let v=sin(6.14*t[g=grd]/350)
plot/t="01-jan-1981:12":"31-dec-1981:12" v
plot/o/t="01-jan-1981:12":"31-dec-1981:12" v[l=@shf:365]
plot/o/t="01-jan-1981:12":"31-dec-1981:12" v[l=@shf:730]
plot/o/t="01-jan-1981:12":"31-dec-1981:12" v[l=@shf:1095]

Regards,
Mick


|!!! Mick.Spillane@noaa.gov !!!|
|__Room 2070 Bldg#3 NOAA/PMEL__|
|____Phone_:_(206)526-6780_____|