pdf 260K - Madagascar

10
76 s f p u t f l o a t (Fo , ”o3” , f t ∗ dt ) ;
77
78 dt2 = dt ∗ dt ;
79
80 /∗ s e t Laplacian c o e f f i c i e n t s ∗/
81 d1 = 1 . 0 / ( d1∗d1 ) ;
82 d2 = 1 . 0 / ( d2∗d2 ) ;
83
84 c11 = 4.0∗ d1 / 3 . 0 ;
85 c12= −d1 / 1 2 . 0 ;
86 c21 = 4.0∗ d2 / 3 . 0 ;
87 c22= −d2 / 1 2 . 0 ;
88 c0 = −2.0 ∗ ( c11+c12+c21+c22 ) ;
89
90 /∗ read wavelet , v e l o c i t y & source p o s i t i o n ∗/
91 ww=s f f l o a t a l l o c ( nt ) ; s f f l o a t r e a d (ww , nt ,Fw) ;
92 vv=s f f l o a t a l l o c 2 ( n1 , n2 ) ; s f f l o a t r e a d ( vv [ 0 ] , n12 , Fv ) ;
93 r r=s f f l o a t a l l o c 2 ( n1 , n2 ) ; s f f l o a t r e a d ( r r [ 0 ] , n12 , Fr ) ;
94
95 /∗ a l l o c a t e temporary arrays ∗/
96 u0=s f f l o a t a l l o c 2 ( n1 , n2 ) ;
97 u1=s f f l o a t a l l o c 2 ( n1 , n2 ) ;
98 u2=s f f l o a t a l l o c 2 ( n1 , n2 ) ;
99 ud=s f f l o a t a l l o c 2 ( n1 , n2 ) ;
100
101 for ( i 2 =0; i2