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

# Convert Files to RSF and update header
Flow('obc_vz_1','obc_vz_1.segy',
     '''segyread read=d |
     put n2=1381 n3=240 o1=0 o2=0 o3=3000
     d2=12.32 d3=25''')
Flow('obc_vz_2','obc_vz_2.segy',
     '''segyread read=d | 
     put n2=1381 n3=240 o1=0 o2=0 o3=11000
     d2=12.32 d3=25
     ''')
# Concatinate Datasets
Flow('vz',['obc_vz_1','obc_vz_2'],'cat ${SOURCES[0:2]} axis=3',stdin=0)

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

Result('zero2','vz',
       '''window $SOURCE  min2=0 max2=0 size2=1 |
       grey  color=I gainpanel=a label2=Position\ X unit2=m
       title=Zero\ Offset\ Data''')

End()