[aspect-devel] The syntax for the function expression

Nan Zhang Nan.Zhang at colorado.edu
Mon Sep 12 23:34:05 PDT 2016


Hi Max,

Could you point me where is your source code that parse the "Viscosity
depth function". I traced back to "viscosity_function.parse_parameters(prm);",
and did not see its difference from the parse_parameters() for
function in "Compositional
initial conditions". Sure, they certainly use the same class ParseFunction.
I want you check where you used "d" as the variable, while "compositional
initial conditions" use "r, theta, phi" as the variables.

Thanks!
Nan

On Tue, Sep 13, 2016 at 12:57 PM, Max Rudolph <maxwellr at gmail.com> wrote:

> Nan,
> Jacky Austermann and I wrote the code for the depth-dependent material
> model at the 2015 hackathon, and we added tests using each of the methods
> to prescribe depth dependence in 2D and 3D, so you can look at these for
> examples. The format of the function is described in the deal.ii manual for
> the class ParsedFunction. If you want a linear increase in the lower
> mantle, I agree that the function may be the easiest way to achieve this.
> It works in 3D spherical geometry and also 2D annulus. Each geometry model
> provides a means to query the depth from the boundary.
>
> Max
>
> On Mon, Sep 12, 2016 at 8:42 PM, Nan Zhang <Nan.Zhang at colorado.edu> wrote:
>
>> Hi Max,
>>
>> Thanks for your comments! Do you think if it works for the spherical
>> geometry? I am going to try.
>>
>> The reason that I prefer using "function" option is the need for a
>> linear-increasing viscosity prefactor for the lower mantle. The "List"
>> option won't work. I guess the "File" option should work too. Is there any
>> example in the manual?
>>
>> My concern is there are many "function expression" within various
>> subsections. It seems that they have different syntax or format. When I
>> checked the source code as I mentioned in my previous email, it is hard to
>> see the difference because they all kind of use "parse_parameters()".
>>
>> How did you find there examples in the user manual? I have read the
>> manual, while certainly forgot it quickly.
>>
>> Cheers,
>> Nan
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20160913/691f8d9c/attachment.html>


More information about the Aspect-devel mailing list