[cig-commits] [commit] master: Add a new test. (af6ce9c)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/8cd10da0b2e753ddd449a2bce24790abdb4af621...3407c2fccc8c98ef3939d3c2d29077b4d899ff4b

>---------------------------------------------------------------

commit af6ce9c287316bc9459577309b0c802eab231782
Author: Wolfgang Bangerth <bangerth at math.tamu.edu>
Date:   Mon May 19 14:59:05 2014 -0500

    Add a new test.


>---------------------------------------------------------------

af6ce9c287316bc9459577309b0c802eab231782
 tests/dynamic_topography.prm | 127 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 127 insertions(+)

diff --git a/tests/dynamic_topography.prm b/tests/dynamic_topography.prm
new file mode 100644
index 0000000..32ddaa1
--- /dev/null
+++ b/tests/dynamic_topography.prm
@@ -0,0 +1,127 @@
+# A test that verifies that the dynamic topography is calculated correctly
+
+set Dimension                              = 2
+set CFL number                             = 1.0
+set End time                               = 1e12
+set Start time                             = 0
+set Adiabatic surface temperature          = 1600
+set Pressure normalization                 = surface
+set Surface pressure                       = 0
+set Use years in output instead of seconds = false  # default: true
+set Nonlinear solver scheme                = IMPES
+
+
+subsection Boundary temperature model
+  set Model name = box
+end
+
+
+
+subsection Gravity model
+  set Model name = vertical
+
+  subsection Vertical
+    set Magnitude = 10
+  end
+end
+
+
+subsection Geometry model
+  set Model name = box
+
+  subsection Box
+    set X extent = 1800000
+    set Y extent = 700000
+  end
+end
+
+
+subsection Initial conditions
+  set Model name = function
+
+  subsection Function
+    set Variable names     = x,y
+    set Function expression = if(sqrt((x-900000)*(x-900000) + (y-300000)*(y-300000))<50000,1700,1600)
+  end
+end
+
+
+subsection Compositional fields
+  set Number of fields = 1
+end
+
+subsection Compositional initial conditions
+  set Model name = function
+
+  subsection Function
+    set Variable names      = x,y
+    set Function expression = if(y < 600000, 0.0, 1.0)
+  end
+end
+
+
+subsection Boundary temperature model
+  set Model name = box
+
+  subsection Box
+    set Bottom temperature  = 1600
+    set Left temperature    = 1600
+    set Right temperature   = 1600
+    set Top temperature     = 1600
+  end
+end
+
+
+subsection Material model
+  set Model name = simple
+
+  subsection Simple model
+    set Reference density             = 3300    # default: 3300
+    set Reference specific heat       = 1250
+    set Reference temperature         = 1600    # default: 293
+    set Thermal conductivity          = 0 # default: 4.7
+    set Thermal expansion coefficient = 3e-4
+    set Viscosity                     = 1e21    # default: 5e24
+    set Composition viscosity prefactor = 100
+    set Thermal viscosity exponent    = 36.84136
+  end
+end
+
+
+subsection Mesh refinement
+  set Initial adaptive refinement        = 3
+  set Initial global refinement          = 6
+  set Refinement fraction                = 0.3
+  set Coarsening fraction                = 0.00
+  set Strategy                           = temperature
+end
+
+
+subsection Model settings
+  set Fixed temperature boundary indicators   = 2,3
+  set Prescribed velocity boundary indicators =
+  set Tangential velocity boundary indicators = 0,1,3
+  set Zero velocity boundary indicators       = 2
+
+  set Include adiabatic heating           = false
+  set Include shear heating               = false
+end
+
+
+subsection Postprocess
+  set List of postprocessors = dynamic topography, visualization
+  
+  subsection Dynamic Topography
+    set Subtract mean of dynamic topography = true
+  end
+ 
+#  subsection Visualization
+#    set Time between graphical output = 0
+#    set List of output variables = density, viscosity, nonadiabatic temperature, dynamic topography
+#  end
+
+end
+
+
+set Output directory = output-dynamic_topography
+