[Fortran] 纯文本查看 复制代码
c***************write data to TECPLOT file
c
open(unit=15,file='tecplot.dat',form='formatted',status='unknown')
c-------three-dimensional
text='VARIABLES = X Y Z U V W P K E VIS'
write(15,90)text
90 format(1x,70a)
write(15,92)ni,nj,nk
92 format(1x,'ZONE T="Zone-One", I= ',i3,' ,J= ',i3,' ,K= ',i3,
. ' ,F=POINT')
do k=1,nk
do j=1,nj
do i=1,ni
write(15,98)xp(i,j,k),yp(i,j,k),zp(i,j,k),phi(i,j,k,u),
&phi(i,j,k,v),phi(i,j,k,w),phi(i,j,k,p),phi(i,j,k,te),
&phi(i,j,k,ed),vis(i,j,k)
end do
end do
end do
close(15)
88 format(1x,8(1pe13.6,2x))
98 format(1x,10(1pe13.6,2x))
return
end