[CIG-SHORT] How to output energy quantities from Pylith

[Charles Williams]

Dear Xiao,

These quantities are not computed within PyLith; however, you can easily compute them from output stresses, strains, and other state variables.  The simplest solution would be to read the state variables output file using h5py, compute the necessary quantities, and then write HDF5 output again using h5py.  You also need to call the PETSc function to generate an associated XDMF file for visualization.  I am attaching a simple example that may have something similar to what you need.

Cheers,
Charles



> On 1/08/2016, at 10:44 AM, Xiao Ma <xiaoma5 at illinois.edu> wrote:
> 
> Hi all, 
> 
> 
> I have a question about how to output energy quantities such as strain energy , kinetic energy, frictional dissipation energy ,etc from Pylith ? 
> 
> 
> Thank you,
> 
> Xiao
> 
> -- 
> -------------------------------------------------------------------------------------
> Xiao Ma
> Graduate Research Assistant
> PhD student in Structural Engineering
> University of Illinois at Urbana-Champaign
> E-mail: xiaoma5 at illinois.edu <mailto:xiaoma5 at illinois.edu>
> --------------------------------------------------------------------------------------
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short

Charles A. Williams
Scientist
GNS Science
1 Fairway Drive, Avalon
PO Box 30368
Lower Hutt  5040
New Zealand
ph (office): 0064-4570-4566
fax (office): 0064-4570-4600
C.Williams at gns.cri.nz <mailto:C.Williams at gns.cri.nz>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20160801/554cc99f/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: stress_invars.py
Type: text/x-python-script
Size: 9925 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20160801/554cc99f/attachment.bin>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20160801/554cc99f/attachment-0001.html>