Commit b34d04c5 authored by Dennis Gläser's avatar Dennis Gläser
Browse files

[mpfa][darcyslaw] initialize rho even for gravity disabled

parent fb08e150
......@@ -84,9 +84,7 @@ public:
const auto& tij = fluxVarsCache.advectionTij(phaseIdx);
// interface density as arithmetic mean of the neighbors (when gravity is on)
Scalar rho;
if (gravity)
rho = interpolateDensity(elemVolVars, scvf, phaseIdx);
Scalar rho = gravity ? interpolateDensity(elemVolVars, scvf, phaseIdx) : 0.0;
// calculate Tij*pj
Scalar flux(0.0);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment