from rsf.proj import *

n1=101
n2=101

Flow('ref.asc',None,
     '''
     echo 0.5 0.5 0.8 1.0 0.8 0.5 0.5
     n1=7 d1=1 o1=0 data_format=ascii_float
     in=$TARGET
     ''')
Flow('ref','ref.asc','dd form=native | remap1 order=2 n1=601 d1=0.01 o1=0')

Flow('syn','ref',
     '''
     kirmod nt=%d dt=0.04 nh=1 
     ns=%d ds=0.05 s0=0.5 vel=1.0 | window
     ''' % (n1,n2))
Result('syn',
       'grey title="Synthetic Data" label1=Time unit1=s label2=Midpoint unit2=km')

Flow('vel','syn','window n2=1 | math output=1')
Flow('mig','syn vel',
     'cosft sign2=1 | gazdag velocity=${SOURCES[1]} | cosft sign2=-1')
Result('mig',
       '''
       grey label1="Migrated Time" unit1=s 
       label2=Midpoint unit2=km 
       title="Phase Shift Migration"
       ''')

End()