Commit 69f99683 authored by Lukas Riedel's avatar Lukas Riedel

removed diagnostic output

parent f723a504
......@@ -433,11 +433,6 @@ public:
else if (BoundaryCondition::isEvaporation(bc))
{
// switches between Neumann and Dirichlet bc
// this choice is kept constant for one time step
RF u_0 = boundary.g(ig.intersection(),it.position(),time);
normal_flux = boundary.j(ig.intersection(),it.position(),time);
// estimate the dirichlet flux
// evaluate gradient of basis functions (we assume Galerkin method lfsu = lfsv)
const auto& gradphi_s = cache[order_s].evaluateJacobian(local_s,lfsu_s.finiteElement().localBasis());
......@@ -473,10 +468,12 @@ public:
// compute numerical flux estimation
const RF numFlux = satCond_s * ( - (gradu_s * normal) + penalty * jump);
// switches between Neumann and Dirichlet bc
// this choice is kept constant for one time step
normal_flux = boundary.j(ig.intersection(),it.position(),time);
if(numFlux < normal_flux) {
bcType = BCType::dirichlet;
bc_type_cache.emplace(it_id,BoundaryCondition::Dirichlet);
std::cout << "- Falling back to Dirichlet BC." << std::endl;
}
else {
bcType = BCType::vonNeumann;
......
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