from rsf.proj import *

Flow('spike',None,'spike n1=30 n2=30 k1=5 k2=1 d1=1 label1= label2= label3=')

planes = []
for nw in (2,3,8,13):
    plane = 'plane%d' % nw
    flt = 'flt%d' % nw
    lag = 'lag%d' % nw

    Flow(plane,'spike',
         'plane n1=100 n2=100 a1=%d b1=1 lag=${TARGETS[1]}' % nw)
    Plot('g'+plane,plane,
         '''
         grey pclip=100 title="Lagrange-%d" wantaxis=n wheretitle=t 
         titlesz=18 crowd=0.88
         ''' % nw)
    Plot('w'+plane,plane,
         '''
         wiggle pclip=100 wanttitle=n wantaxis=n transp=1 yreverse=1 
         poly=y crowd=0.88  titlesz=18
         ''')
    Plot('s'+plane,plane,
         '''
         spectra all=y |  
         graph crowd=0.88 labelsz=16 titlesz=18 wheretitle=t title=Spectrum
         ''')
    Plot(plane,['g'+plane,'w'+plane,'s'+plane],'SideBySideAniso')
    planes.append(plane)

Result('steer-lagrange',planes,'OverUnderAniso')

Flow('hyp','spike','spike n1=30 n2=30 nsp=3 k1=5,10,15  k2=1,1,1 d1=1 | plane a1=7 hyp=1')
Plot('ghyp','hyp',
     '''
     grey pclip=100 wanttitle=n wantaxis=n crowd=0.88
     ''')
Plot('whyp','hyp',
     '''
     wiggle pclip=100 wanttitle=n wantaxis=n transp=1 yreverse=1 
     poly=y crowd=0.88  
     ''')
Result('steer-hyp7','ghyp whyp','SideBySideAniso')

End()