Чтение двоичного файла с помощью Julia

Предположим, что полный код F90 для записи данных выглядит следующим образом:

program SaveGrid
    implicit none

    integer :: nv
    integer :: Nx, Ny, Nzt
    real, dimension(:), allocatable :: xx, yy, zz
    integer :: i

    nv = 3
    Nx = 5
    Ny = 7
    Nzt = 10

    allocate(xx(Nx), yy(Ny), zz(Nzt))

    do i = 1, Nx
        xx(i) = float(i)
    end do
    do i = 1, Ny
        yy(i) = float(i) * 2.0
    end do
    do i = 1, Nzt
        zz(i) = float(i) * 3.0
    end do

    open(unit=11, file = "ring.grid", form = "unformatted")

    write(11) nv
    write(11) Nx, Ny, Nzt
    write(11) xx
    write(11) yy
    write(11) zz

    close(11)

    deallocate(xx, yy, zz)
end program SaveGrid

Чтобы прочитать это в Julia, вам нужно обработать способ записи неформатированных данных F90.

function read_fortran_binary(filename)
    open(filename, "r") do file
        # Read record marker.
        record_size = read(file, Int32)
        # Read data.
        nv = read(file, Int32)
        # Read record marker.
        end_record_size = read(file, Int32)

        # Read record marker.
        read(file, Int32)
        # Read data.
        Nx = read(file, Int32)
        Ny = read(file, Int32)
        Nzt = read(file, Int32)
        # Read record marker.
        read(file, Int32)

        return nv, Nx, Ny, Nzt
    end
end

nv, Nx, Ny, Nzt = read_fortran_binary("ring.grid")

println("nv:  ", nv)
println("Nx:  ", Nx)
println("Ny:  ", Ny)
println("Nzt: ", Nzt)

Маркеры записей имеются как до, так и после каждого фрагмента данных, записанного в файл. Вам тоже нужно прочитать это. Очевидно, вам не нужно игнорировать эти данные, и размеры записей можно использовать, чтобы сделать ваш код Julia более надежным и гибким за счет изменений в неформатированных данных F90.

nv:  3
Nx:  5
Ny:  7
Nzt: 10