From 0598cb9712c52da0213d23d64413f7d57023f545 Mon Sep 17 00:00:00 2001
From: danieldouglas92 <daniel.douglas@student.nmt.edu>
Date: Wed, 6 Nov 2024 16:44:13 -0700
Subject: [PATCH] add test for original.wb

---
 tests/original_wb.prm         |  6 ++++++
 tests/original_wb/original.wb | 31 +++++++++++++++++++++++++++++++
 tests/original_wb/statistics  | 14 ++++++++++++++
 3 files changed, 51 insertions(+)
 create mode 100644 tests/original_wb.prm
 create mode 100644 tests/original_wb/original.wb
 create mode 100644 tests/original_wb/statistics

diff --git a/tests/original_wb.prm b/tests/original_wb.prm
new file mode 100644
index 00000000000..5df3b63ac55
--- /dev/null
+++ b/tests/original_wb.prm
@@ -0,0 +1,6 @@
+# test that output/original.wb is written correctly
+
+# based on world_builder_simple.prm:
+# set World builder file = $ASPECT_SOURCE_DIR/tests/world_builder_simple.wb
+include $ASPECT_SOURCE_DIR/tests/world_builder_simple.prm
+set World builder file                     = $ASPECT_SOURCE_DIR/tests/world_builder_simple.wb
\ No newline at end of file
diff --git a/tests/original_wb/original.wb b/tests/original_wb/original.wb
new file mode 100644
index 00000000000..92a5536332a
--- /dev/null
+++ b/tests/original_wb/original.wb
@@ -0,0 +1,31 @@
+{
+  "version":"1.0",
+  "cross section":[[0,0],[100,0]],
+  "features":
+  [
+    {"model":"oceanic plate", "name":"oceanic plate", "coordinates":[[-1e3,-1e3],[1150e3,-1e3],[1150e3,1e3],[-1e3,1e3]],
+     "temperature models":[{"model":"plate model", "max depth":95e3, "bottom temperature":1600, "spreading velocity":0.005, "ridge coordinates":[[[100e3,-1e3],[100e3,1e3]]]}],
+     "composition models":[{"model":"uniform", "compositions":[0], "max depth":10e3},
+                           {"model":"uniform", "compositions":[1], "min depth":10e3, "max depth":95e3}]},
+
+    {"model":"continental plate", "name":"continental plate", "coordinates":[[1150e3,-1e3],[2001e3,-1e3],[2001e3,1e3],[1150e3,1e3]],
+     "temperature models":[{"model":"linear", "max depth":95e3, "bottom temperature":1600}],
+     "composition models":[{"model":"uniform", "compositions":[2], "max depth":30e3},
+                           {"model":"uniform", "compositions":[3], "min depth":30e3, "max depth":65e3}]},
+
+    {"model":"mantle layer", "name":"upper mantle", "min depth":95e3, "max depth":660e3, "coordinates":[[-1e3,-1e3],[2001e3,-1e3],[2001e3,1e3],[-1e3,1e3]],
+     "temperature models":[{"model":"linear", "min depth":95e3, "max depth":660e3, "top temperature":1600, "bottom temperature":1820}],
+     "composition models":[{"model":"uniform", "compositions":[4]}]},
+
+    {"model":"mantle layer", "name":"lower mantle", "min depth":660e3, "max depth":1160e3, "coordinates":[[-1e3,-1e3],[2001e3,-1e3],[2001e3,1e3],[-1e3,1e3]],
+     "temperature models":[{"model":"linear", "min depth":660e3, "max depth":1160e3, "top temperature":1820, "bottom temperature":2000}],
+     "composition models":[{"model":"uniform", "compositions":[5]}]},
+
+    {"model":"subducting plate", "name":"Subducting plate", "coordinates":[[1150e3,-1e3],[1150e3,1e3]], "dip point":[2000e3,0],
+     "segments":[{"length":200e3, "thickness":[95e3], "angle":[0,45]}, {"length":200e3, "thickness":[95e3], "angle":[45]}, 
+                 {"length":200e3, "thickness":[95e3], "angle":[45,0]}, {"length":100e3, "thickness":[95e3], "angle":[0]}],
+     "temperature models":[{"model":"plate model", "density":3300, "plate velocity":0.01, "adiabatic heating":false}],
+     "composition models":[{"model":"uniform", "compositions":[0], "max distance slab top":10e3},
+                           {"model":"uniform", "compositions":[1], "min distance slab top":10e3, "max distance slab top":95e3 }]}
+  ]
+}
diff --git a/tests/original_wb/statistics b/tests/original_wb/statistics
new file mode 100644
index 00000000000..3d03d497d3c
--- /dev/null
+++ b/tests/original_wb/statistics
@@ -0,0 +1,14 @@
+# 1: Time step number
+# 2: Time (seconds)
+# 3: Time step size (seconds)
+# 4: Number of mesh cells
+# 5: Number of Stokes degrees of freedom
+# 6: Number of temperature degrees of freedom
+# 7: Number of degrees of freedom for all compositions
+# 8: Iterations for temperature solver
+# 9: Iterations for composition solver 1
+# 10: Iterations for Stokes solver
+# 11: Velocity iterations in Stokes preconditioner
+# 12: Schur complement iterations in Stokes preconditioner
+# 13: Visualization file name
+0 0.000000000000e+00 0.000000000000e+00 1024 9539 4225 4225 0 0 36 38 38 output-original_wb/solution/solution-00000