def visualize(slip, faultxyz, horxyz):
fault = geoprobe.swfault('/data/nankai/data/swFaults/jdk_oos_splay_large_area_depth.swf')
azimuth = np.degrees(np.arctan2(*slip[::-1]))
azimuth = 90 - azimuth
mag = np.linalg.norm(slip) / 1000
f = FaultModel(fault, horxyz, origxyz=horxyz, ve=3, azimuth=azimuth,
slip=mag, alpha=data.alpha, calc_fault=faultxyz)
f.configure_traits()
def main():
fault = geoprobe.swfault('/data/nankai/data/swFaults/jdk_oos_splay_large_area_depth.swf')
faultxyz = data.to_world(data.to_xyz(data.fault))
horxyz = data.to_world(data.to_xyz(data.horizons[0]))[::100]
slip = invert_slip(faultxyz, horxyz, alpha=data.alpha)
azimuth = np.degrees(np.arctan2(*slip[::-1]))
azimuth = 90 - azimuth
mag = np.linalg.norm(slip) / 1000
f = FaultModel(fault, horxyz, origxyz=horxyz, ve=3, azimuth=azimuth,
slip=mag, alpha=data.alpha, calc_fault=faultxyz)
f.configure_traits()
def main():
fault = geoprobe.swfault(
'/data/nankai/data/swFaults/jdk_oos_splay_large_area_depth.swf')
faultxyz = data.to_world(data.to_xyz(data.fault))
horxyz = data.to_world(data.to_xyz(data.horizons[0]))[::100]
slip = invert_slip(faultxyz, horxyz, alpha=data.alpha)
azimuth = np.degrees(np.arctan2(*slip[::-1]))
azimuth = 90 - azimuth
mag = np.linalg.norm(slip) / 1000
f = FaultModel(fault,
horxyz,
origxyz=horxyz,
ve=3,
azimuth=azimuth,
slip=mag,
alpha=data.alpha,
calc_fault=faultxyz)
f.configure_traits()
