Fix BCs + default output filename

icepack · Jan 19, 2024 · ff10c2c · ff10c2c
1 parent f93d160
commit ff10c2c
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/demos/kangerd/initialize.py b/demos/kangerd/initialize.py
@@ -16,8 +16,8 @@
 parser.add_argument("--outline")
 parser.add_argument("--num-levels", type=int, default=1)
 parser.add_argument("--degree", type=int, default=1)
-parser.add_argument("--regularization", type=float, default=5e3)
-parser.add_argument("--output")
+parser.add_argument("--regularization", type=float, default=2.5e3)
+parser.add_argument("--output", default="kangerdlugssuaq-initial.h5")
 args = parser.parse_args()
 
 # Fetch the glacier outline, generate a mesh, and create function spaces
@@ -170,7 +170,7 @@ def bed_friction(**kwargs):
 
 flow_model = icepack.models.IceStream(friction=bed_friction)
 opts = {
-    "dirichlet_ids": [1, 2, 4],
+    "dirichlet_ids": [1, 2, 3, 4],
     "diagnostic_solver_type": "petsc",
     "diagnostic_solver_parameters": {
         "snes_type": "newtontr",