pdf 260K - Madagascar

16
24 r i c k e r 1 frequency=10
25 ’ ’ ’ )
26
27 # Source l o c a t i o n
28 Flow ( ’ source ’ , None ,
29 ’ ’ ’
30 s p i k e n1=201 n2=301 d1=0.01 d2=0.01
31 l a b e l 1=x1 unit1=km l a b e l 2=x2 unit2=km
32 k1=101 k2=151
33 ’ ’ ’ )
34
35 # V e l o c i t y model
36 Flow ( ’ v1 ’ , ’ source ’ , ’ math output=1 ’ )
37 Flow ( ’ v2 ’ , ’ source ’ , ’ math output =1.5 ’ )
38
39 # Modeling
40 Flow ( ’ wave ’ , ’ source %s wavelet v1 v2 ’ % exe [ 0 ] ,
41 ’ ’ ’
42 . / ${SOURCES[ 1 ] } wav=${SOURCES[ 2 ] }
43 v=${SOURCES[ 3 ] } vx=${SOURCES[ 4 ] }
44 f t =200 j t =5
45 ’ ’ ’ )
46
47 Plot ( ’ wave ’ , ’ grey gainpanel=a l l t i t l e=Wave ’ , view=1)
48
49 Result ( ’ wave ’ ,
50 ’ ’ ’
51 window n3=1 min3=0.9 |
52 grey t i t l e =Wave s c r e e n h t=8 screenwd=12
53 ’ ’ ’ )
54
55 # Download Hess VTI model
56 #########################
57 zcat = WhereIs ( ’ gzcat ’ ) or WhereIs ( ’ zcat ’ )
58 for case in ( ’ vp ’ , ’ e p s i l o n ’ ) :
59 sgy = ’ t i m o d e l %s . segy ’ % case
60 sgyz = sgy + ’ . gz ’
61 Fetch ( sgyz , d i r=’ Hess VTI ’ ,
62 s e r v e r=’ f t p : / / s o f tware . seg . org ’ ,
63 top=’ pub/ d a t a s e t s /2D ’ )
64 # Uncompress
65 Flow ( sgy , sgyz , zcat + ’ $SOURCE ’ , s t d i n =0)
66 # Convert to RSF format
67 Flow ( case , sgy ,
68 ’ ’ ’