Stress and strain orientation in PyLith

Dear all,

I would like to ask about the orientation of the stress and strain tensor in PyLith result.
As previously, I have mentioned that I run the simulation of 2D subduction and put also a splay fault in the geometry. I would like to calculate the coulomb stress on this splay fault.

My question is, do the stress and strain tensor in Pylith already have orientation thus I can calculate the coulomb stress directly using the stress and strain tensor from the result? Or I should consider the geometry (strike,dip) of the splay fault manually and do post-calculation to retrieve the strain and normal traction?

Thank you in advance.

Best Regards,
Alvina

The traction changes on the fault computed by PyLith are in the local fault coordinate system, and the components correspond to left-lateral shear, updip shear, and fault-normal stress (look in the PyLith manual for the order). These components may then be used to compute the Coulomb stress changes. You may also find it useful to use the orientation info provided in the fault_info file (strike_dir, dip_dir, normal_dir), which provides the orientation of the local fault surface in the global system.