Fix up continuation

Author: danshapero

notebooks/rainier-monolithic.ipynb

@@ -424,30 +424,8 @@
    "id": "4d91d056-05ff-4abb-a702-c1683b6a3718",
    "metadata": {},
    "source": [
-    "First, we'll do an initial solve with a linear rheology."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "b99a1d2c-bdf8-48d7-b450-9155d6276caa",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "n.assign(1.0)\n",
-    "terms = [fn(**fields, **rheology, test_function=ϕ) for fn, ϕ in fns]\n",
-    "H = firedrake.Constant(50.0)\n",
-    "terms[0] = firedrake.replace(terms[0], {h: H})\n",
-    "F = sum(terms)\n",
-    "firedrake.solve(F == 0, z, **pparams, **sparams)"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "c545438e-e0f2-4224-be60-cce046bfa0a6",
-   "metadata": {},
-   "source": [
-    "Now we'll gradually step the flow law exponent from 1.0 up to 3.0."
+    "First, we'll do an initial solve with a linear rheology.\n",
+    "Then we'll gradually step the flow law exponent from 1.0 up to 3.0."
    ]
   },
   {
@@ -457,7 +435,6 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "n.assign(1.1)\n",
     "terms = [fn(**fields, **rheology, test_function=ϕ) for fn, ϕ in fns]\n",
     "terms[0] = firedrake.replace(terms[0], {h: H})\n",
     "F = sum(terms)\n",
@@ -482,7 +459,7 @@
    "outputs": [],
    "source": [
     "num_continuation_steps = 5\n",
-    "for exponent in np.linspace(1.0, 3.0, num_continuation_steps)[1:]:\n",
+    "for exponent in np.linspace(1.0, 3.0, num_continuation_steps):\n",
     "    n.assign(exponent)\n",
     "    momentum_solver.solve()"
    ]