Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
allow ascii data reference profile material model to output seismic v…
…elocities
- Loading branch information
Showing
9 changed files
with
22,114 additions
and
2 deletions.
There are no files selected for viewing
305 changes: 305 additions & 0 deletions
305
data/adiabatic-conditions/ascii-data/example_seismic_velocities.txt
Large diffs are not rendered by default.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,149 @@ | ||
# A cookbook based on the burman_test cookbook, but modified testing | ||
# the output of seismic velocity data generated by burnman and reading | ||
# them in through an ascii data file. | ||
|
||
set Dimension = 2 | ||
set Use years in output instead of seconds = true | ||
set End time = 0 | ||
set Linear solver tolerance = 1e-6 | ||
|
||
subsection Material model | ||
set Model name = ascii reference profile | ||
|
||
subsection Ascii reference profile | ||
set Thermal viscosity exponent = 10.0 | ||
set Viscosity prefactors = 1.0, 0.1, 1.0, 10.0 | ||
|
||
subsection Ascii data model | ||
set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/ | ||
set Data file name = example_seismic_velocities.txt | ||
end | ||
end | ||
|
||
set Material averaging = harmonic average | ||
end | ||
|
||
# The geometry is a spherical shell with the inner and | ||
# outer radius of the Earth's mantle. | ||
subsection Geometry model | ||
set Model name = spherical shell | ||
|
||
subsection Spherical shell | ||
set Inner radius = 3481000 | ||
set Outer radius = 6336000 | ||
end | ||
end | ||
|
||
subsection Adiabatic conditions model | ||
set Model name = ascii data | ||
|
||
subsection Ascii data model | ||
set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/ | ||
set Data file name = example_seismic_velocities.txt | ||
end | ||
end | ||
|
||
|
||
# The gravity model reads an ascii data file that contains a gravity | ||
# profile consistent with the material properties computed using the | ||
# Birch-Murnaghan equation of state. It automatically uses the same file | ||
# specified in the adiabatic conditions model. | ||
subsection Gravity model | ||
set Model name = ascii data | ||
subsection Ascii data model | ||
set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/ | ||
set Data file name = example_seismic_velocities.txt | ||
end | ||
end | ||
|
||
# The model uses the anelastic liquid approximation. | ||
subsection Formulation | ||
set Formulation = anelastic liquid approximation | ||
end | ||
|
||
# The present-day plate velocities are imposed on the upper model | ||
# boundary. | ||
subsection Model settings | ||
set Prescribed velocity boundary indicators = outer:gplates | ||
set Tangential velocity boundary indicators = inner | ||
|
||
set Fixed temperature boundary indicators = inner, outer | ||
end | ||
|
||
|
||
subsection Boundary velocity model | ||
subsection GPlates model | ||
set Data directory = $ASPECT_SOURCE_DIR/data/velocity-boundary-conditions/gplates/ | ||
set Velocity file name = current_day.gpml | ||
set Data file time step = 1e6 | ||
set Point one = 1.5708,4.87 | ||
set Point two = 1.5708,5.24 | ||
end | ||
end | ||
|
||
|
||
subsection Boundary temperature model | ||
set Model name = spherical constant | ||
subsection Spherical constant | ||
set Inner temperature = 3500 | ||
set Outer temperature = 273 | ||
end | ||
end | ||
|
||
# The initial temperature is an adiabatic profile, which is | ||
# taken from the adabatic conditions model. In this case, this is | ||
# is the ascii data model that loads a file containing profiles | ||
# corresponding to the Birch-Murnaghan equation of state. | ||
# In addition to the adiabatic profile, there are thermal boundary | ||
# layers at the sirface and the core-mantle boundary. | ||
subsection Initial temperature model | ||
set List of model names = adiabatic, function | ||
|
||
subsection Adiabatic | ||
set Age top boundary layer = 5e7 | ||
end | ||
|
||
subsection Function | ||
set Variable names = x,z | ||
set Function constants = pi=3.1415926 | ||
set Function expression = if(sqrt(x*x+(z-4500000)*(z-4500000))<1000000,800,0) | ||
end | ||
end | ||
|
||
|
||
subsection Mesh refinement | ||
set Refinement fraction = 0.4 | ||
set Coarsening fraction = 0.05 | ||
set Initial adaptive refinement = 0 | ||
set Initial global refinement = 4 | ||
set Strategy = temperature | ||
set Time steps between mesh refinement = 5 | ||
end | ||
|
||
|
||
subsection Heating model | ||
set List of model names = adiabatic heating, shear heating | ||
|
||
subsection Adiabatic heating | ||
set Use simplified adiabatic heating = true | ||
end | ||
end | ||
|
||
subsection Postprocess | ||
set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics | ||
|
||
subsection Visualization | ||
set Output format = gnuplot | ||
set Time between graphical output = 0 | ||
set Number of grouped files = 0 | ||
set List of output variables = material properties, gravity, adiabat, nonadiabatic temperature, nonadiabatic pressure, named additional outputs | ||
end | ||
|
||
subsection Depth average | ||
set Time between graphical output = 1e6 | ||
end | ||
end | ||
|
||
subsection Checkpointing | ||
set Time between checkpoint = 1800 | ||
end |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,44 @@ | ||
----------------------------------------------------------------------------- | ||
----------------------------------------------------------------------------- | ||
|
||
|
||
Setting up GPlates boundary velocity plugin. | ||
|
||
Input point 1 spherical coordinates: 1.571 4.870 | ||
Input point 1 normalized cartesian coordinates: 0.157 -0.988 0.000 | ||
Input point 1 rotated model coordinates: 0.157 -0.988 0.000 | ||
Input point 2 spherical coordinates: 1.571 5.240 | ||
Input point 2 normalized cartesian coordinates: 0.503 -0.864 0.000 | ||
Input point 2 rotated model coordinates: 0.503 -0.864 0.000 | ||
|
||
Model will be rotated by 0.00 degrees around axis 0.00 0.00 1.00 | ||
The ParaView rotation angles are: 0.00 0.00 0.00 | ||
The inverse ParaView rotation angles are: 0.00 0.00 0.00 | ||
|
||
Loading GPlates data boundary file ASPECT_DIR/data/velocity-boundary-conditions/gplates/current_day.gpml. | ||
|
||
|
||
Loading new velocity file did not succeed. | ||
Assuming constant boundary conditions for rest of model run. | ||
|
||
Number of active cells: 3,072 (on 5 levels) | ||
Number of degrees of freedom: 41,280 (25,344+3,264+12,672) | ||
|
||
*** Timestep 0: t=0 years | ||
Solving temperature system... 0 iterations. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 156+0 iterations. | ||
|
||
Postprocessing: | ||
Writing graphical output: output-burnman_seismic_test/solution/solution-00000 | ||
RMS, max velocity: 0.319 m/year, 1.61 m/year | ||
Temperature min/avg/max: 273 K, 2082 K, 3500 K | ||
Heat fluxes through boundary parts: -1.315e+06 W, 3.033e+06 W | ||
|
||
Termination requested by criterion: end time | ||
|
||
|
||
+---------------------------------------------+------------+------------+ | ||
+---------------------------------+-----------+------------+------------+ | ||
+---------------------------------+-----------+------------+------------+ | ||
|
Oops, something went wrong.