Commit 2a166dec authored by Benjamin Faigle's avatar Benjamin Faigle
Browse files

Tiny tweaks in initialization output:

- uses output routine of derived pressure modules through asImp_()
- Total concentrations are calculated in every initialization step to avoid "nan" in output


git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@10145 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent b6e92074
......@@ -302,7 +302,7 @@ public:
{
std::cout << "Writing debug for current time step\n";
initializationOutputWriter_.beginWrite(problem_.timeManager().time() + pseudoTS);
addOutputVtkFields(initializationOutputWriter_);
asImp_().addOutputVtkFields(initializationOutputWriter_);
#if DUNE_MINIMAL_DEBUG_LEVEL <= 2
int size_ = problem_.gridView().size(0);
......@@ -614,9 +614,10 @@ void FVPressureCompositional<TypeTag>::initialMaterialLaws(bool compositional)
pressure, problem_.spatialParams().porosity(*eIt), temperature_);
}
} //end conc initial condition
cellData.calculateMassConcentration(problem_.spatialParams().porosity(*eIt));
} //end compositional
cellData.calculateMassConcentration(problem_.spatialParams().porosity(*eIt));
problem_.transportModel().totalConcentration(wCompIdx,globalIdx) = cellData.massConcentration(wCompIdx);
problem_.transportModel().totalConcentration(nCompIdx,globalIdx) = cellData.massConcentration(nCompIdx);
......
Supports Markdown
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