1

我有一个 Fortran 例程，它打开很多文本文件，从时间循环中写入数据。此例程open与newunit选项一起使用，此单元存储在一个对象中，以便以后将内容写入文件中。这在大多数情况下都可以正常工作，但是当程序需要同时打开大量 N 个文件时，我收到以下错误：

**forrtl: severe (104): incorrect STATUS= specifier value for connected file, unit -1, file CONOUT$**

指子程序中的第一个open函数createFiles。无论文件是否已存在，都会发生此错误。我不知道这是否有帮助，但在这个阶段应该生成的新单位是-32768.

我包含一个带有“timeSeries”类的最小代码示例，其中包括一个创建两个文件的例程：

  • 第一个文件fileName1在里面写东西后直接打开和关闭
  • 第二个文件fileName2保持打开状态，以便写入稍后在时间循环中计算的内容，并在时间循环结束时关闭

该示例由以下两个文件组成。它因 i=32639 而中断。

main.f90 ：

program writeFiles

use TS

logical :: stat
integer :: i, istep, N, NtimeSteps
character(len=16) :: fileName1, fileName2
character(len=300) :: path
type(timeSeries), dimension(:), allocatable :: myTS

call getcwd( path )

path = trim(path) // '\Output_files'
inquire(directory = trim(path), exist = stat )
if (.not. stat) call system("mkdir " // '"' // trim(path) // '"' )

N = 50000
NtimeSteps = 100
allocate(myTS(N))

do i = 1, N
    write(fileName1,'(a6,i6.6,a4)') 'file1_', i, '.txt'
    write(fileName2,'(a6,i6.6,a4)') 'file2_', i, '.txt'
    call myTS(i)%createFiles(trim(path),fileName1,fileName2)
end do

do istep = 1, NtimeSteps
    #
    #compute stuff
    #
    do i = 1, N 
         write(myTS(i)%fileUnit,*) 'stuff'
    end do
end do

do i = 1, N
    close(myTS(i)%fileUnit)
end do

end program writeFiles

模块.f90：

module TS

    type timeSeries
    integer :: fileUnit

    contains

    procedure :: createFiles       => timeSeries_createFiles

    end type timeSeries

    contains

    subroutine timeSeries_createFiles(this,dir,fileName1,fileName2)
        class(timeSeries) :: this
        character(*) :: dir, fileName1, fileName2 

        open(newunit = this%fileUnit , file = dir // '\' // fileName1, status = 'replace') !error occurs here after multiple function calls
        write(this%fileUnit,*) 'Write stuff'
        close(this%fileUnit)
        open(newunit = this%fileUnit , file = dir // '\' // fileName2, status = 'replace') 
    end subroutine timeSeries_createFiles

end module

关于这个错误的原因的任何想法？同时打开的文件数量有限制吗？它可能与内存问题有关吗？

我正在使用 Intel(R) Visual Fortran 编译器 17.0.4.210

4

1 回答 1

3

Windows 有一个有趣的习惯，即在您关闭文件后短时间内不会释放已关闭文件的所有资源。几十年来，我断断续续地看到过这种问题。我通常的建议是，当您打算在同一文件上再次打开时，在关闭后半秒内调用 SLEEPQQ。但你不是在这里这样做。

这里还有更多令人费解的地方。当打开显式文件并使用 NEWUNIT 时，不应出现引用单元 -1 和 CONOUT$ 的错误消息。在 Intel 的实施中，NEWUNIT 数字从 -129 开始，从那里开始变得更负。单位 -1 用于 PRINT 或 WRITE(*)，CONOUT$ 是控制台。STATUS='REPLACE' 对于连接到控制台的设备无效。newunit 编号为 -32768 说明并建议 Intel 库中 NEWUNIT 的内部限制。

我自己做了一个测试，发现如果你使用 NEWUNIT 并关闭单元，单元号会在循环回到 -129 之前低至 -16384。如果您确实要关闭单元，那没关系，但您永远不会关闭您打开的第二个文件，因此您至少要达到打开的 NEWUNIT 文件的最大数量。我建议找出一种不同的方法来解决不需要打开数千个文件的问题。

于 2018-03-02T16:58:16.030 回答