better threshold

NDF-Poli-USP · Apr 14, 2021 · d8913e5 · d8913e5
1 parent e45081a
commit d8913e5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/spyro/solvers/leapfrog_adjoint_level_set.py b/spyro/solvers/leapfrog_adjoint_level_set.py
@@ -186,7 +186,7 @@ def Leapfrog_adjoint_level_set(
         # analytical background gradient
         # Z, _ = SpatialCoordinate(mesh)
         # c_background = Function(V).interpolate(Min(1.5 + 4.0 * abs(Z), 4.1))
-        indicator = Function(V).interpolate(c > 3.9)  # 1 inside shape, 0 outside
+        indicator = Function(V).interpolate(c > 4.3)  # 1 inside shape, 0 outside
 
         gradc = Function(VF, name="grad_c").interpolate(
             indicator * grad(c_salt) + (1 - indicator) * grad(c_background)