from rsf.proj import *

Flow('vel',None,
        '''
        math output=2.0 n1=100 n2=100 n3=100 
        d1=0.005 d2=0.005 d3=0.005 
        label1=x1 unit1=km label2=x2 unit2=km 
        label3=x3 unit3=km 
        ''')


Flow('snapfd3d','vel','sfTestfd3d fm=20 dt=0.001 nt=250 kt=200 verb=y')
Result('snapfd3d',
       '''
       byte allpos=n gainpanel=all |
       grey3 flat=n frame1=10 frame2=25 frame3=25 color=j
       title="3D FD modeling:kt=250" point1=0.8 point2=0.8
       ''')

End()