[CIG-SHORT] on the description of power-law-coefficient

[Ikuo Cho]

Dear Pylith Developers:

I have a question on a description of the parameter
power-law-coefficient in a spatial database: Does the dimension corresponding to the power-law-coefficient have a
meaning? I'm asking this because I want to solve a problem of spatially variable power-law exponents.
More concretely, is the following description correct or not?

0. 0. 0. 2500. 3000. 5500 1.e-37 3.0
0. 0. 20. 2500. 3000. 5500 1.e-37 2.0

with 

value-names = density vs vp power-law-coefficient power-law-exponent 
value-units= kg/m**3 m/s m/s Pa**(-3.0)/s None 

when the flow law is given as

(strain rate) = At (stress)**3.0, where At=1.e-37 [Pa**(-3.0)/s], at
(x,y,z)=(0,0,0)

and 

(strain rate) = At (stress)**2.0, where At=1.e-37 [Pa**(-2.0)/s], at
(x,y,z)=(0,0,20).


Ikuo Cho
----------------------------------------------------------
Geological Survey of Japan,
National Institute of Advanced Industrial Science and Technology
Tsukuba Central 7, Tsukuba 305-8567 Japan
Tel +81-29-861-3891, Fax +81-29-861-3682
----------------------------------------------------------