Plot3D网格文件的读写

web81 posted @ 2009年1月14日 06:06 in CFD编程应用 with tags Plot3D CFD xyz , 12803 阅读

PLot3D是一种在CFD领域很常见的格式。几乎所有的网格生成软件都支持这种格式的输出，因此搞清楚这种文件的格式的读入和输出对CFD研究人员的工作带来很大的帮助。这里给出的都是读入ASCII格式文件的代码，输出方式类似。

PLot3D的网格文件后缀名为"xyz"，可能是二进制（bin）格式或者文本文件格式(ASCII)。

1.单块网格（Single Grid）

读入3维单块网格

PROGRAM PLOT3D_SINGLE
IMPLICIT NONE
CHARACTER :: BUFFER*30
INTEGER :: I,J,K,IMAX,JMAX,KMAX
REAL,ALLOCATABLE::MESH(:,:,:,:) !MESH存储网格的点坐标

WRITE(*,*)"Please input file name:"
READ(*,*)BUFFER
OPEN(200,FILE=TRIM(buffer))
READ(200,*) IMAX,JMAX,KMAX !获得网格尺寸
ALLOCATE(MESH(3,IMAX,JMAX,KMAX)) !划分网格所需要的内存空间
READ(200,*) (((MESH(1,I,J,K), I=1,IMAX), J=1,JMAX), K=1,KMAX),&
(((MESH(2,I,J,K), I=1,IMAX), J=1,JMAX), K=1,KMAX),&
(((MESH(3,I,J,K), I=1,IMAX), J=1,JMAX), K=1,KMAX)
CLOSE(200)
STOP
END PROGRAM

2.多块网格（Multiple Grids）

读入3维多块网格（似乎GRIDGEN读入网格时需采用这种格式的文件，你可以将其改成输出程序，所输出的网格就可以用GRIDGEN读入了）

PROGRAM PLOT3D_MULTI
IMPLICIT NONE
CHARACTER :: BUFFER*30
INTEGER :: I,J,K,IMAX,JMAX,KMAX,NMESH,N
INTEGER,ALLOCATABLE::IJK(:,:) !IJK存储每块网格的尺寸
REAL,ALLOCATABLE::MESH(:,:,:,:,:) !MESH存储网格的点坐标

WRITE(*,*)"Please input file name:"
READ(*,*)BUFFER
OPEN(200,FILE=TRIM(buffer))
READ(200,*)NMESH !读入网格块数
ALLOCATE(IJK(3,NMESH))
READ(200,*) ((IJK(1,N),IJK(2,N),IJK(3,N)), N=1,NMESH) !获得多块网格尺寸
DO N=1,NMESH
IMAX=MAX(IMAX,IJK(1,N))
JMAX=MAX(IMAX,IJK(2,N))
KMAX=MAX(IMAX,IJK(3,N))
END DO
ALLOCATE(MESH(3,IMAX,JMAX,KMAX,NMESH)) !划分网格所需要的内存空间
READ(200,*) ((((MESH(1,I,J,K,N), I=1,IJK(1,N)), J=1,IJK(2,N)), K=1,IJK(3,N)),&
(((MESH(2,I,J,K,N), I=1,IJK(1,N)), J=1,IJK(2,N)), K=1,IJK(3,N)),&
(((MESH(3,I,J,K,N), I=1,IJK(1,N)), J=1,IJK(2,N)), K=1,IJK(3,N)),N=1,NMESH)
CLOSE(200)
STOP
END PROGRAM

以上代码是我自己编写，如果你有更好的改进方案或者发现错误，请与我联系。