将文件中的数据读入二维数组

Question

我对 fortran 和编程完全陌生（经验不足）。我很难为以下内容编写代码。

我需要使用两个输入文件中的数据将数据存储在二维数组中 然后我需要确定两个文件的大小 然后需要从两个维度中选择一个变量的特定值的变量 然后计算平均值 最后将结果写入一个新的数据文件。

例如：

• 有两个无线电探空仪数据文件，不同数据大小的压力、高度、温度、露点值
• 需要在 1000、850 和 500 的特定压力下获取高度、温度和露点的平均值（两个文件都有这些记录）并写入新文件。

我已经阅读了几本书，但无法掌握我正在寻找的确切内容。你们中的任何人都可以给我一些关于编码的见解吗？我真的很感谢你们。

Answer 1

您的代码将如下所示：

program main
    implicit none

    real::height(100)
    integer::i

    open(unit=16, file='height.dat')

    do i = 1, 100
        read(16, *) height(i)
    enddo

    close(16)

    open(unit=16, file='mean_height.dat')

    write(16, *) sum(height)/100

    close(16)

end program

这假设有一个名为height.dat 的文件，每行有一个高度值，总共100 行。 它读取 100 个高度，然后将平均高度写入文件 mean_height.dat。 您显然必须对此进行重大修改。 例如，您不知道文件中有 100 行，因此您需要为读取的值使用 allocatable 数组：

program main
    implicit none

    real,allocatable::height(:)
    integer::i, N

    !Code that figures out how many lines there are
    !...

    allocate(height(N))

    open(unit=16, file='height.dat')

    do i = 1, N

这应该可以帮助您入门...理想情况下，您可以从比您描述的更简单的输入文件开始学习基础知识。