[cig-commits] [commit] master: Add a test for the existing GPlates plugin. (775ae08)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Sat May 24 14:21:37 PDT 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/e2d2428a203331b4832dbc4e5f4d716934c22694...f68c1b584dc0211691228ef830ea1d2538da258d
>---------------------------------------------------------------
commit 775ae083ef88136b64c5cd45aa9d4082a331a919
Author: Rene Gassmoeller <R.Gassmoeller at mailbox.org>
Date: Fri May 23 21:16:36 2014 -0500
Add a test for the existing GPlates plugin.
>---------------------------------------------------------------
775ae083ef88136b64c5cd45aa9d4082a331a919
cookbooks/gplates-2d.prm => tests/gplates-1-3.prm | 22 +++------
tests/gplates-1-3/screen-output | 59 +++++++++++++++++++++++
tests/gplates-1-3/statistics | 17 +++++++
3 files changed, 82 insertions(+), 16 deletions(-)
diff --git a/cookbooks/gplates-2d.prm b/tests/gplates-1-3.prm
similarity index 76%
copy from cookbooks/gplates-2d.prm
copy to tests/gplates-1-3.prm
index dfb1b1e..cb86c9c 100644
--- a/cookbooks/gplates-2d.prm
+++ b/tests/gplates-1-3.prm
@@ -1,9 +1,9 @@
-# A simple setup for for using the GPlates interface in a 2d shell. See the
-# manual for more information.
+# A simple setup for for using the GPlates interface in a 2d shell.
+# This test uses the grid produced by GPlates 1.3
set Dimension = 2
set Use years in output instead of seconds = true
-set End time = 3e8
+set End time = 0
set Output directory = output
set Adiabatic surface temperature = 1600
@@ -79,23 +79,13 @@ end
subsection Mesh refinement
- set Initial adaptive refinement = 1
- set Initial global refinement = 5
+ set Initial adaptive refinement = 0
+ set Initial global refinement = 3
set Strategy = temperature
set Time steps between mesh refinement = 2
end
subsection Postprocess
- set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics, depth average
-
- subsection Visualization
- set Output format = vtu
- set Time between graphical output = 1e6
- set Number of grouped files = 0
- end
-
- subsection Depth average
- set Time between graphical output = 1e6
- end
+ set List of postprocessors = velocity statistics, temperature statistics, heat flux statistics
end
diff --git a/tests/gplates-1-3/screen-output b/tests/gplates-1-3/screen-output
new file mode 100644
index 0000000..438969c
--- /dev/null
+++ b/tests/gplates-1-3/screen-output
@@ -0,0 +1,59 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
+-- . running in DEBUG mode
+-- . running with 1 MPI process
+-- . using Trilinos
+-----------------------------------------------------------------------------
+
+Number of active cells: 768 (on 4 levels)
+Number of degrees of freedom: 10,656 (6,528+864+3,264)
+
+
+ Set up GPlates boundary velocity module.
+
+ Input point 1 spherical coordinates: 1.571 4.870
+ Input point 1 Cartesian output coordinates: 0.157 -0.988 0.000
+ Input point 2 spherical coordinates: 1.571 5.240
+ Input point 2 Cartesian output coordinates: 0.503 -0.864 0.000
+
+ Model will be rotated by 0.00 degrees around axis 0.00 0.00 1.00
+
+
+ Loading GPlates boundary velocity file /home/rengas/Software/Aspect-Versionen/aspect/data/velocity-boundary-conditions/gplates/current_day.gpml.
+
+
+ Loading GPlates boundary velocity file /home/rengas/Software/Aspect-Versionen/aspect/data/velocity-boundary-conditions/gplates/current_day.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+12 iterations.
+
+ Postprocessing:
+ RMS, max velocity: 0.0257 m/year, 0.101 m/year
+ Temperature min/avg/max: 273 K, 1575 K, 2600 K
+ Heat fluxes through boundary parts: -9.361e+05 W, 1.919e+06 W
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 4.74s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 1 | 1.02s | 21% |
+| Assemble temperature system | 1 | 0.888s | 19% |
+| Build Stokes preconditioner | 1 | 0.713s | 15% |
+| Build temperature preconditioner| 1 | 0.00956s | 0.2% |
+| Solve Stokes system | 1 | 0.409s | 8.6% |
+| Solve temperature system | 1 | 0.00272s | 0.057% |
+| Initialization | 2 | 0.649s | 14% |
+| Postprocessing | 1 | 0.118s | 2.5% |
+| Setup dof systems | 1 | 0.456s | 9.6% |
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/gplates-1-3/statistics b/tests/gplates-1-3/statistics
new file mode 100644
index 0000000..0c1d1db
--- /dev/null
+++ b/tests/gplates-1-3/statistics
@@ -0,0 +1,17 @@
+# 1: Time step number
+# 2: Time (years)
+# 3: Number of mesh cells
+# 4: Number of Stokes degrees of freedom
+# 5: Number of temperature degrees of freedom
+# 6: Iterations for temperature solver
+# 7: Iterations for Stokes solver
+# 8: Time step size (years)
+# 9: RMS velocity (m/year)
+# 10: Max. velocity (m/year)
+# 11: Minimal temperature (K)
+# 12: Average temperature (K)
+# 13: Maximal temperature (K)
+# 14: Average nondimensional temperature (K)
+# 15: Outward heat flux through boundary with indicator 0 (W)
+# 16: Outward heat flux through boundary with indicator 1 (W)
+0 0.0000e+00 768 7392 3264 0 42 1.6661e+06 2.56804950e-02 1.01393828e-01 2.73000000e+02 1.57484593e+03 2.60000000e+03 5.59452482e-01 -9.36071065e+05 1.91866167e+06
More information about the CIG-COMMITS
mailing list