def test_dip_shear():
ds = dip_shear(strike=strike,
dip=dip,
rho=rho,
g=g,
mxx=mxx,
myy=myy,
mxy=mxy,
mzz=mzz,
mxz=mxz,
myz=myz,
txx=txx,
tyy=tyy,
txy=txy,
depth=depth)
ds = np.int(np.round(ds))
assert ds == 924052, 'dip shear stress'
mxy=search_df.mxy,
mzz=search_df.mzz,
mxz=search_df.mxz,
myz=search_df.myz,
txx=search_df.txx,
tyy=search_df.tyy,
txy=search_df.txy,
depth=search_df.depth)
search_df['tau_d'] = scv.dip_shear(strike=search_df.strike,
dip=search_df.dip,
rho=rho,
g=g,
mxx=search_df.mxx,
myy=search_df.myy,
mxy=search_df.mxy,
mzz=search_df.mzz,
mxz=search_df.mxz,
myz=search_df.myz,
txx=search_df.txx,
tyy=search_df.tyy,
txy=search_df.txy,
depth=search_df.depth)
search_df['tau_rake'] = hsp.get_rake_from_shear_components(
strike_shear=search_df.tau_s, dip_shear=search_df.tau_d)
search_df['rake_misfit_rad'] = np.radians(
scv.angle_difference(search_df.slip_rake,
search_df.tau_rake,
return_abs=True))
# OK, now let's do some calculationsi
print('calculating fault stresses')
search_df['tau_s'] = scv.strike_shear(strike = search_df.strike,
dip=search_df.dip, rho=rho, g=g,
mxx=search_df.mxx, myy=search_df.myy,
mxy=search_df.mxy, mzz=search_df.mzz,
mxz=search_df.mxz, myz=search_df.myz,
txx=search_df.txx, tyy=search_df.tyy,
txy=search_df.txy, depth=search_df.depth)
search_df['tau_d'] = scv.dip_shear(strike = search_df.strike,
dip=search_df.dip, rho=rho, g=g,
mxx=search_df.mxx, myy=search_df.myy,
mxy=search_df.mxy, mzz=search_df.mzz,
mxz=search_df.mxz, myz=search_df.myz,
txx=search_df.txx, tyy=search_df.tyy,
txy=search_df.txy, depth=search_df.depth)
print("dropping columns that won't be used more")
search_df.drop(labels=['mxx', 'myy', 'mxy', 'mzz', 'mxz', 'myz', 'depth',
'strike', 'dip'], axis=1, inplace=True)
search_df['tau_rake'] = hsp.get_rake_from_shear_components(strike_shear=
search_df.tau_s,
dip_shear=
search_df.tau_d)
search_df['rake_misfit_rad'] = np.radians(scv.angle_difference(
