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

[nicolas riel]

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".

I attached my modified NonNewtonian.c file.

Has someone already tried to do this kind of modifications ?

Thank you
Nicolas






-- 
Nicolas RIEL
Doctorant / Phd student

Institut des Sciences de la Terre
(ISTerre-UMR CNRS/UJF 5275) Maison des Géosciences
1381, Rue de la Piscine
38400 - Saint Martin d'Hères

http://nicolas.riel.free.fr
Tél.   : 33.(0)4.76.51.40.57
Fax.  : 33.(0)4.76.51.40.58
Mob. : 33.(0)6.65.58.97.99
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20110612/51ee8756/attachment.htm