Hello,
I have completed the strikeslip-2d simulation. After completing the static co-seismic slip, I can use pylith_viz to visualize the deformation. But how can I quantitatively inspect the parameters of the simulation results, such as the deformation velocity field?
I also have another question: if I want to study the viscosity of the upper mantle in a region, how can I perform inversion using PyLith? I would appreciate some guidance or ideas.
Thanks!
The data in the HDF5 files output from PyLith can be loaded into many visualization tools (pylith_viz, ParaView, Visit) and software tools (Matlab, Python, C/C++). You can perform additional post-processing and analysis using these tools.
For quasi-static problems, the solution field only contains displacement, so you will have to post-process the results to compute velocity.
To constrain the viscosity of the upper mantle in a region using PyLith, I would start with a sensitivity study that examines the range of values given constraints from observations, etc. PyLith does not yet have data assimilation and adjoint capabilities to directly constrain material properties; this is a feature we would like to add.