[Transport] Fix Dirichlet mode switch in local operator

parent fb080ea0
...@@ -376,10 +376,10 @@ public: ...@@ -376,10 +376,10 @@ public:
// Convert input to total solute if needed // Convert input to total solute if needed
if (_dirichlet_mode == DirichletMode::TotalSolute) if (_dirichlet_mode == DirichletMode::TotalSolute)
{ {
if (Dune::FloatCmp::gt(water_content_i,0.)) if (Dune::FloatCmp::gt(water_content_i,WaterContent{0.}))
g /=water_content_i; g /=water_content_i;
else else
g = 0. g = 0.;
} }
// Compute the effective hydrodynamic dispersion // Compute the effective hydrodynamic dispersion
......
Markdown is supported
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