将数组变矩阵

li913

   do r=1,6
            STR(k,r)=strain(k+(r-1)*nblock)
   end do
没问题，也可以更简单：
str = reshape(strain,[nblck,6])