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('zvel','zvstr',
     '''
     put d1=0.025 o2=10.025 d2=0.025
     label1=Depth unit1=kft label2=Lateral unit2=kft |
     scale dscale=0.001
     ''')

Result('zvel',
       '''
       window f1=200 |
       grey title=Velocity titlesz=7 color=j
       screenratio=0.3125 screenht=4 labelsz=5
       mean=y
       ''')


# Compute approximate reflectivity
##################################
Flow('zref','zvel',
     '''
     depth2time velocity=$SOURCE nt=2501 dt=0.004 |
     ai2refl | ricker1 frequency=10 |
     time2depth velocity=$SOURCE
     ''')

Result('zref',
       '''
       window f1=200 |
       grey title=Reflectivity titlesz=7
       screenratio=0.3125 screenht=4 labelsz=5
       ''')


End()