from rsf.proj import *

# Download velocity model from the data server
##############################################
vstr = 'sigsbee2a_stratigraphy.sgy'
Fetch(vstr,'sigsbee')
Flow('zvstr',vstr,'segyread read=data')

Flow('vel','zvstr',
     '''
     put d1=0.025 d2=0.025 o1=0 o2=10.025 
     label1=Depth unit1=kft label2=Lateral unit2=kft |
     scale dscale=0.001 
     ''')

Result('sigsbee','vel',
       '''
       grey title=Velocity titlesz=7 color=j
       screenratio=0.375 screenht=4 labelsz=5 scalebar=y
       mean=y bartype=v barlabel="V" barunit="kft/s"
       ''')


Flow('simag1 simag2','vel',
        '''
        gpurtm imag2=${TARGETS[1]} NJ=6 phost=75
        fm=15 dt=0.0008 tdmute=300
        nt=13000 ns=55 ng=301 
        jsx=50 jsz=0 jgx=1 jgz=0 
        sxbeg=150 szbeg=1 gxbeg=0 gzbeg=1 
        vmute=4.72
        ''')

Result('simag1','grey allpos=n title="correlation"')
Result('simag2','grey allpos=n title="normalized_correlation"')

End()