Skip to content
Snippets Groups Projects
Commit 07c9d910 authored by Martin Utz's avatar Martin Utz
Browse files

[swe][flux] Make the LET parameter changeable.

parent ac781f2c
No related branches found
No related tags found
1 merge request!1810Resolve "[shallowwater] Make the borders of the LET limiting changeable"
...@@ -83,10 +83,12 @@ std::array<Scalar,3> riemannProblem(const Scalar waterDepthLeft, ...@@ -83,10 +83,12 @@ std::array<Scalar,3> riemannProblem(const Scalar waterDepthLeft,
const Scalar waterDepthRightReconstructed = max(0.0, waterDepthRight - dzr); const Scalar waterDepthRightReconstructed = max(0.0, waterDepthRight - dzr);
// compute the mobility of the flux with the fluxlimiter // compute the mobility of the flux with the fluxlimiter
static const Scalar upperWaterDepthFluxLimiting = getParam<Scalar>("FluxLimiterLET.UpperWaterDepth", 1e-3);
static const Scalar lowerWaterDepthFluxLimiting = getParam<Scalar>("FluxLimiterLET.LowerWaterDepth", 1e-5);
const Scalar mobility = ShallowWater::fluxLimiterLET(waterDepthLeftReconstructed, const Scalar mobility = ShallowWater::fluxLimiterLET(waterDepthLeftReconstructed,
waterDepthRightReconstructed, waterDepthRightReconstructed,
0.001, upperWaterDepthFluxLimiting,
0.00001); lowerWaterDepthFluxLimiting);
// make rotation of the flux we compute an 1d flux // make rotation of the flux we compute an 1d flux
Scalar tempFlux = velocityXLeft; Scalar tempFlux = velocityXLeft;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment