※3/18追記
- 大前提なはずのデータの宣言を間違っていたので修正．
- 上記に伴い各例修正
- コメントを参考に方法を一つ追加
1次元の配列xに対して，1行にdim個，line_num行出力する．
プログラム骨組み
program main
integer, parameter :: dim = 3, line_num = 51
integer :: x(dim*line_num)
open(11, file='data.dat', status='old')
read(11,*) x
close(11)
!
! 主にここに追記
!
end program main
要素をdim個改行をせずに出力し，その後に改行する方法
二重ループを用いて
do i=0,line_num-1
do l = 1, dim
write(6, '(I4)', advance='no') x(i*dim+l)
end do
write(6, '()')
end do
無理矢理反復を使う方法（これより4つ目の方法を推奨）
反復を用いて
do i=0,line_num-1
write(6,'(1000I4)') x(i*dim+1:(i+1)*dim)
end do
反復の指定をとりあえず大きい値で指定してます．
実行時間的にも一つ目よりこっちのが早い．当然，dimが1000個を超えるとダメ．
書式を気にしないなら
結局
do i=0,line_num-1
write(6,*) x(i*dim+1:(i+1)*dim)
end do
これでいい．dimの制約もないし，こっちのが早い．
※コメントを参考に追記
出力formatをcharacterで定義する方法
formatをcharacterとして定義
character :: cdim*20 !これは宣言に
write(cdim, *) dim
cdim = '('//trim(adjustl(cdim))//'I4)'
write(6, cdim) x
これが一番よさげかな…？
ペーペーにはこんなもんですが，他にいい方法ないですかね？（反復に変数使う方法，調べても出てこない…）