注意! データの並びは「右手系」で!
open(7,form='formatted', file='P3D.xyz')
write(7,*) 2 ! # of blocks (ここでは2個)
write(7,*) nout+1, nrad+1, 2 ! block 1 の次元
write(7,*) nout+1, nrad+1, 2 ! block 2 の次元
! Block 1
write(7,*) &
((( xxxu(i,j), i=0,nout), j=0,nrad), k=1,2), &
((( yyyu(i,j), i=0,nout), j=0,nrad), k=1,2), &
((( zzz(k) , i=0,nout), j=0,nrad), k=1,2)
WRITE(7,*) ' ' ! ブロック毎に区切りの空白をいれておく
! Block 2
write(7,*) &
((( xxxl(i,j), i=0,nout), j=0,nrad), k=1,2), &
((( yyyl(i,j), i=0,nout), j=0,nrad), k=1,2), &
((( zzz(k) , i=0,nout), j=0,nrad), k=1,2)
close(7)