Commit 59d7b337 authored by Lukas Riedel's avatar Lukas Riedel

Store BC value for outflow both as head and flux

This makes the outflow BC usable as Dirichlet BC replacement in the
Richards model
parent 44737794
......@@ -646,7 +646,7 @@ namespace Dorie{
else if(buffer == "outflow"){
bcData[i][j].type.push_back(BoundaryCondition::Outflow);
instream >> buffer;
bcData[i][j].headValue.push_back(0.0);
bcData[i][j].headValue.push_back(stod(buffer));
bcData[i][j].fluxValue.push_back(stod(buffer));
}
else
......
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