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]}')

Result('vel','grey title="velocity model: 3 layers" color=j scalebar=y')

Flow('shots','vel',
        '''
        modeling2d nt=1100 dt=0.001 ns=10 ng=200        
        sxbeg=5 szbeg=2 jsx=20 jsz=0 
        gxbeg=0 gzbeg=3 jgx=1 jgz=0     
        ''')

Result('shots',
        '''
       byte allpos=n gainpanel=all |
       grey3 flat=n frame1=300 frame2=100 frame3=5 
       label1=Time unit1=s 
       label2="Receiver no." label3="Shot no."
       title="Shot records" point1=0.8 point2=0.8
        ''')

Plot('shots','grey title=Shots',view=1)


End()