from rsf.proj import *
# Fetch Files from repository
Fetch("obc_vx_1.segy","marm2")
Fetch("obc_vx_2.segy","marm2")

# Convert Files to RSF and update header
Flow('obc_vx_1','obc_vx_1.segy', '''segyread tape=$SOURCE |
        put n2=1381 n3=320 o1=0 o2=0 o3=3000
        d2=12.32 d3=25''',stdin=0)
Flow('obc_vx_2','obc_vx_2.segy', '''segyread tape=$SOURCE | 
        put n2=1381 n3=160 o1=0 o2=0 o3=11000
        d2=12.32 d3=25''',stdin=0)
# Concatinate Datasets
Flow('vx',['obc_vx_1','obc_vx_2'],'cat ${SOURCES[0:2]} axis=3',stdin=0)

# Plot Data
Result('zero','vx','''window $SOURCE  min2=0 max2=0 size2=1 |
               grey  color=I gainpanel=a label2=Position\ X unit2=m
               pclip=98 title=Zero\ Offset\ Data''')

Plot('movie','vx',
     '''window $SOURCE  j3=20 |
     grey  color=I gainpanel=a label2=Position\ X unit2=m
     title=Zero\ Offset\ Data''',view=1)

End()