Fortran将粒子数据输出到.vtu文件中进行截视

Question

我试图编写一个Fortran子例程来输出(随机定位)粒子数据到.vtu文件中，以便在Paraview中进行后处理。该程序生成一个文件，格式如下：

<?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文件出了问题吗？

对于每个粒子，标量图值设置为1，这可能是不必要的。我只是把它包括在内，以防我希望粒子在某一点上有不同的颜色(可能应该是一个浮子)。

这必须包含单元格数据吗？连通性和偏移量刚刚被设置为粒子/点数，而对于每个粒子/点，类型刚刚被设置为"1“。这会导致问题吗？

Answer 1

结果表明，上面的输出格式没有什么问题。问题在于“.”中使用的格式说明符。组件。其中一些是使用"*“指定的，而”*“必须生成不一致/不兼容的输出文件。通过为坐标指定这些(例如"//TRIM(no_dim)//"(4X，E15.8))，该文件与Paraview一起工作。对于感兴趣的人，还应该编写.pvd文件(即指向.pvu文件的文件)，以便在Paraview中生成动画。