diff --git a/dumux/decoupled/2p/diffusion/fv/fvpressure2padaptive.hh b/dumux/decoupled/2p/diffusion/fv/fvpressure2padaptive.hh
index 36921788bac1075dfb5b435eb691e9b608128cd6..9d2ffcf272b853c3b18fcea60fb313be0d1b84d3 100644
--- a/dumux/decoupled/2p/diffusion/fv/fvpressure2padaptive.hh
+++ b/dumux/decoupled/2p/diffusion/fv/fvpressure2padaptive.hh
@@ -224,15 +224,8 @@ void FVPressure2Padaptive<TypeTag>::getFlux(Dune::FieldVector<Scalar, 2>& entrie
// get face area
Scalar faceArea = intersection.geometry().volume();
- // distance vector between barycenters
- GlobalPosition distVecIJ = globalPosJ - globalPosI;
-
- // compute distance between cell centers
- Scalar distIJ = distVecIJ.two_norm();
-
// Count number of hanging nodes
// not really necessary
- int isIndexJ = intersection.indexInOutside();
int globalIdxK = 0;
ElementPointer elementK = intersection.outside();
// We are looking for two things:
@@ -299,8 +292,6 @@ void FVPressure2Padaptive<TypeTag>::getFlux(Dune::FieldVector<Scalar, 2>& entrie
ng = gravity_ * unitOuterNormal;
- Scalar lambdaWK = cellDataK.mobility(wPhaseIdx);
- Scalar lambdaNWK = cellDataK.mobility(nPhaseIdx);
Scalar fractionalWK = cellDataK.fracFlowFunc(wPhaseIdx);
Scalar fractionalNWK = cellDataK.fracFlowFunc(nPhaseIdx);
diff --git a/dumux/decoupled/2p/diffusion/fv/fvvelocity2padaptive.hh b/dumux/decoupled/2p/diffusion/fv/fvvelocity2padaptive.hh
index 9063b2871ba0f1c110f959aa8f4156a938f74742..872d4943b4ee65a3d75e2de0215fd922ad80d266 100644
--- a/dumux/decoupled/2p/diffusion/fv/fvvelocity2padaptive.hh
+++ b/dumux/decoupled/2p/diffusion/fv/fvvelocity2padaptive.hh
@@ -202,15 +202,6 @@ void FVVelocity2Padaptive<TypeTag>::calculateVelocity(const Intersection& inters
//get face normal
const Dune::FieldVector<Scalar, dim>& unitOuterNormal = intersection.centerUnitOuterNormal();
- // get face area
- Scalar faceArea = intersection.geometry().volume();
-
- // distance vector between barycenters
- GlobalPosition distVecIJ = globalPosJ - globalPosI;
-
- // compute distance between cell centers
- Scalar distIJ = distVecIJ.two_norm();
-
// Count number of hanging nodes
// not really necessary
int isIndexJ = intersection.indexInOutside();
@@ -281,8 +272,6 @@ void FVVelocity2Padaptive<TypeTag>::calculateVelocity(const Intersection& inters
ng = gravity_ * unitOuterNormal;
- Scalar lambdaWK = cellDataK.mobility(wPhaseIdx);
- Scalar lambdaNWK = cellDataK.mobility(nPhaseIdx);
Scalar fractionalWK = cellDataK.fracFlowFunc(wPhaseIdx);
Scalar fractionalNWK = cellDataK.fracFlowFunc(nPhaseIdx);
diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaopressure2p.hh b/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaopressure2p.hh
index 3b9e194d7da575e0686a941cf210c59d1297fd4d..e41e109992b39518cdbfd6053f96a8f5bd45eb3e 100644
--- a/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaopressure2p.hh
+++ b/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaopressure2p.hh
@@ -1102,9 +1102,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1221,10 +1218,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1349,9 +1342,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1483,9 +1473,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1728,9 +1715,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1862,9 +1846,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1918,9 +1899,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -2105,9 +2083,6 @@ void FVMPFAOPressure2P<TypeTag>::assemble()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -2407,8 +2382,6 @@ void FVMPFAOPressure2P<TypeTag>::updateMaterialLaws()
CellData& cellData = problem_.variables().cellData(globalIdx);
- Scalar temperature = problem_.temperature(*eIt);
-
//determine phase saturations from primary saturation variable
Scalar satW = cellData.saturation(wPhaseIdx);
diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaovelocity2p.hh b/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaovelocity2p.hh
index 5e352a76d5412235070b286c6080e635613e3174..3f7088ac8e4aa399472de7d2122de6103687b7c3 100644
--- a/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaovelocity2p.hh
+++ b/dumux/decoupled/2p/diffusion/fvmpfa/fvmpfaovelocity2p.hh
@@ -827,9 +827,6 @@ void FVMPFAOVelocity2P<TypeTag>::calculateVelocity()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -950,9 +947,6 @@ void FVMPFAOVelocity2P<TypeTag>::calculateVelocity()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1089,9 +1083,6 @@ void FVMPFAOVelocity2P<TypeTag>::calculateVelocity()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1238,9 +1229,6 @@ void FVMPFAOVelocity2P<TypeTag>::calculateVelocity()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
@@ -1628,9 +1616,6 @@ void FVMPFAOVelocity2P<TypeTag>::calculateVelocity()
}
}
- Scalar temperature = problem_.temperature(*eIt);
- Scalar referencePressure = problem_.referencePressure(*eIt);
-
Scalar lambdaWBound = 0;
Scalar lambdaNWBound = 0;
diff --git a/dumux/decoupled/2p/diffusion/mimetic/mimeticgroundwater.hh b/dumux/decoupled/2p/diffusion/mimetic/mimeticgroundwater.hh
index 1546c17d3d07795ed6dd6796c647350f7a5d5275..8da2600ea89da5dbce92de6a06fe66eee5f38181 100644
--- a/dumux/decoupled/2p/diffusion/mimetic/mimeticgroundwater.hh
+++ b/dumux/decoupled/2p/diffusion/mimetic/mimeticgroundwater.hh
@@ -419,8 +419,6 @@ private:
if (this->bctype[faceIndex].isNeumann(pressEqIdx))
{
- int globalIdx = problem_.variables().index(element);
-
problem_.neumann(boundValues, *isIt);
Scalar J = (boundValues[wPhaseIdx]/density_[wPhaseIdx] + boundValues[nPhaseIdx]/density_[nPhaseIdx]);
this->b[faceIndex] -= J * isIt->geometry().volume();
diff --git a/dumux/decoupled/2p/diffusion/mimetic/mimeticpressure2p.hh b/dumux/decoupled/2p/diffusion/mimetic/mimeticpressure2p.hh
index 310abcaf7bc9506b855846313294431689d30d7a..c91e5a4d8a499f19b0502c8cfc4c264d97fd262c 100644
--- a/dumux/decoupled/2p/diffusion/mimetic/mimeticpressure2p.hh
+++ b/dumux/decoupled/2p/diffusion/mimetic/mimeticpressure2p.hh
@@ -315,8 +315,6 @@ void MimeticPressure2P<TypeTag>::updateMaterialLaws()
CellData& cellData = problem_.variables().cellData(globalIdx);
- Scalar temperature = problem_.temperature(*eIt);
-
//determine phase saturations from primary saturation variable
Scalar satW = cellData.saturation(wPhaseIdx);