def test_domain_new_zealand(tmpdir): domain.RectangularSphericalSection( min_latitude=-70, max_latitude=-10, min_longitude=140, max_longitude=200, min_depth_in_km=0, max_depth_in_km=1440, boundary_width_in_degree=2.5, rotation_axis=[1.0, 1.0, 0.2], rotation_angle_in_degree=0.0).plot(plot_simulation_domain=True) images_are_identical("domain_new_zealand", str(tmpdir), tol=6.5)
def test_domain_equality(tmpdir): d1 = domain.RectangularSphericalSection(min_latitude=-70, max_latitude=-10, min_longitude=140, max_longitude=200, min_depth_in_km=0, max_depth_in_km=1440, boundary_width_in_degree=2.5, rotation_axis=[1.0, 1.0, 0.2], rotation_angle_in_degree=0.0) d2 = domain.RectangularSphericalSection(min_latitude=-70, max_latitude=-10, min_longitude=140, max_longitude=200, min_depth_in_km=0, max_depth_in_km=1441, boundary_width_in_degree=2.5, rotation_axis=[1.0, 1.0, 0.2], rotation_angle_in_degree=0.0) assert d1 == copy.deepcopy(d1) assert d1 != d2
def test_simple_european_domain(tmpdir): """ Tests the plotting of a simple European domain. """ domain.RectangularSphericalSection( min_latitude=-10, max_latitude=10, min_longitude=-10, max_longitude=10, min_depth_in_km=0, max_depth_in_km=1440, boundary_width_in_degree=2.5, rotation_axis=[1.0, 1.0, 0.2], rotation_angle_in_degree=-65.0).plot(plot_simulation_domain=True) images_are_identical("domain_simple_europe", str(tmpdir), tol=6.5)
def test_plotting_domain_north_america(tmpdir): """ Tests plotting the North American domain. """ domain.RectangularSphericalSection( min_latitude=-25, max_latitude=45, min_longitude=-120, max_longitude=0, min_depth_in_km=0, max_depth_in_km=1440, boundary_width_in_degree=9.0, rotation_axis=[0.766044443118978, 0.6427876096865393, 0], rotation_angle_in_degree=-30.0).plot(plot_simulation_domain=True) images_are_identical("domain_north_america", str(tmpdir), tol=6.5)
def test_plotting_domain_anatolia(tmpdir): """ Tests plotting of a regional scale domain. """ domain.RectangularSphericalSection( min_latitude=34.1, max_latitude=42.9, min_longitude=23.1, max_longitude=42.9, min_depth_in_km=0, max_depth_in_km=471, boundary_width_in_degree=1.46, rotation_axis=[0, 0, 1], rotation_angle_in_degree=0.0).plot(plot_simulation_domain=True) images_are_identical("domain_anatolia", str(tmpdir))
def test_plotting_domain_close_to_north_pole(tmpdir): """ Tests the plotting of a domain including the North pole. """ domain.RectangularSphericalSection( min_latitude=-30, max_latitude=40, min_longitude=-30, max_longitude=60, min_depth_in_km=0, max_depth_in_km=1440, boundary_width_in_degree=2.0, rotation_axis=[0, 1, 0], rotation_angle_in_degree=-57.5).plot(plot_simulation_domain=True) images_are_identical("domain_with_north_pole", str(tmpdir), tol=6.5)
def test_plotting_domain_iberia(tmpdir): """ Tests plotting a domain around Iberia. """ domain.RectangularSphericalSection( min_latitude=-33.0, max_latitude=-3, min_longitude=-20.0, max_longitude=27.0, min_depth_in_km=0, max_depth_in_km=599.0, boundary_width_in_degree=2.0, rotation_axis=[0, 1, 0], rotation_angle_in_degree=-57.5).plot(plot_simulation_domain=False) images_are_identical("domain_iberia", str(tmpdir))
def test_plotting_domain_japan(tmpdir): """ Tests plotting of domain around Japan. """ domain.RectangularSphericalSection( min_latitude=15.5, max_latitude=55, min_longitude=115.5, max_longitude=155.0, min_depth_in_km=0, max_depth_in_km=599.0, boundary_width_in_degree=3.0, rotation_axis=[1, 1, 1], rotation_angle_in_degree=0.0).plot(plot_simulation_domain=True) images_are_identical("domain_japan", str(tmpdir))
def test_plotting_edge_case__domain(tmpdir): """ Tests the plotting of a heavily deformed domain. Not realistic but useful as a test case. """ domain.RectangularSphericalSection( min_latitude=-10, max_latitude=10, min_longitude=-120, max_longitude=170, min_depth_in_km=0, max_depth_in_km=1440, boundary_width_in_degree=2.5, rotation_axis=[1.0, 1.0, 0.2], rotation_angle_in_degree=-75.0).plot(plot_simulation_domain=True) images_are_identical("domain_edge_case", str(tmpdir), tol=8)
def test_domain_new_zealand_regional_dateline_wraparound(tmpdir): """ This is a much smaller domain but still having a dateline wraparound. Make sure this works with basemap. """ domain.RectangularSphericalSection( min_latitude=-55.0, max_latitude=-28.0, min_longitude=155.0, max_longitude=200.0, min_depth_in_km=0, max_depth_in_km=200, boundary_width_in_degree=1.0, rotation_axis=[1.0, 1.0, 1.0], rotation_angle_in_degree=0.0).plot(plot_simulation_domain=False) images_are_identical("domain_new_zealand_regional_dateline_wrap", str(tmpdir))