[CIG-SHORT] Inconsistency Between Linear and Powerlaw Viscoelastic Gravitational Forces

[Jeffrey Thompson]

Dear CIG Developers,

I've been running into difficulties with the magnitude of gravitational
stresses using a PowerLaw3D rheology.  As a consistency check, I ran two
models with the same material parameters, using a MaxwellIsotropic3D
rheology first and a PowerLaw3D rheology (with n=1) second.  I obtained
wildly different values of stress within the two model runs, with the
magnitudes of stresses varying by about 10 orders of magnitude.  I've
attached a tarball with my model parameters for both cases, and was
wondering if you could take a look at this problem.

I'm currently running PyLith 1.5.1 using the tarballed version available off
of the CIG website (i.e. not an svn version).

Thanks,
Jeffrey Thompson
Seismolab
California Institute of Technology
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-short/attachments/20101012/402be80f/attachment-0001.htm 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: RunFiles.tar.gz
Type: application/x-gzip
Size: 872957 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-short/attachments/20101012/402be80f/attachment-0001.bin