Parallel Distributed Simulation of Elastic Waves

Page 66

58 - !!J.thCp [ v(i + 1/2,3/2, k)- v(i + 1/2, 1/2, k)] v(i + 1/2, NY- 1/2, k, n + 1) = v(i + 1/2, NY- 1/2, k) - !!J.thCp [v(i + 1/2, NY- 1/2, k)- v(i + 1/2, NY- 3/2, k)l v(i + 1/2, j + 1/2, 0, N Z- 1/2) = v(i + 1/2, j + 1/2, N Z- 1/2) - !!J.thCp [v(i + 1/2,j + 1/2, N Z- 1/2)- v(i + 1/2,j + 1/2, N Z- 3/2)] w{1/2,j,k+ 1/2,n+ 1) = w{1/2,j,k+ 1/2) - !!J.th Cp [ w(3/2, j, k + 1/2)- w(1/2, j, k + 1/2)] w(NX -1/2,j,k+ 1/2,n+ 1) = w(NX -1/2,j,k+ 1/2) - !!J.thCp [w(N X- 1/2,j, k + 1/2)- w(N X- 3/2, j, k + 1/2)] w(i + 1/2, 0, k + 1/2, n + 1) = w(i + 1/2, 0, k + 1/2) -~Tt c [w (z. + 1/2, 1, k + 1/2)- w(z. + 1/2, 0, k + 1/2)l w(i + 1/2, NY- 1/2, k + 1/2, n +I)= w(i + 1/2, NY- 1/2, k + 1/2) ~tC [ -T w(i + 1/2, NY- 1/2, k + 1/2)- w(i + 1/2, NY- 3/2, k + 1/2)l w(i + l/2,j, 0, NZ- 1/2) = w(i + 1/2,j, NZ- 1/2) ~tC [ -T w(i + 1/2,j, N Z- 1/2)- w(i + l/2,j, N Z- 3/2)l