Made diffusion work again: The density at the integration points has not been...

Made diffusion work again: The density at the integration points has not been calculated in the flux variables and was thus 0 --> no diffusive fluxes.

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@7912 2fb0f335-1f38-0410-981e-8018bf24f1b0

dumux/boxmodels/2p2c/2p2cfluxvariables.hh

@@ -112,6 +112,21 @@ protected:
                           const Element &element,
                           const ElementVolumeVariables &elemVolVars)
     {
+        // calculate densities at the integration points of the face
+        Vector tmp(0.0);
+        for (int idx = 0;
+             idx < fvGeom_.numVertices;
+             idx++) // loop over adjacent vertices
+        {
+            for (int phaseIdx = 0; phaseIdx < numPhases; phaseIdx++)
+            {
+               densityAtIP_[phaseIdx] += elemVolVars[idx].density(phaseIdx)*
+                       face().shapeValue[idx];
+               molarDensityAtIP_[phaseIdx] += elemVolVars[idx].molarDensity(phaseIdx)*
+                       face().shapeValue[idx];
+            }
+        }
+
         calculateGradients_(problem, element, elemVolVars);
         calculateVelocities_(problem, element, elemVolVars);
         calculateDiffCoeffPM_(problem, element, elemVolVars);