如何将gridgen生成的网格输出

pilotlvhy

请教:
如何将gridgen生成的网格输出为tecplot或者fortran认知的那种格式?

ilovecfd

采用gridgen或plot3d格式输出，再用程序转化即可

boling2000

tecplot 格式可以认识plot3d格式，读入数据的时候不要选用Load Data File(s), 而选用 Import ,在从里面选择 Plot3d Loader 就可以读入了。

pilotlvhy

thanks
昨天我编了一个将gridgen格式的转化为fortran格式的程序，效果蛮不错的，呵呵

pilotlvhy

西工大的学长，我输出的格式是ｐｌｏｔ３Ｄ，为什么有ｂｌｏｃｋ的网格还是不能在ｔｅｃｐｌｏｔ里面打开？而没有ｂｌｏｃｋ的二维网格却可以采用ｇｒｉｄｇｅｎ格式输出，然后在ｔｅｃｐｌｏｔ中打开？

boling2000

there are several options in tecplot, I guess you could have a try or read the help file, actually I am not familar with that kind of format , you know, I just like to write the standard format for tecplot.
but one of my friend liked to use plot3d format.

pilotlvhy

ok,i get it

bjs30

pilotlvhy 请您指教如何编程将gridgen格式的转化为fortran格式的程序？
谢谢！

sillybear

谢谢！pilotlvhy
转化为fortran格式的程序能不能提供给我们学习呀。

npuzx

hao a

pilotlvhy

好的

zhmeww

3x

pilotlvhy

real,allocatable::x(:,:,,y(:,:,,z(:,:,
real,allocatable::x1(:),y1(:),z1(:)
write(*,*) 'grid file name'
read(*,*) gridfile
ind=index(gridfile,' ')
open(1,file=gridfile(:ind-1),STATUS='UNKNOWN')
read(1,*)  n_block
      read(1,*) io,jo,ko
close(1)  
nijk=io*jo*ko
m=int(nijk/4)
mmod=nijk-4*m
allocate(x(io,jo,ko),y(io,jo,ko),z(io,jo,ko))
allocate(x1(nijk),y1(nijk),z1(nijk))
open(1,file=gridfile(:ind-1),STATUS='UNKNOWN')
read(1,*) n_block
      read(1,*) io,jo,ko
do 111 i=0,m-1
111read(1,*) x1(4*i+1),x1(4*i+2),x1(4*i+3),x1(4*i+4)
if(mmod.ne.0) then
do i=4*m+1,nijk
read(1,*) x1(i)
enddo
endif
do 112 i=0,m-1
112read(1,*) y1(4*i+1),y1(4*i+2),y1(4*i+3),y1(4*i+4)
if(mmod.ne.0) then
do i=4*m+1,nijk
read(1,*) y1(i)   
enddo
endif
do 113 i=0,m-1
113read(1,*) z1(4*i+1),z1(4*i+2),z1(4*i+3),z1(4*i+4)
if(mmod.ne.0) then
do i=4*m+1,nijk
read(1,*) z1(i)   
enddo
endif
close(1)
mm=0
do 11 k=ko,1,-1
do 11 j=jo,1,-1
do 11 i=io,1,-1
mm=1+mm
x(i,j,k)=x1(mm)
y(i,j,k)=y1(mm)
z(i,j,k)=z1(mm)
11continue
OPEN(12,FILE=gridfile(:ind-5)//'_ijk.plt',STATUS='UNKNOWN')
      write(12,*)'zone i=',io,' j=',jo,' k=',ko,' f=point'
do 114 k=1,ko
do 114 j=1,jo
do 114 i=1,io
write(12,*) x(i,j,k),y(i,j,k),z(i,j,k)
114continue
      CLOSE(12)
end

pilotlvhy

批评指正，其实也是没必要的，直接用gridgen输出就行了

zhujianguo

提供给大家一块学习学习把