[problem] Pass more paramters to neumann method that it is possibly soldependent

dumux/implicit/box/localresidual.hh

@@ -188,7 +188,7 @@ protected:
         // temporary vector to store the neumann boundary fluxes
         PrimaryVariables flux(0);
 
-        auto neumannFluxes = this->problem().neumann(element, scvf);
+        auto neumannFluxes = this->problem().neumann(element, fvGeometry, elemVolVars, scvf);
 
         // multiply neumann fluxes with the area and the extrusion factor
         neumannFluxes *= scvf.area()*elemVolVars[insideScv].extrusionFactor();

dumux/implicit/cellcentered/localresidual.hh

@@ -176,7 +176,7 @@ protected:
         // temporary vector to store the neumann boundary fluxes
         PrimaryVariables flux(0);
 
-        auto neumannFluxes = this->problem().neumann(element, scvf);
+        auto neumannFluxes = this->problem().neumann(element, fvGeometry, elemVolVars, scvf);
 
         // multiply neumann fluxes with the area and the extrusion factor
         auto&& scv = fvGeometry.scv(scvf.insideScvIdx());

dumux/implicit/problem.hh

@@ -273,7 +273,10 @@ public:
      * in normal direction of each phase. Negative values mean influx.
      * E.g. for the mass balance that would the mass flux in \f$ [ kg / (m^2 \cdot s)] \f$.
      */
-    PrimaryVariables neumann(const Element &element, const SubControlVolumeFace &scvFace) const
+    PrimaryVariables neumann(const Element& element,
+                             const FVElementGeometry& fvGeometry,
+                             const ElementVolumeVariables& elemVolvars,
+                             const SubControlVolumeFace& scvFace) const
     {
         // forward it to the interface with only the global position
         return asImp_().neumannAtPos(scvFace.center());