[CIG-SHORT] Output normal and shear stress on a fault plane in .h5 format

[Brad Aagaard]

Josimar,

"slip", "slip_rate", and "traction" are all vectors. You request output 
by naming the field and get all components. For fields defined on the 
fault, such as slip, slip_rate, and traction, the output is given in 
fault system coordinates, so the components are shear and normal in 2-D 
and shear-horiz, shear-up-dip, and normal in 3-D. The sign convention 
for traction is the same as slip (positive shear traction is left-lateral).

Regards,
Brad


On 02/12/2018 12:37 PM, Josimar Alves da Silva wrote:
> Hi Brad,
> 
> Could you please let me know if PyLith can output directly normal and 
> shear stress on a fault plane ? Or do we have to post process PyLith 
> outputs to get these values ?
> 
> I have the impression that when I set the fault output for 
> vertex_data_field = [traction]  I get X,Y,Z values that do not 
> correspond to the fault normal, for example (Please confirm).
> 
> Note that I tried something like vertex_data_field = [traction_normal] 
> but it seems that this is not available (see error below),
> 
> thank you
> Josimar
> 
> 
> #### This outputs X,Y,Z tractions on a fault plane
> [pylithapp.problem.interfaces.Thrust_fault.output]
> ##vertex_data_fields = [traction, slip ]
> vertex_data_fields = [traction, slip]
> skip = 0
> writer = pylith.meshio.DataWriterHDF5
> writer.filename = output/Thrust_fault.h5
> 
> 
> 
> #### Error message when I set vertex_data_field=[traction_normal]
> 
> ValueError: Requested fields not available for output.
> 
> Data provider: 'faultcohesivedyn'
> 
> Field type: 'vertex'
> 
> Data type: 'data'
> 
> Available fields:  'slip' 'slip_rate' 'traction' 'cumulative_slip' 
> 'previous_slip'
> 
> Fields not available:  'traction_normal'
> 
>