def test_mendez_distances(self): dist = Distances(None, None, None, None, None, None, None, None, None, None, None) rupture_centroid_lat = asarray((-31.0, -33, -36)) rupture_centroid_lon = asarray((116.0, 118, 222)) site_lat = asarray((-31, -32, -33, -34)) site_lon = asarray((116.0, 116.0, 116.0, 116.0)) distance1 = dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, 'Epicentral', projection) distance2 = dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, 'Obsolete_Mendez_epicentral', projection, trace_start_lat=rupture_centroid_lat, trace_start_lon=rupture_centroid_lon, rupture_centroid_x=0.0, rupture_centroid_y=0.0) assert allclose(distance1, distance2)
def test_Epicentral2(self): dist = Distances(None,None,None,None,None,None,None,None,None,None,None) distance_type='Epicentral' rupture_centroid_lat=asarray((-31.0,-33,-36)) rupture_centroid_lon=asarray((116.0,118,222)) site_lat=asarray((-31,-32,-33,-34)) site_lon=asarray((116.0,116.0,116.0,116.0)) distance=dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, distance_type, projection) d = As_The_Cockey_Flies(rupture_centroid_lat, rupture_centroid_lon, site_lat, site_lon) assert allclose(d,distance,rtol=0.1)
def test_Epicentral(self): dist = Distances(None,None,None,None,None,None,None,None,None,None,None) distance_type='Epicentral' rupture_centroid_lat=asarray((-31.0)) rupture_centroid_lon=asarray((116.0)) site_lat=asarray((-31,-32,-33,-34)) site_lon=asarray((116.0,116.0,116.0,116.0)) distance=dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, distance_type, projection) d=asarray((0,1,2,3))*(1.852*60) d=d[:,newaxis] d2=As_The_Cockey_Flies(rupture_centroid_lat, rupture_centroid_lon, site_lat, site_lon) assert allclose(d,distance,rtol=0.001) assert allclose(d,d2)
def test_mendez_distances(self): dist = Distances(None,None,None,None,None,None,None,None,None,None,None) rupture_centroid_lat=asarray((-31.0,-33,-36)) rupture_centroid_lon=asarray((116.0,118,222)) site_lat=asarray((-31,-32,-33,-34)) site_lon=asarray((116.0,116.0,116.0,116.0)) distance1=dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, 'Epicentral', projection) distance2=dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, 'Obsolete_Mendez_epicentral', projection, trace_start_lat=rupture_centroid_lat, trace_start_lon=rupture_centroid_lon, rupture_centroid_x=0.0, rupture_centroid_y=0.0) assert allclose(distance1,distance2)
def test_Hypocentral(self): dist = Distances(None,None,None,None,None,None,None,None,None,None,None) distance_type='Hypocentral' depths=10.0 rupture_centroid_lat=asarray((-31.0,-33,-36)) rupture_centroid_lon=asarray((116.0,118,222)) site_lat=asarray((-31,-32,-33,-34)) site_lon=asarray((116.0,116.0,116.0,116.0)) d=dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, 'Epicentral', projection) distance=dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, distance_type, projection) d=sqrt(d*d+depths*depths) assert allclose(d,distance)
def test_Hypocentral(self): dist = Distances(None, None, None, None, None, None, None, None, None, None, None) distance_type = 'Hypocentral' depths = 10.0 rupture_centroid_lat = asarray((-31.0, -33, -36)) rupture_centroid_lon = asarray((116.0, 118, 222)) site_lat = asarray((-31, -32, -33, -34)) site_lon = asarray((116.0, 116.0, 116.0, 116.0)) d = dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, 'Epicentral', projection) distance = dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, distance_type, projection) d = sqrt(d * d + depths * depths) assert allclose(d, distance)
def test_Epicentral2(self): dist = Distances(None, None, None, None, None, None, None, None, None, None, None) distance_type = 'Epicentral' rupture_centroid_lat = asarray((-31.0, -33, -36)) rupture_centroid_lon = asarray((116.0, 118, 222)) site_lat = asarray((-31, -32, -33, -34)) site_lon = asarray((116.0, 116.0, 116.0, 116.0)) distance = dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, distance_type, projection) d = As_The_Cockey_Flies(rupture_centroid_lat, rupture_centroid_lon, site_lat, site_lon) assert allclose(d, distance, rtol=0.1)
def test_Epicentral(self): dist = Distances(None, None, None, None, None, None, None, None, None, None, None) distance_type = 'Epicentral' rupture_centroid_lat = asarray((-31.0)) rupture_centroid_lon = asarray((116.0)) site_lat = asarray((-31, -32, -33, -34)) site_lon = asarray((116.0, 116.0, 116.0, 116.0)) distance = dist.raw_distances(site_lat, site_lon, rupture_centroid_lat, rupture_centroid_lon, lengths, azimuths, widths, dips, depths, depths_to_top, distance_type, projection) d = asarray((0, 1, 2, 3)) * (1.852 * 60) d = d[:, newaxis] d2 = As_The_Cockey_Flies(rupture_centroid_lat, rupture_centroid_lon, site_lat, site_lon) assert allclose(d, distance, rtol=0.001) assert allclose(d, d2)