Fix to the flores indices

This commit is contained in:
David Landa Marban
2023-02-07 00:41:21 +01:00
parent bdec8e1a6b
commit 94b8f5a941

View File

@@ -377,7 +377,8 @@ public:
darcyFlux = pressureDifference *
(Toolbox::value(up.mobility(phaseIdx, facedir)) * Toolbox::value(transMult) * (-trans / faceArea));
}
darcy[phaseIdx] = darcyFlux.value() * faceArea; // For the FLORES fluxes
unsigned activeCompIdx = Indices::canonicalToActiveComponentIndex(FluidSystem::solventComponentIndex(phaseIdx));
darcy[conti0EqIdx + activeCompIdx] = darcyFlux.value() * faceArea; // For the FLORES fluxes
unsigned pvtRegionIdx = up.pvtRegionIndex();
// if (upIdx == globalFocusDofIdx){