diff --git a/dumux/decoupled/2p2c/fvpressure2p2cmultiphysics.hh b/dumux/decoupled/2p2c/fvpressure2p2cmultiphysics.hh index 5067f72619dd11cd358127975420466028bb05e4..d7dc58268e1f3246490fe50584ce175685816742 100644 --- a/dumux/decoupled/2p2c/fvpressure2p2cmultiphysics.hh +++ b/dumux/decoupled/2p2c/fvpressure2p2cmultiphysics.hh @@ -528,8 +528,6 @@ void FVPressure2P2CMultiPhysics<TypeTag>::get1pFlux(Dune::FieldVector<Scalar, 2> // due to "safety cell" around subdomain, both cells I and J // have single-phase conditions, although one is in 2p domain. int phaseIdx = std::min(cellDataI.subdomain(), cellDataJ.subdomain()); - // determine respective equation Idx - int contiEqIdx = (phaseIdx == wPhaseIdx) ? contiWEqIdx : contiNEqIdx; Scalar rhoMean = 0.5 * (cellDataI.density(phaseIdx) + cellDataJ.density(phaseIdx)); //Scalar density = 0;