From 2a166dec6dcf558f75e29b79b27ca42a2bd86481 Mon Sep 17 00:00:00 2001
From: Benjamin Faigle <benjamin.faigle@posteo.de>
Date: Fri, 1 Feb 2013 13:49:53 +0000
Subject: [PATCH] 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
---
dumux/decoupled/2p2c/fvpressurecompositional.hh | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/dumux/decoupled/2p2c/fvpressurecompositional.hh b/dumux/decoupled/2p2c/fvpressurecompositional.hh
index 2de2d7ca5e..6be4a3592b 100644
--- a/dumux/decoupled/2p2c/fvpressurecompositional.hh
+++ b/dumux/decoupled/2p2c/fvpressurecompositional.hh
@@ -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);
--
GitLab