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))
