from rsf.proj import *

# download cmp2.rsf from the server
Fetch('cmp2.rsf','cant12')

# convert cmp2.rsf to native format
Flow('cmp','cmp2','dd form=native')

# create cmp.vpl file
Plot('cmp','grey title=CMP')

wvel=1500

# NMO with water velocity
Flow('nmo1','cmp','nmostretch half=n v0=%g'%wvel)

# mute the stretched amplitudes for shallow events
Flow('nmo','nmo1','mutter v0=1500 x0=400 abs=n')

# create nmo.vpl
Plot('nmo','grey title=NMO')

# create nmo0.vpl
Plot('nmo1','grey title=NMO without mute')

# create nmo0-nmo.vpl file under Fig directory
# nmo0.vpl and nmo.vpl created earlier using Plot
# command will be ploted side by side
Result('nmo1-nmo','nmo1 nmo','SideBySideAniso')

#---------------------
# ADD CODE
#---------------------
# Your task is to add the necessary code to attenuate
# multiples for this CMP.
# The same work flow used in the SConstruct under demo
# directory should work here with only changes to 
# - x0
# - t0
# - v0
# where it says CHANGE ME in the comments


End()