You can do that using reshape and shape intrinsics. Something like: INTEGER, DIMENSION(3, 3) :: array array = reshape((/ 1, 2, 3, 4, 5, 6, 7, 8, 9 /), shape(array)) But remember the column-major order. The array will be 1 4 7 2 5 8 3 6 9 after reshaping. So to get: 1 2 … Read more
Assumed shape dummy arguments (those with (:)) require explicit interface to the procedure to be available at the call site. That means the calling code must know how exactly the subroutine header looks like. See also Module calling an external procedure with implicit interface That explicit interface can be provided in several ways 1. preferred … Read more
A real which isn’t marked as double precision will be assumed to be single precision. Just because sometime later you assign it to a double precision variable, or convert it to double precision, that doesn’t mean that the value will ‘magically’ be double precision. It doesn’t look ahead to see how the value will be … Read more
Modules are always the right thing to use 😉 If you have a very simple F90 program you can include functions and subroutines in the ‘contains’ block: program simple implicit none integer :: x, y x = … y = myfunc(x) contains function myfunc(x) result(y) implicit none integer, intent(in) :: x integer :: y … … Read more