| 查看: 4838 | 回复: 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 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
25