Shrink input file.

tests/composition_active_nans.prm

@@ -6,127 +6,12 @@
 
 # MPI: 2
 
-set Dimension                              = 2
-set Start time                             = 0
-set End time                               = 0.1
-set Use years in output instead of seconds = false
-
-
-
-subsection Geometry model
-  set Model name = box
-
-  subsection Box
-    set X extent = 2
-    set Y extent = 1
-  end
-end
-
-
-# The parameters below this comment were created by the update script
-# as replacement for the old 'Model settings' subsection. They can be
-# safely merged with any existing subsections with the same name.
-
-subsection Boundary temperature model
-  set Fixed temperature boundary indicators   = 2, 3
-end
-
-subsection Boundary velocity model
-  set Tangential velocity boundary indicators = 0, 1, 2
-end
-
-subsection Boundary velocity model
-  set Prescribed velocity boundary indicators = 3: function
-end
-
-subsection Heating model
-  set List of model names = shear heating
-end
-
-subsection Boundary temperature model
-  set List of model names = box
-
-  subsection Box
-    set Bottom temperature = 1
-    set Top temperature    = 0
-  end
-end
-
-
-subsection Boundary velocity model
-  subsection Function
-    set Variable names      = x,z,t
-    set Function constants  = pi=3.1415926
-    set Function expression = if(x>1+sin(0.5*pi*t), 1, -1); 0
-  end
-end
-
-
-subsection Gravity model
-  set Model name = vertical
-end
-
-
-subsection Initial temperature model
-  set Model name = function
-
-  subsection Function
-    set Variable names      = x,z
-    set Function expression = (1-z)
-  end
-end
-
-
-# Compared to the passive material model, we here make
-# the density composition dependent by letting it depend
-# linearly on the value of the first compositional field.
-subsection Material model
-  set Model name = simple
-
-  subsection Simple model
-    set Thermal conductivity                           = 1e-6
-    set Thermal expansion coefficient                  = 0.01
-    set Viscosity                                      = 1
-    set Reference density                              = 1
-    set Reference temperature                          = 0
-    set Density differential for compositional field 1 = 100
-  end
-end
-
-
-subsection Mesh refinement
-  set Initial adaptive refinement        = 0
-  set Initial global refinement          = 3
-  set Time steps between mesh refinement = 0
-end
+include $ASPECT_SOURCE_DIR/tests/composition_active.prm
 
+set End time                               = 0.1
 
 subsection Postprocess
-  set List of postprocessors = visualization, temperature statistics, composition statistics
-
   subsection Visualization
-    set Interpolate output = false
     set List of output variables = density, artificial viscosity
-    set Time between graphical output = 0.1
-  end
-end
-
-
-# This is the new part: We declare that there will
-# be two compositional fields that will be
-# advected along. Their initial conditions are given by
-# a function that is one for the lowermost 0.2 height
-# units of the domain and zero otherwise in the first case,
-# and one in the top most 0.2 height units in the latter.
-subsection Compositional fields
-  set Number of fields = 2
-end
-
-subsection Initial composition model
-  set Model name = function
-
-  subsection Function
-    set Variable names      = x,y
-    set Function expression = if(y<0.2, 1, 0) ; if(y>0.8, 1, 0)
   end
 end