Skip to content

S
SWE
Commit a83ba300 authored by Sebastian Rettenberger's avatar Sebastian Rettenberger
Browse files
Options

Change vtk writer to behave like netcdf writer

parent 2746fc63
Hide whitespace changes
Inline Side-by-side
Showing with 11 additions and 11 deletions
src/writer/VtkWriter.cpp
View file @ a83ba300
...@@ -63,7 +63,7 @@ void io::VtkWriter::writeTimeStep( ...@@ -63,7 +63,7 @@ void io::VtkWriter::writeTimeStep(
std::ofstream vtkFile(generateFileName().c_str()); std::ofstream vtkFile(generateFileName().c_str());
assert(vtkFile.good()); assert(vtkFile.good());
// VTK HEADER // VTK header
vtkFile << "<?xml version=\"1.0\"?>" << std::endl vtkFile << "<?xml version=\"1.0\"?>" << std::endl
<< "<VTKFile type=\"StructuredGrid\">" << std::endl << "<VTKFile type=\"StructuredGrid\">" << std::endl
<< "<StructuredGrid WholeExtent=\"" << offsetX << " " << offsetX+nX << "<StructuredGrid WholeExtent=\"" << offsetX << " " << offsetX+nX
...@@ -74,7 +74,7 @@ void io::VtkWriter::writeTimeStep( ...@@ -74,7 +74,7 @@ void io::VtkWriter::writeTimeStep(
vtkFile << "<Points>" << std::endl vtkFile << "<Points>" << std::endl
<< "<DataArray NumberOfComponents=\"3\" type=\"Float32\" format=\"ascii\">" << std::endl; << "<DataArray NumberOfComponents=\"3\" type=\"Float32\" format=\"ascii\">" << std::endl;
//GITTER PUNKTE //Grid points
for (int j=0; j < nX+1; j++) for (int j=0; j < nX+1; j++)
for (int i=0; i < nX+1; i++) for (int i=0; i < nX+1; i++)
vtkFile << (offsetX+i)*dX << " " << (offsetY+j)*dY <<" 0" << std::endl; vtkFile << (offsetX+i)*dX << " " << (offsetY+j)*dY <<" 0" << std::endl;
...@@ -84,28 +84,28 @@ void io::VtkWriter::writeTimeStep( ...@@ -84,28 +84,28 @@ void io::VtkWriter::writeTimeStep(
vtkFile << "<CellData>" << std::endl; vtkFile << "<CellData>" << std::endl;
// Water surface height (h+b) // Water surface height h
vtkFile << "<DataArray Name=\"H\" type=\"Float32\" format=\"ascii\">" << std::endl; vtkFile << "<DataArray Name=\"h\" type=\"Float32\" format=\"ascii\">" << std::endl;
for (int j=1; j < nY+1; j++) for (int j=1; j < nY+1; j++)
for (int i=1; i < nX+1; i++) for (int i=1; i < nX+1; i++)
vtkFile << i_h[i][j]+b[i][j] << std::endl; vtkFile << i_h[i][j] << std::endl;
vtkFile << "</DataArray>" << std::endl; vtkFile << "</DataArray>" << std::endl;
// Velocities // Momentums
vtkFile << "<DataArray Name=\"U\" type=\"Float32\" format=\"ascii\">" << std::endl; vtkFile << "<DataArray Name=\"hu\" type=\"Float32\" format=\"ascii\">" << std::endl;
for (int j=1; j < nY+1; j++) for (int j=1; j < nY+1; j++)
for (int i=1; i < nX+1; i++) for (int i=1; i < nX+1; i++)
vtkFile << ((i_h[i][j]>0) ? i_hu[i][j]/i_h[i][j] : 0.0 ) << std::endl; vtkFile << i_hu[i][j] << std::endl;
vtkFile << "</DataArray>" << std::endl; vtkFile << "</DataArray>" << std::endl;
vtkFile << "<DataArray Name=\"V\" type=\"Float32\" format=\"ascii\">" << std::endl; vtkFile << "<DataArray Name=\"hv\" type=\"Float32\" format=\"ascii\">" << std::endl;
for (int j=1; j < nY+1; j++) for (int j=1; j < nY+1; j++)
for (int i=1; i<nX+1; i++) for (int i=1; i<nX+1; i++)
vtkFile << ((i_h[i][j]>0) ? i_hv[i][j]/i_h[i][j] : 0.0 ) << std::endl; vtkFile << i_hv[i][j] << std::endl;
vtkFile << "</DataArray>" << std::endl; vtkFile << "</DataArray>" << std::endl;
// Bathymetry // Bathymetry
vtkFile << "<DataArray Name=\"B\" type=\"Float32\" format=\"ascii\">" << std::endl; vtkFile << "<DataArray Name=\"b\" type=\"Float32\" format=\"ascii\">" << std::endl;
for (int j=1; j<nY+1; j++) for (int j=1; j<nY+1; j++)
for (int i=1; i<nX+1; i++) for (int i=1; i<nX+1; i++)
vtkFile << b[i][j] << std::endl; vtkFile << b[i][j] << std::endl;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment