Hello there,
I am currently using version 2.2.0 of ASPECT and am trying to include phase transitions in the model. When working with only one transition, everything works fine. But upon including two transition depths I keep getting error messages that the phase transition depths are not provided in the correct format, although I provided them as comma separated lists, just as ASPECT demands. See the code and error message below. Any help would be much appreciated. Thank you all in advance.
##
# Setup
# - with unstable initial conditions
# - simple material model
set Dimension = 2
set Start time = 0
set End time = 1e9
set Use years in output instead of seconds = true
set Output directory = NGP-phase-2layer
# Geometry definition
# Rectangular box of 3000 x 3000 km
subsection Geometry model
set Model name = box
subsection Box
set X extent = 3000000
set Y extent = 3000000
end
end
############### Boundary conditions
# Set only tangential velocities at the boundaries (free slip)
subsection Boundary velocity model
set Tangential velocity boundary indicators = 0,1,2,3
set Prescribed velocity boundary indicators =
end
# Set temperature at top to 1600 K and 4000 K at bottom
subsection Boundary temperature model
# Fix temperatures at top and bottom
set Fixed temperature boundary indicators = 2,3
set Model name = box
subsection Box
set Bottom temperature = 4000
set Top temperature = 1600
end
end
subsection Initial temperature model
set Model name = function
subsection Function
set Variable names = x,z
# Initial temperature is roughly adiabatic + a big hot blob at the bottom
set Function expression = 3100 - z*0.5e-3 + if (x>1.5e6 & x<2.5e6 & z<1e6,500,0)
end
end
# Material properties
subsection Material model
set Model name = latent heat
subsection Latent heat
set Thermal expansion coefficient = 4e-5
set Viscosity = 1e22
set Thermal conductivity = 4.7
set Compressibility = 4e-12
## Add bogus phase transitions - this is nonsense!
set Phase transition Clapeyron slopes = 2e6, -3e6
set Phase transition temperatures = 1700, 2000
set Phase transition depths = 500e3, 1500e3
set Phase transition density jumps = 250, 400
set Phase transition widths = 1000, 1000
set Viscosity prefactors = 1.0, 1.0, 1.0
set Corresponding phase for density jump = 0, 0
end
end
# Gravity is constant
subsection Gravity model
set Model name = vertical
subsection Vertical
set Magnitude = 10.0
end
end
#without adiabatic heating and constant heating
subsection Heating model
set List of model names = shear heating, latent heat
end
subsection Mesh refinement
set Initial adaptive refinement = 0
set Initial global refinement = 5
set Time steps between mesh refinement = 10
set Strategy = temperature,density
end
subsection Postprocess
set List of postprocessors = visualization
subsection Visualization
set Number of grouped files = 0
set Output format = vtu
set Time between graphical output = 1e6
set List of output variables = density,heating
set Interpolate output = true
end
end
The error message:
----------------------------------------------------
Exception 'ExcMessage ("The required format for property <" + property_name + "> was not found. Specify a comma separated " + "list of `<double>' or `<key1> : <double>|<double>|..., " + "<key2> : <double>|... , ... '.")' on rank 0 on processing:
An error occurred in line <198> of file </home/ubuntu/aspect/source/utilities.cc> in function
std::multimap<std::__cxx11::basic_string<char>, double> aspect::Utilities::{anonymous}::parse_string_to_map(const string&, const std::vector<std::__cxx11::basic_string<char> >&, const string&)
The violated condition was:
false
Additional information:
The required format for property <Phase transition depths> was not found. Specify a comma separated list of `<double>' or `<key1> : <double>|<double>|..., <key2> : <double>|... , ... '.