いわて駐在研究日誌2。

NEVER STAND BEHIND ME

Fortranでのデータ書き出し(マルチブロック、kの次元2としてあるのは、OpenFoamで2D計算の予定のためです)

注意! データの並びは「右手系」で!

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)