bound.for

2D有限差分法程序,未加边界条件,适合一般情况下的波动方程模拟

!     copyright:wangaifeng 
!     DATE     :2007-07-25
!     Company  :GUCAS
      subroutine bound(u,nx,ny,ighost)
	implicit real (a-h,o-z)
      real u(-ighost:nx+ighost,-ighost:ny+ighost,0:3)

      do k=0,3
	do j=-ighost,ny+ighost 
        do i=-ighost+1,-ighost,-1
            u(i,j,k)=u(i+1,j,k)
        enddo 
          do i=nx+1,nx+ighost
            u(i,j,k)=u(i-1,j,k)
          enddo  
      enddo
        do i=-ighost,nx+ighost
          do j=-ighost+1,-ighost,-1
            u(i,j,k)=u(i,j+1,k)
          enddo 
          do j=ny+1,ny+ighost
            u(i,j,k)=u(i,j-1,k)
          enddo  
        enddo  
	end do

	end