geodynamics / vq

Virtual Quake is a boundary element code designed to investigate long term fault system behavior and interactions between faults through stress transfer.
Other
12 stars 24 forks source link

Dip/Strike fix, Greens function fix, added PyVQ options --levels, --traces #75

Closed kwschultz closed 9 years ago

kwschultz commented 9 years ago

SimElement: Fixed the handling of Strike angle to match convention that 0 < strike < 2*pi

SimElement: Fixed the order of the cross product to find an element's normal vector to match the VQ mesher convention that faults dip to the left of the strike direction. This change will affect the computed dip angles for sim elements.

QuakeLibOkada: Fixed the sign on the maximum depth for the fault. Green's functions now match those presented in Okubo 1992, with the caveat that there's still some small location error that is fixed when increasing mesh resolution. e.g. gravity field for a single thrust fault does not match expected pattern, but meshing the fault into 10+ elements fixes the field and it agrees with the expected pattern.

--traces option for pyvq: Will plot fault traces given a fault model, can be used in conjunction with --use_sections to highlight specific sections

--levels option for pyvq: Used to draw contours on field plots (e.g. --levels -50 -20 -10 0 10 20 50)

eheien commented 9 years ago

Thanks!