from rsfproj import *
import string

n=128
tower = 0.5
spike = []
for i2 in range(n):
    ramp =i2/(n-0.5)
    tower = tower + ramp
    if tower > 1:
        tower = tower - 1
        spike.append(i2+1)

Flow('top',None,
     'spike n2=%d d2=0.01 nsp=%d k2=%s k1=33 n1=64 label1=" " label2=" " ' %
     (n,len(spike),string.join(map(str,spike),',')))
Plot('top','grey pclip=100 wanttitle=n')

Flow('bot','top',
     'cosft sign2=1 | gazdag vel=1.5 inv=1 | cosft sign2=-1')
Plot('bot','grey pclip=99 wanttitle=n')

Result('stormhole','top bot','OverUnderAniso',vppen='txscale=1.5')

Flow('diag',None,
     '''
     spike n1=129 n2=131
     d1=0.004 d2=0.004
     nsp=6 
     k1=10,30,50,70,90,110 
     k2=10,30,50,70,90,110
     ''')
Plot('diag','grey pclip=100 wanttitle=n')

Flow('modl','diag',
     'cosft sign2=1 | gazdag vel=1.5 inv=1 | cosft sign2=-1')
Plot('modl','grey pclip=99.8 wanttitle=n')

Result('diag','diag modl','SideBySideAniso',vppen='txscale=1.5')

End()