output.for

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

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

      open(unit=1,file='result.plt') 
      write(1,*)'TITLE = "EXAMPLE: 3D GEOMETRIES" '
      write(1,*)'VARIABLES = "X ", "Y  ", "rou" "u" "v" "p" "E"'
      write(1,*)'ZONE T="Floor", I=',nx,' J=',ny,' F=POINT '
	do i=1,nx
	  do j=1,ny
	   rou=u(i,j,0)
         uu=u(i,j,1)/rou
	   vv=u(i,j,2)/rou
	   p=(gama-1)*(u(i,j,3)-0.5*rou*(uu*uu+vv*vv))
        write(1,*) i*dx,j*dy,rou,uu,vv,p,u(i,j,3)
	enddo
	enddo
      close(1)
	
	
	
	end