from rsf.proj import *


Flow('vp',None,
        '''
        math output=2 n1=200 n2=200 d1=0.005 d2=0.005
        label1=x1 unit1=km label2=x2 unit2=km 
        ''')
Flow('vs','vp',
        '''
        math output="0.707*input" n1=200 n2=200 d1=0.005 d2=0.005
        label1=x1 unit1=km label2=x2 unit2=km 
        ''')
Flow('rho',None,
        '''
        math output=1 n1=200 n2=200 d1=0.005 d2=0.005
        label1=x1 unit1=km label2=x2 unit2=km 
        ''')

Flow('wavz wavx','vp vs rho',
        '''
        Testelastic2d  vs=${SOURCES[1]} rho=${SOURCES[2]} wavx=${TARGETS[1]} 
        nb=30 fm=25 nt=300 kt=270 dt=0.001
        ''')

Plot('wavx',
     '''
     grey title="Component-x" screenratio=1
     label1=Depth unit1=km label2=Lateral unit2=km color=i 
     ''')
Plot('wavz',
     '''
     grey title="Component-z" screenratio=1
     label1=Depth unit1=km label2=Lateral unit2=km color=i 
     ''')
Result('elasticxz','wavx wavz', 'SideBySideIso')

End()