Как правильно вернуть из функции динамический массив/матрицу

По православному нужно так:

integer, allocatable :: m(:)

Но как у тебя тоже всё должно работать, просто более старый стандарт.

У тебя кстати вся функция как-то неправильно оформлена. И вообще это не функция, а процедура, похоже.

Я немного поменял:

module two_d
  use utility
  implicit none
  type map
     real, dimension(:), allocatable :: xs, hs
  end type map


contains
  function get_map() result(m)
    implicit none
    type(map) ::  m
    integer :: N, i
    character(len=10) :: file_name

    file_name = 'field2.txt'

    N = get_num_rows(file_name)

    allocate(m%xs(N))
    allocate(m%hs(N))

    open(unit = 1, file = file_name, status = 'old')
    do i = 1, N
       read (1, *) m%xs(i), m%hs(i)
    end do

    close(1)
    return
  end function get_map

end module

Если убрать dimension, то он начинает ругаться. Он правда и сейчас ругается на вызов и я еще разбираюсь почему:

gfortran -Wall -Wpedantic -Wextra utility.f90 two_d.f90 main.f90 
main.f90:4:20:

    4 |   print *, get_map()
      |                    1
Error: Data transfer element at (1) cannot have ALLOCATABLE components unless it is processed by a defined input/output procedure

У тебя кстати вся функция как-то неправильно оформлена

Примеры смотрю тут http://www.mathcs.emory.edu/~cheung/Courses/561/Syllabus/6-Fortran/

И вообще это не функция, а процедура, похоже.

Возможно, я еще разбираюсь чем оно отличается. В Raku они отличаются тем, что function умеет возвращать значение, а subroutine - нет

Он правда и сейчас ругается на вызов

Все, понял, оно же возвращает составной тип и не знает как его распечатать:

  type(map) :: x

  x = get_map()
  print *, x%xs
  print *, x%hs

ну так надо intent(out) как бы в объявление переменной. Там где type(map), intent(out) :: m

Еще если функция, а не subroutine, то нифига напечатать не получится.

И книгу сейчас еще вспомню посоветую. Хотя изначально кажется достаточно простой и красивый синтаксис, тонкостей много. Лучше книга (есть таковая).

Книга такая: Илья Чернов «Фортран сегодня». Очень, очень рекомендую к прочтению, потому как тонкостей весьма много и они повсюду.

https://sites.google.com/site/fortrantoday/books/book

Спасибо за книгу, обязательно прочитаю! Она как раз небольшая.

Да, тонкостей много, по примерам в интернете на stackoverflow не научишься)

Вижу многочисленные ошибки я. Используй процедуру, а не функцию. Не делай так

type map
     real, dimension(:), allocatable :: xs, hs
end type map

.....

allocate(m%xs(N))
allocate(m%hs(N))

Делай так

type map
     real :: xs, hs
end type map

.....

type(map), allocatable :: m(:)

allocate(m(N))

Мнение эксперта по фортрану с 20-ти летним стажем, так ничего в нём и не понимающего:

module two_d
  use utility
  implicit none
  type map
     real :: xs, hs
  end type map


contains
  subroutine get_map(m)
    type(map), allocatable, intent(out) ::  m(:)
    integer :: N, i
    character(len=10) :: file_name

    file_name = 'field2.txt'

    N = get_num_rows(file_name)

    allocate(m(N))

    open(unit = 1, file = file_name, status = 'old')
    do i = 1, N
       read (1, *) m(i)%xs, m(i)%hs
    end do

  end subroutine get_map

end module

Попробуй, должно вроде работать.

Но ведь если я буду делать allocate(m(N)) он не будет знать куда это N деть: компилятор не соврет

two_d.f90:20:13:

   20 |     allocate(m(N))
      |             1
Error: Allocate-object at (1) is neither a data pointer nor an allocatable variable

Ну либо я не знаю какой-то хак.

type(map), allocatable, intent(out) ::  m(:)
integer :: N

N=10 ! или чего там ещё

allocate(m(N))

УМВР

Спасибо за пример, сейчас гляну.

Да, у меня было в обьявлении типа real, allocatable :: xs(:), hs(:). Сделал как у вас и теперь все компилируется

Финкция в фортране – это синтакцический сахар. Та же процедура (subroutine) с ограничениями.

Главное, чтобы считалось правильно.