reading matrices in FORTRAN