From 17ce44f68eb20090c986d59e21f56b902a295ec3 Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Sun, 2 Feb 2025 23:01:32 -0500 Subject: [PATCH] update parameters --- doc/parameter_view/parameters.xml | 622 +++++++++++---------- doc/sphinx/parameters/Material_20model.md | 6 +- doc/sphinx/parameters/Mesh_20refinement.md | 4 +- 3 files changed, 318 insertions(+), 314 deletions(-) diff --git a/doc/parameter_view/parameters.xml b/doc/parameter_view/parameters.xml index 6a290651f06..93c4dcb6c37 100644 --- a/doc/parameter_view/parameters.xml +++ b/doc/parameter_view/parameters.xml @@ -83,7 +83,7 @@ The number of space dimensions you want to run this program in. ASPECT can run i The end time of the simulation. The default value is a number so that when converted from years to seconds it is approximately equal to the largest number representable in floating point arithmetic. For all practical purposes, this equals infinity. Units: Years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise. -436 +438 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -11867,16 +11867,20 @@ The value of the constant viscosity. Units: \si{\pascal\second}. - - + +all:1 + + +all:1 + -A list of prefactors for the viscosity for each phase. The reference viscosity will be multiplied by this factor to get the corresponding viscosity for each phase. List must have one more entry than Phase transition depths. Units: non-dimensional. +A list of prefactors for the viscosity for each phase. The ``Viscosity'' parameter (modified by the ``Composition viscosity prefactor'', depending on composition) will be multiplied by this factor to get the corresponding viscosity for each phase. List must have the same number of entries as there are phases, that is one more than ``Phase transition depths'' for each composition that is used in the model. Units: non-dimensional. 711 -[List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] +[Anything] @@ -18644,7 +18648,7 @@ $ASPECT_SOURCE_DIR/data/geometry-model/initial-topography-model/ascii-data/test/ The name of a directory that contains the model data. This path may either be absolute (if starting with a `/') or relative to the current directory. The path may also include the special text `$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -133 +127 [DirectoryName] @@ -18661,7 +18665,7 @@ box_3d_%s.0.txt The file name of the model data. -134 +128 [Anything] @@ -18678,7 +18682,7 @@ The file name of the model data. Scalar factor, which is applied to the model data. You might want to use this to scale the input to a reference model. Another way to use this factor is to convert units of the input files. For instance, if you provide velocities in cm/yr set this factor to 0.01. -135 +129 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -18695,7 +18699,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -132 +135 [Anything] @@ -18714,7 +18718,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -131 +134 [Anything] @@ -18731,7 +18735,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -130 +133 [Anything] @@ -18750,7 +18754,7 @@ The names of the variables as they will be used in the function, separated by co The hillslope transport coefficient $\kappa$ used to diffuse the free surface, either as a stabilization step or to mimic erosional and depositional processes. Units: $\si{m^2/s}$. -127 +130 [Double 0...MAX_DOUBLE (inclusive)] @@ -18767,7 +18771,7 @@ The hillslope transport coefficient $\kappa$ used to diffuse the free surface, e The number of time steps between each application of diffusion. -128 +131 [Integer range 0...2147483647 (inclusive)] @@ -18803,7 +18807,7 @@ normal After each time step the free surface must be advected in the direction of the velocity field. Mass conservation requires that the mesh velocity is in the normal direction of the surface. However, for steep topography or large curvature, advection in the normal direction can become ill-conditioned, and instabilities in the mesh can form. Projection of the mesh velocity onto the local vertical direction can preserve the mesh quality better, but at the cost of slightly poorer mass conservation of the domain. -129 +132 [Selection normal|vertical ] @@ -19074,7 +19078,7 @@ This refinement criterion computes the gradient of the compositional field at qu On the other hand, for discontinuous finite elements (see the `Use discontinuous composition discretization' parameter in the `Discretization' section), the gradient at quadrature points does not include the contribution of jumps in the compositional field between cells, and consequently will not be an accurate approximation of the true gradient. As an extreme example, consider the case of using piecewise constant finite elements for compositional fields; in that case, the gradient of the solution at quadrature points inside each cell will always be exactly zero, even if the finite element solution is different from each cell to the next. Consequently, the current refinement criterion will likely not be useful in this situation. That said, the `composition approximate gradient' refinement criterion exists for exactly this purpose. -`composition threshold': A mesh refinement criterion that computes refinement indicators from the compositional fields. If any field exceeds the threshold given in the input file, the cell is marked for refinement. +`composition threshold': A mesh refinement criterion that computes refinement indicators from the compositional fields. One threshold per compositional is given in the input file, and if any field exceeds its threshold, the cell is marked for refinement. `density': A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the density, $\rho$. Because this quantity may not be a continuous function ($\rho$ and $C_p$ may be discontinuous functions along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1+d/2}$ where $h_K$ is the diameter of each cell and $d$ is the dimension. This scaling ensures that the error indicators converge to zero as $h_K\rightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$. @@ -19164,7 +19168,7 @@ A list of scaling factors by which every individual compositional field will be If the list of scaling factors given in this parameter is empty, then this indicates that they should all be chosen equal to 0. If the list is not empty then it needs to have as many entries as there are compositional fields. -430 +431 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -19181,7 +19185,7 @@ If the list of scaling factors given in this parameter is empty, then this indic A scaling factor for the artificial viscosity of the temperature equation. Use 0.0 to disable. -429 +430 [Double 0...MAX_DOUBLE (inclusive)] @@ -19198,7 +19202,7 @@ A comma separated list of names denoting those boundaries where there should be The names of the boundaries listed here can either be numbers (in which case they correspond to the numerical boundary indicators assigned by the geometry object), or they can correspond to any of the symbolic names the geometry object may have provided for each part of the boundary. You may want to compare this with the documentation of the geometry model you use in your model. -431 +432 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -19217,7 +19221,7 @@ The names of the boundaries listed here can either be numbers (in which case the The desired ratio between compaction length and size of the mesh cells, or, in other words, how many cells the mesh should (at least) have per compaction length. Every cell where this ratio is smaller than the value specified by this parameter (in places with fewer mesh cells per compaction length) is marked for refinement. -432 +424 [Double 0...MAX_DOUBLE (inclusive)] @@ -19234,7 +19238,7 @@ A list of scaling factors by which every individual compositional field will be If the list of scaling factors given in this parameter is empty, then this indicates that they should all be chosen equal to one. If the list is not empty then it needs to have as many entries as there are compositional fields. -420 +425 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -19251,7 +19255,7 @@ A list of scaling factors by which every individual compositional field gradient If the list of scaling factors given in this parameter is empty, then this indicates that they should all be chosen equal to one. If the list is not empty then it needs to have as many entries as there are compositional fields. -421 +426 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -19268,7 +19272,7 @@ A list of scaling factors by which every individual compositional field gradient If the list of scaling factors given in this parameter is empty, then this indicates that they should all be chosen equal to one. If the list is not empty then it needs to have as many entries as there are compositional fields. -422 +427 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -19280,10 +19284,10 @@ If the list of scaling factors given in this parameter is empty, then this indic -A list of thresholds that every individual compositional field will be evaluated against. +A list of thresholds, one for each compositional field to be evaluated against. -423 +428 [List of <[Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -19300,7 +19304,7 @@ A list of isosurfaces separated by semi-colons (;). Each isosurface entry consis The first two entries for each isosurface, describing the minimum and maximum grid levels, can be two numbers or contain one of the key values 'min' and 'max'. This indicates the key will be replaced with the global minimum and maximum refinement levels. The 'min' and 'max' keys also accept adding values to be added or subtracted from them respectively. This is done by adding a '+' or '-' and a number behind them (e.g. min+2 or max-1). Note that you can't subtract a value from a minimum value or add a value to the maximum value. If, for example, `max-4` drops below the minimum or `min+4` goes above the maximum, it will simply use the global minimum and maximum values respectively. The same holds for any mesh refinement level below the global minimum or above the global maximum. -424 +429 [Anything] @@ -19319,7 +19323,7 @@ depth A selection that determines the assumed coordinate system for the function variables. Allowed values are `depth', `cartesian' and `spherical'. `depth' will create a function, in which only the first variable is non-zero, which is interpreted to be the depth of the point. `spherical' coordinates are interpreted as r,phi or r,phi,theta in 2d/3d respectively with theta being the polar angle. -425 +414 [Selection depth|cartesian|spherical ] @@ -19334,7 +19338,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -428 +417 [Anything] @@ -19353,7 +19357,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -427 +416 [Anything] @@ -19370,7 +19374,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -426 +415 [Anything] @@ -19389,7 +19393,7 @@ depth A selection that determines the assumed coordinate system for the function variables. Allowed values are `depth', `cartesian' and `spherical'. `depth' will create a function, in which only the first variable is non-zero, which is interpreted to be the depth of the point. `spherical' coordinates are interpreted as r,phi or r,phi,theta in 2d/3d respectively with theta being the polar angle. -414 +418 [Selection depth|cartesian|spherical ] @@ -19404,7 +19408,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -417 +421 [Anything] @@ -19423,7 +19427,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -416 +420 [Anything] @@ -19440,7 +19444,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -415 +419 [Anything] @@ -19459,7 +19463,7 @@ absolute value What type of temperature anomaly should be considered when evaluating against the threshold: Only negative anomalies (negative only), only positive anomalies (positive only) or the absolute value of the nonadiabatic temperature. -419 +423 [Selection negative only|positive only|absolute value ] @@ -19476,7 +19480,7 @@ What type of temperature anomaly should be considered when evaluating against th A threshold that the nonadiabatic temperature will be evaluated against. Units: \si{\kelvin} -418 +422 [Double 0...MAX_DOUBLE (inclusive)] @@ -19539,7 +19543,7 @@ false By default, every cell needs to contain particles to use this interpolator plugin. If this parameter is set to true, cells are allowed to have no particles. In case both the current cell and its neighbors are empty, the interpolator will return 0 for the current cell's properties. -249 +254 [Bool] @@ -19568,7 +19572,7 @@ Select one of the following models: `rk4': Runge Kutta fourth order integrator, where $y_{n+1} = y_n + \frac{1}{6} k_1 + \frac{1}{3} k_2 + \frac{1}{3} k_3 + \frac{1}{6} k_4$ and $k_1$, $k_2$, $k_3$, $k_4$ are defined as usual. -242 +245 [Selection euler|rk2|rk4 ] @@ -19597,7 +19601,7 @@ Select one of the following models: `quadratic least squares': Interpolates particle properties onto a vector of points using a quadratic least squares method. Note that deal.II must be configured with BLAS/LAPACK. -244 +247 [Selection bilinear least squares|cell average|distance weighted average|harmonic average|nearest neighbor|quadratic least squares ] @@ -19650,7 +19654,7 @@ The following properties are available: `viscoplastic strain invariants': A plugin that calculates the finite strain invariant a particle has experienced and assigns it to either the plastic and/or viscous strain field based on whether the material is plastically yielding, or the total strain field used in the visco plastic material model. The implementation of this property is equivalent to the implementation for compositional fields that is located in the plugin in \texttt{benchmarks/buiter\_et\_al\_2008\_jgr/plugin/},and is effectively the same as what the visco plastic material model uses for compositional fields. -252 +255 [MultipleSelection composition|cpo bingham average|cpo elastic tensor|crystal preferred orientation|elastic stress|elastic tensor decomposition|function|grain size|initial composition|initial position|integrated strain|integrated strain invariant|melt particle|pT path|position|reference position|strain rate|velocity|viscoplastic strain invariants ] @@ -19667,7 +19671,7 @@ repartition Strategy that is used to balance the computational load across processors for adaptive meshes. -205 +208 [MultipleSelection none|remove particles|add particles|remove and add particles|repartition ] @@ -19684,7 +19688,7 @@ Strategy that is used to balance the computational load across processors for ad Upper limit for particle number per cell. This limit is useful for adaptive meshes to prevent coarse cells from slowing down the whole model. It will be checked and enforced after mesh refinement, after MPI transfer of particles and after particle movement. If there are \texttt{n\_number\_of\_particles} $>$ \texttt{max\_particles\_per\_cell} particles in one cell then \texttt{n\_number\_of\_particles} - \texttt{max\_particles\_per\_cell} particles in this cell are randomly chosen and destroyed. -207 +210 [Integer range 0...2147483647 (inclusive)] @@ -19701,7 +19705,7 @@ Upper limit for particle number per cell. This limit is useful for adaptive mesh Lower limit for particle number per cell. This limit is useful for adaptive meshes to prevent fine cells from being empty of particles. It will be checked and enforced after mesh refinement and after particle movement. If there are \texttt{n\_number\_of\_particles} $<$ \texttt{min\_particles\_per\_cell} particles in one cell then \texttt{min\_particles\_per\_cell} - \texttt{n\_number\_of\_particles} particles are generated and randomly placed in this cell. If the particles carry properties the individual property plugins control how the properties of the new particles are initialized. -206 +209 [Integer range 0...2147483647 (inclusive)] @@ -19749,7 +19753,7 @@ Select one of the following models: `uniform radial': Generate a uniform distribution of particles over a spherical domain in 2d or 3d. Uniform here means the particles will be generated with an equal spacing in each spherical spatial dimension, i.e., the particles are created at positions that increase linearly with equal spacing in radius, colatitude and longitude around a certain center point. Note that in order to produce a regular distribution the number of generated particles might not exactly match the one specified in the input file. -210 +213 [Selection ascii file|probability density function|quadrature points|random uniform|reference cell|uniform box|uniform radial ] @@ -19766,7 +19770,7 @@ Select one of the following models: Weight that is associated with the computational load of a single particle. The sum of particle weights will be added to the sum of cell weights to determine the partitioning of the mesh if the `repartition' particle load balancing strategy is selected. The optimal weight depends on the used integrator and particle properties. In general for a more expensive integrator and more expensive properties a larger particle weight is recommended. Before adding the weights of particles, each cell already carries a weight of 1000 to account for the cost of field-based computations. -208 +211 [Integer range 0...2147483647 (inclusive)] @@ -19783,7 +19787,7 @@ true Some particle interpolation algorithms require knowledge about particles in neighboring cells. To allow this, particles in ghost cells need to be exchanged between the processes neighboring this cell. This parameter determines whether this transport is happening. This parameter is deprecated and will be removed in the future. Ghost particle updates are always performed. Please set the parameter to `true'. -209 +212 [Bool] @@ -19801,7 +19805,7 @@ Some particle interpolation algorithms require knowledge about particles in neig This determines how many samples are taken when using the random draw volume averaging. Setting it to zero means that the number of samples is set to be equal to the number of grains. -255 +262 [Double 0...MAX_DOUBLE (inclusive)] @@ -19818,7 +19822,7 @@ This determines how many samples are taken when using the random draw volume ave The seed used to generate random numbers. This will make sure that results are reproducible as long as the problem is run with the same amount of MPI processes. It is implemented as final seed = Random number seed + MPI Rank. -254 +261 [Integer range 0...2147483647 (inclusive)] @@ -19837,7 +19841,7 @@ Spin tensor Options: Spin tensor -261 +268 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -19854,7 +19858,7 @@ Options: Spin tensor The number of grains of each different mineral each particle contains. -257 +264 [Integer range 1...2147483647 (inclusive)] @@ -19871,7 +19875,7 @@ The number of grains of each different mineral each particle contains. The Backward Euler property advection method involve internal iterations. This option allows for setting the maximum number of iterations. Note that when the iteration is ended by the max iteration amount an assert is thrown. -260 +267 [Integer range 0...2147483647 (inclusive)] @@ -19888,7 +19892,7 @@ Backward Euler Options: Forward Euler, Backward Euler -258 +265 [Anything] @@ -19905,7 +19909,7 @@ Options: Forward Euler, Backward Euler The Backward Euler property advection method involve internal iterations. This option allows for setting a tolerance. When the norm of tensor new - tensor old is smaller than this tolerance, the iteration is stopped. -259 +266 [Double 0...MAX_DOUBLE (inclusive)] @@ -19922,7 +19926,7 @@ The Backward Euler property advection method involve internal iterations. This o The seed used to generate random numbers. This will make sure that results are reproducible as long as the problem is run with the same number of MPI processes. It is implemented as final seed = user seed + MPI Rank. -256 +263 [Integer range 0...2147483647 (inclusive)] @@ -19940,7 +19944,7 @@ The seed used to generate random numbers. This will make sure that results are r This is exponent p as defined in equation 11 of Kaminski et al., 2004. -268 +275 [Double 0...MAX_DOUBLE (inclusive)] @@ -19957,7 +19961,7 @@ This is exponent p as defined in equation 11 of Kaminski et al., 2004. The dimensionless intrinsic grain boundary mobility for both olivine and enstatite. -265 +272 [Double 0...MAX_DOUBLE (inclusive)] @@ -19974,7 +19978,7 @@ The dimensionless intrinsic grain boundary mobility for both olivine and enstati This is the dimensionless nucleation rate as defined in equation 8 of Kaminski et al., 2004. -269 +276 [Double 0...MAX_DOUBLE (inclusive)] @@ -19991,7 +19995,7 @@ This is the dimensionless nucleation rate as defined in equation 8 of Kaminski e This is the power law exponent that characterizes the rheology of the slip systems. It is used in equation 11 of Kaminski et al., 2004. -267 +274 [Double 0...MAX_DOUBLE (inclusive)] @@ -20008,7 +20012,7 @@ This is the power law exponent that characterizes the rheology of the slip syste The Dimensionless Grain Boundary Sliding (GBS) threshold. This is a grain size threshold below which grain deform by GBS and become strain-free grains. -270 +277 [Double 0...MAX_DOUBLE (inclusive)] @@ -20025,7 +20029,7 @@ The Dimensionless Grain Boundary Sliding (GBS) threshold. This is a grain size t The volume fraction for the different minerals. There need to be the same amount of values as there are minerals -266 +273 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -20044,7 +20048,7 @@ Olivine: Karato 2008, Enstatite This determines what minerals and fabrics or fabric selectors are used used for the LPO/CPO calculation. The options are Olivine: Passive, A-fabric, Olivine: B-fabric, Olivine: C-fabric, Olivine: D-fabric, Olivine: E-fabric, Olivine: Karato 2008 or Enstatite. Passive sets all RRSS entries to the maximum. The Karato 2008 selector selects a fabric based on stress and water content as defined in figure 4 of the Karato 2008 review paper (doi: 10.1146/annurev.earth.36.031207.124120). -263 +270 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -20061,7 +20065,7 @@ Uniform grains and random uniform rotations The model used to initialize the CPO for all particles. Currently 'Uniform grains and random uniform rotations' and 'World Builder' are the only valid option. -262 +269 [Anything] @@ -20078,7 +20082,7 @@ The model used to initialize the CPO for all particles. Currently 'Uniform The volume fractions for the different minerals. There need to be the same number of values as there are minerals.Note that the currently implemented scheme is incompressible and does not allow chemical interaction or the formation of new phases -264 +271 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -20096,7 +20100,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -274 +259 [Anything] @@ -20115,7 +20119,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -273 +258 [Anything] @@ -20132,7 +20136,7 @@ If the function you are describing represents a vector-valued function with mult The number of function components where each component is described by a function expression delimited by a ';'. -271 +256 [Integer range 0...2147483647 (inclusive)] @@ -20149,7 +20153,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -272 +257 [Anything] @@ -20169,7 +20173,7 @@ $ASPECT_SOURCE_DIR/data/particle/generator/ascii/ The name of a directory that contains the particle data. This path may either be absolute (if starting with a '/') or relative to the current directory. The path may also include the special text '$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -233 +236 [DirectoryName] @@ -20186,7 +20190,7 @@ particle.dat The name of the particle file. -234 +237 [Anything] @@ -20203,7 +20207,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -237 +240 [Anything] @@ -20220,7 +20224,7 @@ A typical example would be to set this runtime parameter to `pi=3.1415926536&apo The formula that denotes the spatially variable probability density function. This expression may contain any of the usual operations such as addition or multiplication, as well as all of the common functions such as `sin' or `cos'. In addition, it may contain expressions like `if(x>0, 1, 0)' where the expression evaluates to the second argument if the first argument is true, and to the third argument otherwise; this example would result in no particles at all in that part of the domain where $x==0$, and a constant particle density in the rest of the domain. For a full overview of possible expressions accepted see the documentation of the muparser library at http://muparser.beltoforion.de/. Note that the function has to be non-negative everywhere in the domain, and needs to be positive in at least some parts of the domain. -238 +241 [Anything] @@ -20237,7 +20241,7 @@ The formula that denotes the spatially variable probability density function. Th Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -239 +242 [Double 0...MAX_DOUBLE (inclusive)] @@ -20254,7 +20258,7 @@ true If true, particle numbers per cell are calculated randomly according to their respective probability density. This means particle numbers per cell can deviate statistically from the integral of the probability density. If false, first determine how many particles each cell should have based on the integral of the density over each of the cells, and then once we know how many particles we want on each cell, choose their locations randomly within each cell. -240 +243 [Bool] @@ -20271,7 +20275,7 @@ If true, particle numbers per cell are calculated randomly according to their re The seed for the random number generator that controls the particle generation. Keep constant to generate identical particle distributions in subsequent model runs. Change to get a different distribution. In parallel computations the seed is further modified on each process to ensure different particle patterns on different processes. Note that the number of particles per processor is not affected by the seed. -241 +244 [Integer range 0...2147483647 (inclusive)] @@ -20288,7 +20292,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -235 +238 [Anything] @@ -20307,7 +20311,7 @@ The names of the variables as they will be used in the function, separated by co Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -211 +214 [Double 0...MAX_DOUBLE (inclusive)] @@ -20324,7 +20328,7 @@ true If true, particle numbers per cell are calculated randomly according to their respective probability density. This means particle numbers per cell can deviate statistically from the integral of the probability density. If false, first determine how many particles each cell should have based on the integral of the density over each of the cells, and then once we know how many particles we want on each cell, choose their locations randomly within each cell. -212 +215 [Bool] @@ -20341,7 +20345,7 @@ If true, particle numbers per cell are calculated randomly according to their re The seed for the random number generator that controls the particle generation. Keep constant to generate identical particle distributions in subsequent model runs. Change to get a different distribution. In parallel computations the seed is further modified on each process to ensure different particle patterns on different processes. Note that the number of particles per processor is not affected by the seed. -213 +216 [Integer range 0...2147483647 (inclusive)] @@ -20360,7 +20364,7 @@ The seed for the random number generator that controls the particle generation. List of number of particles to create per cell and spatial dimension. The size of the list is the number of spatial dimensions. If only one value is given, then each spatial dimension is set to the same value. The list of numbers are parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -214 +217 [List of <[Integer range 1...2147483647 (inclusive)]> of length 0...4294967295 (inclusive)] @@ -20379,7 +20383,7 @@ List of number of particles to create per cell and spatial dimension. The size o Maximum x coordinate for the region of particles. -217 +220 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20396,7 +20400,7 @@ Maximum x coordinate for the region of particles. Maximum y coordinate for the region of particles. -219 +222 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20413,7 +20417,7 @@ Maximum y coordinate for the region of particles. Maximum z coordinate for the region of particles. -221 +224 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20430,7 +20434,7 @@ Maximum z coordinate for the region of particles. Minimum x coordinate for the region of particles. -216 +219 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20447,7 +20451,7 @@ Minimum x coordinate for the region of particles. Minimum y coordinate for the region of particles. -218 +221 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20464,7 +20468,7 @@ Minimum y coordinate for the region of particles. Minimum z coordinate for the region of particles. -220 +223 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20481,7 +20485,7 @@ Minimum z coordinate for the region of particles. Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -215 +218 [Double 0...MAX_DOUBLE (inclusive)] @@ -20500,7 +20504,7 @@ Total number of particles to create (not per processor or per element). The numb x coordinate for the center of the spherical region, where particles are generated. -223 +226 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20517,7 +20521,7 @@ x coordinate for the center of the spherical region, where particles are generat y coordinate for the center of the spherical region, where particles are generated. -224 +227 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20534,7 +20538,7 @@ y coordinate for the center of the spherical region, where particles are generat z coordinate for the center of the spherical region, where particles are generated. -225 +228 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20551,7 +20555,7 @@ z coordinate for the center of the spherical region, where particles are generat Maximum latitude coordinate for the region of particles in degrees. Measured from the center position, and from the north pole. -231 +234 [Double 0...180 (inclusive)] @@ -20568,7 +20572,7 @@ Maximum latitude coordinate for the region of particles in degrees. Measured fro Maximum longitude coordinate for the region of particles in degrees. Measured from the center position. -229 +232 [Double -180...360 (inclusive)] @@ -20585,7 +20589,7 @@ Maximum longitude coordinate for the region of particles in degrees. Measured fr Maximum radial coordinate for the region of particles. Measured from the center position. -227 +230 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -20602,7 +20606,7 @@ Maximum radial coordinate for the region of particles. Measured from the center Minimum latitude coordinate for the region of particles in degrees. Measured from the center position, and from the north pole. -230 +233 [Double 0...180 (inclusive)] @@ -20619,7 +20623,7 @@ Minimum latitude coordinate for the region of particles in degrees. Measured fro Minimum longitude coordinate for the region of particles in degrees. Measured from the center position. -228 +231 [Double -180...360 (inclusive)] @@ -20636,7 +20640,7 @@ Minimum longitude coordinate for the region of particles in degrees. Measured fr Minimum radial coordinate for the region of particles. Measured from the center position. -226 +229 [Double 0...MAX_DOUBLE (inclusive)] @@ -20653,7 +20657,7 @@ Minimum radial coordinate for the region of particles. Measured from the center Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -222 +225 [Double 0...MAX_DOUBLE (inclusive)] @@ -20670,7 +20674,7 @@ Total number of particles to create (not per processor or per element). The numb The number of radial shells of particles that will be generated around the central point. -232 +235 [Integer range 1...2147483647 (inclusive)] @@ -20691,7 +20695,7 @@ true Whether to correctly evaluate old and current velocity solution to reach higher-order accuracy in time. If set to 'false' only the old velocity solution is evaluated to simulate a first order method in time. This is only recommended for benchmark purposes. -243 +246 [Bool] @@ -20712,7 +20716,7 @@ false Extends the range used by 'Use linear least squares limiter' by linearly interpolating values at cell boundaries from neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. Enabling 'Use boundary extrapolation' requires enabling 'Use linear least squares limiter'. -246 +251 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -20729,7 +20733,7 @@ true Limit the interpolation of particle properties onto the cell, so that the value of each property is no smaller than its minimum and no larger than its maximum on the particles of each cell, and the average of neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. -245 +250 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -20748,7 +20752,7 @@ false Extends the range used by 'Use quadratic least squares limiter' by linearly interpolating values at cell boundaries from neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. Enabling 'Use boundary extrapolation' requires enabling 'Use quadratic least squares limiter'. -251 +249 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -20765,7 +20769,7 @@ true Limit the interpolation of particle properties onto the cell, so that the value of each property is no smaller than its minimum and no larger than its maximum on the particles of each cell, and the average of neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. -250 +248 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -20785,7 +20789,7 @@ Limit the interpolation of particle properties onto the cell, so that the value The minimum porosity that has to be present at the position of a particle for it to be considered a melt particle (in the sense that the melt presence property is set to 1). -253 +260 [Double 0...1 (inclusive)] @@ -20805,7 +20809,7 @@ false By default, every cell needs to contain particles to use this interpolator plugin. If this parameter is set to true, cells are allowed to have no particles. In case both the current cell and its neighbors are empty, the interpolator will return 0 for the current cell's properties. -319 +324 [Bool] @@ -20834,7 +20838,7 @@ Select one of the following models: `rk4': Runge Kutta fourth order integrator, where $y_{n+1} = y_n + \frac{1}{6} k_1 + \frac{1}{3} k_2 + \frac{1}{3} k_3 + \frac{1}{6} k_4$ and $k_1$, $k_2$, $k_3$, $k_4$ are defined as usual. -312 +315 [Selection euler|rk2|rk4 ] @@ -20863,7 +20867,7 @@ Select one of the following models: `quadratic least squares': Interpolates particle properties onto a vector of points using a quadratic least squares method. Note that deal.II must be configured with BLAS/LAPACK. -314 +317 [Selection bilinear least squares|cell average|distance weighted average|harmonic average|nearest neighbor|quadratic least squares ] @@ -20916,7 +20920,7 @@ The following properties are available: `viscoplastic strain invariants': A plugin that calculates the finite strain invariant a particle has experienced and assigns it to either the plastic and/or viscous strain field based on whether the material is plastically yielding, or the total strain field used in the visco plastic material model. The implementation of this property is equivalent to the implementation for compositional fields that is located in the plugin in \texttt{benchmarks/buiter\_et\_al\_2008\_jgr/plugin/},and is effectively the same as what the visco plastic material model uses for compositional fields. -322 +325 [MultipleSelection composition|cpo bingham average|cpo elastic tensor|crystal preferred orientation|elastic stress|elastic tensor decomposition|function|grain size|initial composition|initial position|integrated strain|integrated strain invariant|melt particle|pT path|position|reference position|strain rate|velocity|viscoplastic strain invariants ] @@ -20933,7 +20937,7 @@ repartition Strategy that is used to balance the computational load across processors for adaptive meshes. -275 +278 [MultipleSelection none|remove particles|add particles|remove and add particles|repartition ] @@ -20950,7 +20954,7 @@ Strategy that is used to balance the computational load across processors for ad Upper limit for particle number per cell. This limit is useful for adaptive meshes to prevent coarse cells from slowing down the whole model. It will be checked and enforced after mesh refinement, after MPI transfer of particles and after particle movement. If there are \texttt{n\_number\_of\_particles} $>$ \texttt{max\_particles\_per\_cell} particles in one cell then \texttt{n\_number\_of\_particles} - \texttt{max\_particles\_per\_cell} particles in this cell are randomly chosen and destroyed. -277 +280 [Integer range 0...2147483647 (inclusive)] @@ -20967,7 +20971,7 @@ Upper limit for particle number per cell. This limit is useful for adaptive mesh Lower limit for particle number per cell. This limit is useful for adaptive meshes to prevent fine cells from being empty of particles. It will be checked and enforced after mesh refinement and after particle movement. If there are \texttt{n\_number\_of\_particles} $<$ \texttt{min\_particles\_per\_cell} particles in one cell then \texttt{min\_particles\_per\_cell} - \texttt{n\_number\_of\_particles} particles are generated and randomly placed in this cell. If the particles carry properties the individual property plugins control how the properties of the new particles are initialized. -276 +279 [Integer range 0...2147483647 (inclusive)] @@ -20998,7 +21002,7 @@ Select one of the following models: `uniform radial': Generate a uniform distribution of particles over a spherical domain in 2d or 3d. Uniform here means the particles will be generated with an equal spacing in each spherical spatial dimension, i.e., the particles are created at positions that increase linearly with equal spacing in radius, colatitude and longitude around a certain center point. Note that in order to produce a regular distribution the number of generated particles might not exactly match the one specified in the input file. -280 +283 [Selection ascii file|probability density function|quadrature points|random uniform|reference cell|uniform box|uniform radial ] @@ -21015,7 +21019,7 @@ Select one of the following models: Weight that is associated with the computational load of a single particle. The sum of particle weights will be added to the sum of cell weights to determine the partitioning of the mesh if the `repartition' particle load balancing strategy is selected. The optimal weight depends on the used integrator and particle properties. In general for a more expensive integrator and more expensive properties a larger particle weight is recommended. Before adding the weights of particles, each cell already carries a weight of 1000 to account for the cost of field-based computations. -278 +281 [Integer range 0...2147483647 (inclusive)] @@ -21032,7 +21036,7 @@ true Some particle interpolation algorithms require knowledge about particles in neighboring cells. To allow this, particles in ghost cells need to be exchanged between the processes neighboring this cell. This parameter determines whether this transport is happening. This parameter is deprecated and will be removed in the future. Ghost particle updates are always performed. Please set the parameter to `true'. -279 +282 [Bool] @@ -21050,7 +21054,7 @@ Some particle interpolation algorithms require knowledge about particles in neig This determines how many samples are taken when using the random draw volume averaging. Setting it to zero means that the number of samples is set to be equal to the number of grains. -325 +332 [Double 0...MAX_DOUBLE (inclusive)] @@ -21067,7 +21071,7 @@ This determines how many samples are taken when using the random draw volume ave The seed used to generate random numbers. This will make sure that results are reproducible as long as the problem is run with the same amount of MPI processes. It is implemented as final seed = Random number seed + MPI Rank. -324 +331 [Integer range 0...2147483647 (inclusive)] @@ -21086,7 +21090,7 @@ Spin tensor Options: Spin tensor -331 +338 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -21103,7 +21107,7 @@ Options: Spin tensor The number of grains of each different mineral each particle contains. -327 +334 [Integer range 1...2147483647 (inclusive)] @@ -21120,7 +21124,7 @@ The number of grains of each different mineral each particle contains. The Backward Euler property advection method involve internal iterations. This option allows for setting the maximum number of iterations. Note that when the iteration is ended by the max iteration amount an assert is thrown. -330 +337 [Integer range 0...2147483647 (inclusive)] @@ -21137,7 +21141,7 @@ Backward Euler Options: Forward Euler, Backward Euler -328 +335 [Anything] @@ -21154,7 +21158,7 @@ Options: Forward Euler, Backward Euler The Backward Euler property advection method involve internal iterations. This option allows for setting a tolerance. When the norm of tensor new - tensor old is smaller than this tolerance, the iteration is stopped. -329 +336 [Double 0...MAX_DOUBLE (inclusive)] @@ -21171,7 +21175,7 @@ The Backward Euler property advection method involve internal iterations. This o The seed used to generate random numbers. This will make sure that results are reproducible as long as the problem is run with the same number of MPI processes. It is implemented as final seed = user seed + MPI Rank. -326 +333 [Integer range 0...2147483647 (inclusive)] @@ -21189,7 +21193,7 @@ The seed used to generate random numbers. This will make sure that results are r This is exponent p as defined in equation 11 of Kaminski et al., 2004. -338 +345 [Double 0...MAX_DOUBLE (inclusive)] @@ -21206,7 +21210,7 @@ This is exponent p as defined in equation 11 of Kaminski et al., 2004. The dimensionless intrinsic grain boundary mobility for both olivine and enstatite. -335 +342 [Double 0...MAX_DOUBLE (inclusive)] @@ -21223,7 +21227,7 @@ The dimensionless intrinsic grain boundary mobility for both olivine and enstati This is the dimensionless nucleation rate as defined in equation 8 of Kaminski et al., 2004. -339 +346 [Double 0...MAX_DOUBLE (inclusive)] @@ -21240,7 +21244,7 @@ This is the dimensionless nucleation rate as defined in equation 8 of Kaminski e This is the power law exponent that characterizes the rheology of the slip systems. It is used in equation 11 of Kaminski et al., 2004. -337 +344 [Double 0...MAX_DOUBLE (inclusive)] @@ -21257,7 +21261,7 @@ This is the power law exponent that characterizes the rheology of the slip syste The Dimensionless Grain Boundary Sliding (GBS) threshold. This is a grain size threshold below which grain deform by GBS and become strain-free grains. -340 +347 [Double 0...MAX_DOUBLE (inclusive)] @@ -21274,7 +21278,7 @@ The Dimensionless Grain Boundary Sliding (GBS) threshold. This is a grain size t The volume fraction for the different minerals. There need to be the same amount of values as there are minerals -336 +343 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -21293,7 +21297,7 @@ Olivine: Karato 2008, Enstatite This determines what minerals and fabrics or fabric selectors are used used for the LPO/CPO calculation. The options are Olivine: Passive, A-fabric, Olivine: B-fabric, Olivine: C-fabric, Olivine: D-fabric, Olivine: E-fabric, Olivine: Karato 2008 or Enstatite. Passive sets all RRSS entries to the maximum. The Karato 2008 selector selects a fabric based on stress and water content as defined in figure 4 of the Karato 2008 review paper (doi: 10.1146/annurev.earth.36.031207.124120). -333 +340 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -21310,7 +21314,7 @@ Uniform grains and random uniform rotations The model used to initialize the CPO for all particles. Currently 'Uniform grains and random uniform rotations' and 'World Builder' are the only valid option. -332 +339 [Anything] @@ -21327,7 +21331,7 @@ The model used to initialize the CPO for all particles. Currently 'Uniform The volume fractions for the different minerals. There need to be the same number of values as there are minerals.Note that the currently implemented scheme is incompressible and does not allow chemical interaction or the formation of new phases -334 +341 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -21345,7 +21349,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -344 +329 [Anything] @@ -21364,7 +21368,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -343 +328 [Anything] @@ -21381,7 +21385,7 @@ If the function you are describing represents a vector-valued function with mult The number of function components where each component is described by a function expression delimited by a ';'. -341 +326 [Integer range 0...2147483647 (inclusive)] @@ -21398,7 +21402,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -342 +327 [Anything] @@ -21418,7 +21422,7 @@ $ASPECT_SOURCE_DIR/data/particle/generator/ascii/ The name of a directory that contains the particle data. This path may either be absolute (if starting with a '/') or relative to the current directory. The path may also include the special text '$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -303 +306 [DirectoryName] @@ -21435,7 +21439,7 @@ particle.dat The name of the particle file. -304 +307 [Anything] @@ -21452,7 +21456,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -307 +310 [Anything] @@ -21469,7 +21473,7 @@ A typical example would be to set this runtime parameter to `pi=3.1415926536&apo The formula that denotes the spatially variable probability density function. This expression may contain any of the usual operations such as addition or multiplication, as well as all of the common functions such as `sin' or `cos'. In addition, it may contain expressions like `if(x>0, 1, 0)' where the expression evaluates to the second argument if the first argument is true, and to the third argument otherwise; this example would result in no particles at all in that part of the domain where $x==0$, and a constant particle density in the rest of the domain. For a full overview of possible expressions accepted see the documentation of the muparser library at http://muparser.beltoforion.de/. Note that the function has to be non-negative everywhere in the domain, and needs to be positive in at least some parts of the domain. -308 +311 [Anything] @@ -21486,7 +21490,7 @@ The formula that denotes the spatially variable probability density function. Th Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -309 +312 [Double 0...MAX_DOUBLE (inclusive)] @@ -21503,7 +21507,7 @@ true If true, particle numbers per cell are calculated randomly according to their respective probability density. This means particle numbers per cell can deviate statistically from the integral of the probability density. If false, first determine how many particles each cell should have based on the integral of the density over each of the cells, and then once we know how many particles we want on each cell, choose their locations randomly within each cell. -310 +313 [Bool] @@ -21520,7 +21524,7 @@ If true, particle numbers per cell are calculated randomly according to their re The seed for the random number generator that controls the particle generation. Keep constant to generate identical particle distributions in subsequent model runs. Change to get a different distribution. In parallel computations the seed is further modified on each process to ensure different particle patterns on different processes. Note that the number of particles per processor is not affected by the seed. -311 +314 [Integer range 0...2147483647 (inclusive)] @@ -21537,7 +21541,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -305 +308 [Anything] @@ -21556,7 +21560,7 @@ The names of the variables as they will be used in the function, separated by co Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -281 +284 [Double 0...MAX_DOUBLE (inclusive)] @@ -21573,7 +21577,7 @@ true If true, particle numbers per cell are calculated randomly according to their respective probability density. This means particle numbers per cell can deviate statistically from the integral of the probability density. If false, first determine how many particles each cell should have based on the integral of the density over each of the cells, and then once we know how many particles we want on each cell, choose their locations randomly within each cell. -282 +285 [Bool] @@ -21590,7 +21594,7 @@ If true, particle numbers per cell are calculated randomly according to their re The seed for the random number generator that controls the particle generation. Keep constant to generate identical particle distributions in subsequent model runs. Change to get a different distribution. In parallel computations the seed is further modified on each process to ensure different particle patterns on different processes. Note that the number of particles per processor is not affected by the seed. -283 +286 [Integer range 0...2147483647 (inclusive)] @@ -21609,7 +21613,7 @@ The seed for the random number generator that controls the particle generation. List of number of particles to create per cell and spatial dimension. The size of the list is the number of spatial dimensions. If only one value is given, then each spatial dimension is set to the same value. The list of numbers are parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -284 +287 [List of <[Integer range 1...2147483647 (inclusive)]> of length 0...4294967295 (inclusive)] @@ -21628,7 +21632,7 @@ List of number of particles to create per cell and spatial dimension. The size o Maximum x coordinate for the region of particles. -287 +290 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21645,7 +21649,7 @@ Maximum x coordinate for the region of particles. Maximum y coordinate for the region of particles. -289 +292 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21662,7 +21666,7 @@ Maximum y coordinate for the region of particles. Maximum z coordinate for the region of particles. -291 +294 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21679,7 +21683,7 @@ Maximum z coordinate for the region of particles. Minimum x coordinate for the region of particles. -286 +289 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21696,7 +21700,7 @@ Minimum x coordinate for the region of particles. Minimum y coordinate for the region of particles. -288 +291 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21713,7 +21717,7 @@ Minimum y coordinate for the region of particles. Minimum z coordinate for the region of particles. -290 +293 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21730,7 +21734,7 @@ Minimum z coordinate for the region of particles. Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -285 +288 [Double 0...MAX_DOUBLE (inclusive)] @@ -21749,7 +21753,7 @@ Total number of particles to create (not per processor or per element). The numb x coordinate for the center of the spherical region, where particles are generated. -293 +296 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21766,7 +21770,7 @@ x coordinate for the center of the spherical region, where particles are generat y coordinate for the center of the spherical region, where particles are generated. -294 +297 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21783,7 +21787,7 @@ y coordinate for the center of the spherical region, where particles are generat z coordinate for the center of the spherical region, where particles are generated. -295 +298 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21800,7 +21804,7 @@ z coordinate for the center of the spherical region, where particles are generat Maximum latitude coordinate for the region of particles in degrees. Measured from the center position, and from the north pole. -301 +304 [Double 0...180 (inclusive)] @@ -21817,7 +21821,7 @@ Maximum latitude coordinate for the region of particles in degrees. Measured fro Maximum longitude coordinate for the region of particles in degrees. Measured from the center position. -299 +302 [Double -180...360 (inclusive)] @@ -21834,7 +21838,7 @@ Maximum longitude coordinate for the region of particles in degrees. Measured fr Maximum radial coordinate for the region of particles. Measured from the center position. -297 +300 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -21851,7 +21855,7 @@ Maximum radial coordinate for the region of particles. Measured from the center Minimum latitude coordinate for the region of particles in degrees. Measured from the center position, and from the north pole. -300 +303 [Double 0...180 (inclusive)] @@ -21868,7 +21872,7 @@ Minimum latitude coordinate for the region of particles in degrees. Measured fro Minimum longitude coordinate for the region of particles in degrees. Measured from the center position. -298 +301 [Double -180...360 (inclusive)] @@ -21885,7 +21889,7 @@ Minimum longitude coordinate for the region of particles in degrees. Measured fr Minimum radial coordinate for the region of particles. Measured from the center position. -296 +299 [Double 0...MAX_DOUBLE (inclusive)] @@ -21902,7 +21906,7 @@ Minimum radial coordinate for the region of particles. Measured from the center Total number of particles to create (not per processor or per element). The number is parsed as a floating point number (so that one can specify, for example, '1e4' particles) but it is interpreted as an integer, of course. -292 +295 [Double 0...MAX_DOUBLE (inclusive)] @@ -21919,7 +21923,7 @@ Total number of particles to create (not per processor or per element). The numb The number of radial shells of particles that will be generated around the central point. -302 +305 [Integer range 1...2147483647 (inclusive)] @@ -21940,7 +21944,7 @@ true Whether to correctly evaluate old and current velocity solution to reach higher-order accuracy in time. If set to 'false' only the old velocity solution is evaluated to simulate a first order method in time. This is only recommended for benchmark purposes. -313 +316 [Bool] @@ -21961,7 +21965,7 @@ false Extends the range used by 'Use linear least squares limiter' by linearly interpolating values at cell boundaries from neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. Enabling 'Use boundary extrapolation' requires enabling 'Use linear least squares limiter'. -316 +321 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -21978,7 +21982,7 @@ true Limit the interpolation of particle properties onto the cell, so that the value of each property is no smaller than its minimum and no larger than its maximum on the particles of each cell, and the average of neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. -315 +320 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -21997,7 +22001,7 @@ false Extends the range used by 'Use quadratic least squares limiter' by linearly interpolating values at cell boundaries from neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. Enabling 'Use boundary extrapolation' requires enabling 'Use quadratic least squares limiter'. -321 +319 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -22014,7 +22018,7 @@ true Limit the interpolation of particle properties onto the cell, so that the value of each property is no smaller than its minimum and no larger than its maximum on the particles of each cell, and the average of neighboring cells. If more than one value is given, it will be treated as a list with one component per particle property. -320 +318 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -22034,7 +22038,7 @@ Limit the interpolation of particle properties onto the cell, so that the value The minimum porosity that has to be present at the position of a particle for it to be considered a melt particle (in the sense that the melt presence property is set to 1). -323 +330 [Double 0...1 (inclusive)] @@ -22209,7 +22213,7 @@ $ASPECT_SOURCE_DIR/data/postprocess/boundary-strain-rate-residual/ The name of a directory that contains the ascii data. This path may either be absolute (if starting with a `/') or relative to the current directory. The path may also include the special text `$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -396 +397 [DirectoryName] @@ -22226,7 +22230,7 @@ box_3d_boundary_strain_rate.txt The file name of the input surface strain rate an ascii data. The file has one column in addition to the coordinate columns corresponding to the second invariant of strain rate. -397 +398 [Anything] @@ -22243,7 +22247,7 @@ The file name of the input surface strain rate an ascii data. The file has one c Scalar factor, which is applied to the model data. You might want to use this to scale the input to a reference model. -398 +399 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -22262,7 +22266,7 @@ $ASPECT_SOURCE_DIR/data/boundary-velocity/gplates/ The name of a directory that contains the GPlates model or the ascii data. This path may either be absolute (if starting with a `/') or relative to the current directory. The path may also include the special text `$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -399 +400 [DirectoryName] @@ -22279,7 +22283,7 @@ current_day.gpml The file name of the input velocity as a GPlates model or an ascii data. For the GPlates model, provide file in the same format as described in the 'gplates' boundary velocity plugin. For the ascii data, provide file in the same format as described in 'ascii data' initial composition plugin. -400 +401 [Anything] @@ -22296,7 +22300,7 @@ The file name of the input velocity as a GPlates model or an ascii data. For the Scalar factor, which is applied to the model data. You might want to use this to scale the input to a reference model. Another way to use this factor is to convert units of the input files. For instance, if you provide velocities in cm/year set this factor to 0.01. -401 +402 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -22313,7 +22317,7 @@ false Use ascii data files (e.g., GPS) for computing residual velocities instead of GPlates data. -403 +404 [Bool] @@ -22330,7 +22334,7 @@ false Specify velocity as r, phi, and theta components instead of x, y, and z. Positive velocities point up, east, and north (in 3d) or out and clockwise (in 2d). This setting only makes sense for spherical geometries.GPlates data is always interpreted to be in east and north directions and is not affected by this parameter. -402 +403 [Bool] @@ -22345,7 +22349,7 @@ Specify velocity as r, phi, and theta components instead of x, y, and z. Positiv Command to execute. -406 +407 [Anything] @@ -22362,7 +22366,7 @@ false Whether to run command from all processes (true), or only on process 0 (false). -405 +406 [Bool] @@ -22379,7 +22383,7 @@ false Select whether \aspect{} should terminate if the command returns a non-zero exit status. -404 +405 [Bool] @@ -22394,7 +22398,7 @@ Select whether \aspect{} should terminate if the command returns a non-zero exit A list of names for each of the compositional fields that you want to compute the combined RMS velocity for. -407 +408 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -22413,7 +22417,7 @@ true Whether to compress the raw and weighted cpo data output files with zlib. -354 +374 [Bool] @@ -22430,7 +22434,7 @@ Whether to compress the raw and weighted cpo data output files with zlib. The seed used to generate random numbers. This will make sure that results are reproducible as long as the problem is run with the same amount of MPI processes. It is implemented as final seed = random number seed + MPI Rank. -349 +369 [Integer range 0...2147483647 (inclusive)] @@ -22443,7 +22447,7 @@ The seed used to generate random numbers. This will make sure that results are r On large clusters it can be advantageous to first write the output to a temporary file on a local file system and later move this file to a network file system. If this variable is set to a non-empty string it will be interpreted as a temporary storage location. -351 +371 [Anything] @@ -22462,7 +22466,7 @@ The time interval between each generation of output files. A value of zero indic Units: years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise. -348 +368 [Double 0...MAX_DOUBLE (inclusive)] @@ -22479,7 +22483,7 @@ false File operations can potentially take a long time, blocking the progress of the rest of the model run. Setting this variable to `true' moves this process into background threads, while the rest of the model continues. -350 +370 [Bool] @@ -22497,7 +22501,7 @@ A list containing the what part of the random draw volume weighted particle cpo Note that the rotation matrix and Euler angles both contain the same information, but in a different format. Euler angles are recommended over the rotation matrix since they only require to write 3 values instead of 9. If the list is empty, this file will not be written. Furthermore, the entries will be written out in the order given, and if entries are entered multiple times, they will be written out multiple times. -353 +373 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -22515,7 +22519,7 @@ A list containing what particle cpo data needs to be written out after the parti Note that the rotation matrix and Euler angles both contain the same information, but in a different format. Euler angles are recommended over the rotation matrix since they only require to write 3 values instead of 9. If the list is empty, this file will not be written.Furthermore, the entries will be written out in the order given, and if entries are entered multiple times, they will be written out multiple times. -352 +372 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -22530,7 +22534,7 @@ Note that the rotation matrix and Euler angles both contain the same information The depth boundaries of zones within which we are to compute averages. By default this list is empty and we subdivide the entire domain into equidistant depth zones and compute averages within each of these zones. If this list is not empty it has to contain one more entry than the 'Number of zones' parameter, representing the upper and lower depth boundary of each zone. It is not necessary to cover the whole depth-range (i.e. you can select to only average in a single layer by choosing 2 arbitrary depths as the boundaries of that layer). -357 +377 [List of <[Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -22552,7 +22556,7 @@ List of options: all|temperature|composition|adiabatic temperature|adiabatic pressure|adiabatic density|adiabatic density derivative|velocity magnitude|sinking velocity|rising velocity|Vs|Vp|log viscosity|viscosity|vertical heat flux|vertical mass flux|composition mass -359 +379 [MultipleSelection all|temperature|composition|adiabatic temperature|adiabatic pressure|adiabatic density|adiabatic density derivative|velocity magnitude|sinking velocity|rising velocity|Vs|Vp|log viscosity|viscosity|vertical heat flux|vertical mass flux|composition mass ] @@ -22569,7 +22573,7 @@ all|temperature|composition|adiabatic temperature|adiabatic pressure|adiabatic d The number of zones in depth direction within which we are to compute averages. By default, we subdivide the entire domain into 10 depth zones and compute temperature and other averages within each of these zones. However, if you have a very coarse mesh, it may not make much sense to subdivide the domain into so many zones and you may wish to choose less than this default. It may also make computations slightly faster. On the other hand, if you have an extremely highly resolved mesh, choosing more zones might also make sense. -356 +376 [Integer range 1...2147483647 (inclusive)] @@ -22586,7 +22590,7 @@ gnuplot, txt A list of formats in which the output shall be produced. The format in which the output is generated also determines the extension of the file into which data is written. The list of possible output formats that can be given here is documented in the appendix of the manual where the current parameter is described. By default the output is written as gnuplot file (for plotting), and as a simple text file. -358 +378 [MultipleSelection none|dx|ucd|gnuplot|povray|eps|gmv|tecplot|vtk|vtu|hdf5|svg|deal.II intermediate|txt ] @@ -22603,7 +22607,7 @@ A list of formats in which the output shall be produced. The format in which the 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. -355 +375 [Double 0...MAX_DOUBLE (inclusive)] @@ -22622,7 +22626,7 @@ false Output the excess entropy only instead the each entropy terms. -408 +367 [Bool] @@ -22641,7 +22645,7 @@ Output the excess entropy only instead the each entropy terms. Dynamic topography is calculated as the excess or lack of mass that is supported by mantle flow. This value depends on the density of material that is moved up or down, i.e. crustal rock, and the density of the material that is displaced (generally water or air). While the density of crustal rock is part of the material model, this parameter `Density above' allows the user to specify the density value of material that is displaced above the solid surface. By default this material is assumed to be air, with a density of 0. Units: \si{\kilogram\per\meter\cubed}. -360 +380 [Double 0...MAX_DOUBLE (inclusive)] @@ -22658,7 +22662,7 @@ Dynamic topography is calculated as the excess or lack of mass that is supported Dynamic topography is calculated as the excess or lack of mass that is supported by mantle flow. This value depends on the density of material that is moved up or down, i.e. mantle above CMB, and the density of the material that is displaced (generally outer core material). While the density of mantle rock is part of the material model, this parameter `Density below' allows the user to specify the density value of material that is displaced below the solid surface. By default this material is assumed to be outer core material with a density of 9900. Units: \si{\kilogram\per\meter\cubed}. -361 +381 [Double 0...MAX_DOUBLE (inclusive)] @@ -22675,7 +22679,7 @@ true Whether to output a file containing the bottom (i.e., CMB) dynamic topography. -363 +383 [Bool] @@ -22692,7 +22696,7 @@ true Whether to output a file containing the surface dynamic topography. -362 +382 [Bool] @@ -22751,7 +22755,7 @@ false The density value above the surface boundary. -369 +389 [Double 0...MAX_DOUBLE (inclusive)] @@ -22768,7 +22772,7 @@ The density value above the surface boundary. The density value below the CMB boundary. -370 +390 [Double 0...MAX_DOUBLE (inclusive)] @@ -22785,7 +22789,7 @@ true Option to include the contribution from CMB topography on geoid. The default is true. -365 +385 [Bool] @@ -22802,7 +22806,7 @@ true Option to include the contribution from surface topography on geoid. The default is true. -364 +384 [Bool] @@ -22819,7 +22823,7 @@ Option to include the contribution from surface topography on geoid. The default This parameter can be a random positive integer. However, the value normally should not exceed the maximum degree of the initial perturbed temperature field. For example, if the initial temperature uses S40RTS, the maximum degree should not be larger than 40. -366 +386 [Integer range 0...2147483647 (inclusive)] @@ -22836,7 +22840,7 @@ This parameter can be a random positive integer. However, the value normally sho This parameter normally is set to 2 since the perturbed gravitational potential at degree 1 always vanishes in a reference frame with the planetary center of mass same as the center of figure. -367 +387 [Integer range 0...2147483647 (inclusive)] @@ -22853,7 +22857,7 @@ false Option to output the spherical harmonic coefficients of the CMB topography contribution to the maximum degree. The default is false. -373 +393 [Bool] @@ -22870,7 +22874,7 @@ false Option to output the geoid anomaly in geographical coordinates (latitude and longitude). The default is false, so the postprocessor will output the data in geocentric coordinates (x,y,z) as normally. -368 +388 [Bool] @@ -22887,7 +22891,7 @@ false Option to output the spherical harmonic coefficients of the density anomaly contribution to the maximum degree. The default is false. -374 +394 [Bool] @@ -22904,7 +22908,7 @@ false Option to output the spherical harmonic coefficients of the geoid anomaly up to the maximum degree. The default is false, so the postprocessor will only output the geoid anomaly in grid format. -371 +391 [Bool] @@ -22921,7 +22925,7 @@ false Option to output the free-air gravity anomaly up to the maximum degree. The unit of the output is in SI, hence $m/s^2$ ($1mgal = 10^-5 m/s^2$). The default is false. -375 +395 [Bool] @@ -22938,7 +22942,7 @@ false Option to output the spherical harmonic coefficients of the surface topography contribution to the maximum degree. The default is false. -372 +392 [Bool] @@ -22957,7 +22961,7 @@ false Whether to put every nonlinear iteration into a separate line in the statistics file (if true), or to output only one line per time step that contains the total number of iterations of the Stokes and advection linear system solver. -376 +396 [Bool] @@ -22972,7 +22976,7 @@ Whether to put every nonlinear iteration into a separate line in the statistics Parameter for the list of points sampling scheme: List of satellite latitude coordinates. -393 +364 [List of <[Double -90...90 (inclusive)]> of length 0...4294967295 (inclusive)] @@ -22985,7 +22989,7 @@ Parameter for the list of points sampling scheme: List of satellite latitude coo Parameter for the list of points sampling scheme: List of satellite longitude coordinates. -392 +363 [List of <[Double -180...180 (inclusive)]> of length 0...4294967295 (inclusive)] @@ -22998,7 +23002,7 @@ Parameter for the list of points sampling scheme: List of satellite longitude co Parameter for the list of points sampling scheme: List of satellite radius coordinates. Just specify one radius if all points values have the same radius. If not, make sure there are as many radius as longitude and latitude -391 +362 [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] @@ -23015,7 +23019,7 @@ Parameter for the list of points sampling scheme: List of satellite radius coord Parameter for the uniform distribution sampling scheme: Gravity may be calculated for a sets of points along the latitude between a minimum and maximum latitude. -388 +359 [Double -90...90 (inclusive)] @@ -23032,7 +23036,7 @@ Parameter for the uniform distribution sampling scheme: Gravity may be calculate Parameter for the uniform distribution sampling scheme: Gravity may be calculated for a sets of points along the longitude between a minimum and maximum longitude. -387 +358 [Double -180...180 (inclusive)] @@ -23049,7 +23053,7 @@ Parameter for the uniform distribution sampling scheme: Gravity may be calculate Parameter for the map sampling scheme: Maximum radius can be defined in or outside the model. -384 +355 [Double 0...MAX_DOUBLE (inclusive)] @@ -23066,7 +23070,7 @@ Parameter for the map sampling scheme: Maximum radius can be defined in or outsi Parameter for the uniform distribution sampling scheme: Gravity may be calculated for a sets of points along the latitude between a minimum and maximum latitude. -386 +357 [Double -90...90 (inclusive)] @@ -23083,7 +23087,7 @@ Parameter for the uniform distribution sampling scheme: Gravity may be calculate Parameter for the uniform distribution sampling scheme: Gravity may be calculated for a sets of points along the longitude between a minimum and maximum longitude. -385 +356 [Double -180...180 (inclusive)] @@ -23100,7 +23104,7 @@ Parameter for the uniform distribution sampling scheme: Gravity may be calculate Parameter for the map sampling scheme: Minimum radius may be defined in or outside the model. Prescribe a minimum radius for a sampling coverage at a specific height. -383 +354 [Double 0...MAX_DOUBLE (inclusive)] @@ -23117,7 +23121,7 @@ Parameter for the map sampling scheme: Minimum radius may be defined in or outsi Parameter for the fibonacci spiral sampling scheme: This specifies the desired number of satellites per radius layer. The default value is 200. Note that sampling becomes more uniform with increasing number of satellites -378 +349 [Integer range 0...2147483647 (inclusive)] @@ -23134,7 +23138,7 @@ Parameter for the fibonacci spiral sampling scheme: This specifies the desired n Parameter for the map sampling scheme: This specifies the number of points along the latitude (e.g. gravity map) between a minimum and maximum latitude. -382 +353 [Integer range 0...2147483647 (inclusive)] @@ -23151,7 +23155,7 @@ Parameter for the map sampling scheme: This specifies the number of points along Parameter for the map sampling scheme: This specifies the number of points along the longitude (e.g. gravity map) between a minimum and maximum longitude. -381 +352 [Integer range 0...2147483647 (inclusive)] @@ -23168,7 +23172,7 @@ Parameter for the map sampling scheme: This specifies the number of points along Parameter for the map sampling scheme: This specifies the number of points along the radius (e.g. depth profile) between a minimum and maximum radius. -380 +351 [Integer range 0...2147483647 (inclusive)] @@ -23185,7 +23189,7 @@ Parameter for the map sampling scheme: This specifies the number of points along Set the precision of gravity acceleration, potential and gradients in the gravity output and statistics file. -390 +361 [Integer range 1...2147483647 (inclusive)] @@ -23202,7 +23206,7 @@ Set the precision of gravity acceleration, potential and gradients in the gravit Quadrature degree increase over the velocity element degree may be required when gravity is calculated near the surface or inside the model. An increase in the quadrature element adds accuracy to the gravity solution from noise due to the model grid. -379 +350 [Integer range -1...2147483647 (inclusive)] @@ -23219,7 +23223,7 @@ Quadrature degree increase over the velocity element degree may be required when Gravity anomalies may be computed using density anomalies relative to a reference density. -389 +360 [Double 0...MAX_DOUBLE (inclusive)] @@ -23236,7 +23240,7 @@ map Choose the sampling scheme. By default, the map produces a grid of equally angled points between a minimum and maximum radius, longitude, and latitude. A list of points contains the specific coordinates of the satellites. The fibonacci spiral sampling scheme produces a uniformly distributed map on the surface of sphere defined by a minimum and/or maximum radius. -377 +348 [Selection map|list|list of points|fibonacci spiral ] @@ -23253,7 +23257,7 @@ Choose the sampling scheme. By default, the map produces a grid of equally angle The time interval between each generation of gravity output files. A value of 0 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. -394 +365 [Double 0...MAX_DOUBLE (inclusive)] @@ -23270,7 +23274,7 @@ The time interval between each generation of gravity output files. A value of 0 The maximum number of time steps between each generation of gravity output files. -395 +366 [Integer range 0...2147483647 (inclusive)] @@ -23289,7 +23293,7 @@ true If set to 'true', also output the peak virtual memory usage (computed as the maximum over all processors). -198 +201 [Bool] @@ -23308,7 +23312,7 @@ vtu A comma separated list of file formats to be used for graphical output. The list of possible output formats that can be given here is documented in the appendix of the manual where the current parameter is described. -200 +203 [MultipleSelection none|dx|ucd|gnuplot|povray|eps|gmv|tecplot|vtk|vtu|hdf5|svg|deal.II intermediate|ascii ] @@ -23321,7 +23325,7 @@ A comma separated list of file formats to be used for graphical output. The list A comma separated list of particle properties that should \textit{not} be output. If this list contains the entry `all', only the id of particles will be provided in graphical output files. -204 +207 [Anything] @@ -23338,7 +23342,7 @@ A comma separated list of particle properties that should \textit{not} be output VTU file output supports grouping files from several CPUs into a given number of files 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. A value of 1 will generate one big file containing the whole solution, while a larger value will create that many files (at most as many as there are MPI ranks). -201 +204 [Integer range 0...2147483647 (inclusive)] @@ -23351,7 +23355,7 @@ VTU file output supports grouping files from several CPUs into a given number of On large clusters it can be advantageous to first write the output to a temporary file on a local file system and later move this file to a network file system. If this variable is set to a non-empty string it will be interpreted as a temporary storage location. -203 +206 [Anything] @@ -23370,7 +23374,7 @@ The time interval between each generation of output files. A value of zero indic Units: years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise. -199 +202 [Double 0...MAX_DOUBLE (inclusive)] @@ -23387,7 +23391,7 @@ false File operations can potentially take a long time, blocking the progress of the rest of the model run. Setting this variable to `true' moves this process into a background thread, while the rest of the model continues. -202 +205 [Bool] @@ -23402,7 +23406,7 @@ File operations can potentially take a long time, blocking the progress of the r The list of points at which the solution should be evaluated. Points need to be separated by semicolons, and coordinates of each point need to be separated by commas. -346 +187 [List of <[List of <[Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)]> of length 2...2 (inclusive)]> of length 0...4294967295 (inclusive) separated by <;>] @@ -23419,7 +23423,7 @@ The list of points at which the solution should be evaluated. Points need to be The time interval between each generation of point values output. 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. -345 +186 [Double 0...MAX_DOUBLE (inclusive)] @@ -23436,7 +23440,7 @@ false Whether or not the Evaluation points are specified in the natural coordinates of the geometry model, e.g. radius, lon, lat for the chunk model. Currently, natural coordinates for the spherical shell and sphere geometries are not supported. -347 +188 [Bool] @@ -23455,7 +23459,7 @@ false Whether to write the full moment of inertia tensor into the statistics output instead of its norm for the current rotation axis. This is a second-order symmetric tensor with 6 components in 3d. In 2d this option has no effect, because the rotation axis is fixed and thus the moment of inertia is always a scalar. -187 +190 [Bool] @@ -23472,7 +23476,7 @@ false Whether to use a constant density of one for the computation of the angular momentum and moment of inertia. This is an approximation that assumes that the 'volumetric' rotation is equal to the 'mass' rotation. If this parameter is true this postprocessor computes 'net rotation' instead of 'angular momentum'. -186 +189 [Bool] @@ -23491,7 +23495,7 @@ $ASPECT_SOURCE_DIR/data/postprocess/sea-level/ The name of a directory that contains the ice height [m] ascii data. This path may either be absolute (if starting with a `/') or relative to the current directory. The path may also include the special text `$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -192 +195 [DirectoryName] @@ -23508,7 +23512,7 @@ $ASPECT_SOURCE_DIR/data/postprocess/sea-level/ The name of a directory that contains the topography ascii data. This path may either be absolute (if starting with a `/') or relative to the current directory. The path may also include the special text `$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -190 +193 [DirectoryName] @@ -23525,7 +23529,7 @@ shell_3d_ice_top.0.txt The file name of the ice height ascii data. For the ascii data, provide file in the same format as described in 'ascii data' initial composition plugin. -193 +196 [Anything] @@ -23542,7 +23546,7 @@ shell_3d_topo_top.0.txt The file name of the topography ascii data. For the ascii data, provide file in the same format as described in 'ascii data' initial composition plugin. -191 +194 [Anything] @@ -23559,7 +23563,7 @@ The file name of the topography ascii data. For the ascii data, provide file in The density of ice [kg/m3] -188 +191 [Double 0...MAX_DOUBLE (inclusive)] @@ -23576,7 +23580,7 @@ false Whether or not to write sea level to a text file named named 'sea_level.NNNNN' in the output directory -194 +197 [List of <[Bool]> of length 0...4294967295 (inclusive)] @@ -23593,7 +23597,7 @@ Whether or not to write sea level to a text file named named 'sea_level.NNN The time interval between each generation of text 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. -195 +198 [Double 0...MAX_DOUBLE (inclusive)] @@ -23610,7 +23614,7 @@ The time interval between each generation of text output files. A value of zero The density of water [kg/m3] -189 +192 [Double 0...MAX_DOUBLE (inclusive)] @@ -23629,7 +23633,7 @@ false Whether or not to write topography to a text file named named 'topography.NNNNN' in the output directory -196 +199 [Bool] @@ -23646,7 +23650,7 @@ Whether or not to write topography to a text file named named 'topography.N The time interval between each generation of text 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. -197 +200 [Double 0...MAX_DOUBLE (inclusive)] @@ -24190,7 +24194,7 @@ false A boolean flag that controls whether to output the heat flux map as a point wise value, or as a cell-wise averaged value. The point wise output is more accurate, but it currently omits prescribed heat flux values at boundaries and advective heat flux that is caused by velocities non-tangential to boundaries. If you do not use these two features it is recommended to switch this setting on to benefit from the increased output resolution. -162 +180 [Bool] @@ -24211,7 +24215,7 @@ A comma separated list of material properties that should be written whenever wr viscosity|density|thermal expansivity|specific heat|thermal conductivity|thermal diffusivity|compressibility|entropy derivative temperature|entropy derivative pressure|reaction terms|melt fraction -163 +181 [MultipleSelection viscosity|density|thermal expansivity|specific heat|thermal conductivity|thermal diffusivity|compressibility|entropy derivative temperature|entropy derivative pressure|reaction terms|melt fraction ] @@ -24230,7 +24234,7 @@ viscosity|density|thermal expansivity|specific heat|thermal conductivity|thermal Constant parameter in the quadratic function that approximates the solidus of peridotite. Units: \si{\degreeCelsius}. -165 +157 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24247,7 +24251,7 @@ Constant parameter in the quadratic function that approximates the solidus of pe Prefactor of the linear pressure term in the quadratic function that approximates the solidus of peridotite. \si{\degreeCelsius\per\pascal}. -166 +158 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24264,7 +24268,7 @@ Prefactor of the linear pressure term in the quadratic function that approximate Prefactor of the quadratic pressure term in the quadratic function that approximates the solidus of peridotite. \si{\degreeCelsius\per\pascal\squared}. -167 +159 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24281,7 +24285,7 @@ Prefactor of the quadratic pressure term in the quadratic function that approxim Constant parameter in the quadratic function that approximates the lherzolite liquidus used for calculating the fraction of peridotite-derived melt. Units: \si{\degreeCelsius}. -168 +160 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24298,7 +24302,7 @@ Constant parameter in the quadratic function that approximates the lherzolite li Prefactor of the linear pressure term in the quadratic function that approximates the lherzolite liquidus used for calculating the fraction of peridotite-derived melt. \si{\degreeCelsius\per\pascal}. -169 +161 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24315,7 +24319,7 @@ Prefactor of the linear pressure term in the quadratic function that approximate Prefactor of the quadratic pressure term in the quadratic function that approximates the lherzolite liquidus used for calculating the fraction of peridotite-derived melt. \si{\degreeCelsius\per\pascal\squared}. -170 +162 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24332,7 +24336,7 @@ Prefactor of the quadratic pressure term in the quadratic function that approxim Constant parameter in the quadratic function that approximates the liquidus of peridotite. Units: \si{\degreeCelsius}. -171 +163 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24349,7 +24353,7 @@ Constant parameter in the quadratic function that approximates the liquidus of p Prefactor of the linear pressure term in the quadratic function that approximates the liquidus of peridotite. \si{\degreeCelsius\per\pascal}. -172 +164 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24366,7 +24370,7 @@ Prefactor of the linear pressure term in the quadratic function that approximate Prefactor of the quadratic pressure term in the quadratic function that approximates the liquidus of peridotite. \si{\degreeCelsius\per\pascal\squared}. -173 +165 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24383,7 +24387,7 @@ Prefactor of the quadratic pressure term in the quadratic function that approxim Constant parameter in the quadratic function that approximates the solidus of pyroxenite. Units: \si{\degreeCelsius}. -178 +170 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24400,7 +24404,7 @@ Constant parameter in the quadratic function that approximates the solidus of py Prefactor of the linear pressure term in the quadratic function that approximates the solidus of pyroxenite. Note that this factor is different from the value given in Sobolev, 2011, because they use the potential temperature whereas we use the absolute temperature. \si{\degreeCelsius\per\pascal}. -179 +171 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24417,7 +24421,7 @@ Prefactor of the linear pressure term in the quadratic function that approximate Prefactor of the quadratic pressure term in the quadratic function that approximates the solidus of pyroxenite. \si{\degreeCelsius\per\pascal\squared}. -180 +172 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24434,7 +24438,7 @@ Prefactor of the quadratic pressure term in the quadratic function that approxim Prefactor of the linear depletion term in the quadratic function that approximates the melt fraction of pyroxenite. \si{\degreeCelsius\per\pascal}. -181 +173 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24451,7 +24455,7 @@ Prefactor of the linear depletion term in the quadratic function that approximat Prefactor of the quadratic depletion term in the quadratic function that approximates the melt fraction of pyroxenite. \si{\degreeCelsius\per\pascal\squared}. -182 +174 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24468,7 +24472,7 @@ Prefactor of the quadratic depletion term in the quadratic function that approxi Mass fraction of clinopyroxene in the peridotite to be molten. Units: non-dimensional. -177 +169 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24485,7 +24489,7 @@ Mass fraction of clinopyroxene in the peridotite to be molten. Units: non-dimens Exponent of the melting temperature in the melt fraction calculation. Units: non-dimensional. -176 +168 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24502,7 +24506,7 @@ Exponent of the melting temperature in the melt fraction calculation. Units: non Constant in the linear function that approximates the clinopyroxene reaction coefficient. Units: non-dimensional. -174 +166 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24519,7 +24523,7 @@ Constant in the linear function that approximates the clinopyroxene reaction coe Prefactor of the linear pressure term in the linear function that approximates the clinopyroxene reaction coefficient. Units: \si{\per\pascal}. -175 +167 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24540,7 +24544,7 @@ A comma separated list of melt properties that should be written whenever writin compaction viscosity|fluid viscosity|permeability|fluid density|fluid density gradient|is melt cell|darcy coefficient|darcy coefficient no cutoff|compaction length -164 +182 [MultipleSelection compaction viscosity|fluid viscosity|permeability|fluid density|fluid density gradient|is melt cell|darcy coefficient|darcy coefficient no cutoff|compaction length ] @@ -24559,7 +24563,7 @@ false Whether to use the deviatoric stress tensor instead of the full stress tensor to compute principal stress directions and values. -157 +175 [Bool] @@ -24667,7 +24671,7 @@ reference profile Scheme to compute the average velocity-depth profile. The reference profile option evaluates the conditions along the reference adiabat according to the material model. The lateral average option instead calculates a lateral average from subdivision of the mesh. The lateral average option may produce spurious results where there are sharp velocity changes. -160 +178 [Selection reference profile|lateral average ] @@ -24684,7 +24688,7 @@ Scheme to compute the average velocity-depth profile. The reference profile opti Number of depth slices used to define average seismic compressional wave velocities from which anomalies are calculated. Units: non-dimensional. -161 +179 [Integer range 1...2147483647 (inclusive)] @@ -24703,7 +24707,7 @@ reference profile Scheme to compute the average velocity-depth profile. The reference profile option evaluates the conditions along the reference adiabat according to the material model. The lateral average option instead calculates a lateral average from subdivision of the mesh. The lateral average option may produce spurious results where there are sharp velocity changes. -158 +176 [Selection reference profile|lateral average ] @@ -24720,7 +24724,7 @@ Scheme to compute the average velocity-depth profile. The reference profile opti Number of depth slices used to define average seismic shear wave velocities from which anomalies are calculated. Units: non-dimensional. -159 +177 [Integer range 1...2147483647 (inclusive)] @@ -24765,7 +24769,7 @@ $ASPECT_SOURCE_DIR/data/prescribed-stokes-solution/ The name of a directory that contains the model data. This path may either be absolute (if starting with a `/') or relative to the current directory. The path may also include the special text `$ASPECT_SOURCE_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the `data/' subdirectory of ASPECT. -1325 +1310 [DirectoryName] @@ -24782,7 +24786,7 @@ box_2d.txt The file name of the model data. -1326 +1311 [Anything] @@ -24799,7 +24803,7 @@ The file name of the model data. Point that determines the plane in which the 2d slice lies in. This variable is only used if 'Slice dataset in 2d plane' is true. The slice will go through this point, the point defined by the parameter 'Second point on slice', and the center of the model domain. After the rotation, this first point will lie along the (0,1,0) axis of the coordinate system. The coordinates of the point have to be given in Cartesian coordinates. -1329 +1314 [Anything] @@ -24816,7 +24820,7 @@ Point that determines the plane in which the 2d slice lies in. This variable is Scalar factor, which is applied to the model data. You might want to use this to scale the input to a reference model. Another way to use this factor is to convert units of the input files. For instance, if you provide velocities in cm/yr set this factor to 0.01. -1327 +1312 [Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)] @@ -24833,7 +24837,7 @@ Scalar factor, which is applied to the model data. You might want to use this to Second point that determines the plane in which the 2d slice lies in. This variable is only used if 'Slice dataset in 2d plane' is true. The slice will go through this point, the point defined by the parameter 'First point on slice', and the center of the model domain. The coordinates of the point have to be given in Cartesian coordinates. -1330 +1315 [Anything] @@ -24850,7 +24854,7 @@ false Whether to use a 2d data slice of a 3d data file or the entire data file. Slicing a 3d dataset is only supported for 2d models. -1328 +1313 [Bool] @@ -24867,7 +24871,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -1321 +1327 [Anything] @@ -24886,7 +24890,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -1320 +1326 [Anything] @@ -24903,7 +24907,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -1319 +1325 [Anything] @@ -24920,7 +24924,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -1318 +1324 [Anything] @@ -24939,7 +24943,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -1317 +1323 [Anything] @@ -24956,7 +24960,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -1316 +1322 [Anything] @@ -24973,7 +24977,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -1324 +1330 [Anything] @@ -24992,7 +24996,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -1323 +1329 [Anything] @@ -25009,7 +25013,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -1322 +1328 [Anything] @@ -25026,7 +25030,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -1315 +1321 [Anything] @@ -25045,7 +25049,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -1314 +1320 [Anything] @@ -25062,7 +25066,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -1313 +1319 [Anything] @@ -25079,7 +25083,7 @@ Sometimes it is convenient to use symbolic constants in the expression that desc A typical example would be to set this runtime parameter to `pi=3.1415926536' and then use `pi' in the expression of the actual formula. (That said, for convenience this class actually defines both `pi' and `Pi' by default, but you get the idea.) -1312 +1318 [Anything] @@ -25098,7 +25102,7 @@ The formula that denotes the function you want to evaluate for particular values If the function you are describing represents a vector-valued function with multiple components, then separate the expressions for individual components by a semicolon. -1311 +1317 [Anything] @@ -25115,7 +25119,7 @@ x,y,t The names of the variables as they will be used in the function, separated by commas. By default, the names of variables at which the function will be evaluated are `x' (in 1d), `x,y' (in 2d) or `x,y,z' (in 3d) for spatial coordinates and `t' for time. You can then use these variable names in your function expression and they will be replaced by the values of these variables at which the function is currently evaluated. However, you can also choose a different set of names for the independent variables at which to evaluate your function expression. For example, if you work in spherical coordinates, you may wish to set this input parameter to `r,phi,theta,t' and then use these variable names in your function expression. -1310 +1316 [Anything] @@ -25829,7 +25833,7 @@ Whether to checkpoint the simulation right before termination. Terminate the simulation once the specified timestep has been reached. -435 +437 [Integer range 0...2147483647 (inclusive)] @@ -25879,7 +25883,7 @@ The criterion considers the total heat flux over all boundaries listed by their The wall time of the simulation. Unit: hours. -437 +439 [Double 0...MAX_DOUBLE (inclusive)] @@ -25895,7 +25899,7 @@ A comma separated list of names denoting those boundaries that should be taken i The names of the boundaries listed here can either be numbers (in which case they correspond to the numerical boundary indicators assigned by the geometry object), or they can correspond to any of the symbolic names the geometry object may have provided for each part of the boundary. You may want to compare this with the documentation of the geometry model you use in your model. -440 +442 [List of <[Anything]> of length 0...4294967295 (inclusive)] @@ -25912,7 +25916,7 @@ The names of the boundaries listed here can either be numbers (in which case the The maximum relative deviation of the heat flux in recent simulation time for the system to be considered in steady state. If the actual deviation is smaller than this number, then the simulation will be terminated. -438 +440 [Double 0...MAX_DOUBLE (inclusive)] @@ -25929,7 +25933,7 @@ The maximum relative deviation of the heat flux in recent simulation time for th The minimum length of simulation time that the system should be in steady state before termination. Note that if the time step size is similar to or larger than this value, the termination criterion will only have very few (in the most extreme case, just two) heat flux values to check. To ensure that a larger number of time steps are included in the check for steady state, this value should be much larger than the time step size. Units: years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise. -439 +441 [Double 0...MAX_DOUBLE (inclusive)] @@ -25948,7 +25952,7 @@ The minimum length of simulation time that the system should be in steady state The maximum relative deviation of the temperature in recent simulation time for the system to be considered in steady state. If the actual deviation is smaller than this number, then the simulation will be terminated. -443 +434 [Double 0...MAX_DOUBLE (inclusive)] @@ -25965,7 +25969,7 @@ The maximum relative deviation of the temperature in recent simulation time for The minimum length of simulation time that the system should be in steady state before termination.Units: years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise. -444 +435 [Double 0...MAX_DOUBLE (inclusive)] @@ -25984,7 +25988,7 @@ The minimum length of simulation time that the system should be in steady state The maximum relative deviation of the RMS in recent simulation time for the system to be considered in steady state. If the actual deviation is smaller than this number, then the simulation will be terminated. -441 +443 [Double 0...MAX_DOUBLE (inclusive)] @@ -26001,7 +26005,7 @@ The maximum relative deviation of the RMS in recent simulation time for the syst The minimum length of simulation time that the system should be in steady state before termination.Units: years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise. -442 +444 [Double 0...MAX_DOUBLE (inclusive)] @@ -26020,7 +26024,7 @@ terminate-aspect The name of a file that, if it exists in the output directory (whose name is also specified in the input file) will lead to termination of the simulation. The file's location is chosen to be in the output directory, rather than in a generic location such as the ASPECT directory, so that one can run multiple simulations at the same time (which presumably write to different output directories) and can selectively terminate a particular one. -434 +436 [FileName (Type: input)] diff --git a/doc/sphinx/parameters/Material_20model.md b/doc/sphinx/parameters/Material_20model.md index dde22d0e446..252df554cc5 100644 --- a/doc/sphinx/parameters/Material_20model.md +++ b/doc/sphinx/parameters/Material_20model.md @@ -1776,11 +1776,11 @@ Units: \si{\pascal\second} (parameters:Material_20model/Latent_20heat/Viscosity_20prefactors)= ### __Parameter name:__ Viscosity prefactors -**Default value:** +**Default value:** all:1 -**Pattern:** [List of <[Double 0...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] +**Pattern:** [Anything] -**Documentation:** A list of prefactors for the viscosity for each phase. The reference viscosity will be multiplied by this factor to get the corresponding viscosity for each phase. List must have one more entry than Phase transition depths. Units: non-dimensional. +**Documentation:** A list of prefactors for the viscosity for each phase. The “Viscosity” parameter (modified by the “Composition viscosity prefactor”, depending on composition) will be multiplied by this factor to get the corresponding viscosity for each phase. List must have the same number of entries as there are phases, that is one more than “Phase transition depths” for each composition that is used in the model. Units: non-dimensional. (parameters:Material_20model/Latent_20heat_20melt)= ## **Subsection:** Material model / Latent heat melt diff --git a/doc/sphinx/parameters/Mesh_20refinement.md b/doc/sphinx/parameters/Mesh_20refinement.md index d8bb20912c9..fd4934886da 100644 --- a/doc/sphinx/parameters/Mesh_20refinement.md +++ b/doc/sphinx/parameters/Mesh_20refinement.md @@ -156,7 +156,7 @@ This refinement criterion computes the gradient of the compositional field at qu On the other hand, for discontinuous finite elements (see the ‘Use discontinuous composition discretization’ parameter in the ‘Discretization’ section), the gradient at quadrature points does not include the contribution of jumps in the compositional field between cells, and consequently will not be an accurate approximation of the true gradient. As an extreme example, consider the case of using piecewise constant finite elements for compositional fields; in that case, the gradient of the solution at quadrature points inside each cell will always be exactly zero, even if the finite element solution is different from each cell to the next. Consequently, the current refinement criterion will likely not be useful in this situation. That said, the ‘composition approximate gradient’ refinement criterion exists for exactly this purpose. -‘composition threshold’: A mesh refinement criterion that computes refinement indicators from the compositional fields. If any field exceeds the threshold given in the input file, the cell is marked for refinement. +‘composition threshold’: A mesh refinement criterion that computes refinement indicators from the compositional fields. One threshold per compositional is given in the input file, and if any field exceeds its threshold, the cell is marked for refinement. ‘density’: A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the density, $\rho$. Because this quantity may not be a continuous function ($\rho$ and $C_p$ may be discontinuous functions along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1+d/2}$ where $h_K$ is the diameter of each cell and $d$ is the dimension. This scaling ensures that the error indicators converge to zero as $h_K\rightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$. @@ -306,7 +306,7 @@ If the list of scaling factors given in this parameter is empty, then this indic **Pattern:** [List of <[Double -MAX_DOUBLE...MAX_DOUBLE (inclusive)]> of length 0...4294967295 (inclusive)] -**Documentation:** A list of thresholds that every individual compositional field will be evaluated against. +**Documentation:** A list of thresholds, one for each compositional field to be evaluated against. (parameters:Mesh_20refinement/Isosurfaces)= ## **Subsection:** Mesh refinement / Isosurfaces