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(