马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
A lof of software accept the PLOT3D format, and here it is a small introduction of plot3d format.
==============================================================================
Grid Files
The grid file (XYZ file) defines the coordinates of the grid points in the structured mesh. The points are stored in IJK order with all of the X components stored first, followed by the Y coordinates and the Z coordinates. If the file contains IBlank information, it is stored after the Z coordinates.
XYZ Single Grid
READ(IUNIT) IDIM, JDIM, KDIM
READ(IUNIT) (((X(I,J,K), I=1,IDIM), J=1,JDIM),K=1,KDIM),
C (((Y(I,J,K), I=1,IDIM), J=1,JDIM),K=1,KDIM),
C (((Z(I,J,K), I=1,IDIM), J=1,JDIM),K=1,KDIM)
XYZ with IBlank
READ(IUNIT) IDIM, JDIM, KDIM
READ(IUNIT) (((X(I,J,K), I=1,IDIM), J=1,JDIM),K=1,KDIM),
C (((Y(I,J,K), I=1,IDIM), J=1,JDIM),K=1,KDIM),
C (((Z(I,J,K), I=1,IDIM), J=1,JDIM),K=1,KDIM),
C (((IBLANK(I,J,K), I=1,IDIM), J=1,JDIM),K=1,KDIM)
2D XYZ File Single Grid
READ(IUNIT) IDIM, JDIM
READ(IUNIT) (((X(I,J,K), I=1,IDIM), J=1,JDIM),
C (((Y(I,J,K), I=1,IDIM), J=1,JDIM)
XYZ File with Multiple Grids
READ(IUNIT) NGRID
READ(IUNIT) (IDIM(IGRID), JDIM(IGRID), KDIM(IGRID), IGRID=1,NGRID)
DO 10 IGRID=1,NGRID
READ(IUNIT)
C (((X(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),K=1,KDIM(IGRID)),
C (((Y(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),K=1,KDIM(IGRID)),
C (((Z(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),K=1,KDIM(IGRID))
10CONTINUE
XYZ File with Multiple Grids and IBlanks
READ(IUNIT) NGRID
READ(IUNIT) (IDIM(IGRID), JDIM(IGRID), KDIM(IGRID), IGRID=1,NGRID)
DO 10 IGRID=1,NGRID
READ(IUNIT)
C (((X(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),K=1,KDIM(IGRID)),
C (((Y(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),K=1,KDIM(IGRID)),
C (((Z(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),K=1,KDIM(IGRID)),
C (((IBLANK(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),K=1,KDIM(IGRID))
10CONTINUE
2D XYZ File Multiple Grid
READ(IUNIT) NGRID
READ(IUNIT) (IDIM(IGRID), JDIM(IGRID), GRID=1,NGRID)
DO 10 IGRID=1,NGRID
READ(IUNIT)
C (((X(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)),
C (((Y(I,J,K), I=1,IDIM(IGRID)), J=1,JDIM(IGRID))
10CONTINUE
Q Files
The Plot3D Q file contains flow quantities for the grid points defined in the XYZ file. In order to properly use the Plot3D functions, the variables in the Q file must be properly normalized. The free-stream density is assumed to be 1. In other words, the Density (Q1) is assumed to have been divided by the free-stream density. Similarly, the free-stream speed of sound is assumed to be 1. Thus Momentum (Q2, Q3, Q4) must be normalized based upon these assumptions (it is assumed to have been divided by free-stream density and free-stream speed of sound).
The free-stream velocity magnitude is given by the free-stream mach number, FSMACH. The direction of the free stream velocity is computed from the angle of attack, ALPHA (in degrees). (The angle of attack is currently only used in the computation of the perturbation velocity vectors).
The other two constants: RE and TIME are not currently used.
In computing the Plot3D functions, the fluid is assumed to be air and behave as a perfect gas. The ratio of specific heats is assumed to be 1.4 and the gas constant is assumed to be 1.
Q File Single Grid
READ(IUNIT) IDIM,JDIM,KDIM
READ(IUNIT) FSMACH, ALPHA, RE, TIME
READ(IUNIT) ((((Q(I,J,K,NX), I=1,IDIM), J=1,JDIM), K=1,KDIM), NX=1,5)
2D Q File
READ(IUNIT) IDIM,JDIM
READ(IUNIT) FSMACH, ALPHA, RE, TIME
READ(IUNIT) (((Q(I,J, NX), I=1,IDIM), J=1,JDIM), NX=1,4)
Q File Multiple Grid
READ(IUNIT) NGRID
READ(IUNIT) (IDIM(IGRID),JDIM(IGRID),KDIM(IGRID), IGRID=1,NGRID)
DO 10 IGRID=1,NGRID
READ(IUNIT) FSMACH, ALPHA, RE, TIME
READ(IUNIT) ((((Q(I,J,K,NX), I=1,IDIM(IGRID)), J=1,JDIM(IGRID)), K=1,KDIM(IGRID)), NX=1,5)
10 CONTINUE
|