Born modeling in Madagascar#

1
from rsf.proj import *
2

3
# Prepare source and velocity
4
#############################
5

6
# Source wavelet
7
Flow('wavelet',None,
8
     '''
9
     spike n1=1500 d1=0.001 k1=201 |
10
     ricker1 frequency=10
11
     ''')
12
# Source spectra
13
Result('spectra','wavelet',
14
    '''
15
    spectra |window n1=50 |
16
    graph title="Source Spectra"
17
    ''')
18

19
# Source location
20
Flow('source',None,
21
     '''
22
     spike n1=501 n2=501 d1=0.01 d2=0.01
23
     label1=x1 unit1=km label2=x2 unit2=km
24
     k1=201 k2=251
25
     ''')
26

27
# Background velocity
28
Flow('v0','source','math output=2.')
29

30
# Velocity perturbation
31
Flow('detv',None,
32
    '''
33
    spike n1=501 d1=0.01 o1=0 k1=300,301,302 mag=0.02 |
34
    spray axis=2 n=501 d=0.01 o=0 |
35
    costaper nw2=25
36
    ''')
37

38
# True velocity
39
Flow('vel','v0 detv','add ${SOURCES[1]}')
40
Result('vel',
41
    '''
42
    grey scalebar=y barreverse=y
43
    color=j bias=2.01 title="Velocity"
44
    ''')