将未知大小的矩阵读取到Fortran

Question

我想知道如何在矩阵尺寸未知的情况下将文件中的矩阵读取到Fortran 90。我知道如何获取矩阵的行数，但不知道如何知道列数以在程序内分配矩阵。 可以说矩阵就是这样

1 2 3 4 50
2 50 2 50 1
21 3 15 515 7 

我得到这个是为了获取外部文件的行数：

! GETTING NUMBER OF ROWS
open(unit=11, file=esmfile)
n=0
DO
  READ(11, *, IOSTAT=iostatus) dummy
  IF (iostatus < 0) EXIT
   n = n + 1
END DO