From 3a8d1a6d9320c3a77fb154cab89ef78d4e51463a Mon Sep 17 00:00:00 2001
From: Melanie Darcis <mdarcis@gmx.de>
Date: Thu, 11 Nov 2010 14:49:20 +0000
Subject: [PATCH] doc
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@4657 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
tutorial/tutorialproblem_coupled.hh | 23 ++++++++++++++++++++---
1 file changed, 20 insertions(+), 3 deletions(-)
diff --git a/tutorial/tutorialproblem_coupled.hh b/tutorial/tutorialproblem_coupled.hh
index 83fcb90930..5439a48946 100644
--- a/tutorial/tutorialproblem_coupled.hh
+++ b/tutorial/tutorialproblem_coupled.hh
@@ -114,8 +114,6 @@ public:
const GridView &gridView)
: ParentType(timeManager, gridView)
{
- // initialize the tables of the fluid system
- GET_PROP_TYPE(TypeTag, PTAG(FluidSystem))::init();
}
// Specified the problem name. This is used as a prefix for files
@@ -123,6 +121,24 @@ public:
const char *name() const
{ return "tutorial_coupled"; }
+ //! Returns true if a restart file should be written.
+ /* The default behaviour is to write no restart file.
+ */
+ bool shouldWriteRestartFile() const /*@\label{tutorial-coupled:restart}@*/
+ {
+ return false;
+ }
+
+ //! Returns true if the current solution should be written to disk (i.e. as a VTK file)
+ /*! The default behaviour is to write out the solution for
+ * every time step. Else, change divisor.
+ */
+ bool shouldWriteOutput() const /*@\label{tutorial-coupled:output}@*/
+ {
+ return this->timeManager().timeStepIndex() > 0 &&
+ (this->timeManager().timeStepIndex() % 1 == 0);
+ }
+
// Return the temperature within a finite volume. We use constant
// 10 degrees Celsius.
Scalar temperature(const Element &element,
@@ -164,7 +180,8 @@ public:
const GlobalPosition &pos =
fvElemGeom.boundaryFace[boundaryFaceIdx].ipGlobal;
Scalar right = this->bboxMax()[0];
- if (pos[0] > right - eps_) {
+ // extraction of oil on the right boundary for approx. 1.e6 seconds
+ if (pos[0] > right - eps_ && this->timeManager().time() <= 1.e6) {
// oil outflux of 0.3 g/(m * s) on the right boundary.
values[Indices::contiWEqIdx] = 0;
values[Indices::contiNEqIdx] = 3e-4;
--
GitLab