[Paraview] Error in MultiBlockPLOT3DReader
Richard GRENON
richard.grenon at onera.fr
Wed Dec 19 09:33:50 EST 2007
Hello,
I discovered Paraview 3.2.1 recently, Nice software.
I am using the MultibockPLOT3DReader for my data and I have found an
error when reading Fortran Unformatted Plot3D files, that is reading
with the options "Binary File", "Multi Grid" and "Has Byte Count" set
on: when reading MultiGrid files with 5 grids, I have found that
Paraview computed wrong min and max values of density, momentum
components and stagnation energy. In particular, the min value of
density was set to zero, athough I knew it should be 0.5857 for my data,
and computing other Plot3d functions failed.
My Fortran Unformatted PLOT3D files were tested without problem in
Fieldview, so I don't think they are corrupted. And also, I have
verified with a C program that the true min and max values were not the
values computed by Paraview. But all became OK in Paraview after I
converted my data in pure Binary files without Byte Count, so the
problem in Paraview is really when reading Fortran Unformatted PLOT3D
files that require "Binary File", "Multi Grid" and "Has Byte Count" to
be set on together.
I looked at the source and I have found that the MultiBlockPLOT3DReader
skips the four leading "Byte Count" bytes before reading all density,
momentum and stagnation energy values in a single record, but it does
not skip again the trailing "Byte Count" bytes after reading the data
record as it should do! That does not matter if you have a single grid
file, but when you have several grids, data are shifted in the arrays
each you read a new grid. The problems is only when reading the "q" file
in the loop on i from 0 to numBlocks starting at line 943 in the file
"vtkMultiBlockPLOT3DReader.cxx". I think that you should add
"this->SkipByteCount(qFp)" at the end of the loop just after line 1059.
I hope that developers can check if I am right, and maybe include the
corrected reader in a future Paraview release?
Best regards.
--
Richard GRENON
ONERA
Departement d'Aerodynamique Appliquee - DAAP/ACI
8 rue des Vertugadins
92190 MEUDON - FRANCE
phone : +33 1 46 73 42 17
fax : +33 1 46 73 41 46
mailto:Richard.Grenon at onera.fr
http://www.onera.fr
More information about the ParaView
mailing list