[cig-commits] [commit] master: Add test for gplates-3d. (9bf63b5)

[cig_noreply at geodynamics.org]

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