对于你这样的问题，fortran77中较好的解决办法是用数据块子程序，通过公用块的方式实现。fortran的数组不同于其他语言的数组，它是列优先的，很容易对矩阵进行水平向扩展。
a(:,:4)=reshape([1,2,3,4,5,6,7,8],[2,4])
