Clefフォーマットのファイルをfortranで読み込む

open(20,file='fort.20',form='unformatted')
do
  read(20,iostat=ists) n,time
  if(ists /= 0) exit
  read(20) (((u(i,j,k),i=1,imax),j=1,jmax),k=1,kmax)
  read(20) (((v(i,j,k),i=1,imax),j=1,jmax),k=1,kmax)
  read(20) (((w(i,j,k),i=1,imax),j=1,jmax),k=1,kmax)
  read(20) (((p(i,j,k),i=1,imax),j=1,jmax),k=1,kmax)
! --- ここに，処理したい文を入れる．例えば，圧力を2倍したり，平均をとったりして，fort.21に書き出すなど．
end do
close(20)