关闭输入文件时出现Fortran分段错误

时间:2018-07-11 19:32:38

标签: io fortran modeling numerical

希望有人可以帮助我解决这个问题。将不胜感激。

背景。我正在为我们的研究小组扩展一些传统代码,该代码可以进行一些热传递,热对流等类型的数值建模。我正在尝试做的是读取温度输入文件,以便可以在模型随时间推移时调整模型的温度。虽然该信息基本上无关紧要。尝试关闭我打开并正在读取的文件时,我不断遇到分段错误。我已在下面附上了代码以及典型的温度输入文件的外观(但由于将有成千上万个数据点而大大简化了)。在消除分段错误和/或更正代码以使其更有效方面的任何帮助将不胜感激。

subroutine Read_Temp_Input(Temp_in)!(dt)!(filename)!,dt,Temp_out)
implicit none

character*80                                :: filename
integer, parameter                          :: dp = selected_real_kind(15)
real(kind=dp), allocatable, dimension(:,:)  :: Temp_in, Temp_out
integer                                     :: i,j,ierror,n
real (kind=dp)                              :: in_time, in_temp, dt, inter_time, inter_temp, max_time
character*80                                :: line


!Variable definitions
!in_time is the time given from the input file
!in_temp is the temp given from the input file
!inter_time is the interpolated time using dt, max_time, and linear interpolation
!inter_temp is the interpolated temp using dt, max_time, and linear interpolation
!max_time is the largest time value in the input file
!dt is the time step utilized in the simulation, provided by main routine

!read in the temperature input file
filename = "temps.txt"
Open(UNIT=1,FILE=filename,ACTION="Read")

!Determine the size of the necessary allocation
n = 0
do
    Read(1,*,iostat=ierror) in_time, in_temp
    if (ierror.ne.0) then
        exit
    else
        print*, in_time, in_temp
        n = n + 1
    endif
enddo

!Allocate the Temp_in array to a n x 2 array.
allocate(Temp_in(n,2))

close(unit=1)
Open(UNIT=1,FILE=filename,ACTION="Read")

n = 0

do
    Read(1,*,iostat=ierror) in_time, in_temp
    if (n.ne.0) then
        if (ierror.ne.0) then
            exit
        else
            Temp_in(n-1,0) = in_time
            Temp_in(n-1,1) = in_temp
        endif
    endif
    n = n + 1
enddo

dt = 0.5_dp
print*, 'is the fault before here?'
close(1)
print*, 'is the fault here?'
end subroutine

温度文件

 Time  Temp
 1     300
 2     400
 3     500
 5     600

附加也是我一直得到的输出。我会及时摆脱随机打印语句,只是将它们用作测试代码中断位置的一种方式。

enter image description here

非常感谢您的帮助。

1 个答案:

答案 0 :(得分:4)

默认情况下,Fortran数组是基于1的,因此示例中的Temp_in也是如此。这样的线

Temp_in(n-1,0) = in_time
Temp_in(n-1,1) = in_temp

将超出范围(对于n=1,两者都超出范围,并且第一个始终超出范围)。

您可以使用--check=bounds标志使编译器为您检查该标志,如果尝试以这种方式访问​​数组,则它将引发运行时错误。

因此,将Temp_in分配为

 allocate(Temp_in(0:n-1,0:1))

或使用基于1的索引作为解决方案。


更新

此外,在确定条目数时,读取示例temps.txt文件的第一行将失败,因为它尝试读取两个real,但会发现其他内容。因此,您需要在打开文件以在第一行中读取后进行虚拟读取(或使用与第二个循环中相同的附加检查)

Open(UNIT=1,FILE=filename,ACTION="Read")
! read in the first line and throw it away
read(1,*)
! go on