pdf 260K - Madagascar

12
20 integer n1 , n2
21
22 n1 = size ( uin , 1 )
23 n2 = size ( uin , 2 )
24
25 uout ( 3 : n1 −2 ,3: n2−2) = &
26 c11 ∗( uin ( 2 : n1 −3 ,3: n2−2) + uin ( 4 : n1 −1 ,3: n2 −2)) + &
27 c12 ∗( uin ( 1 : n1 −4 ,3: n2−2) + uin ( 5 : n1 , 3 : n2 −2)) + &
28 c21 ∗( uin ( 3 : n1 −2 ,2: n2−3) + uin ( 3 : n1 −2 ,4: n2 −1)) + &
29 c22 ∗( uin ( 3 : n1 −2 ,1: n2−4) + uin ( 3 : n1 −2 ,5: n2 ) ) + &
30 c0 ∗ uin ( 3 : n1 −2 ,3: n2−2)
31 end subroutine l a p l a c i a n
32 end module l a p l a c e
33
34 program Wave
35 use r s f
36 use l a p l a c e
37 implicit none
38
39 integer : : i t , nt , ft , jt , n1 , n2
40 real : : dt , d1 , d2 , dt2
41
42 real , dimension ( : ) , allocatable : : ww
43 real , dimension ( : , : ) , allocatable : : vv , r r
44 real , dimension ( : , : ) , allocatable : : u0 , u1 , u2 , ud
45
46 type ( f i l e ) : : Fw, Fv , Fr , Fo ! I /O f i l e s
47
48 c a l l s f i n i t ( ) ! i n i t i a l i z e Madagascar
49
50 ! setup I /O f i l e s
51 Fr = r s f i n p u t ( ” in ” ) ! source p o s i t i o n
52 Fo = r s f o u t p u t ( ” out ” ) ! output w a v e f i e l d
53
54 Fw = r s f i n p u t ( ”wav” ) ! source w a v e l e t
55 Fv = r s f i n p u t ( ”v” ) ! v e l o c i t y
56
57 ! Read/Write axes
58 c a l l from par ( Fr , ”n1” , n1 )
59 c a l l from par ( Fr , ”n2” , n2 )
60 c a l l from par ( Fr , ”d1” , d1 )
61 c a l l from par ( Fr , ”d2” , d2 )
62
63 c a l l from par (Fw, ”n1” , nt )
64 c a l l from par (Fw, ”d1” , dt )