Fortran-默认情况下以列形式写入磁盘?

我正在编写一个简单的代码,将一些大型矩阵输出到磁盘,然后在matlab中读取。

1. 磁盘写入效率(寻找速度不太慢的东西)
2. 易于在matlab中读取

一

PROGRAM WriteDisk

character(80)  :: filename = ' '    
INTEGER :: indt
INTEGER :: ind1, n1 = 161
INTEGER :: ind2, n2=20
INTEGER :: ind3, n3=2
INTEGER :: ind4, n4=2
INTEGER :: ind5, n5=21
INTEGER :: ind6, n6=20
INTEGER :: ind7, n7=2
INTEGER :: ind8, n8=2

INTEGER :: dummy    
REAL, ALLOCATABLE :: m1(:,:,:,:,:,:,:,:,:)
ALLOCATE(m1(2,n1,n2,n3,n4,n5,n6,n7,n8))   


dummy = 1

do ind8 = 1,n8
do ind7 = 1,n7
do ind6 = 1,n6    
do ind5 = 1,n5
do ind4 = 1,n4        
do ind3 = 1,n3        
do ind2 = 1,n2      
do ind1=1,n1    

    m1(2,ind1,ind2,ind3,ind4,ind5,ind6,ind7,ind8) = dummy
    dummy = dummy + 1

end do
end do
end do
end do
end do
end do
end do
end do

indt = 1
write(filename,'(a,i0,a)')'PF_m1_',indt,'.txt'    
OPEN(UNIT=25,FILE=filename,STATUS='replace',ACTION='write')
WRITE(25, *)  m1(2,:,:,:,:,:,:,:,:)
CLOSE(UNIT=25)



END PROGRAM

上面的程序写入矩阵 m1 作为一个 4327680 x 5 。这使得在matlab中对其进行整形很麻烦(尽管完全可能),因为在matlab中,我需要执行以下操作:

n1 = 161;
n2 = 20;
n3 = 2;
n4 = 2;
n5 = 21;
n6 = 20;
n7 = 2;
n8 = 2;

m1 = load('PF_m1_1.txt'); %This is a two dimensional matrix that needs to be transposed and reshaped TWICE to get the original matrix
m1 = m1';
m1 = m1(:);
m1 = reshape(m1, n1,n2,n3,n4,n5,n6,n7,n8)

有没有办法把它写成一个向量,一个元素一个元素 m1(2,1,1,1,1,1,1,1,1) m1(2,2,1,1,1,1,1,1,1) 作为第二个元素, m1(2,end,end,end,end,end,end,end,end) 作为最后一个元素,等等?

.mat 文件