Work around DirichletBC state lock

jrmaddison · Nov 10, 2023 · 84470c2 · 84470c2
1 parent 0f9d3db
commit 84470c2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/icepack/solvers/flow_solver.py b/src/icepack/solvers/flow_solver.py
@@ -247,7 +247,7 @@ def setup(self, **kwargs):
         # Create homogeneous BCs for the Dirichlet part of the boundary
         u = self._fields["velocity"]
         V = u.function_space()
-        bcs = firedrake.DirichletBC(V, u, self._dirichlet_ids)
+        bcs = firedrake.DirichletBC(V, u.copy(deepcopy=True), self._dirichlet_ids)
         if not self._dirichlet_ids:
             bcs = None