Fix tests for Aug 2023 updated remote datasets (#2636)

Co-authored-by: Wei Ji <[email protected]>

Author: seisman

pygmt/tests/test_datasets_earth_age.py

@@ -40,8 +40,8 @@ def test_earth_age_01d():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), 0.167381, rtol=1e-5)
-    npt.assert_allclose(data.max(), 338.0274, rtol=1e-5)
+    npt.assert_allclose(data.min(), 0.17, atol=0.01)
+    npt.assert_allclose(data.max(), 338.02, atol=0.01)
 
 
 def test_earth_age_01d_with_region():
@@ -53,8 +53,8 @@ def test_earth_age_01d_with_region():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), 11.293945, rtol=1e-5)
-    npt.assert_allclose(data.max(), 125.1189, rtol=1e-5)
+    npt.assert_allclose(data.min(), 11.29, atol=0.01)
+    npt.assert_allclose(data.max(), 125.12, atol=0.01)
 
 
 def test_earth_age_01m_without_region():
@@ -86,5 +86,5 @@ def test_earth_age_01m_default_registration():
     assert data.coords["lat"].data.max() == 5.0
     assert data.coords["lon"].data.min() == -10.0
     assert data.coords["lon"].data.max() == -9.0
-    npt.assert_allclose(data.min(), 88.63)
-    npt.assert_allclose(data.max(), 125.25)
+    npt.assert_allclose(data.min(), 88.63, atol=0.01)
+    npt.assert_allclose(data.max(), 125.25, atol=0.01)

pygmt/tests/test_datasets_earth_free_air_anomaly.py

@@ -40,8 +40,8 @@ def test_earth_faa_01d():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -275.75)
-    npt.assert_allclose(data.max(), 308.375)
+    npt.assert_allclose(data.min(), -275.85, atol=0.025)
+    npt.assert_allclose(data.max(), 308.35, atol=0.025)
 
 
 def test_earth_faa_01d_with_region():
@@ -53,8 +53,8 @@ def test_earth_faa_01d_with_region():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), -58.75)
-    npt.assert_allclose(data.max(), 69.524994)
+    npt.assert_allclose(data.min(), -58.475, atol=0.025)
+    npt.assert_allclose(data.max(), 69.975, atol=0.025)
 
 
 def test_earth_faa_01m_without_region():
@@ -78,5 +78,5 @@ def test_earth_faa_01m_default_registration():
     npt.assert_allclose(data.coords["lat"].data.max(), 4.991666666)
     npt.assert_allclose(data.coords["lon"].data.min(), -9.99166666)
     npt.assert_allclose(data.coords["lon"].data.max(), -9.00833333)
-    npt.assert_allclose(data.min(), -51)
-    npt.assert_allclose(data.max(), 113.675)
+    npt.assert_allclose(data.min(), -49.225, atol=0.025)
+    npt.assert_allclose(data.max(), 115.0, atol=0.025)

pygmt/tests/test_datasets_earth_geoid.py

@@ -40,8 +40,8 @@ def test_earth_geoid_01d():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -106.45)
-    npt.assert_allclose(data.max(), 83.619995)
+    npt.assert_allclose(data.min(), -106.45, atol=0.01)
+    npt.assert_allclose(data.max(), 83.62, atol=0.01)
 
 
 def test_earth_geoid_01d_with_region():
@@ -53,8 +53,8 @@ def test_earth_geoid_01d_with_region():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), 4.87)
-    npt.assert_allclose(data.max(), 29.89)
+    npt.assert_allclose(data.min(), 4.87, atol=0.01)
+    npt.assert_allclose(data.max(), 29.89, atol=0.01)
 
 
 def test_earth_geoid_01m_without_region():
@@ -86,5 +86,5 @@ def test_earth_geoid_01m_default_registration():
     assert data.coords["lat"].data.max() == 5.0
     assert data.coords["lon"].data.min() == -10.0
     assert data.coords["lon"].data.max() == -9.0
-    npt.assert_allclose(data.min(), 20.34)
-    npt.assert_allclose(data.max(), 30.039999)
+    npt.assert_allclose(data.min(), 20.34, atol=0.01)
+    npt.assert_allclose(data.max(), 30.04, atol=0.01)

pygmt/tests/test_datasets_earth_magnetic_anomaly.py

@@ -40,8 +40,8 @@ def test_earth_mag_01d():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -384)
-    npt.assert_allclose(data.max(), 1057.2)
+    npt.assert_allclose(data.min(), -384.0, atol=0.2)
+    npt.assert_allclose(data.max(), 1057.2, atol=0.2)
 
 
 def test_earth_mag_01d_with_region():
@@ -53,8 +53,8 @@ def test_earth_mag_01d_with_region():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), -180.40002, rtol=1e-5)
-    npt.assert_allclose(data.max(), 127.39996, rtol=1e-5)
+    npt.assert_allclose(data.min(), -180.4, atol=0.2)
+    npt.assert_allclose(data.max(), 127.4, atol=0.2)
 
 
 def test_earth_mag_02m_without_region():
@@ -92,8 +92,8 @@ def test_earth_mag_02m_default_registration():
     npt.assert_allclose(data.coords["lat"].data.max(), 4.983333333)
     npt.assert_allclose(data.coords["lon"].data.min(), -9.98333333)
     npt.assert_allclose(data.coords["lon"].data.max(), -9.01666667)
-    npt.assert_allclose(data.min(), -231)
-    npt.assert_allclose(data.max(), 131.79999)
+    npt.assert_allclose(data.min(), -231.0, atol=0.2)
+    npt.assert_allclose(data.max(), 131.8, atol=0.2)
 
 
 def test_earth_mag4km_01d():
@@ -109,8 +109,8 @@ def test_earth_mag4km_01d():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -799.19995)
-    npt.assert_allclose(data.max(), 3226.4)
+    npt.assert_allclose(data.min(), -799.2, atol=0.2)
+    npt.assert_allclose(data.max(), 3226.4, atol=0.2)
 
 
 def test_earth_mag4km_01d_with_region():
@@ -126,8 +126,8 @@ def test_earth_mag4km_01d_with_region():
     assert data.shape == (11, 21)
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), -153.19995, rtol=1e-5)
-    npt.assert_allclose(data.max(), 113.59985, rtol=1e-5)
+    npt.assert_allclose(data.min(), -153.2, atol=0.2)
+    npt.assert_allclose(data.max(), 113.6, atol=0.2)
 
 
 def test_earth_mag4km_02m_default_registration():
@@ -146,8 +146,8 @@ def test_earth_mag4km_02m_default_registration():
     npt.assert_allclose(data.coords["lat"].data.max(), 5.98333333)
     npt.assert_allclose(data.coords["lon"].data.min(), -114.98333333)
     npt.assert_allclose(data.coords["lon"].data.max(), -112.01666667)
-    npt.assert_allclose(data.min(), -132.80005, rtol=1e-5)
-    npt.assert_allclose(data.max(), 79.59985, rtol=1e-5)
+    npt.assert_allclose(data.min(), -132.8, atol=0.2)
+    npt.assert_allclose(data.max(), 79.6, atol=0.2)
 
 
 def test_earth_mag_01d_wdmam():
@@ -164,8 +164,8 @@ def test_earth_mag_01d_wdmam():
     assert data.shape == (181, 361)
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -773.5)
-    npt.assert_allclose(data.max(), 1751.3)
+    npt.assert_allclose(data.min(), -794.0, atol=0.2)
+    npt.assert_allclose(data.max(), 2169.8, atol=0.2)
 
 
 def test_earth_mag_01d_wdmam_with_region():
@@ -181,8 +181,8 @@ def test_earth_mag_01d_wdmam_with_region():
     assert data.shape == (11, 21)
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), -103.900024, rtol=1e-5)
-    npt.assert_allclose(data.max(), 102.19995, rtol=1e-5)
+    npt.assert_allclose(data.min(), -145.6, atol=0.2)
+    npt.assert_allclose(data.max(), 107.6, atol=0.2)
 
 
 def test_earth_mag_03m_wdmam_with_region():
@@ -198,8 +198,8 @@ def test_earth_mag_03m_wdmam_with_region():
     assert data.lat.max() == -58
     assert data.lon.min() == 10
     assert data.lon.max() == 13
-    npt.assert_allclose(data.min(), -639.7001)
-    npt.assert_allclose(data.max(), 629.6)
+    npt.assert_allclose(data.min(), -790.2, atol=0.2)
+    npt.assert_allclose(data.max(), 528.0, atol=0.2)
 
 
 def test_earth_mag_03m_wdmam_without_region():

pygmt/tests/test_datasets_earth_relief.py

@@ -37,8 +37,8 @@ def test_earth_relief_01d_igpp_synbath(data_source):
     assert data.shape == (181, 361)
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -8600.5)
-    npt.assert_allclose(data.max(), 5559.0)
+    npt.assert_allclose(data.min(), -8600.5, atol=0.5)
+    npt.assert_allclose(data.max(), 5559.0, atol=0.5)
 
 
 @pytest.mark.parametrize("data_source", ["gebco", "gebcosi"])
@@ -48,12 +48,16 @@ def test_earth_relief_01d_gebco(data_source):
     data.
     """
     data = load_earth_relief(resolution="01d", data_source=data_source)
+    assert data.attrs["units"] == "meters"
+    assert data.attrs["long_name"] == "Earth elevation relative to the geoid"
+    assert data.attrs["vertical_datum"] == "EGM96"
+    assert data.attrs["horizontal_datum"] == "WGS84"
     assert data.shape == (181, 361)
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -8598)
-    npt.assert_allclose(data.max(), 5559.0)
+    npt.assert_allclose(data.min(), -8597.0, atol=1.0)
+    npt.assert_allclose(data.max(), 5559.0, atol=1.0)
 
 
 def test_earth_relief_01d_with_region_srtm():
@@ -69,8 +73,8 @@ def test_earth_relief_01d_with_region_srtm():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), -5154)
-    npt.assert_allclose(data.max(), 805.5)
+    npt.assert_allclose(data.min(), -5151, atol=0.5)
+    npt.assert_allclose(data.max(), 805.5, atol=0.5)
 
 
 def test_earth_relief_01d_with_region_gebco():
@@ -86,8 +90,8 @@ def test_earth_relief_01d_with_region_gebco():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), -5146)
-    npt.assert_allclose(data.max(), 806)
+    npt.assert_allclose(data.min(), -5151.0, atol=1.0)
+    npt.assert_allclose(data.max(), 806.0, atol=1.0)
 
 
 def test_earth_relief_30m():
@@ -99,8 +103,8 @@ def test_earth_relief_30m():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 90.5, 0.5))
     npt.assert_allclose(data.lon, np.arange(-180, 180.5, 0.5))
-    npt.assert_allclose(data.min(), -9454.5)
-    npt.assert_allclose(data.max(), 5887.5)
+    npt.assert_allclose(data.min(), -9454.5, atol=0.5)
+    npt.assert_allclose(data.max(), 5887.5, atol=0.5)
 
 
 def test_earth_gebcosi_15m_with_region():
@@ -117,8 +121,8 @@ def test_earth_gebcosi_15m_with_region():
     assert data.gmt.registration == 1
     npt.assert_allclose(data.lat, np.arange(-87.875, -84, 0.25))
     npt.assert_allclose(data.lon, np.arange(85.125, 87, 0.25))
-    npt.assert_allclose(data.min(), -531)
-    npt.assert_allclose(data.max(), 474)
+    npt.assert_allclose(data.min(), -531, atol=1.0)
+    npt.assert_allclose(data.max(), 474, atol=1.0)
 
 
 def test_earth_relief_30s_synbath():
@@ -132,8 +136,8 @@ def test_earth_relief_30s_synbath():
         data_source="synbath",
     )
     assert data.shape == (60, 120)
-    npt.assert_allclose(data.min(), -3552.5)
-    npt.assert_allclose(data.max(), -2154)
+    npt.assert_allclose(data.min(), -3552.5, atol=0.5)
+    npt.assert_allclose(data.max(), -2257.5, atol=0.5)
 
 
 def test_earth_relief_01m_without_region():
@@ -244,8 +248,8 @@ def test_earth_relief_15s_default_registration():
     npt.assert_allclose(data.coords["lat"].data.max(), 4.997917)
     npt.assert_allclose(data.coords["lon"].data.min(), -9.997917)
     npt.assert_allclose(data.coords["lon"].data.max(), -9.502083)
-    npt.assert_allclose(data.min(), -3897)
-    npt.assert_allclose(data.max(), -74)
+    npt.assert_allclose(data.min(), -3897, atol=0.5)
+    npt.assert_allclose(data.max(), -71, atol=0.5)
 
 
 def test_earth_relief_03s_default_registration():
@@ -260,5 +264,5 @@ def test_earth_relief_03s_default_registration():
     npt.assert_allclose(data.coords["lat"].data.max(), 5)
     npt.assert_allclose(data.coords["lon"].data.min(), -10)
     npt.assert_allclose(data.coords["lon"].data.max(), -9.8)
-    npt.assert_allclose(data.min(), -2069.996)
-    npt.assert_allclose(data.max(), -924.0801)
+    npt.assert_allclose(data.min(), -2070.0, atol=0.5)
+    npt.assert_allclose(data.max(), -924.5, atol=0.5)

pygmt/tests/test_datasets_earth_vertical_gravity_gradient.py

@@ -42,8 +42,8 @@ def test_earth_vertical_gravity_gradient_01d():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-90, 91, 1))
     npt.assert_allclose(data.lon, np.arange(-180, 181, 1))
-    npt.assert_allclose(data.min(), -136.34375)
-    npt.assert_allclose(data.max(), 104.59375)
+    npt.assert_allclose(data.min(), -137.125, atol=1 / 32)
+    npt.assert_allclose(data.max(), 104.59375, atol=1 / 32)
     assert data[1, 1].isnull()
 
 
@@ -58,8 +58,8 @@ def test_earth_vertical_gravity_gradient_01d_with_region():
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))
     npt.assert_allclose(data.lon, np.arange(-10, 11, 1))
-    npt.assert_allclose(data.min(), -16.34375)
-    npt.assert_allclose(data.max(), 19.78125)
+    npt.assert_allclose(data.min(), -15.6875, atol=1 / 32)
+    npt.assert_allclose(data.max(), 19.875, atol=1 / 32)
 
 
 def test_earth_vertical_gravity_gradient_01m_without_region():
@@ -95,5 +95,5 @@ def test_earth_vertical_gravity_gradient_01m_default_registration():
     npt.assert_allclose(data.coords["lat"].data.max(), 4.991666666)
     npt.assert_allclose(data.coords["lon"].data.min(), -9.99166666)
     npt.assert_allclose(data.coords["lon"].data.max(), -9.00833333)
-    npt.assert_allclose(data.min(), -40.25)
-    npt.assert_allclose(data.max(), 81.75)
+    npt.assert_allclose(data.min(), -37.5625, atol=1 / 32)
+    npt.assert_allclose(data.max(), 82.59375, atol=1 / 32)