A Problem when using “Initial temperature model/Spherical gaussian perturbation” option

I set the .prm file parameters as follows:

subsection Geometry model
  set Model name =  spherical shell
  subsection Spherical shell
    set Outer radius  = 1740000
    set Inner radius  = 400000
    set Opening angle = 360

subsection Boundary temperature model
  set Fixed temperature boundary indicators = bottom, top
  set List of model names =  spherical constant 
  subsection Spherical constant
    set Inner temperature = 1843
    set Outer temperature = 273

subsection Initial temperature model
  set List of model names =  ascii profile ,  spherical gaussian perturbation
  set List of model operators = add
  subsection Ascii profile
    set Data directory  = data/
    set Data file name  = Temperature.txt

and the error information is:

An error occurred in line <301> of file </THL7/home/maqz2/.spack/stage/spack-stage-aspect-develop-p2hfes7wazqttrpllhwuax4iohoxsxap/spack-src/source/initial_temperature/spherical_shell.cc> in function
    void aspect::InitialTemperature::SphericalGaussianPerturbation<dim>::parse_parameters(dealii::ParameterHandler&) [with int dim = 2]
The violated condition was: 
Additional information: 
    This initial condition can only be used if a boundary temperature is prescribed.

As mentioned above, I have set the boundary temperature as Fixed temperature boundary, so I couldn’t find out how to set the perturbation next. I’d be grateful to you for answering my question.

Ah, this is unfortunate. It’s a bug in ASPECT. I’m not sure how this ever worked, because we are initializing the initial temperature models before the boundary temperature models, and as a consequence, when we get to the point where you get the error, the temperature model has not yet been set up. What happens if you choose only the spherical gaussian perturbation as initial temperature, i.e., not add it together with the ascii profile? Does that work?


1 Like

Hi, Wolfgang
Sorry for replying you late. I’ve tested four .prm files in “Initial temperature model” option: (1)My .prm file with only gaussian perturbation, (2)My .prm file with only hexagonal perturbation, (3) “shell_simple_2d.prm” in cookbooks (with only hexagonal perturbation), (4) “shell_simple_2d.prm” with only gaussian perturbation.
And only case1, i.e. My .prm file with only gaussian perturbation can’t work , the error message is the same as before. Is this means spherical gaussian perturbation doesn’t work in ASPECT? Maybe I should use spherical hexagonal perturbation next.