Add test case

test/tests/userobjects/element_subdomain_modifier/adaptivity.i

@@ -0,0 +1,99 @@
+[Problem]
+  solve = false
+[]
+
+[Mesh]
+  [gen]
+    type = GeneratedMeshGenerator
+    dim = 2
+    xmin = -1
+    xmax = 1
+    ymin = -1
+    ymax = 1
+    nx = 16
+    ny = 16
+  []
+  [left]
+    type = SubdomainBoundingBoxGenerator
+    input = 'gen'
+    block_id = 1
+    bottom_left = '-1 -1 0'
+    top_right = '0 1 1'
+  []
+  [right]
+    type = SubdomainBoundingBoxGenerator
+    input = 'left'
+    block_id = 2
+    bottom_left = '0 -1 0'
+    top_right = '1 1 1'
+  []
+  [moving_boundary]
+    type = SideSetsAroundSubdomainGenerator
+    input = 'right'
+    block = 1
+    new_boundary = 'moving_boundary'
+    normal = '1 0 0'
+  []
+[]
+
+[UserObjects]
+  [moving_circle]
+    type = CoupledVarThresholdElementSubdomainModifier
+    coupled_var = 'phi'
+    block = 2
+    criterion_type = ABOVE
+    threshold = 0.5
+    subdomain_id = 1
+    moving_boundary_name = moving_boundary
+    execute_on = 'TIMESTEP_BEGIN'
+  []
+[]
+
+[Functions]
+  [moving_gauss]
+    type = ParsedFunction
+    value = 'exp(-((x+0.5-t)^2+(y)^2)/0.25)'
+  []
+[]
+
+[AuxVariables]
+  [phi]
+  []
+[]
+
+[AuxKernels]
+  [phi]
+    type = FunctionAux
+    variable = phi
+    function = moving_gauss
+    execute_on = 'INITIAL TIMESTEP_BEGIN TIMESTEP_END'
+  []
+[]
+
+[Adaptivity]
+  steps = 1
+  marker = marker
+  initial_marker = marker
+  max_h_level = 1
+  [Indicators/indicator]
+    type = GradientJumpIndicator
+    variable = phi
+  []
+  [Markers/marker]
+    type = ErrorFractionMarker
+    indicator = indicator
+    coarsen = 0.2
+    refine = 0.5
+    check_subdomain_consistent_for_coarsen = true
+  []
+[]
+
+[Executioner]
+  type = Transient
+  dt = 0.1
+  num_steps = 10
+[]
+
+[Outputs]
+  exodus = true
+[]