From 541e49e3a4836c2cfc88c73a46be80f042cc8b1e Mon Sep 17 00:00:00 2001
From: Johannes Hommel <johannes.hommel@iws.uni-stuttgart.de>
Date: Mon, 20 Nov 2017 11:50:31 +0100
Subject: [PATCH] [2pncmin][next] ported the .cc and problem files according to
their 2pnc versions
---
.../2pncmin/implicit/test_cc2pncmin.cc | 12 +++++-------
1 file changed, 5 insertions(+), 7 deletions(-)
diff --git a/test/porousmediumflow/2pncmin/implicit/test_cc2pncmin.cc b/test/porousmediumflow/2pncmin/implicit/test_cc2pncmin.cc
index fdbafc6593..0a7f57a387 100644
--- a/test/porousmediumflow/2pncmin/implicit/test_cc2pncmin.cc
+++ b/test/porousmediumflow/2pncmin/implicit/test_cc2pncmin.cc
@@ -113,9 +113,8 @@ int main(int argc, char** argv) try
auto problem = std::make_shared<Problem>(fvGridGeometry);
// the solution vector
- using GridView = typename GET_PROP_TYPE(TypeTag, GridView);
using SolutionVector = typename GET_PROP_TYPE(TypeTag, SolutionVector);
- SolutionVector x(leafGridView.size(GridView::dimension));
+ SolutionVector x(leafGridView.size(0));
problem->applyInitialSolution(x);
auto xOld = x;
@@ -141,6 +140,9 @@ int main(int argc, char** argv) try
VtkOutputModule<TypeTag> vtkWriter(*problem, *fvGridGeometry, *gridVariables, x, problem->name());
VtkOutputFields::init(vtkWriter); //! Add model specific output fields
//add specific output
+ vtkWriter.addField(problem->getCurrentDensity(), "currentDensity [A/cm^2]");
+ vtkWriter.addField(problem->getReactionSourceH2O(), "reactionSourceH2O [mol/(sm^2)]");
+ vtkWriter.addField(problem->getReactionSourceO2(), "reactionSourceO2 [mol/(sm^2)]");
vtkWriter.addField(problem->getKxx(), "Kxx");
vtkWriter.addField(problem->getKyy(), "Kyy");
vtkWriter.write(0.0);
@@ -166,10 +168,6 @@ int main(int argc, char** argv) try
// time loop
timeLoop->start(); do
{
- // set time for problem for implicit Euler scheme
- problem->setTime( timeLoop->time() + timeLoop->timeStepSize() );
- problem->setTimeStepSize( timeLoop->timeStepSize() );
-
// set previous solution for storage evaluations
assembler->setPreviousSolution(xOld);
@@ -251,4 +249,4 @@ catch (...)
{
std::cerr << "Unknown exception thrown! ---> Abort!" << std::endl;
return 4;
-}
\ No newline at end of file
+}
--
GitLab