[CIG-LONG] Gale: NonNewtonian viscosity

[Walter Landry]

<Guillaume.Duclaux at csiro.au> wrote:
> Hi all,
>  
> Few month ago I pointed out a small error in the encoding of the NonNewtonian
> viscosity to Walter. I just paste below a part of Walter's answer to my
> message.
> 
>> the line
> 
> viscosity=exp(T_0/(2*n*T))/(2*A);
> 
> has an extraneous factor of 2, so it should be:
> 
> viscosity=exp(T_0/(n*T))/(2*A);
> 
> Also, to correct this small feature, one can multiply T_0 with a factor of 2
> in the input file.

For some reason I never committed the fix for this.  Done (rev 15938).

> I found today that the minViscosity and maxViscosity properties are inverted.
> (Never used these properties before...)
> Indeed, minViscosity applies an upper threshold (max viscosity) to the
> effective Viscosity, and maxViscosity a lower one (min Viscosity). I guess it
> should the other way around. 

Are you sure?  I am looking at the code, and it seems correct.  Do you
have an input file?

Cheers,
Walter Landry
walter at geodynamics.org