[cig-commits] [commit] master: Add test for gplates-3d. (9bf63b5)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Oct 29 11:29:23 PDT 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/aad39993908bfd5cc020de4548464f6027d37173...aed4d3ee69f4a126db2602d3d0bb6fcc4a49633a
>---------------------------------------------------------------
commit 9bf63b550ce123e7772f51b3e4338fb363fb5160
Author: Rene Gassmoeller <R.Gassmoeller at mailbox.org>
Date: Wed Oct 29 11:27:33 2014 +0100
Add test for gplates-3d.
>---------------------------------------------------------------
9bf63b550ce123e7772f51b3e4338fb363fb5160
tests/{gplates-1-3.prm => gplates-3d.prm} | 11 ++++-----
tests/gplates-3d/screen-output | 35 ++++++++++++++++++++++++++++
tests/{gplates-1-3 => gplates-3d}/statistics | 6 ++++-
3 files changed, 45 insertions(+), 7 deletions(-)
diff --git a/tests/gplates-1-3.prm b/tests/gplates-3d.prm
similarity index 86%
copy from tests/gplates-1-3.prm
copy to tests/gplates-3d.prm
index cb86c9c..dda8ede 100644
--- a/tests/gplates-1-3.prm
+++ b/tests/gplates-3d.prm
@@ -1,7 +1,6 @@
-# A simple setup for for using the GPlates interface in a 2d shell.
-# This test uses the grid produced by GPlates 1.3
+# A simple setup for for using the GPlates interface in a 3d shell.
-set Dimension = 2
+set Dimension = 3
set Use years in output instead of seconds = true
set End time = 0
set Output directory = output
@@ -42,7 +41,7 @@ end
subsection Boundary velocity model
subsection GPlates model
set Data directory = $ASPECT_SOURCE_DIR/data/velocity-boundary-conditions/gplates/
- set Velocity file name = current_day.gpml
+ set Velocity file name = current_day_1.4.gpml
set Time step = 1e6
set Point one = 1.5708,4.87
set Point two = 1.5708,5.24
@@ -80,12 +79,12 @@ end
subsection Mesh refinement
set Initial adaptive refinement = 0
- set Initial global refinement = 3
+ set Initial global refinement = 1
set Strategy = temperature
set Time steps between mesh refinement = 2
end
subsection Postprocess
- set List of postprocessors = velocity statistics, temperature statistics, heat flux statistics
+ set List of postprocessors = velocity statistics, temperature statistics, heat flux statistics, velocity boundary statistics
end
diff --git a/tests/gplates-3d/screen-output b/tests/gplates-3d/screen-output
new file mode 100644
index 0000000..b188e0b
--- /dev/null
+++ b/tests/gplates-3d/screen-output
@@ -0,0 +1,35 @@
+
+Number of active cells: 768 (on 2 levels)
+Number of degrees of freedom: 28,690 (20,790+970+6,930)
+
+
+ Setting up GPlates boundary velocity plugin.
+
+
+ Loading GPlates boundary velocity file ASPECT_DIR/data/velocity-boundary-conditions/gplates/current_day_1.4.gpml.
+
+
+ Loading GPlates boundary velocity file ASPECT_DIR/data/velocity-boundary-conditions/gplates/current_day_1.4.gpml.
+
+
+ Loading new velocity file did not succeed.
+ Assuming constant boundary conditions for rest of model run.
+
+*** Timestep 0: t=0 years
+ Solving temperature system... 0 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 30+34 iterations.
+
+ Postprocessing:
+ RMS, max velocity: 0.0466 m/year, 0.152 m/year
+ Temperature min/avg/max: 273 K, 1456 K, 2600 K
+ Heat fluxes through boundary parts: -8.33e+12 W, 7.445e+12 W
+ Max and min velocity along boundary parts: 0.1011 m/yr, 0.001509 m/yr, 0.1176 m/yr, 0.0007523 m/yr
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
++---------------------------------+-----------+------------+------------+
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/gplates-1-3/statistics b/tests/gplates-3d/statistics
similarity index 55%
copy from tests/gplates-1-3/statistics
copy to tests/gplates-3d/statistics
index 3d59c62..efd01f8 100644
--- a/tests/gplates-1-3/statistics
+++ b/tests/gplates-3d/statistics
@@ -16,4 +16,8 @@
# 16: Average nondimensional temperature (K)
# 17: Outward heat flux through boundary with indicator 0 ("inner") (W)
# 18: Outward heat flux through boundary with indicator 1 ("outer") (W)
-0 0.0000e+00 768 7392 3264 0 45 108 64 1.6661e+06 2.56804983e-02 1.01393827e-01 2.73000000e+02 1.57484593e+03 2.60000000e+03 5.59452482e-01 -9.36071065e+05 1.91866167e+06
+# 19: Maximum velocity magnitude on boundary with indicator 0 ("inner") (m/yr)
+# 20: Minimum velocity magnitude on boundary with indicator 0 ("inner") (m/yr)
+# 21: Maximum velocity magnitude on boundary with indicator 1 ("outer") (m/yr)
+# 22: Minimum velocity magnitude on boundary with indicator 1 ("outer") (m/yr)
+0 0.0000e+00 768 21760 6930 0 64 196 210 8.5232e+05 4.66045988e-02 1.51878113e-01 2.73000000e+02 1.45591131e+03 2.60000000e+03 5.08341777e-01 -8.32986869e+12 7.44494632e+12 1.01104202e-01 1.50888268e-03 1.17586379e-01 7.52285535e-04
More information about the CIG-COMMITS
mailing list