[aspect-devel] Depletion density change

[Juliane Dannberg]

Hi Mohamed,

sorry for the late reply.
You are right, we should allow negative density changes (and the 
documentation of the parameter already describes what happens in this 
case).
I have created a pull request here: 
https://github.com/geodynamics/aspect/pull/1281

Best,
Juliane


On 11/24/2016 06:35 AM, Mohamed Gouiza wrote:
>
> Hi,
>
> In the Material model/Melt simple, the Depletion density change 
> parameter takes only positive values (0 <= v <= MAX_DOUBLE).
>
> It should take negative values too, as depleted material usually has 
> lower densities.
>
> The same parameter in the Material model/Melt global accepts negative 
> value: -MAX_DOUBLE <= v <= MAX_DOUBLE
>
> Cheers,
>
> -------------------------------------------------
> Mohamed Gouiza, Research Fellow
> Basin Structure Group, Institute of Applied Geosciences
> University of Leeds, School of Earth and Environment
>
> Leeds,  LS2 9JT, UK
>
> M.Gouiza at leeds.ac.uk <mailto:M.Gouiza at leeds.ac.uk>
> +44 7985 782073
> -------------------------------------------------
>
>
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161130/8ea8ac69/attachment-0001.html>