up [pdf]
from rsf.proj import *

Flow('vel',None,
        '''
        math output=2.0 n1=200 n2=200 d1=0.005 d2=0.005
        label1=x1 unit1=km label2=x2 unit2=km 
        ''')

Flow('wav pz px','vel','sfTestspml pz=${TARGETS[1]} px=${TARGETS[2]} nb=30 nt=400 dt=0.001 verb=y kt=250')
Plot('wav','grey gainpanel=all title="Forward (8th order)" scalebar=y color=j scalebar=y',view=1)
Result('snapspml','wav',
        '''
        window n3=1 min3=0.35 |grey title="SPML(8th order):t=0.35s" label1=Depth 
        unit1=km label2=Lateral unit2=km screenratio=1 color=j 
        ''')

Result('pz','grey title="Component-z" label1=Depth unit1=km label2=Lateral unit2=km screenratio=1 color=j')
Result('px','grey title="Component-x" label1=Depth unit1=km label2=Lateral unit2=km screenratio=1 color=j')

End()

sfmath
sfTestspml
sfgrey
sfwindow