Fortran Coder

查看: 11228|回复: 8
打印 上一主题 下一主题

[文件读写] 无法保存dat文件

[复制链接]

35

帖子

17

主题

0

精华

熟手

F 币
136 元
贡献
240 点
楼主
发表于 2017-11-14 16:17:34 | 显示全部楼层
我测试过你的代码,本身写入.dat文件是不存在问题的,是可以生成输出文件的。
主要问题应该出现在你的数组中,可以试着分别将每个数组输出来看一下。

[Fortran] 纯文本查看 复制代码
    !< write data to TECPLOT file
    program printDat
    implicit none
    integer :: i,j,k
    integer :: ni
    integer :: nj
    integer :: nk
    real(kind=8) :: xp(5,5,5) 
    character :: text
  
    ni = 5
    nj = 5
    nk = 5
    
    open(unit=15,file='tecplot.dat',form='formatted',status='unknown')
    !< three-dimensional
    text='VARIABLES = X Y Z U V W P K E VIS'
    write(15,90) text
90  format(1x,70a)
      
    write(15,92) ni,nj,nk
92  format(1x,'ZONE T="Zone-One", I= ',i3,' ,J= ',i3,' ,K= ',i3,' ,F=POINT')
    
    do k=1,nk
        do j=1,nj
            do i=1,ni
!                write(15,98)xp(i,j,k),yp(i,j,k),zp(i,j,k),phi(i,j,k,u),
!                    &phi(i,j,k,v),phi(i,j,k,w),phi(i,j,k,p),phi(i,j,k,te),
!                    &phi(i,j,k,ed),vis(i,j,k)
                write(15,108) xp(i,j,k)
            enddo
        enddo
    enddo
98  format(1x,10(1pe13.6,2x))
108 format(1x,1(1pe13.6,2x))    
    
    close(15)
    endprogram
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-3 15:40

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表