Drop flotation scaling

icepack · Oct 3, 2023 · ce2bb25 · ce2bb25
1 parent 124edf9
commit ce2bb25
Showing 1 changed file with 3 additions and 6 deletions.
diff --git a/demos/helheim/helheim.ipynb b/demos/helheim/helheim.ipynb
@@ -177,7 +177,7 @@
     "T = firedrake.VectorFunctionSpace(mesh, \"CG\", 1)\n",
     "Z = V * Σ * T\n",
     "z = firedrake.Function(Z)\n",
-    "z.sub(0).assign(u)"
+    "z.sub(0).assign(u);"
    ]
   },
   {
@@ -192,10 +192,7 @@
     "\n",
     "ϕ = Constant(0.9)\n",
     "τ_f = firedrake.interpolate(ϕ * sqrt(inner(τ_d, τ_d)), Q)\n",
-    "area = assemble(Constant(1) * dx(mesh))\n",
-    "u_rms = np.sqrt(assemble(inner(u, u) * dx) / area)\n",
-    "print(u_rms)\n",
-    "u_f = Constant(u_rms)"
+    "u_f = firedrake.interpolate(sqrt(inner(u, u)), Q)"
    ]
   },
   {
@@ -284,7 +281,7 @@
     "h_min = Constant(10.0)\n",
     "p_I = ρ_I * g * max_value(h_min, h)\n",
     "p_W = -ρ_W * g * min_value(0, s - h)\n",
-    "f = (1 - p_W / p_I) ** m\n",
+    "f = Constant(1.0)\n",
     "\n",
     "fields = {\n",
     "    \"velocity\": u,\n",