| 查看: 4849 | 回复: 3 | ||
[求助]
fortran程序中输出tecplot二进制文件
|
|
下面这段程序是我在某大神德基础上修改的binary二进制文件的输出程序，得到的文件会提示“invalid header in the binary datafile！” （见附件“捕获.JPG”） 程序改了好多遍，都没解决问题。希望各位前辈，师兄指点一下。 "捕获2.JPG"是我想输出的三维图形，就是一个网格，没有其他的东西。 f90文件是我的程序节选。 open(99,file="tecbin.plt",form="binary"  !----1.1 Magic number, Version number write(99) "#!TDV101" !---1.2. Integer value of 1.---------------------------------------------------------- ! +-----------+ ! | INT32 | This is used to determine the byte order ! +-----------+ of the reader relative to the writer. write(99) 1 !---1.3. Title and variable names.------------------------------------------------- !---1.3.1. The TITLE. (See note 1.) write(99) ICHAR('A') write(99) 0 !---1.3.2 Number of variables (NumVar) in the datafile. write(99) 3 !----1.3.3 Variable names. N = L[1] + L[2] + .... L[NumVar] write(99) ichar('X') ! name of variables (NumVar) in the datafile write(99) 0 write(99) ichar('Y') write(99) 0 write(99) ichar('Z') write(99) 0 !---1.4. Zones------------------------------------------------------------------- !---1.4.1 Zone 1+++++++++++++++++++++++++++++++++++++++++++ write(99) 299.0 ! Zone marker. Value = 299.0 write(99) ichar('Z') ! Zone name. (See note 1.) N = length of zone name + 1. write(99) ichar('O') write(99) ichar('N') write(99) ichar('E') write(99) ichar('0') write(99) ichar('1') write(99) 0 write(99) -1 !-------Zone Color (set to -1 if you want tecplot todetermine). write(99) 0 !-------ZoneType 0=ORDERED,1=FELINESEG,2=FETRIANGLE,3=FEQUADRILATERAL,4=FETETRAHEDRON,5=FEBRICK write(99) 1 !-------DataPacking 0=Block, 1=Point write(99) 0 !-------Specify Var Location. 0 = Don't specify, all data is located at the nodes. 1 = Specify write(99) 0 !-------Number of user defined face neighbor connections (value >= 0) write(99) IMax !-------IMax,JMax,KMax write(99) JMax write(99) KMax write(99) 0 ! 1=Auxiliary name/value pair to follow 0=No more Auxiliar name/value pairs ! 0=No more Auxiliar name/value pairs !---1.4.2 Zone 2+++++++++++++++++++++++++++++++++++++++++++ write(99) 299.0 ! Zone marker. Value = 299.0 write(99) ichar('Z') ! Zone name. (See note 1.) N = length of zone name + 1. write(99) ichar('O') write(99) ichar('N') write(99) ichar('E') write(99) ichar('0') write(99) ichar('1') write(99) 0 write(99) -1 !-------Zone Color (set to -1 if you want tecplot todetermine). write(99) 0 !-------ZoneType 0=ORDERED,1=FELINESEG,2=FETRIANGLE,3=FEQUADRILATERAL,4=FETETRAHEDRON,5=FEBRICK write(99) 1 !-------DataPacking 0=Block, 1=Point write(99) 0 !-------Specify Var Location. 0 = Don't specify, all data is located at the nodes. 1 = Specify write(99) 0 !-------Number of user defined face neighbor connections (value >= 0) write(99) IMax !-------IMax=10,JMax=10,KMax=30. the same value as follows write(99) JMax write(99) KMax write(99) 0 ! 1=Auxiliary name/value pair to follow 0=No more Auxiliar name/value pairs ! 0=No more Auxiliar name/value pairs !---1.4.3 Zone 3+++++++++++++++++++++++++++++++++++++++++++ write(99) 299.0 ! Zone marker. Value = 299.0 write(99) ichar('Z') ! Zone name. (See note 1.) N = length of zone name + 1. write(99) ichar('O') write(99) ichar('N') write(99) ichar('E') write(99) ichar('0') write(99) ichar('1') write(99) 0 write(99) -1 !-------Zone Color (set to -1 if you want tecplot todetermine). write(99) 0 !-------ZoneType 0=ORDERED,1=FELINESEG,2=FETRIANGLE,3=FEQUADRILATERAL,4=FETETRAHEDRON,5=FEBRICK write(99) 1 !-------DataPacking 0=Block, 1=Point write(99) 0 !-------Specify Var Location. 0 = Don't specify, all data is located at the nodes. 1 = Specify write(99) 0 !-------Number of user defined face neighbor connections (value >= 0) write(99) IMax !-------IMax,JMax,KMax write(99) JMax write(99) KMax write(99) 0 ! 1=Auxiliary name/value pair to follow 0=No more Auxiliar name/value pairs ! 0=No more Auxiliar name/value pairs !---2---------data section----------------------------- write(99) 357.0 ! EOHMARKER, value=357.0 !---2.1------zone 1++++++++++++++++++ write(99) 299.0 ! Zone marker Value = 299.0 write(99) 2 ! variable data format, N=Total number of vars write(99) 2 ! 1=Float, 2=Double, 3=LongInt, 4=ShortInt, 5=Byte, 6=Bit write(99) 2 write(99) 0 ! Has passive variables: 0 = no, 1 = yes. write(99) -1 ! Zero based zone number to share connectivity list with (-1 = no sharing) do k=1,kmax do j=1,jmax do i=1,imax write(99) xg(i,j,k,1),yg(i,j,k,1),zg(i,j,k,1) !--------Zone Data. Each variable is in data format asspecified above. end do end do end do !---2.2------zone 2 +++++++++++++++++ write(99) 299.0 ! Zone marker Value = 299.0 write(99) 2 ! variable data format, N=Total number of vars write(99) 2 ! 1=Float, 2=Double, 3=LongInt, 4=ShortInt, 5=Byte, 6=Bit write(99) 2 write(99) 0 ! Has passive variables: 0 = no, 1 = yes. write(99) -1 ! Zero based zone number to share connectivity list with (-1 = no sharing) do k=1,kmax do j=1,jmax do i=1,imax write(99) xg(i,j,k,2),yg(i,j,k,2),zg(i,j,k,2) !--------Zone Data. Each variable is in data format asspecified above. end do end do end do !---2.3------zone 3 +++++++++++++ write(99) 299.0 ! Zone marker Value = 299.0 write(99) 2 ! variable data format, N=Total number of vars write(99) 2 ! 1=Float, 2=Double, 3=LongInt, 4=ShortInt, 5=Byte, 6=Bit write(99) 2 write(99) 0 ! Has passive variables: 0 = no, 1 = yes. write(99) -1 ! Zero based zone number to share connectivity list with (-1 = no sharing) do k=1,kmax do j=1,jmax do i=1,imax write(99) xg(i,j,k,3),yg(i,j,k,3),zg(i,j,k,3) !--------Zone Data. Each variable is in data format asspecified above. end do end do end do close(99)  捕获.JPG  捕获2.JPG |
回复此楼» 猜你喜欢
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fortran中最大整数问题
已经有7人回复
- fortran程序中write输出格式有问题
已经有5人回复
- FORTRAN中不定个数传递参数
已经有5人回复
- Fortran输出二进制数据，用matlab读出，结果出错，求助问题何在？
已经有4人回复
- fortran程序运行结果输出及绘图
已经有9人回复
- fortran编译中$符号的使用
已经有4人回复
- fortran中对变量取值
已经有3人回复
- 自编CFD程序怎样做后处理
已经有22人回复
- fortran中运行中出现无穷大的数怎么处理
已经有9人回复
- 用Fortran编写程序输出的数据一般用什么程序作图比较好？
已经有28人回复
- fortran文件操作，接着写不覆盖
已经有3人回复
- FORTRAN中如何实现重复读取一行？
已经有24人回复
- fortran输出求助
已经有6人回复
- fortran中到底该如何比较两个实数是否相等
已经有8人回复
- 请教关于fortran的ex后缀可执行文件，谢谢
已经有3人回复
- fortran两个问题，文件中的空格，屏幕同一行覆盖输出
已经有23人回复
- 【求助】FORTRAN写文件时每行开头都多了一个空格？
已经有6人回复
- 【求助】Fortran编译目标文件后的链接问题
已经有4人回复
- 【求助】fortran中怎么存储不确定长度的数据【已完结】
已经有6人回复
- 【求助】Fortran程序，如何让exe文件可以单独运行？
已经有6人回复
- 【求助】matlab读取fortran输出的txt文件【已解决】
已经有3人回复
2楼2014-07-21 23:44:43
3楼2014-07-22 08:56:47
4楼2014-07-22 11:30:17
投票:
30