diff --git a/dumux/porousmediumflow/2pncmin/implicit/model.hh b/dumux/porousmediumflow/2pncmin/implicit/model.hh index 59884fc075e526bd44d9c6c5be3b26b1763e1bf6..3f20df868d2ed564d5aebfa1c26e913ef0ba7775 100644 --- a/dumux/porousmediumflow/2pncmin/implicit/model.hh +++ b/dumux/porousmediumflow/2pncmin/implicit/model.hh @@ -444,7 +444,7 @@ protected: else if (phasePresence == nPhaseOnly) { Scalar sumxw = 0; - //Calculate sum of mole fractions (water and air) in the hypothetical wetting phase + //Calculate sum of mole fractions in the hypothetical wetting phase for (int compIdx = 0; compIdx < numComponents; compIdx++) { sumxw += volVars.moleFraction(wPhaseIdx, compIdx); diff --git a/dumux/porousmediumflow/2pncmin/implicit/volumevariables.hh b/dumux/porousmediumflow/2pncmin/implicit/volumevariables.hh index d1e8f39617b97cb89467b91337dcf73a10e716ef..b83bfd9003d6c98320094513674e97855a41776c 100644 --- a/dumux/porousmediumflow/2pncmin/implicit/volumevariables.hh +++ b/dumux/porousmediumflow/2pncmin/implicit/volumevariables.hh @@ -352,13 +352,13 @@ public: moleFrac[compIdx] = priVars[compIdx]; } moleFrac[nCompIdx] = priVars[switchIdx]; - Scalar sumMoleFracNotWater = 0; + Scalar sumMoleFracOtherComponents = 0; for (int compIdx=numMajorComponents; compIdx<numComponents; ++compIdx) { - sumMoleFracNotWater+=moleFrac[compIdx]; + sumMoleFracOtherComponents+=moleFrac[compIdx]; } - sumMoleFracNotWater += moleFrac[nCompIdx]; - moleFrac[wCompIdx] = 1 -sumMoleFracNotWater; + sumMoleFracOtherComponents += moleFrac[nCompIdx]; + moleFrac[wCompIdx] = 1 -sumMoleFracOtherComponents; // convert mass to mole fractions and set the fluid state for (int compIdx=0; compIdx<numComponents; ++compIdx)