Commit 9a67c93e authored by Kilian Weishaupt's avatar Kilian Weishaupt Committed by Martin Schneider
Browse files

[test][md][stokesdarcy1p] Use new indices

parent 12f29881
......@@ -94,9 +94,9 @@ int main(int argc, char** argv)
auto couplingManager = std::make_shared<CouplingManager>(stokesFvGridGeometry, darcyFvGridGeometry);
// the indices
constexpr auto stokesCellCenterIdx = CouplingManager::stokesCellCenterIdx;
constexpr auto stokesFaceIdx = CouplingManager::stokesFaceIdx;
constexpr auto darcyIdx = CouplingManager::darcyIdx;
constexpr auto stokesCellCenterIdx = CouplingManager::freeFlowCellCenterIdx;
constexpr auto stokesFaceIdx = CouplingManager::freeFlowFaceIdx;
constexpr auto darcyIdx = CouplingManager::porousMediumIdx;
// the problem (initial and boundary conditions)
using StokesProblem = GetPropType<StokesTypeTag, Properties::Problem>;
......
......@@ -105,7 +105,7 @@ public:
BoundaryTypes values;
values.setAllNeumann();
if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
if (couplingManager().isCoupledEntity(CouplingManager::porousMediumIdx, scvf))
values.setAllCouplingNeumann();
if (verticalFlow_)
......@@ -150,8 +150,8 @@ public:
{
NumEqVector values(0.0);
if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
values[Indices::conti0EqIdx] = couplingManager().couplingData().massCouplingCondition(element, fvGeometry, elemVolVars, scvf);
if (couplingManager().isCoupledEntity(CouplingManager::porousMediumIdx, scvf))
values[Indices::conti0EqIdx] = couplingManager().couplingData().massCouplingCondition(element, fvGeometry, elemVolVars, elemFluxVarsCache, scvf);
return values;
}
......
......@@ -156,7 +156,7 @@ public:
}
if(couplingManager().isCoupledEntity(CouplingManager::stokesIdx, scvf))
if(couplingManager().isCoupledEntity(CouplingManager::freeFlowIdx, scvf))
{
values.setCouplingNeumann(Indices::conti0EqIdx);
values.setCouplingNeumann(Indices::momentumYBalanceIdx);
......@@ -192,7 +192,7 @@ public:
{
NumEqVector values(0.0);
if(couplingManager().isCoupledEntity(CouplingManager::stokesIdx, scvf))
if(couplingManager().isCoupledEntity(CouplingManager::freeFlowIdx, scvf))
{
values[Indices::conti0EqIdx] = couplingManager().couplingData().massCouplingCondition(element, fvGeometry, elemVolVars, elemFaceVars, scvf);
values[Indices::momentumYBalanceIdx] = couplingManager().couplingData().momentumCouplingCondition(element, fvGeometry, elemVolVars, elemFaceVars, scvf);
......@@ -250,7 +250,7 @@ public:
*/
Scalar alphaBJ(const SubControlVolumeFace& scvf) const
{
return couplingManager().problem(CouplingManager::darcyIdx).spatialParams().beaversJosephCoeffAtPos(scvf.center());
return couplingManager().problem(CouplingManager::porousMediumIdx).spatialParams().beaversJosephCoeffAtPos(scvf.center());
}
// \}
......
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