ASPECT
Changes after the latest release (v2.0.0)

This is the list of changes made after the release of Aspect version 2.0.0. All entries are signed with the names of the author.

  1. Fixed: The global entropy variation that is used in the artificial viscosity computation was always computed with a quadrature of the polynomial degree of the temperature (even for compositional fields). This is fixed now, and causes minor changes in the stabilization for compositional fields if they use a different element order than the temperature (usually not the case).
    (Rene Gassmoeller, 2018/08/31)

  2. New: There is now a 'material statistics' postprocessor that computes the average density, average viscosity, and total mass in a model.
    (Rene Gassmoeller, 2018/08/22)

  3. New: ASPECT has two visualization postprocessors which calculate and output the grain lag angle and the infinite strain axis (ISA) rotation timescale, respectively. These two quantities can be used to calculate the grain orientation lag parameter of Kaminski and Ribe (G3, 2002)
    (Hannah Mark, 2018/08/12)

  4. Fixed: Fixed a bug in the 'Named additional outputs' postprocessor, which would cause a segmentation fault if more than one named additional material output object was active in the same model.
    (Rene Gassmoeller, 2018/08/06)

  5. Fixed: Fixed a bug in the 'Steinberger' material model that caused ASPECT to use wrong compressibilities when only a single compositional field, and a single material data file were prescribed. The correct behavior would have been to ignore the compositional field and use the data file as background material, but instead the compositional field was interpreted as volume fraction and multiplied with the compressibility calculated from the data file. This is fixed now.
    (Rene Gassmoeller, 2018/07/29)

  6. New: ASPECT now outputs a dynamically generated URL based on used features to ask people to cite appropriate papers.
    (Timo Heister, 2018/07/25)

  7. New: Added a postprocessor for the temperature anomaly, defined as temperature - (depth average of temperature).
    (Max Rudolph, 2018/07/11)

  8. New: ASPECT now has a Patch on S40RTS initial temperature model which combines S40RTS and an ascii grid containing an upper mantle tomography. Data is input in Vs and converted to temperature within the function (using the same method as S40RTS_perturbation).
    (Sophie Coulson, 2018/07/09)

  9. New: Analytical benchmark for 2D Cartesian incompressible Stokes flow following the FEM book by Donea and Huerta.
    (Cedric Thieulot, 2018/07/02)

  10. New: The benchmarks now include a viscoplastic thermal convection benchmark as described and performed in Tosi et al. (G3 2015).
    (Anne Glerum, 2018/06/27)

  11. Impoved: The assembly speed for the Newton solver has been significantly improved by caching various variables.
    (Menno Fraters and Rene Gassmoeller, 2018/06/26)

  12. Fixed: The default latitude and longitude limits for the 'uniform radial' particle generator were set to (0, pi) even though they are measured in degrees. The new default limits correspond to the entire surface of the sphere, and the manual clarifies that latitude is in (0, 180) and not (-90, 90).
    (Bart Niday, 2018/06/25)

  13. New: ASPECT can now call PerpleX to calculate material properties, phase amounts and compositions on-the-fly. A simple material model (PerpleX lookup) is provided which uses PerpleX to provide densities, compressibilities, thermal expansivities and heat capacities on a cell-by-cell basis. This model is provided as a proof-of-concept; more efficient procedures will be required for production runs.
    (Bob Myhill, 2018/06/25)

  14. New: We now have unit_tests in the unit_test/ folder using the Catch testing framework. They should be used instead of writing a test in tests/ if they don't require an ASPECT run.
    (Timo Heister, 2018/06/24)

  15. New: The ascii data and function boundary velocity plugins now allow velocities to be specified along spherical (up, north, east) unit vectors.
    (Bart Niday, 2018/06/24)
    New: Added a visualization plugin that directly outputs the strain rate tensor.
    (Bart Niday, 2018/06/24)

  16. New: ASPECT can now read in a depth-dependent initial temperature from ASCII data. This is intended to provide a background temperature profile on which to add perturbations for composite initial temperature models.
    (Marie Kajan, 2018/06/24)

  17. New: ASPECT can now read in a depth-dependent vs to density conversion file, which can be used for S40RTS and SAVANI.
    (Jacky Austermann, 2018/06/24)

  18. New: Function GeometryModel::Interface::cartesian_to_other_coordinates converts between defined coordinate systems. This might break some plugins, which can be fixed by including aspect/geometry_model/interface.h
    (Bart Niday, 2018/06/23)

  19. Changed: We now only set time step number and timestep sizes after initialization is finished.

    Right now, there is no way to find out, for example in the hook that is called from the function that builds constraints, whether we are already in the process of time stepping, or only in the initialization phase. So move the initialization of the time step number and time step variables to the very end of the initialization phase. Before this, these variables have invalid values.
    (Wolfgang Bangerth, 2018/06/20)

  20. New: There is now a visualization postprocessor that outputs the compaction length, the characteristic length scale of melt transport.
    (Joe Schools, 2018/06/19)

  21. Fixed: When gnuplot was used as the output format for the 'depth average' postprocessor plugin, the first two columns were labeled 'x' 'y', even though 'depth' and 'time' are always the first two variables. This is now fixed provided deal.II 9.0 or newer is used.
    (Bob Myhill, 2018/06/19)

  22. New: Compositional fields can now be advected with the melt velocity (as opposed to the solid velocity), and this option can be used by selecting the 'melt field' advection method in the input file.
    (Juliane Dannberg, 2018/06/19)

  23. New: ASPECT now does some basic checks that essential parameters (such as the number of compositional fields) are the same before creating a checkpoint and after restarting from it.
    (Wolfgang Bangerth, 2018/06/18)

  24. Fixed: The 'compositional heating' heating plugin had a parameter 'Use compositional field for heat production averaging' that was used inconsistently with its description. Its first entry did not correspond to the background field, but to the first compositional field, and the last value was ignored. This is fixed now, the first entry is used for the background field, and all following values determine whether to include the compositional fields.
    (Rene Gassmoeller, 2018/05/21)

  25. Fixed: The nullspace removal had a minor bug that caused a small remaining rotation in three dimensional models even if the net rotation or angular momentum should have been removed. The bug was fixed, and a new postprocessor called 'rotation statistics' was added that can be used to calculate net rotation and angular momentum.
    (Rene Gassmoeller, Shangxin Liu, 2018/05/03)

  26. New: ASPECT now stores the steps for the continuous integration setup inside the repository in form of a Jenkins pipeline file. This file documents each step in our testing procedure, and allows to easily setup new Jenkins server that reproduce the official ASPECT tester.
    (Rene Gassmoeller, Timo Heister, Tyler Esser, 2018/05/02)