[aspect-devel] questions

Timo Heister heister at math.tamu.edu
Wed Apr 17 07:50:59 PDT 2013


This is not exposed in the paper, but assumed to be = 2.0 (the only
options are 1.0 and 2.0 and 1.0 uses a slightly different definition
of the entropy). I recommend keeping it set at 2.0. I will add a
comment about this in the code.

On Wed, Apr 17, 2013 at 9:33 AM, Thieulot, C. (Cedric) <c.thieulot at uu.nl> wrote:
>>
>>> - turning now to the advection stabilisation parameters, and looking for instance at box.prm, I see three parameters: alpha, beta and cR.
>>> Now lookig at the GJI article about Aspect (kronbichler et al, 2012), section 3.2.6 of the article details the algorithm.
>>> As it turns out the alpha_max of the article seems to be beta in the code (a bit misleading ?)
>>> However, what are the equivalent of alpha and cR in the article ?
>>
>> Yeah, the naming is somewhat misleading:
>> parameters.stabilization_beta = alpha_max in the paper
>> parameters.stabilization_c_R = alpha_E in the paper
>
> and so what does parameters.stabilization_alpha correspond to in the paper ?
>
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel



-- 
Timo Heister
http://www.math.tamu.edu/~heister/


More information about the Aspect-devel mailing list