from rsf.proj import *

Flow('in',None,
     '''
     spike n1=101 k1=50 d1=1 o1=0 n2=1 o2=0 mag=10
     label2=iter unit2='' label1=sample unit1='' label3=amp unit3=''
     ''')
Flow('parab0.asc',None,
     'echo -0.5 1 -0.5 n1=3 esize=0 data_format=ascii_float in=$TARGET')
Flow('parab','parab0.asc','dd form=native')

for step in ('cg','cd'):
    miss = step+'miss'
    wigg = step+'wigg'
    thpl = step+'thpl'
    cont = step+'cont'
    Flow(miss,'in parab',
         'miss1 step=%s filtin=${SOURCES[1]} niter=450 diter=30' % step)
    Plot(wigg,miss,'wiggle wherexlabel=t labelsz=10 pclip=100 title=" " ')
    Plot(thpl+'0',miss,'thplot labelsz=10 title=" " ')
    Plot(thpl,thpl+'0','SideBySideIso')
    Plot(cont,miss,
         '''
         contour nc=5 dc=2 c0=0 labelsz=10 title=" " allpos=n
         transp=n yreverse=n
         ''')
    Result('dirm'+step,[wigg,thpl,cont],'SideBySideIso')

End()