Чтение из файла двумерного массива вещественных чисел

 program prog
 
      implicit none
 
      ! Variables
      INTEGER N,I,J,temp 
 
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: A !Массив, размерность которого заранее неизвестна
          
      ! Body
      write (*,*) "------------Program Starts-----------"
 
     OPEN(5,FILE='input.txt')
     READ(5,*)N
 
    ALLOCATE(A(N,N));
 
     DO 10 I=1,N
        READ(5,20)(A(I,J),J=1,N)
       10 CONTINUE
       20 FORMAT(D9.2) !Хочется считывать числа, скорее всего они будут небольшие - до 9 знаков точно
     CLOSE(5)
    
      read (*,*) temp
      
      end program prog

program read_from_file
! описание переменных
implicit none
integer j,k !индексы элементов массива
integer N ! размер массива
double precision, allocatable :: array(:,:) ! двумерный массив чисел двойной точности с автоматическим выделение памяти под него 
! тело программы
write(*,*) 'The Program starts.'
open( UNIT=1, file= 'indata.dat', STATUS='old', ACTION='read' ) ! прикрепление файла к устройству
read(1,*) N ! чтение размера массива
allocate(array(N,N)) ! выделение памяти под массив
do j = 1, N, 1  ! пробегаем по строкам
    read(1,*,end=200) (array(j,k),k=1,N) ! пробегаем по стобцам 
end do ! for j
200 j=j-1 ! уличная магия для устранения проблем при чтении последней строки в файле 
! вывод на экран для проверки
do j = 1, N, 1
    do k = 1, N, 1
        write(*,002) array(j,k)
    end do !for k
    write(*,*)
end do ! for j
deallocate(array) ! освобождение памяти занимаемой массивом
close(UNIT=1, STATUS='keep') ! открепление файла от устройства 1
002 format(f12.2,\) !ПРИ ТАКОМ ФОРМАТЕ ВЫВОДА ОН ОКРУГЛЯЕТ ДЕСЯТИЧНУЮ ЧАСТЬ!
!-----------------------------------------------------------------------
write(*,*) 'The Program is complete.'
stop
end program read_from_file