Add the boundary term, seems to go

Author: danshapero

notebooks/frontal-ablation.ipynb

@@ -13,7 +13,7 @@
     "from IPython.display import HTML\n",
     "from tqdm.notebook import trange, tqdm\n",
     "import firedrake\n",
-    "from firedrake import Constant, max_value, inner, grad, jump, dx, dS\n",
+    "from firedrake import Constant, max_value, inner, grad, jump, dx, ds, dS\n",
     "import irksome\n",
     "from irksome import Dt\n",
     "from icepack2 import model"
@@ -119,7 +119,11 @@
     "\n",
     "ν = firedrake.FacetNormal(mesh)\n",
     "F_c = max_value(0, inner(f_c, ν))\n",
-    "F += -inner(f_c, grad(ϕ)) * dx + jump(F_c) * jump(ϕ) * dS"
+    "F += (\n",
+    "    -inner(f_c, grad(ϕ)) * dx\n",
+    "    + jump(F_c) * jump(ϕ) * dS\n",
+    "    + F_c * ϕ * ds\n",
+    ")"
    ]
   },
   {