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','vel','sfTestfd2d nb=30 nt=400 dt=0.001')
Plot('wav','grey gainpanel=all title="Forward (4th order)" scalebar=y color=j scalebar=y',view=1)
Result('snapfd2d','wav',
        '''
        window n3=1 min3=0.35 |grey title="Sponge ABC(4th order):t=0.35s" label1=Depth 
        unit1=km label2=Lateral unit2=km screenratio=1 color=j 
        ''')

End()