diff --git a/dumux/boxmodels/1p2c/1p2cfluxvariables.hh b/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
index 1490d4a2bac10ff8aa02cb1471a71bd13c5eef5f..6a56063ce064381144e6750fe97894a9311ae178 100644
--- a/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
+++ b/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
@@ -90,6 +90,9 @@ public:
     {
         scvfIdx_ = scvfIdx;
 
+        viscosityAtIP_ = Scalar(0);
+        molarDensityAtIP_ = Scalar(0);
+
         calculateGradients_(problem, element, elemDat);
         calculateK_(problem, element, elemDat);
         calculateDiffCoeffPM_(problem, element, elemDat);
@@ -139,6 +142,12 @@ public:
         return diffCoeffPM_;
     };
 
+    Scalar viscosityAtIP() const
+    { return viscosityAtIP_;}
+
+    Scalar molarDensityAtIP() const
+    { return molarDensityAtIP_; }
+
 
     /*!
      * \brief Given the intrinisc permeability times the pressure
@@ -202,6 +211,12 @@ protected:
                 tmp = feGrad;
                 tmp *= elemDat[idx].concentration(1);
                 concentrationGrad_ += tmp;
+
+                // phase viscosity
+                viscosityAtIP_ += elemDat[idx].viscosity()*face().shapeValue[idx];
+
+                //phase moledensity
+                molarDensityAtIP_ += elemDat[idx].molarDensity()*face().shapeValue[idx];
             }
         }
         else {
@@ -334,6 +349,12 @@ protected:
 
     //! the intrinsic permeability tensor
     Tensor K_;
+
+    //! viscosity of the fluid at the integration point
+    Scalar viscosityAtIP_;
+
+    //! molar densities of the fluid at the integration point
+    Scalar molarDensityAtIP_;
 };
 
 } // end namepace