我列一个详细的tecplot的规矩给你看看吧。
打开open(unit=70,file='field.dat')
首先是结构网格的输出:
二维时:
write(70,*)'VARIABLES="X","Y","U"'
write(70,*)'ZONE F=POINT, I=',nx,', J=',ny
do j=1,ny
do i=1,nx
write(70,*)xc(i),yc(j),u(i,j)
enddo
enddo
三维类似:
write(70,*)'VARIABLES="X","Y","Z","U"'
write(70,*)'ZONE F=POINT, I=',nx,', J=',ny,' K=',nz
do k=1,nz
do j=1,ny
do i=1,nx
write(70,*)xc(i),yc(j),zc(k),u(i,j,k)
enddo
enddo
enddo
非结构网格:
二维4边形单元
write(70,*)'variables= "x" "y" "u" '
write(70,*)'zone f=fepoint et=quadrilateral n=',np,' e=',Nequd
do i=1,np
write(70,*)x(i),y(i),u(i)
enddo
do ee=1,Nequd
e=Equd(ee)
write(70,*)(ep(e,j),j=1,4)
enddo
三角形单元:
write(70,*)'zone f=fepoint et=triangle n=',np,' e=',Netri
do i=1,np
write(70,*)x(i),y(i),u(i)
enddo
do ee=1,Netri
e=Etri(ee)
write(70,*)(ep(e,j),j=1,3)
enddo
其他情况类似。多google、baidu和看手册,很容易找到答案。主要是动手!!!
|