Paraview で後処理するために、(ランダムに配置された) 粒子データを .vtu ファイルに出力する Fortran サブルーチンを作成しようとしています。プログラムは、次のようにフォーマットされたファイルを生成します。
<?xml version="1.0"?>
<VTKFile type= "UnstructuredGrid" version= "0.1" byte_order= "BigEndian">
<UnstructuredGrid>
<Piece NumberOfPoints="6664" NumberOfCells="0">
<PointData Scalars="Pressure" Vectors="Velocity">
<DataArray type="Float32" Name="Pressures" format="ascii">
...
</DataArray>
<DataArray type="Float32" Name="Density" format="ascii">
...
</DataArray>
<DataArray type="Float32" Name="Mass" format="ascii">
...
</DataArray>
<DataArray type="Float32" Name="Scalarplot" format="ascii">
...
</DataArray>
<DataArray type="Float32" Name="Velocity" NumberOfComponents="3" format="ascii">
... ... ...
</DataArray>
</PointData>
<Points>
<DataArray type="Float32" NumberOfComponents="3" format="ascii">
... ... ...
</DataArray>
</Points>
<Cells>
<DataArray type="Int32" Name="connectivity" format="ascii">
...
</DataArray>
<DataArray type="Int32" Name="offsets" format="ascii">
...
</DataArray>
<DataArray type="Int32" Name="types" format="ascii">
...
</DataArray>
</Cells>
</Piece>
</UnstructuredGrid>
</VTKFile>
コードは、 SPHYSICSプログラムのサブルーチンとファイル形式のドキュメントに基づいています。
Paraview はファイルを読み込みますが、座標またはフィールド変数を使用して粒子を視覚化できません。これは .vtu ファイルの問題が原因でしょうか?
scalarplot 値は各粒子に対して 1 に設定されており、おそらく不要です。ある時点でパーティクルの色を変えたい場合に備えて、これを含めました (おそらく Float にする必要があります)。
これにはセルデータが含まれている必要がありますか? 接続性とオフセットがパーティクル/ポイント番号に設定され、タイプが各パーティクル/ポイントの「1」に設定されました。これは問題につながる可能性がありますか?