solveall.for

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

!     copyright:wangaifeng 
!     DATE     :2007-07-25
!     Company  :GUCAS
      subroutine solveall(u,nx,ny,dt,dx,dy,ighost,gama)
	implicit real(a-h,o-z)
	real u(-ighost:nx+ighost,-ighost:ny+ighost,0:3)
	real gt(-ighost:nx+ighost,-ighost:ny+ighost,0:3)
	real ft(-ighost:nx+ighost,-ighost:ny+ighost,0:3)

      call Solvex(u,ft,nx,ny,ighost,gama)
      call Solvey(u,gt,nx,ny,ighost,gama)
	
	r1=dt/dx
	r2=dt/dy

	do i=-ighost+2,nx+ighost-2
	do j=-ighost+2,ny+ighost-2 
	do k=0,3
	u(i,j,k)=u(i,j,k)-r1*(ft(i,j,k)-ft(i-1,j,k))
     &	-r2*(gt(i,j,k)-gt(i,j-1,k))
	enddo
	enddo
	enddo

	end