[cig-commits] r1362 - in trunk/aspect/tests: . diffusion
bangerth at dealii.org
bangerth at dealii.org
Wed Nov 14 06:24:18 PST 2012
Author: bangerth
Date: 2012-11-14 07:24:18 -0700 (Wed, 14 Nov 2012)
New Revision: 1362
Added:
trunk/aspect/tests/diffusion.prm
trunk/aspect/tests/diffusion/
trunk/aspect/tests/diffusion/screen-output
trunk/aspect/tests/diffusion/statistics
Log:
New test by Juliane Dannberg checking a half-space diffusion model.
Added: trunk/aspect/tests/diffusion/screen-output
===================================================================
--- trunk/aspect/tests/diffusion/screen-output (rev 0)
+++ trunk/aspect/tests/diffusion/screen-output 2012-11-14 14:24:18 UTC (rev 1362)
@@ -0,0 +1,651 @@
+Number of active cells: 1,024 (on 6 levels)
+Number of degrees of freedom: 13,764 (8,450+1,089+4,225)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 2 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: -7.101e-17 K, 0.005208 K, 1 K
+ Heat fluxes through boundary parts: -6.648e-19 W, -1.105e-18 W, 2.605e-46 W, -9.6e-05 W
+
+*** Timestep 1: t=2.31396e+06 seconds
+ Solving temperature system... 47 iterations.
+ Solving Stokes system... 9 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.04311 K, 1 K
+ Heat fluxes through boundary parts: 3.527e-12 W, -4.018e-12 W, 3.578e-15 W, -2.17e-05 W
+
+*** Timestep 2: t=4.05617e+06 seconds
+ Solving temperature system... 33 iterations.
+ Solving Stokes system... 5 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.06177 K, 1 K
+ Heat fluxes through boundary parts: 1.209e-11 W, 1.608e-11 W, 9.273e-15 W, -1.04e-05 W
+
+*** Timestep 3: t=6.52448e+06 seconds
+ Solving temperature system... 36 iterations.
+ Solving Stokes system... 5 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.08069 K, 1 K
+ Heat fluxes through boundary parts: -4.783e-11 W, -1.152e-11 W, 3.687e-14 W, -7.295e-06 W
+
+*** Timestep 4: t=8.9178e+06 seconds
+ Solving temperature system... 35 iterations.
+ Solving Stokes system... 5 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.09499 K, 1 K
+ Heat fluxes through boundary parts: -2.584e-11 W, -1.09e-11 W, 1.399e-13 W, -6.424e-06 W
+
+*** Timestep 5: t=1.12201e+07 seconds
+ Solving temperature system... 33 iterations.
+ Solving Stokes system... 5 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1067 K, 1 K
+ Heat fluxes through boundary parts: -3.782e-12 W, -2.833e-12 W, 4.631e-13 W, -5.854e-06 W
+
+Number of active cells: 472 (on 6 levels)
+Number of degrees of freedom: 6,575 (4,032+527+2,016)
+
+*** Timestep 6: t=1.48115e+07 seconds
+ Solving temperature system... 37 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 4 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1228 K, 1 K
+ Heat fluxes through boundary parts: 4.044e-12 W, -1.253e-12 W, 2.781e-12 W, -5.138e-06 W
+
+*** Timestep 7: t=1.58214e+07 seconds
+ Solving temperature system... 22 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1269 K, 1 K
+ Heat fluxes through boundary parts: -8.753e-13 W, 2.516e-11 W, 4.105e-12 W, -4.979e-06 W
+
+*** Timestep 8: t=1.6642e+07 seconds
+ Solving temperature system... 17 iterations.
+ Solving Stokes system... 4 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1302 K, 1 K
+ Heat fluxes through boundary parts: -4.629e-12 W, 4.116e-11 W, 5.533e-12 W, -4.862e-06 W
+
+*** Timestep 9: t=1.77375e+07 seconds
+ Solving temperature system... 19 iterations.
+ Solving Stokes system... 3 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1344 K, 1 K
+ Heat fluxes through boundary parts: -9.886e-12 W, 4.666e-11 W, 8.183e-12 W, -4.719e-06 W
+
+*** Timestep 10: t=1.89012e+07 seconds
+ Solving temperature system... 20 iterations.
+ Solving Stokes system... 5 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1387 K, 1 K
+ Heat fluxes through boundary parts: 8.277e-12 W, 1.772e-11 W, 1.232e-11 W, -4.583e-06 W
+
+Number of active cells: 349 (on 6 levels)
+Number of degrees of freedom: 4,965 (3,044+399+1,522)
+
+*** Timestep 11: t=2.06973e+07 seconds
+ Solving temperature system... 24 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 3 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1451 K, 1 K
+ Heat fluxes through boundary parts: 3.92e-12 W, -6.661e-11 W, 2.345e-12 W, -4.388e-06 W
+
+*** Timestep 12: t=2.11623e+07 seconds
+ Solving temperature system... 13 iterations.
+ Solving Stokes system... 4 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1468 K, 1 K
+ Heat fluxes through boundary parts: 5.186e-12 W, -9.828e-11 W, 3.451e-12 W, -4.339e-06 W
+
+*** Timestep 13: t=2.21044e+07 seconds
+ Solving temperature system... 17 iterations.
+ Solving Stokes system... 4 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.15 K, 1 K
+ Heat fluxes through boundary parts: -3.958e-11 W, 8.964e-11 W, 6.86e-12 W, -4.243e-06 W
+
+*** Timestep 14: t=2.27163e+07 seconds
+ Solving temperature system... 15 iterations.
+ Solving Stokes system... 3 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.152 K, 1 K
+ Heat fluxes through boundary parts: -2.328e-11 W, 7.578e-11 W, 9.77e-12 W, -4.191e-06 W
+
+*** Timestep 15: t=2.31166e+07 seconds
+ Solving temperature system... 11 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1533 K, 1 K
+ Heat fluxes through boundary parts: -9.689e-12 W, 1.087e-10 W, 1.197e-11 W, -4.159e-06 W
+
+Number of active cells: 331 (on 6 levels)
+Number of degrees of freedom: 4,729 (2,900+379+1,450)
+
+*** Timestep 16: t=2.35849e+07 seconds
+ Solving temperature system... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1548 K, 1 K
+ Heat fluxes through boundary parts: -1.107e-11 W, 3.135e-11 W, -5.589e-10 W, -4.123e-06 W
+
+*** Timestep 17: t=2.41236e+07 seconds
+ Solving temperature system... 13 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1565 K, 1 K
+ Heat fluxes through boundary parts: 2.215e-11 W, -5.733e-11 W, -7.561e-10 W, -4.083e-06 W
+
+*** Timestep 18: t=2.44991e+07 seconds
+ Solving temperature system... 11 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1578 K, 1 K
+ Heat fluxes through boundary parts: 5.16e-11 W, -8.072e-11 W, -8.875e-10 W, -4.049e-06 W
+
+*** Timestep 19: t=2.47079e+07 seconds
+ Solving temperature system... 8 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1586 K, 1 K
+ Heat fluxes through boundary parts: 6.404e-11 W, -1.005e-10 W, -9.563e-10 W, -4.027e-06 W
+
+*** Timestep 20: t=2.48998e+07 seconds
+ Solving temperature system... 7 iterations.
+ Solving Stokes system... 2 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1593 K, 1 K
+ Heat fluxes through boundary parts: 7.482e-11 W, -9.522e-11 W, -1.017e-09 W, -4.006e-06 W
+
+Number of active cells: 331 (on 6 levels)
+Number of degrees of freedom: 4,739 (2,906+380+1,453)
+
+*** Timestep 21: t=2.52335e+07 seconds
+ Solving temperature system... 10 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 2 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1603 K, 1 K
+ Heat fluxes through boundary parts: -2.889e-11 W, 8.296e-11 W, -1.11e-09 W, -3.981e-06 W
+
+*** Timestep 22: t=2.54305e+07 seconds
+ Solving temperature system... 8 iterations.
+ Solving Stokes system... 3 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1608 K, 1 K
+ Heat fluxes through boundary parts: -3.7e-11 W, 7.275e-11 W, -1.159e-09 W, -3.97e-06 W
+
+*** Timestep 23: t=2.61959e+07 seconds
+ Solving temperature system... 15 iterations.
+ Solving Stokes system... 4 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1636 K, 1 K
+ Heat fluxes through boundary parts: 1.237e-10 W, -1.94e-10 W, -1.328e-09 W, -3.893e-06 W
+
+*** Timestep 24: t=2.65784e+07 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1649 K, 1 K
+ Heat fluxes through boundary parts: 9.246e-11 W, -1.672e-10 W, -1.399e-09 W, -3.862e-06 W
+
+*** Timestep 25: t=2.68215e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1656 K, 1 K
+ Heat fluxes through boundary parts: 7.001e-11 W, -1.392e-10 W, -1.439e-09 W, -3.843e-06 W
+
+Number of active cells: 373 (on 6 levels)
+Number of degrees of freedom: 5,276 (3,236+422+1,618)
+
+*** Timestep 26: t=2.70176e+07 seconds
+ Solving temperature system... 8 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1662 K, 1 K
+ Heat fluxes through boundary parts: 9.094e-11 W, -1.045e-10 W, -1.158e-09 W, -3.828e-06 W
+
+*** Timestep 27: t=2.72129e+07 seconds
+ Solving temperature system... 7 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1668 K, 1 K
+ Heat fluxes through boundary parts: 5.069e-11 W, -4.534e-11 W, -8.414e-10 W, -3.812e-06 W
+
+*** Timestep 28: t=2.74692e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1677 K, 1 K
+ Heat fluxes through boundary parts: 3.45e-11 W, 1.071e-11 W, -4.936e-10 W, -3.792e-06 W
+
+*** Timestep 29: t=2.77761e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1687 K, 1 K
+ Heat fluxes through boundary parts: -2.689e-11 W, -6.563e-12 W, -2.065e-10 W, -3.769e-06 W
+
+*** Timestep 30: t=2.81109e+07 seconds
+ Solving temperature system... 10 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1698 K, 1 K
+ Heat fluxes through boundary parts: -3.105e-12 W, -3.251e-11 W, -3.599e-11 W, -3.746e-06 W
+
+Number of active cells: 388 (on 6 levels)
+Number of degrees of freedom: 5,468 (3,354+437+1,677)
+
+*** Timestep 31: t=2.8466e+07 seconds
+ Solving temperature system... 10 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 2 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1708 K, 1 K
+ Heat fluxes through boundary parts: -4.218e-12 W, -6.936e-11 W, 1.75e-11 W, -3.725e-06 W
+
+*** Timestep 32: t=2.89023e+07 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.172 K, 1 K
+ Heat fluxes through boundary parts: 1.397e-11 W, 1.006e-11 W, -4.152e-11 W, -3.7e-06 W
+
+*** Timestep 33: t=2.94151e+07 seconds
+ Solving temperature system... 13 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1735 K, 1 K
+ Heat fluxes through boundary parts: -4.137e-11 W, -1.147e-11 W, -2.224e-10 W, -3.669e-06 W
+
+*** Timestep 34: t=2.96703e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1743 K, 1 K
+ Heat fluxes through boundary parts: -1.226e-11 W, 6.512e-12 W, -3.446e-10 W, -3.652e-06 W
+
+*** Timestep 35: t=2.99188e+07 seconds
+ Solving temperature system... 8 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1751 K, 1 K
+ Heat fluxes through boundary parts: 7.331e-13 W, 9.708e-13 W, -4.758e-10 W, -3.635e-06 W
+
+Number of active cells: 388 (on 6 levels)
+Number of degrees of freedom: 5,468 (3,354+437+1,677)
+
+*** Timestep 36: t=3.02023e+07 seconds
+ Solving temperature system... 9 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1759 K, 1 K
+ Heat fluxes through boundary parts: 1.123e-11 W, 1.969e-11 W, -6.336e-10 W, -3.617e-06 W
+
+*** Timestep 37: t=3.06634e+07 seconds
+ Solving temperature system... 11 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1773 K, 1 K
+ Heat fluxes through boundary parts: 8.687e-12 W, -4.713e-12 W, -8.966e-10 W, -3.589e-06 W
+
+*** Timestep 38: t=3.11404e+07 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1787 K, 1 K
+ Heat fluxes through boundary parts: -3.517e-12 W, -4.53e-11 W, -1.167e-09 W, -3.563e-06 W
+
+*** Timestep 39: t=3.16027e+07 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1799 K, 1 K
+ Heat fluxes through boundary parts: -1.194e-11 W, -4.006e-11 W, -1.422e-09 W, -3.539e-06 W
+
+*** Timestep 40: t=3.18594e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1806 K, 1 K
+ Heat fluxes through boundary parts: -9.595e-12 W, -4.615e-11 W, -1.559e-09 W, -3.525e-06 W
+
+Number of active cells: 388 (on 6 levels)
+Number of degrees of freedom: 5,468 (3,354+437+1,677)
+
+*** Timestep 41: t=3.20634e+07 seconds
+ Solving temperature system... 7 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1812 K, 1 K
+ Heat fluxes through boundary parts: -7.164e-12 W, -4.927e-11 W, -1.665e-09 W, -3.514e-06 W
+
+*** Timestep 42: t=3.22387e+07 seconds
+ Solving temperature system... 6 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1816 K, 1 K
+ Heat fluxes through boundary parts: -5.048e-12 W, -4.978e-11 W, -1.754e-09 W, -3.505e-06 W
+
+*** Timestep 43: t=3.24369e+07 seconds
+ Solving temperature system... 7 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1822 K, 1 K
+ Heat fluxes through boundary parts: -7.134e-12 W, -3.524e-11 W, -1.853e-09 W, -3.494e-06 W
+
+*** Timestep 44: t=3.27187e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.183 K, 1 K
+ Heat fluxes through boundary parts: -5.764e-14 W, -2.675e-11 W, -1.989e-09 W, -3.479e-06 W
+
+*** Timestep 45: t=3.31596e+07 seconds
+ Solving temperature system... 11 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1843 K, 1 K
+ Heat fluxes through boundary parts: 1.476e-11 W, -8.241e-12 W, -2.194e-09 W, -3.453e-06 W
+
+Number of active cells: 436 (on 6 levels)
+Number of degrees of freedom: 6,099 (3,742+486+1,871)
+
+*** Timestep 46: t=3.34635e+07 seconds
+ Solving temperature system... 10 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1852 K, 1 K
+ Heat fluxes through boundary parts: 5.864e-12 W, 1.244e-11 W, -2.328e-09 W, -3.437e-06 W
+
+*** Timestep 47: t=3.37081e+07 seconds
+ Solving temperature system... 8 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1859 K, 1 K
+ Heat fluxes through boundary parts: 1.248e-12 W, 8.698e-12 W, -2.433e-09 W, -3.425e-06 W
+
+*** Timestep 48: t=3.38708e+07 seconds
+ Solving temperature system... 7 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1863 K, 1 K
+ Heat fluxes through boundary parts: -1.511e-12 W, 3.092e-12 W, -2.5e-09 W, -3.417e-06 W
+
+*** Timestep 49: t=3.40038e+07 seconds
+ Solving temperature system... 6 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1867 K, 1 K
+ Heat fluxes through boundary parts: -3.987e-12 W, -2.014e-12 W, -2.554e-09 W, -3.411e-06 W
+
+*** Snapshot created!
+
+*** Timestep 50: t=3.41488e+07 seconds
+ Solving temperature system... 6 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1871 K, 1 K
+ Heat fluxes through boundary parts: -5.881e-12 W, 9.169e-12 W, -2.611e-09 W, -3.404e-06 W
+
+Number of active cells: 427 (on 6 levels)
+Number of degrees of freedom: 5,985 (3,672+477+1,836)
+
+*** Timestep 51: t=3.43202e+07 seconds
+ Solving temperature system... 6 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1875 K, 1 K
+ Heat fluxes through boundary parts: -1.28e-12 W, 1.334e-11 W, -2.678e-09 W, -3.396e-06 W
+
+*** Timestep 52: t=3.45347e+07 seconds
+ Solving temperature system... 7 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1881 K, 1 K
+ Heat fluxes through boundary parts: 7.812e-12 W, 7.985e-12 W, -2.759e-09 W, -3.386e-06 W
+
+*** Timestep 53: t=3.48202e+07 seconds
+ Solving temperature system... 8 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1888 K, 1 K
+ Heat fluxes through boundary parts: 2.25e-11 W, -2.839e-12 W, -2.863e-09 W, -3.372e-06 W
+
+*** Timestep 54: t=3.51163e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1896 K, 1 K
+ Heat fluxes through boundary parts: 4.027e-11 W, -1.664e-11 W, -2.966e-09 W, -3.357e-06 W
+
+*** Timestep 55: t=3.54302e+07 seconds
+ Solving temperature system... 10 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1905 K, 1 K
+ Heat fluxes through boundary parts: 2.721e-11 W, -2.833e-11 W, -3.072e-09 W, -3.342e-06 W
+
+Number of active cells: 418 (on 6 levels)
+Number of degrees of freedom: 5,881 (3,608+469+1,804)
+
+*** Timestep 56: t=3.56478e+07 seconds
+ Solving temperature system... 8 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1911 K, 1 K
+ Heat fluxes through boundary parts: 2.148e-11 W, -2.648e-11 W, -3.142e-09 W, -3.331e-06 W
+
+*** Timestep 57: t=3.5802e+07 seconds
+ Solving temperature system... 6 iterations.
+ Solving Stokes system... 2 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1915 K, 1 K
+ Heat fluxes through boundary parts: 1.229e-11 W, -2.653e-11 W, -3.191e-09 W, -3.323e-06 W
+
+*** Timestep 58: t=3.6275e+07 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 3 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1927 K, 1 K
+ Heat fluxes through boundary parts: 1.569e-11 W, 2.512e-11 W, -3.333e-09 W, -3.304e-06 W
+
+*** Timestep 59: t=3.64394e+07 seconds
+ Solving temperature system... 8 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1931 K, 1 K
+ Heat fluxes through boundary parts: -3.872e-12 W, 2.758e-11 W, -3.379e-09 W, -3.298e-06 W
+
+*** Timestep 60: t=3.65954e+07 seconds
+ Solving temperature system... 6 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1935 K, 1 K
+ Heat fluxes through boundary parts: 8.103e-12 W, -1.668e-13 W, -3.422e-09 W, -3.291e-06 W
+
+Number of active cells: 418 (on 6 levels)
+Number of degrees of freedom: 5,878 (3,606+469+1,803)
+
+*** Timestep 61: t=3.67968e+07 seconds
+ Solving temperature system... 7 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.194 K, 1 K
+ Heat fluxes through boundary parts: 9.254e-12 W, 7.277e-12 W, -3.476e-09 W, -3.283e-06 W
+
+*** Timestep 62: t=3.72663e+07 seconds
+ Solving temperature system... 11 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1953 K, 1 K
+ Heat fluxes through boundary parts: 2.777e-11 W, -5.52e-12 W, -3.594e-09 W, -3.263e-06 W
+
+*** Timestep 63: t=3.75835e+07 seconds
+ Solving temperature system... 10 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1961 K, 1 K
+ Heat fluxes through boundary parts: 3.007e-11 W, -4.195e-12 W, -3.669e-09 W, -3.249e-06 W
+
+*** Timestep 64: t=3.77341e+07 seconds
+ Solving temperature system... 6 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1965 K, 1 K
+ Heat fluxes through boundary parts: 3.168e-11 W, -3.057e-12 W, -3.703e-09 W, -3.242e-06 W
+
+*** Timestep 65: t=3.79365e+07 seconds
+ Solving temperature system... 7 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.197 K, 1 K
+ Heat fluxes through boundary parts: 1.395e-11 W, -2.357e-12 W, -3.748e-09 W, -3.234e-06 W
+
+Number of active cells: 418 (on 6 levels)
+Number of degrees of freedom: 5,868 (3,600+468+1,800)
+
+*** Timestep 66: t=3.83739e+07 seconds
+ Solving temperature system... 11 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1981 K, 1 K
+ Heat fluxes through boundary parts: -5.606e-11 W, 1.854e-11 W, -3.836e-09 W, -3.217e-06 W
+
+*** Timestep 67: t=3.86669e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1988 K, 1 K
+ Heat fluxes through boundary parts: -6.481e-11 W, 1.163e-11 W, -3.888e-09 W, -3.205e-06 W
+
+*** Timestep 68: t=3.89735e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.1996 K, 1 K
+ Heat fluxes through boundary parts: -4.698e-11 W, 9.371e-12 W, -3.938e-09 W, -3.192e-06 W
+
+*** Timestep 69: t=3.92538e+07 seconds
+ Solving temperature system... 9 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.2003 K, 1 K
+ Heat fluxes through boundary parts: -1.409e-11 W, -5.32e-12 W, -3.981e-09 W, -3.179e-06 W
+
+*** Timestep 70: t=3.95933e+07 seconds
+ Solving temperature system... 10 iterations.
+ Solving Stokes system... 1 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 0 K, 0.2013 K, 1 K
+ Heat fluxes through boundary parts: 2.95e-12 W, -3.217e-12 W, -4.03e-09 W, -3.164e-06 W
+
+Number of active cells: 418 (on 6 levels)
+Number of degrees of freedom: 5,881 (3,608+469+1,804)
+
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 62s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 71 | 4.52s | 7.3% |
+| Assemble temperature system | 71 | 33.8s | 55% |
+| Build Stokes preconditioner | 15 | 2.49s | 4% |
+| Build temperature preconditioner| 71 | 0.573s | 0.92% |
+| Solve Stokes system | 71 | 0.96s | 1.5% |
+| Solve temperature system | 71 | 0.553s | 0.89% |
+| Create snapshot | 1 | 0.0344s | 0.056% |
+| Initialization | 2 | 0.0886s | 0.14% |
+| Postprocessing | 71 | 1.78s | 2.9% |
+| Refine mesh structure, part 1 | 14 | 9.8s | 16% |
+| Refine mesh structure, part 2 | 14 | 0.0916s | 0.15% |
+| Setup dof systems | 15 | 5.56s | 9% |
++---------------------------------+-----------+------------+------------+
+
Added: trunk/aspect/tests/diffusion/statistics
===================================================================
--- trunk/aspect/tests/diffusion/statistics (rev 0)
+++ trunk/aspect/tests/diffusion/statistics 2012-11-14 14:24:18 UTC (rev 1362)
@@ -0,0 +1,87 @@
+# 1: Time step number
+# 2: Time (seconds)
+# 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 (seconds)
+# 9: Minimal temperature (K)
+# 10: Average temperature (K)
+# 11: Maximal temperature (K)
+# 12: Average nondimensional temperature (K)
+# 13: Outward heat flux through boundary with indicator 0 (W)
+# 14: Outward heat flux through boundary with indicator 1 (W)
+# 15: Outward heat flux through boundary with indicator 2 (W)
+# 16: Outward heat flux through boundary with indicator 3 (W)
+ 0 0.0000e+00 1024 9539 4225 12 2 2.3140e+06 -7.10110744e-17 5.20833333e-03 1.00000000e+00 5.20833333e-03 -6.64806357e-19 -1.10492141e-18 2.60509697e-46 -9.60000000e-05
+ 1 2.3140e+06 1024 9539 4225 47 9 1.7422e+06 0.00000000e+00 4.31149878e-02 1.00000000e+00 4.31149878e-02 3.52693478e-12 -4.01793624e-12 3.57775604e-15 -2.17006871e-05
+ 2 4.0562e+06 1024 9539 4225 33 5 2.4683e+06 0.00000000e+00 6.17703343e-02 1.00000000e+00 6.17703343e-02 1.20880539e-11 1.60754348e-11 9.27345900e-15 -1.04012461e-05
+ 3 6.5245e+06 1024 9539 4225 36 5 2.3933e+06 0.00000000e+00 8.06935516e-02 1.00000000e+00 8.06935516e-02 -4.78342778e-11 -1.15246765e-11 3.68680191e-14 -7.29514129e-06
+ 4 8.9178e+06 1024 9539 4225 35 5 2.3023e+06 0.00000000e+00 9.49853955e-02 1.00000000e+00 9.49853955e-02 -2.58368113e-11 -1.08964241e-11 1.39912560e-13 -6.42448739e-06
+ 5 1.1220e+07 1024 9539 4225 33 5 3.5914e+06 0.00000000e+00 1.06730226e-01 1.00000000e+00 1.06730226e-01 -3.78209099e-12 -2.83293645e-12 4.63127411e-13 -5.85439473e-06
+ 6 1.4812e+07 472 4559 2016 37 4 1.0098e+06 0.00000000e+00 1.22794054e-01 1.00000000e+00 1.22794054e-01 4.04389415e-12 -1.25270959e-12 2.78091579e-12 -5.13780519e-06
+ 7 1.5821e+07 472 4559 2016 22 0 8.2066e+05 0.00000000e+00 1.26922845e-01 1.00000000e+00 1.26922845e-01 -8.75276754e-13 2.51607031e-11 4.10463087e-12 -4.97897145e-06
+ 8 1.6642e+07 472 4559 2016 17 4 1.0955e+06 0.00000000e+00 1.30181531e-01 1.00000000e+00 1.30181531e-01 -4.62855593e-12 4.11596778e-11 5.53301923e-12 -4.86210739e-06
+ 9 1.7738e+07 472 4559 2016 19 3 1.1636e+06 0.00000000e+00 1.34397641e-01 1.00000000e+00 1.34397641e-01 -9.88629196e-12 4.66591258e-11 8.18336802e-12 -4.71904006e-06
+10 1.8901e+07 472 4559 2016 20 5 1.7962e+06 0.00000000e+00 1.38656437e-01 1.00000000e+00 1.38656437e-01 8.27681635e-12 1.77218708e-11 1.23199534e-11 -4.58271527e-06
+11 2.0697e+07 349 3443 1522 24 3 4.6491e+05 0.00000000e+00 1.45123032e-01 1.00000000e+00 1.45123032e-01 3.91978292e-12 -6.66122535e-11 2.34536230e-12 -4.38833175e-06
+12 2.1162e+07 349 3443 1522 13 4 9.4215e+05 0.00000000e+00 1.46797324e-01 1.00000000e+00 1.46797324e-01 5.18644337e-12 -9.82807707e-11 3.45108709e-12 -4.33940032e-06
+13 2.2104e+07 349 3443 1522 17 4 6.1190e+05 0.00000000e+00 1.50037528e-01 1.00000000e+00 1.50037528e-01 -3.95767106e-11 8.96378126e-11 6.86001630e-12 -4.24287093e-06
+14 2.2716e+07 349 3443 1522 15 3 4.0031e+05 0.00000000e+00 1.52031427e-01 1.00000000e+00 1.52031427e-01 -2.32808743e-11 7.57825726e-11 9.77020817e-12 -4.19089097e-06
+15 2.3117e+07 349 3443 1522 11 1 4.6826e+05 0.00000000e+00 1.53294367e-01 1.00000000e+00 1.53294367e-01 -9.68880409e-12 1.08706686e-10 1.19712962e-11 -4.15893727e-06
+16 2.3585e+07 331 3279 1450 12 0 5.3869e+05 0.00000000e+00 1.54814689e-01 1.00000000e+00 1.54814689e-01 -1.10747173e-11 3.13491743e-11 -5.58893101e-10 -4.12310194e-06
+17 2.4124e+07 331 3279 1450 13 1 3.7555e+05 0.00000000e+00 1.56545646e-01 1.00000000e+00 1.56545646e-01 2.21531167e-11 -5.73334533e-11 -7.56127139e-10 -4.08342871e-06
+18 2.4499e+07 331 3279 1450 11 0 2.0878e+05 0.00000000e+00 1.57835075e-01 1.00000000e+00 1.57835075e-01 5.15976822e-11 -8.07236280e-11 -8.87544533e-10 -4.04871688e-06
+19 2.4708e+07 331 3279 1450 8 1 1.9192e+05 0.00000000e+00 1.58576352e-01 1.00000000e+00 1.58576352e-01 6.40419425e-11 -1.00480679e-10 -9.56272708e-10 -4.02744676e-06
+20 2.4900e+07 331 3279 1450 7 2 3.3363e+05 0.00000000e+00 1.59272020e-01 1.00000000e+00 1.59272020e-01 7.48239328e-11 -9.52240428e-11 -1.01669641e-09 -4.00641344e-06
+21 2.5233e+07 331 3286 1453 10 2 1.9708e+05 0.00000000e+00 1.60306262e-01 1.00000000e+00 1.60306262e-01 -2.88948329e-11 8.29578364e-11 -1.11030639e-09 -3.98054917e-06
+22 2.5431e+07 331 3286 1453 8 3 7.6541e+05 0.00000000e+00 1.60849181e-01 1.00000000e+00 1.60849181e-01 -3.69982333e-11 7.27455227e-11 -1.15915103e-09 -3.97029360e-06
+23 2.6196e+07 331 3286 1453 15 4 3.8249e+05 0.00000000e+00 1.63622307e-01 1.00000000e+00 1.63622307e-01 1.23660664e-10 -1.94027998e-10 -1.32804181e-09 -3.89325919e-06
+24 2.6578e+07 331 3286 1453 12 0 2.4306e+05 0.00000000e+00 1.64851723e-01 1.00000000e+00 1.64851723e-01 9.24612145e-11 -1.67162629e-10 -1.39876333e-09 -3.86172038e-06
+25 2.6821e+07 331 3286 1453 9 0 1.9607e+05 0.00000000e+00 1.65606486e-01 1.00000000e+00 1.65606486e-01 7.00087879e-11 -1.39177410e-10 -1.43856325e-09 -3.84269365e-06
+26 2.7018e+07 373 3658 1618 8 1 1.9538e+05 0.00000000e+00 1.66208423e-01 1.00000000e+00 1.66208423e-01 9.09401800e-11 -1.04458479e-10 -1.15838040e-09 -3.82762364e-06
+27 2.7213e+07 373 3658 1618 7 1 2.5625e+05 0.00000000e+00 1.66826731e-01 1.00000000e+00 1.66826731e-01 5.06854992e-11 -4.53351083e-11 -8.41361746e-10 -3.81240368e-06
+28 2.7469e+07 373 3658 1618 9 1 3.0689e+05 0.00000000e+00 1.67679753e-01 1.00000000e+00 1.67679753e-01 3.45020590e-11 1.07086746e-11 -4.93585759e-10 -3.79182960e-06
+29 2.7776e+07 373 3658 1618 9 1 3.3479e+05 0.00000000e+00 1.68701430e-01 1.00000000e+00 1.68701430e-01 -2.68892574e-11 -6.56260019e-12 -2.06462223e-10 -3.76856141e-06
+30 2.8111e+07 373 3658 1618 10 1 3.5516e+05 0.00000000e+00 1.69769562e-01 1.00000000e+00 1.69769562e-01 -3.10492113e-12 -3.25146828e-11 -3.59931102e-11 -3.74594370e-06
+31 2.8466e+07 388 3791 1677 10 2 4.3629e+05 0.00000000e+00 1.70821277e-01 1.00000000e+00 1.70821277e-01 -4.21824076e-12 -6.93559308e-11 1.74980838e-11 -3.72507827e-06
+32 2.8902e+07 388 3791 1677 12 1 5.1274e+05 0.00000000e+00 1.71998084e-01 1.00000000e+00 1.71998084e-01 1.39673228e-11 1.00581208e-11 -4.15197399e-11 -3.70036745e-06
+33 2.9415e+07 388 3791 1677 13 1 2.5527e+05 0.00000000e+00 1.73501086e-01 1.00000000e+00 1.73501086e-01 -4.13743187e-11 -1.14681533e-11 -2.22360443e-10 -3.66879929e-06
+34 2.9670e+07 388 3791 1677 9 1 2.4853e+05 0.00000000e+00 1.74287307e-01 1.00000000e+00 1.74287307e-01 -1.22616735e-11 6.51234497e-12 -3.44554210e-10 -3.65180086e-06
+35 2.9919e+07 388 3791 1677 8 1 2.8340e+05 0.00000000e+00 1.75065362e-01 1.00000000e+00 1.75065362e-01 7.33056340e-13 9.70834654e-13 -4.75799457e-10 -3.63515252e-06
+36 3.0202e+07 388 3791 1677 9 1 4.6112e+05 0.00000000e+00 1.75936875e-01 1.00000000e+00 1.75936875e-01 1.12309989e-11 1.96930306e-11 -6.33601745e-10 -3.61694822e-06
+37 3.0663e+07 388 3791 1677 11 0 4.7698e+05 0.00000000e+00 1.77303874e-01 1.00000000e+00 1.77303874e-01 8.68725694e-12 -4.71310290e-12 -8.96649169e-10 -3.58930780e-06
+38 3.1140e+07 388 3791 1677 12 1 4.6234e+05 0.00000000e+00 1.78659418e-01 1.00000000e+00 1.78659418e-01 -3.51732666e-12 -4.52989590e-11 -1.16719276e-09 -3.56296251e-06
+39 3.1603e+07 388 3791 1677 12 1 2.5667e+05 0.00000000e+00 1.79922290e-01 1.00000000e+00 1.79922290e-01 -1.19425959e-11 -4.00573823e-11 -1.42169805e-09 -3.53861151e-06
+40 3.1859e+07 388 3791 1677 9 0 2.0399e+05 0.00000000e+00 1.80618182e-01 1.00000000e+00 1.80618182e-01 -9.59545385e-12 -4.61477830e-11 -1.55864859e-09 -3.52505265e-06
+41 3.2063e+07 388 3791 1677 7 0 1.7538e+05 0.00000000e+00 1.81167103e-01 1.00000000e+00 1.81167103e-01 -7.16403886e-12 -4.92741764e-11 -1.66485886e-09 -3.51430190e-06
+42 3.2239e+07 388 3791 1677 6 1 1.9817e+05 0.00000000e+00 1.81635663e-01 1.00000000e+00 1.81635663e-01 -5.04783323e-12 -4.97806628e-11 -1.75413015e-09 -3.50509229e-06
+43 3.2437e+07 388 3791 1677 7 1 2.8177e+05 0.00000000e+00 1.82180205e-01 1.00000000e+00 1.82180205e-01 -7.13405146e-12 -3.52430498e-11 -1.85273229e-09 -3.49438069e-06
+44 3.2719e+07 388 3791 1677 9 0 4.4088e+05 0.00000000e+00 1.82981632e-01 1.00000000e+00 1.82981632e-01 -5.76380601e-14 -2.67541708e-11 -1.98899313e-09 -3.47879975e-06
+45 3.3160e+07 388 3791 1677 11 1 3.0398e+05 0.00000000e+00 1.84293277e-01 1.00000000e+00 1.84293277e-01 1.47641042e-11 -8.24090436e-12 -2.19364606e-09 -3.45344212e-06
+46 3.3464e+07 436 4228 1871 10 1 2.4461e+05 0.00000000e+00 1.85184997e-01 1.00000000e+00 1.85184997e-01 5.86419262e-12 1.24438344e-11 -2.32831301e-09 -3.43687790e-06
+47 3.3708e+07 436 4228 1871 8 0 1.6268e+05 0.00000000e+00 1.85872482e-01 1.00000000e+00 1.85872482e-01 1.24824660e-12 8.69824524e-12 -2.43270236e-09 -3.42462974e-06
+48 3.3871e+07 436 4228 1871 7 0 1.3298e+05 0.00000000e+00 1.86318284e-01 1.00000000e+00 1.86318284e-01 -1.51060752e-12 3.09200526e-12 -2.50004116e-09 -3.41685237e-06
+49 3.4004e+07 436 4228 1871 6 1 1.4498e+05 0.00000000e+00 1.86676875e-01 1.00000000e+00 1.86676875e-01 -3.98740187e-12 -2.01404008e-12 -2.55380520e-09 -3.41068627e-06
+50 3.4149e+07 436 4228 1871 6 1 1.7142e+05 0.00000000e+00 1.87065481e-01 1.00000000e+00 1.87065481e-01 -5.88080320e-12 9.16880273e-12 -2.61122210e-09 -3.40392725e-06
+51 3.4320e+07 427 4149 1836 6 0 2.1450e+05 0.00000000e+00 1.87523939e-01 1.00000000e+00 1.87523939e-01 -1.28030341e-12 1.33353450e-11 -2.67766186e-09 -3.39579683e-06
+52 3.4535e+07 427 4149 1836 7 0 2.8553e+05 0.00000000e+00 1.88094943e-01 1.00000000e+00 1.88094943e-01 7.81161288e-12 7.98456080e-12 -2.75868494e-09 -3.38551574e-06
+53 3.4820e+07 427 4149 1836 8 0 2.9603e+05 0.00000000e+00 1.88849906e-01 1.00000000e+00 1.88849906e-01 2.24982689e-11 -2.83873099e-12 -2.86287123e-09 -3.37169314e-06
+54 3.5116e+07 427 4149 1836 9 1 3.1390e+05 0.00000000e+00 1.89626214e-01 1.00000000e+00 1.89626214e-01 4.02689202e-11 -1.66358212e-11 -2.96644506e-09 -3.35725718e-06
+55 3.5430e+07 427 4149 1836 10 1 2.1763e+05 0.00000000e+00 1.90473484e-01 1.00000000e+00 1.90473484e-01 2.72072108e-11 -2.83295451e-11 -3.07150094e-09 -3.34185372e-06
+56 3.5648e+07 418 4077 1804 8 0 1.5425e+05 0.00000000e+00 1.91081244e-01 1.00000000e+00 1.91081244e-01 2.14818964e-11 -2.64783971e-11 -3.14181042e-09 -3.33105589e-06
+57 3.5802e+07 418 4077 1804 6 2 4.7292e+05 0.00000000e+00 1.91520990e-01 1.00000000e+00 1.91520990e-01 1.22851774e-11 -2.65322165e-11 -3.19061906e-09 -3.32326824e-06
+58 3.6275e+07 418 4077 1804 12 3 1.6445e+05 0.00000000e+00 1.92714263e-01 1.00000000e+00 1.92714263e-01 1.56904651e-11 2.51200923e-11 -3.33280266e-09 -3.30427112e-06
+59 3.6439e+07 418 4077 1804 8 1 1.5601e+05 0.00000000e+00 1.93125654e-01 1.00000000e+00 1.93125654e-01 -3.87234296e-12 2.75805023e-11 -3.37945288e-09 -3.29775566e-06
+60 3.6595e+07 418 4077 1804 6 1 2.0137e+05 0.00000000e+00 1.93518857e-01 1.00000000e+00 1.93518857e-01 8.10348680e-12 -1.66809927e-13 -3.42222032e-09 -3.29131541e-06
+61 3.6797e+07 418 4075 1803 7 1 4.6945e+05 0.00000000e+00 1.94036627e-01 1.00000000e+00 1.94036627e-01 9.25355248e-12 7.27717113e-12 -3.47587902e-09 -3.28288998e-06
+62 3.7266e+07 418 4075 1803 11 1 3.1725e+05 0.00000000e+00 1.95271718e-01 1.00000000e+00 1.95271718e-01 2.77707461e-11 -5.51994984e-12 -3.59424237e-09 -3.26263519e-06
+63 3.7584e+07 418 4075 1803 10 0 1.5055e+05 0.00000000e+00 1.96114975e-01 1.00000000e+00 1.96114975e-01 3.00735955e-11 -4.19463987e-12 -3.66902844e-09 -3.24883999e-06
+64 3.7734e+07 418 4075 1803 6 1 2.0245e+05 0.00000000e+00 1.96518443e-01 1.00000000e+00 1.96518443e-01 3.16766626e-11 -3.05660975e-12 -3.70319061e-09 -3.24226525e-06
+65 3.7937e+07 418 4075 1803 7 1 4.3739e+05 0.00000000e+00 1.97032804e-01 1.00000000e+00 1.97032804e-01 1.39466944e-11 -2.35725260e-12 -3.74751679e-09 -3.23395022e-06
+66 3.8374e+07 418 4068 1800 11 1 2.9298e+05 0.00000000e+00 1.98092584e-01 1.00000000e+00 1.98092584e-01 -5.60620547e-11 1.85406115e-11 -3.83554976e-09 -3.21665220e-06
+67 3.8667e+07 418 4068 1800 9 1 3.0664e+05 0.00000000e+00 1.98804045e-01 1.00000000e+00 1.98804045e-01 -6.48090517e-11 1.16259586e-11 -3.88794779e-09 -3.20483395e-06
+68 3.8974e+07 418 4068 1800 9 1 2.8032e+05 0.00000000e+00 1.99588015e-01 1.00000000e+00 1.99588015e-01 -4.69817892e-11 9.37088483e-12 -3.93793564e-09 -3.19177869e-06
+69 3.9254e+07 418 4068 1800 9 1 3.3946e+05 0.00000000e+00 2.00338377e-01 1.00000000e+00 2.00338377e-01 -1.40850778e-11 -5.31981828e-12 -3.98076308e-09 -3.17935667e-06
+70 3.9593e+07 418 4068 1800 10 1 4.5430e+05 0.00000000e+00 2.01257971e-01 1.00000000e+00 2.01257971e-01 2.95024361e-12 -3.21681160e-12 -4.02950633e-09 -3.16446314e-06
Added: trunk/aspect/tests/diffusion.prm
===================================================================
--- trunk/aspect/tests/diffusion.prm (rev 0)
+++ trunk/aspect/tests/diffusion.prm 2012-11-14 14:24:18 UTC (rev 1362)
@@ -0,0 +1,403 @@
+# A test by Juliane Dannberg
+#
+# This test verifies whether physical diffusion works as expected.
+# It sets up a halfspace heated from above and cooled from below
+# for which we know that (if infinite) the temperature should behave
+# as a time dependent function like
+# T(z,t) = (T_halfspace - T_0) * erf(z / (2*sqrt(kappa * t))) + T_0
+# where
+# T_0 = surface temperature (here =1)
+# T_halfspace = temperature at infinity (here =0)
+# kappa = thermal diffusivity (= k / (rho * c_p))
+# z = distance to surface
+# t = time
+#
+# The results produced by this test have been verified to be correct.
+
+##################################################################
+
+
+
+# Listing of Parameters
+# ---------------------
+
+# At the top, state that the program should be run for a two-dimensional
+# simulation
+set Dimension = 2
+
+
+# In computations, the time step $k$ is chosen according to $k = c \min_K
+# \frac{h_K}{\|u\|_{\infty,K} p_T}$ where $h_K$ is the diameter of cell $K$,
+# and the denominator is the maximal magnitude of the velocity on cell $K$
+# times the polynomial degree $p_T$ of the temperature discretization. The
+# dimensionless constant $c$ is called the CFL number in this program. For
+# time discretizations that have explicit components, $c$ must be less than a
+# constant that depends on the details of the time discretization and that is
+# no larger than one. On the other hand, for implicit discretizations such as
+# the one chosen here, one can choose the time step as large as one wants (in
+# particular, one can choose $c>1$) though a CFL number significantly larger
+# than one will yield rather diffusive solutions. Units: None.
+set CFL number = 0.01
+
+# The end time of the simulation. Units: years if the 'Use years in output
+# instead of seconds' parameter is set; seconds otherwise.
+set End time = 4e7
+
+# The name of the directory into which all output files should be placed. This
+# may be an absolute or a relative path.
+set Output directory = output-diffusion
+
+# A flag indicating whether the computation should be resumed from a
+# previously saved state (if true) or start from scratch (if false).
+set Resume computation = false
+
+# The start time of the simulation. Units: years if the 'Use years in output
+# instead of seconds' parameter is set; seconds otherwise.
+set Start time = 0
+
+# In order to make the problem in the first time step easier to solve, we need
+# a reasonable guess for the temperature and pressure. To obtain it, we use an
+# adiabatic pressure and temperature field. This parameter describes what the
+# `adiabatic' temperature would be at the surface of the domain (i.e. at depth
+# zero). Note that this value need not coincide with the boundary condition
+# posed at this point. Rather, the boundary condition may differ significantly
+# from the adiabatic value, and then typically induce a thermal boundary
+# layer.
+# For more information, see the section in the manual that discusses the
+# general mathematical model.
+set Adiabatic surface temperature = 0
+
+# The mathematical equations that describe thermal convection only determine
+# the pressure up to an arbitrary constant. On the other hand, for comparison
+# and for looking up material parameters it is important that the pressure be
+# normalized somehow. We do this by enforcing a particular average pressure
+# value at the surface of the domain, where the geometry model determines
+# where the surface is. This parameter describes what this average surface
+# pressure value is supposed to be. By default, it is set to zero, but one may
+# want to choose a different value for example for simulating only the volume
+# of the mantle below the lithosphere, in which case the surface pressure
+# should be the lithostatic pressure at the bottom of the lithosphere.
+# For more information, see the section in the manual that discusses the
+# general mathematical model.
+set Surface pressure = 0
+
+# When computing results for mantle convection simulations, it is often
+# difficult to judge the order of magnitude of results when they are stated in
+# MKS units involving seconds. Rather, some kinds of results such as
+# velocities are often stated in terms of meters per year (or, sometimes,
+# centimeters per year). On the other hand, for non-dimensional computations,
+# one wants results in their natural unit system as used inside the code. If
+# this flag is set to 'true' conversion to years happens; if it is 'false', no
+# such conversion happens.
+set Use years in output instead of seconds = false # default: true
+
+# Select one of the following solver schemes:
+#
+# 'IMPES': This is the IMplicit Pressure Explicit Saturation scheme. In a first
+# step the variables velocity and pressure are computed from equations without
+# time derivatives. In a second step the variable temperature is updated by
+# using an explicit time step.
+set Nonlinear solver scheme = IMPES
+
+
+subsection Boundary temperature model
+ # Select one of the following models:
+ #
+ # `box': A model in which the temperature is chosen constant on the left and
+ # right sides of a box.
+ #
+ # `spherical constant': A model in which the temperature is chosen constant
+ # on the inner and outer boundaries of a spherical shell. Parameters are
+ # read from subsection 'Sherical constant'.
+ set Model name = box
+
+ subsection Box
+ set Bottom temperature = 0
+ set Top temperature = 1
+ end
+end
+
+subsection Discretization
+ # The polynomial degree to use for the velocity variables in the Stokes
+ # system. Units: None.
+ set Stokes velocity polynomial degree = 2
+
+ # The polynomial degree to use for the temperature variable. Units: None.
+ set Temperature polynomial degree = 2
+
+ # Whether to use a Stokes discretization that is locally conservative at the
+ # expense of a larger number of degrees of freedom (true), or to go with a
+ # cheaper discretization that does not locally conserve mass, although it is
+ # globally conservative (false).
+ set Use locally conservative discretization = false
+
+ subsection Stabilization parameters
+ # The exponent $\alpha$ in the entropy viscosity stabilization. Units:
+ # None.
+ set alpha = 2
+
+ # The $\beta$ factor in the artificial viscosity stabilization. An
+ # appropriate value for 2d is 0.052 and 0.078 for 3d. Units: None.
+ set beta = 0.078
+
+ # The $c_R$ factor in the entropy viscosity stabilization. Units: None.
+ set cR = 0.5 # default: 0.11
+ end
+
+end
+
+
+subsection Geometry model
+ # Select one of the following models:
+ #
+ # `box': A box geometry parallel to the coordinate directions. The extent of
+ # the box in each coordinate direction is set in the parameter file.
+ #
+ # `spherical shell': A geometry representing a spherical shell or a pice of
+ # it. Inner and outer radii are read from the parameter file in subsection
+ # 'Spherical shell'.
+ set Model name = box
+
+ subsection Box
+ # Extent of the box in x-direction. Units: m.
+ set X extent = 1 # default: 1
+
+ # Extent of the box in y-direction. Units: m.
+ set Y extent = 1
+
+ # Extent of the box in z-direction. This value is ignored if the
+ # simulation is in 2d Units: m.
+ set Z extent = 1
+ end
+end
+
+
+subsection Gravity model
+ # Select one of the following models:
+ #
+ # `radial constant': A gravity model in which the gravity direction is
+ # radially inward and at constant magnitude. The magnitude is read from the
+ # parameter file in subsection 'Radial constant'.
+ #
+ # `radial earth-like': A gravity model in which the gravity direction is
+ # radially inward and with a magnitude that matches that of the earth at the
+ # core-mantle boundary as well as at the surface and in between is
+ # physically correct under the assumption of a constant density.
+ #
+ # `vertical': A gravity model in which the gravity direction is vertically
+ # downward and at constant magnitude.
+ set Model name = vertical
+
+end
+
+
+subsection Initial conditions
+ # Select one of the following models:
+ #
+ # `perturbed box': An initial temperature field in which the temperature is
+ # perturbed slightly from an otherwise constant value equal to one. The
+ # perturbation is chosen in such a way that the initial temperature is
+ # constant to one along the entire boundary.
+ #
+ # `spherical hexagonal perturbation': An initial temperature field in which
+ # the temperature is perturbed following a six-fold pattern in angular
+ # direction from an otherwise spherically symmetric state.
+ #
+ # `spherical gaussian perturbation': An initial temperature field in which
+ # the temperature is perturbed by a single Gaussian added to an otherwise
+ # spherically symmetric state. Additional parameters are read from the
+ # parameter file in subsection 'Spherical gaussian perturbation'.
+ set Model name = function
+
+ subsection Function
+ set Variable names = x,y
+ set Function expression = 0.0
+ end
+
+end
+
+
+subsection Material model
+ # Select one of the following models:
+ #
+ # `simple': A simple material model that has constant values for all
+ # coefficients but the density. This model uses the formulation that assumes
+ # an incompressible medium despite the fact that the density follows the law
+ # $\rho(T)=\rho_0(1-\beta(T-T_{\text{ref}})$. The value for the components
+ # of this formula and additional parameters are read from the parameter file
+ # in subsection 'Simple model'.
+ #
+ # `Steinberger': lookup from the paper of Steinberger/Calderwood
+ #
+ # `table': A material model that reads tables of pressure and temperature
+ # dependent material coefficients from files.
+ set Model name = simple
+
+ subsection Simple model
+ # Reference density $\rho_0$. Units: $kg/m^3$.
+ set Reference density = 1 # default: 3300
+
+ # The value of the specific heat $cp$. Units: $JG/kgK$.
+ set Reference specific heat = 1250
+
+ # The reference temperature $T_0$. Units: $K$.
+ set Reference temperature = 0 # default: 293
+
+ # The value of the thermal conductivity $k$. Units: $W/m/K$.
+ set Thermal conductivity = 1e-6 # default: 4.7
+
+ # The value of the thermal expansion coefficient $\beta$. Units: $1/K$.
+ set Thermal expansion coefficient = 2e-5
+
+ # The value of the constant viscosity. Units: $kg/m/s$.
+ set Viscosity = 1 # default: 5e24
+ end
+
+end
+
+
+subsection Mesh refinement
+ # A list of times so that if the end time of a time step is beyond this
+ # time, an additional round of mesh refinement is triggered. This is mostly
+ # useful to make sure we can get through the initial transient phase of a
+ # simulation on a relatively coarse mesh, and then refine again when we are
+ # in a time range that we are interested in and where we would like to use a
+ # finer mesh. Units: each element of the list has units years if the 'Use
+ # years in output instead of seconds' parameter is set; seconds otherwise.
+ set Additional refinement times =
+
+ # The number of adaptive refinement steps performed after initial global
+ # refinement but while still within the first time step.
+ set Initial adaptive refinement = 0 # default: 2
+
+ # The number of global refinement steps performed on the initial coarse
+ # mesh, before the problem is first solved there.
+ set Initial global refinement = 5 # default: 2
+
+ # The fraction of cells with the largest error that should be flagged for
+ # refinement.
+ set Refinement fraction = 0.3
+
+ # The fraction of cells with the smallest error that should be flagged for
+ # coarsening.
+ set Coarsening fraction = 0.05
+
+ # The method used to determine which cells to refine and which to coarsen.
+ set Strategy = Density c_p temperature
+
+ # The number of time steps after which the mesh is to be adapted again based
+ # on computed error indicators.
+ set Time steps between mesh refinement = 5 # default: 10
+end
+
+
+subsection Model settings
+ # Whether to include adiabatic heating into the model or not. From a
+ # physical viewpoint, adiabatic heating should always be used but may be
+ # undesirable when comparing results with known benchmarks that do not
+ # include this term in the temperature equation.
+ set Include adiabatic heating = false
+
+ # Whether to include shear heating into the model or not. From a physical
+ # viewpoint, shear heating should always be used but may be undesirable when
+ # comparing results with known benchmarks that do not include this term in
+ # the temperature equation.
+ set Include shear heating = false # default: true
+
+ # H0
+ set Radiogenic heating rate = 0
+
+ # A comma separated list of integers denoting those boundaries on which the
+ # temperature is fixed and described by the boundary temperature object
+ # selected in its own section of this input file. All boundary indicators
+ # used by the geometry but not explicitly listed here will end up with
+ # no-flux (insulating) boundary conditions.
+ set Fixed temperature boundary indicators = 2, 3
+
+ # A comma separated list of integers denoting those boundaries on which the
+ # velocity is tangential but prescribed, i.e., where external forces act to
+ # prescribe a particular velocity. This is often used to prescribe a
+ # velocity that equals that of overlying plates.
+ set Prescribed velocity boundary indicators =
+
+ # A comma separated list of integers denoting those boundaries on which the
+ # velocity is tangential and unrestrained, i.e., where no external forces
+ # act to prescribe a particular tangential velocity (although there is a
+ # force that requires the flow to be tangential).
+ set Tangential velocity boundary indicators = 0,1,2,3
+
+ # A comma separated list of integers denoting those boundaries on which the
+ # velocity is zero.
+ set Zero velocity boundary indicators =
+end
+
+subsection Checkpointing
+ # The wall time between performing checkpoints. If 0, will use the
+ # checkpoint step frequency instead.
+ set Time between checkpoint = 0
+ # The number of timesteps between performing checkpoints. If 0 and time
+ # between checkpoint is not specified, checkpointing will not be performed.
+ set Steps between checkpoint = 50
+end
+
+subsection Postprocess
+ # A comma separated list of postprocessor objects that should be run at the
+ # end of each time step. Some of these postprocessors will declare their own
+ # parameters which may, for example, include that they will actually do
+ # something only every so many time steps or years. Alternatively, the text
+ # 'all' indicates that all available postprocessors should be run after each
+ # time step.
+ #
+ # The following postprocessors are available:
+ #
+ # `depth average': A postprocessor that computes depth averaged quantities
+ # and writes them out.
+ #
+ # `heat flux statistics': A postprocessor that computes some statistics
+ # about the heat flux across boundaries.
+ #
+ # `heat flux statistics for the table model': A postprocessor that computes
+ # some statistics about the heat flux across boundaries.
+ #
+ # `velocity statistics for the table model': A postprocessor that computes
+ # some statistics about the velocity field.
+ #
+ # `temperature statistics': A postprocessor that computes some statistics
+ # about the temperature field.
+ #
+ # `velocity statistics': A postprocessor that computes some statistics about
+ # the velocity field.
+ #
+ # `visualization': A postprocessor that takes the solution and writes it
+ # into files that can be read by a graphical visualization program.
+ # Additional run time parameters are read from the parameter subsection
+ # 'Visualization'.
+ set List of postprocessors = temperature statistics, heat flux statistics
+
+ subsection Depth average
+ # The time interval between each generation of graphical output files. A
+ # value of zero indicates that output should be generated in each time
+ # step. Units: years if the 'Use years in output instead of seconds'
+ # parameter is set; seconds otherwise.
+ set Time between graphical output = 1
+ end
+
+ subsection Visualization
+ # VTU file output supports grouping files from several CPUs into one file
+ # using MPI I/O when writing on a parallel filesystem. Select 0 for no
+ # grouping. This will disable parallel file output and instead write one
+ # file per processor in a background thread. A value of 1 will generate
+ # one big file containing the whole solution.
+ set Number of grouped files = 0
+
+ # The file format to be used for graphical output.
+ set Output format = vtu
+ set List of output variables = density,viscosity
+
+ # The time interval between each generation of graphical output files. A
+ # value of zero indicates that output should be generated in each time
+ # step. Units: years if the 'Use years in output instead of seconds'
+ # parameter is set; seconds otherwise.
+ set Time between graphical output = 0 # default: 1e8
+ end
+end
More information about the CIG-COMMITS
mailing list