def disambigue_sdr(self): r2d = 180. / math.pi d2r = 1. / r2d sdr = 'strike', 'dip', 'slip-rake' strike, dip, rake = [self._params[k] for k in sdr] (alpha, beta, gamma) = (dip * d2r, strike * d2r, -rake * d2r) (alpha, beta, gamma) = moment_tensor.unique_euler(alpha, beta, gamma) self._params['dip'] = alpha * r2d self._params['strike'] = beta * r2d self._params['slip-rake'] = -gamma * r2d
def disambigue_sdr(self): r2d = 180./math.pi d2r = 1./r2d sdr = 'strike', 'dip', 'slip-rake' strike, dip, rake = [ self._params[k] for k in sdr ] (alpha, beta, gamma) = (dip*d2r, strike*d2r, -rake*d2r) (alpha, beta, gamma) = moment_tensor.unique_euler(alpha,beta,gamma) self._params['dip'] = alpha*r2d self._params['strike'] = beta*r2d self._params['slip-rake'] = -gamma*r2d
def beachball(source, filename, **conf_overrides): if 'strike' in source.keys() and 'dip' in source.keys() and 'slip-rake' in source.keys(): r2d = 180./math.pi d2r = 1./r2d strike, dip, rake = source['strike'], source['dip'], source['slip-rake'] dip, strike, mrake = [r2d*xx for xx in moment_tensor.unique_euler(dip*d2r, strike*d2r, -rake*d2r)] rake = -mrake return beachball_mt(None, filename, sdr=(strike,dip,rake), **conf_overrides) else: return []