Fortran 代码仅在运行 python 时才会出错

Question

我有一个 python 代码（我们称之为'pycode.py'），它创建一个.txt 文件（'fortc.txt'）。 该文件需要通过fortran90代码生成的a.exe（'f90code.exe'）读取。 fortran 代码然后执行一些操作并将结果写入另一个.txt 文件（'res_tc.txt'），之前通过 pycode.py 创建，pycode.py 最终将读取和打印。 这是 python 代码：

import csv
import sys, os
import numpy as np

with open('fortc.txt', 'w') as g:
    writer = csv.writer(g, delimiter=" ", quoting=csv.QUOTE_NONE, escapechar=' ')
    writer.writerow(['1', '2', '3', '4'])

    f = open('res_tc.txt','w+')
    os.system('./f90code.exe')
    readres = np.loadtxt('res_tc.txt')
    tc = float(readres)
    print(tc)
    f.close()

如果我手动创建文件“fortc.txt”，然后只运行 f90 代码，而不通过 python，一切正常，我得到了预期的结果。 但是，如果我通过 python 运行所有内容，则会正确创建“fortc.txt”，但显然 fortran 无法正确读取它。 事实上，“fortc.txt”中的所有内容都被读取为“0.0000000000000”，执行操作时会产生以下错误：

Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DENORMAL

显然我得到了一个“NaN”。 此外，看起来 fortran 没有正确运行 do 循环，因为如果我尝试在 do 循环内打印某些内容，我不会打印出任何内容。 以下是读取文件的 f90 代码部分：

program tc

      real*8, allocatable:: ne(:), kt(:), abund(:), z(:)
      integer:: i, j, ndati, n
      data ab/0.,0.1,0.31623,1./

      open(96,file='fortc.txt',status='old')

      n=0
      do
        read(96,*,end=320)
        n=n+1
      end do

      320 continue
      ndat=n

      allocate(ne(ndat), kt(ndat), abund(ndat), z(ndat))

      rewind(96)

      do j=1, ndat
         read(96,*) kt(j), abund(j), ne(j), z(j)
         print*, kt(j)
      end do
end program

f90 结束时的打印调用没有结果，只有上面显示的错误和“NaN”（源自 python 代码中的打印）。 另一方面，如果将 print 语句带到 do 循环之外，我会打印一些东西（'0.00000000'，如上所述）。 同样，只有在通过 python 运行编译的 fortran 代码时才会发生这种情况，而如果我只是从终端运行.exe，它会按预期工作。

do循环的处理方式有什么问题吗？ .txt 文件中的数字没有正确读取，并且没有执行打印，这表明我可能存在问题。 但是，我无法修复它。

Answer 1

当您调用 fortran 代码时，您的输入文件为空：

with open('fortc.txt', 'w') as g:

默认情况下，Python 将刷新写入此文件的数据，并在退出with正文时将其关闭。 但是您在其中包含了对 fortran 代码的调用。

如果您想使用with进行文件操作，请以这种方式重新排列您的代码：

import csv
import sys, os
import numpy as np

with open('fortc.txt', 'w') as g:
    writer = csv.writer(g, delimiter=" ", quoting=csv.QUOTE_NONE, escapechar=' ')
    writer.writerow(['1', '2', '3', '4'])

f = open('res_tc.txt','w+')   # not sure why the frotran program would output to "f" just because it is open.  If it does, ok, it will wrk.
os.system('./f90code.exe')
readres = np.loadtxt('res_tc.txt')
tc = float(readres)
print(tc)
f.close()  # you can use explicit open and close calls for fortc.txt as well