Commit 0510e10e authored by Lukas Riedel's avatar Lukas Riedel

Omit regular flux output for Richards FV solver

parent 98ca6f65
......@@ -288,22 +288,30 @@ void RichardsSimulation<Traits>::write_data () const
if (inifile.get<bool>("output.vertexData")) {
vtkwriter->addVertexData(head,"head");
vtkwriter->addVertexData(wflux,"flux");
vtkwriter->addVertexData(cond,"K_0");
vtkwriter->addVertexData(wc,"theta_w");
vtkwriter->addVertexData(sat,"Theta");
if constexpr (order > 0) {
vtkwriter->addVertexData(wflux, "flux");
}
if constexpr (enable_rt_engine)
if (enable_fluxrc) {
auto wfluxr = get_water_flux_reconstructed();
auto RT_name = "flux_RT" + std::to_string(flux_order);
vtkwriter->addVertexData(wfluxr,RT_name);
}
} else {
}
// cell data
else {
vtkwriter->addCellData(head,"head");
vtkwriter->addCellData(wflux,"flux");
vtkwriter->addCellData(cond,"K_0");
vtkwriter->addCellData(wc,"theta_w");
vtkwriter->addCellData(sat,"Theta");
if constexpr (order > 0) {
vtkwriter->addCellData(wflux, "flux");
}
if constexpr (enable_rt_engine)
if (enable_fluxrc) {
auto wfluxr = get_water_flux_reconstructed();
......
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