pdf 260K - Madagascar

14
110 c a l l r s f w r i t e (Fo , u1 )
111 end i f
112 end do
113 write ( 0 , ∗ )
114
115 c a l l exit ( 0 )
116 end program Wave
Your task is to modify the code to implement elliptically-anisotropic wave propagation
according to equation
∂ 2 P
∂t 2 = V 2
1 (x) ∂2 P
∂x 2 1
+ V 2
2 (x) ∂2 P
∂x 2 2
+ F (x, t) (17)
You can test your implementation using a constant velocity example shown in Figure
2.
Figure 2: Wavefield snapshot for
propagation from a point-source
in a homogeneous medium. Modify
the code to make wave propagation
anisotropic.
1. Change directory to geo384w/hw1/code
2. Run
scons wave.vpl
to compile and run the program and to observe a propagating wave on your
screen.
3. Open the file wave.c in your favorite editor and modify it to implement the
wave operator from equation (17).
4. Run