非常感谢白文老师的回复,真的是太感谢了!如果方便的话麻烦您帮我看看我写的这个程序那里有问题呢?死活就是读不出来,急死了!
dimension imax(35),jmax(35),kmax(35)
dimension x(300,110,200),y(300,110,110),z(300,110,110),
1 iblank(300,110,110)
open(1,file='grid.in',form='unformatted')
open(2,file='mesh.dat')
write(2,*)'variables=x,y,z,iblank'
iunit = 1
read(iunit)nblock
read(iunit)(imax(n),jmax(n),kmax(n),n=1,nblock)
do 10 n=1,nblock
read(iunit) (((x(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))
1 ,(((y(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))
2 ,(((z(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))
3 ,(((iblank(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))
write(2,*)'zone i=',imax(n),' j=',jmax(n),' k=',kmax(n)
do k = 1,kmax(n)
do j = 1,jmax(n)
do i = 1,imax(n)
write(2,*)x(i,j,k),y(i,j,k),z(i,j,k),iblank(i,j,k)
end do
end do
end do
10 continue
close(1)
close(2)
end |