from rsf.proj import *

Flow('vel0',None,
        '''
        math output=1.6 n1=50 n2=200 d1=0.005 d2=0.005
        label1=x1 unit1=km label2=x2 unit2=km 
        ''')

Flow('vel1',None,
        '''
        math output=1.8 n1=50 n2=200 d1=0.005 d2=0.005
        label1=x1 unit1=km label2=x2 unit2=km 
        ''')
Flow('vel2',None,
        '''
        math output=2.0 n1=100 n2=200 d1=0.005 d2=0.005
        label1=x1 unit1=km label2=x2 unit2=km 
        ''')
Flow('vel',['vel0','vel1','vel2'], 'cat axis=1 ${SOURCES[1:3]}')


Flow('shot check1 check2','vel',
        '''     
        sffbrec2d check1=${TARGETS[1]} check2=${TARGETS[2]} 
        csdgather=n fm=15 dt=0.001 ns=1 ng=200 nt=1100 ng=200 kt=550    
        sxbeg=100 szbeg=2 jsx=37 jsz=0 
        gxbeg=0 gzbeg=1 jgx=1 jgz=0
        ''')

Result('shot','grey gainpanel=all title=shot')


Flow('diff','check1 check2','sfadd ${SOURCES[1]} scale=1,-1')


Plot('check1','grey title="snapshot forward" scalebar=y')
Plot('check2','grey title="snapshot backward" scalebar=y')
Plot('diff','grey title="difference" scalebar=y')

Result('check','check1 check2 diff','SideBySideIso')

End()