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

# Convert Files to RSF and update header
Flow('obc_div_v_1','obc_div_v_1.segy', '''segyread tape=$SOURCE |
        put n2=1381 n3=320 o1=0 o2=0 o3=3000
        d2=12.32 d3=25 label1=Depth\ Z label2=Distance\ X label3=Shot-Cord
        unit1=s unit2=m unit2=m''')

Flow('obc_div_v_2','obc_div_v_2.segy', '''segyread tape=$SOURCE | 
        put n2=1381 n3=160 o1=0 o2=0 o3=11025
        d2=12.32 d3=25 label1=Depth\ Z label2=Distance\ X label2=Shot-Cord
        unit1=s unit2=m unit3=m''')
# Concatinate Datasets
Flow('div',['obc_div_v_1','obc_div_v_2'],'cat ${SOURCES[0:2]} axis=3',stdin=0)

# Plot Data
Result('movie','div','''window $SOURCE  
               min3=4250 max3=4250 n3=1 |
               grey color=I gainpanel=a 
               title=OBC\ Div\ Shot\ 50''')

End()