在 Fortran 中读取二进制文件

Question

我正在尝试读取一个由带符号的 16 位整数组成的二进制文件，其中正好有 51840000 个。 C 中完成此操作的代码如下所示：

#include<stdio.h>
#include<stdlib.h>
#include<stdint.h>

int main()
{
    int16_t *arr = malloc(51840000*sizeof(int16_t));
    FILE *fp;
    fp = fopen("LDEM_45N_400M.IMG", "rb");
    if(fp == NULL)
    {
        printf("Error opening file\n");
        exit(1);
    }
    printf("Testing fread() function: \n\n");
    fread(arr, sizeof(*arr), 51840000, fp);
    fclose(fp);
    printf("%d \n", arr[51840000-2]);
    free(arr);
    return 0;
}

我如何在 Fortran 中读取这样的文件？ Fortran中的IO对我来说一直很神秘

Answer 1

使用access="stream"，同样可以阅读。将您的整数声明为 integer(int16) 并使用 iso_fortran_env 模块。

use iso_fortran_env
integer :: ierr, n = 51840000
integer(int16) :: arr(n)

open(newunit=iu,file="LDEM_45N_400M.IMG", access="stream", status="old", action="read",iostat=ierr)
if (ierr/=0) stop "Error opening the file."
read(iu, iostat=ierr) arr
if (ierr/=0) stop "Error reading the array."
close(iu)

在 Fortran 中读取二进制文件

reading a binary file in Fortran

fortran

gfortran

fortran90

intel-fortran