Always check if gridline is the default registration for 01d resolution

Author: seisman

pygmt/tests/test_datasets_earth_age.py

@@ -48,9 +48,7 @@ def test_earth_age_01d_with_region():
     """
     Test loading low-resolution earth age with 'region'.
     """
-    data = load_earth_age(
-        resolution="01d", region=[-10, 10, -5, 5], registration="gridline"
-    )
+    data = load_earth_age(resolution="01d", region=[-10, 10, -5, 5])
     assert data.shape == (11, 21)
     assert data.gmt.registration == 0
     npt.assert_allclose(data.lat, np.arange(-5, 6, 1))

pygmt/tests/test_datasets_earth_free_air_anomaly.py

@@ -31,12 +31,13 @@ def test_earth_faa_01d():
     """
     Test some properties of the free air anomaly 01d data.
     """
-    data = load_earth_free_air_anomaly(resolution="01d", registration="gridline")
+    data = load_earth_free_air_anomaly(resolution="01d")
     assert data.name == "free_air_anomaly"
     assert data.attrs["long_name"] == "IGPP Global Earth Free-Air Anomaly"
     assert data.attrs["units"] == "mGal"
     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(), -275.75)
@@ -47,10 +48,9 @@ def test_earth_faa_01d_with_region():
     """
     Test loading low-resolution earth free air anomaly with 'region'.
     """
-    data = load_earth_free_air_anomaly(
-        resolution="01d", region=[-10, 10, -5, 5], registration="gridline"
-    )
+    data = load_earth_free_air_anomaly(resolution="01d", region=[-10, 10, -5, 5])
     assert data.shape == (11, 21)
+    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)

pygmt/tests/test_datasets_earth_geoid.py

@@ -31,12 +31,13 @@ def test_earth_geoid_01d():
     """
     Test some properties of the earth geoid 01d data.
     """
-    data = load_earth_geoid(resolution="01d", registration="gridline")
+    data = load_earth_geoid(resolution="01d")
     assert data.name == "earth_geoid"
     assert data.attrs["units"] == "m"
     assert data.attrs["long_name"] == "EGM2008 Global Earth Geoid"
     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(), -106.45)
@@ -47,10 +48,9 @@ def test_earth_geoid_01d_with_region():
     """
     Test loading low-resolution earth geoid with 'region'.
     """
-    data = load_earth_geoid(
-        resolution="01d", region=[-10, 10, -5, 5], registration="gridline"
-    )
+    data = load_earth_geoid(resolution="01d", region=[-10, 10, -5, 5])
     assert data.shape == (11, 21)
+    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)

pygmt/tests/test_datasets_earth_magnetic_anomaly.py

@@ -31,12 +31,13 @@ def test_earth_mag_01d():
     """
     Test some properties of the magnetic anomaly 01d data.
     """
-    data = load_earth_magnetic_anomaly(resolution="01d", registration="gridline")
+    data = load_earth_magnetic_anomaly(resolution="01d")
     assert data.name == "magnetic_anomaly"
     assert data.attrs["long_name"] == "Earth magnetic anomaly"
     assert data.attrs["units"] == "nT"
     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(), -384)
@@ -47,10 +48,9 @@ def test_earth_mag_01d_with_region():
     """
     Test loading low-resolution earth magnetic anomaly with 'region'.
     """
-    data = load_earth_magnetic_anomaly(
-        resolution="01d", region=[-10, 10, -5, 5], registration="gridline"
-    )
+    data = load_earth_magnetic_anomaly(resolution="01d", region=[-10, 10, -5, 5])
     assert data.shape == (11, 21)
+    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)
@@ -81,14 +81,13 @@ def test_earth_mag4km_01d():
     """
     Test some properties of the magnetic anomaly 4km 01d data.
     """
-    data = load_earth_magnetic_anomaly(
-        resolution="01d", registration="gridline", mag4km=True
-    )
+    data = load_earth_magnetic_anomaly(resolution="01d", mag4km=True)
     assert data.name == "magnetic_anomaly"
     assert data.attrs["long_name"] == "Earth magnetic anomaly"
     assert data.attrs["units"] == "nT"
     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(), -799.19995)

pygmt/tests/test_datasets_earth_relief.py

@@ -47,10 +47,9 @@ def test_earth_relief_01d_gebco(data_source):
     Test some properties of the earth relief 01d data with GEBCO and GEBOCSI
     data.
     """
-    data = load_earth_relief(
-        resolution="01d", registration="gridline", data_source=data_source
-    )
+    data = load_earth_relief(resolution="01d", data_source=data_source)
     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)
@@ -64,10 +63,10 @@ def test_earth_relief_01d_with_region_srtm():
     data = load_earth_relief(
         resolution="01d",
         region=[-10, 10, -5, 5],
-        registration="gridline",
         data_source="igpp",
     )
     assert data.shape == (11, 21)
+    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)
@@ -81,10 +80,10 @@ def test_earth_relief_01d_with_region_gebco():
     data = load_earth_relief(
         resolution="01d",
         region=[-10, 10, -5, 5],
-        registration="gridline",
         data_source="gebco",
     )
     assert data.shape == (11, 21)
+    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)
@@ -95,8 +94,9 @@ def test_earth_relief_30m():
     """
     Test some properties of the earth relief 30m data.
     """
-    data = load_earth_relief(resolution="30m", registration="gridline")
+    data = load_earth_relief(resolution="30m")
     assert data.shape == (361, 721)
+    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)

pygmt/tests/test_datasets_earth_vertical_gravity_gradient.py

@@ -33,14 +33,13 @@ def test_earth_vertical_gravity_gradient_01d():
     """
     Test some properties of the earth vgg 01d data.
     """
-    data = load_earth_vertical_gravity_gradient(
-        resolution="01d", registration="gridline"
-    )
+    data = load_earth_vertical_gravity_gradient(resolution="01d")
     assert data.name == "earth_vgg"
     assert data.attrs["units"] == "Eotvos"
     assert data.attrs["long_name"] == "IGPP Global Earth Vertical Gravity Gradient"
     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(), -136.34375)
@@ -53,9 +52,10 @@ def test_earth_vertical_gravity_gradient_01d_with_region():
     Test loading low-resolution earth vgg with 'region'.
     """
     data = load_earth_vertical_gravity_gradient(
-        resolution="01d", region=[-10, 10, -5, 5], registration="gridline"
+        resolution="01d", region=[-10, 10, -5, 5]
     )
     assert data.shape == (11, 21)
+    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)