[CIG-LONG] Non-Newtonienne viscosity with depth dependency

[Walter Landry]

nicolas riel <a.nicolas.riel at gmail.com> wrote:
> Hi everyone,
> 
> I'm trying to modify the NonNewtonian rheology to add a pressure dependency.
> 
> I would like the Viscosity to be multiplied by 10 for when the pressure >
> 2.4e10 pa.
> I tried to looked in the DruckerPrager rheology to see how the pressure
> dependency is implemented, but for how the pressure dependency do not work and
> send me the following error message: "Error in NonNewtonian: T, n, and A must
> all be non-zero".

This sounds more like an error in your input file.  Does your input
file work with the old NonNewtonian rheology?  Can you send that input
file?

> I attached my modified NonNewtonian.c file.

I do not see any correction for HydrostaticTerm.  You will have to add
that in to get a correct pressure.

Cheers,
Walter Landry
walter at geodynamics.org