Quickfix min elevation for new SRTM15+V2.1 earth relief grids (#401)

Upstream GMT has updated its earth relief grids to SRTM15+V2.1 to fill in NaN values (see GenericMappingTools/gmt#2920). This commit fixes the test failures due to incorrect minimum elevation values, and minor changes to the contour and grid plots.

Author: weiji14

pygmt/clib/session.py

@@ -111,7 +111,7 @@ class Session:
     ...             )
     ...             # Read the contents of the temp file before it's deleted.
     ...             print(fout.read().strip())
-    -180 180 -90 90 -8596 5559 1 1 361 181
+    -180 180 -90 90 -8592 5559 1 1 361 181
     """
 
     # The minimum version of GMT required
@@ -1218,15 +1218,15 @@ def virtualfile_from_grid(self, grid):
         >>> print(data.lat.values.min(), data.lat.values.max())
         -90.0 90.0
         >>> print(data.values.min(), data.values.max())
-        -8596.0 5559.0
+        -8592.0 5559.0
         >>> with Session() as ses:
         ...     with ses.virtualfile_from_grid(data) as fin:
         ...         # Send the output to a file so that we can read it
         ...         with GMTTempFile() as fout:
         ...             args = '{} -L0 -Cn ->{}'.format(fin, fout.name)
         ...             ses.call_module('grdinfo', args)
         ...             print(fout.read().strip())
-        -180 180 -90 90 -8596 5559 1 1 361 181
+        -180 180 -90 90 -8592 5559 1 1 361 181
         >>> # The output is: w e s n z0 z1 dx dy n_columns n_rows
 
         """

pygmt/tests/baseline/test_grdcontour.png

@@ -62,7 +62,7 @@ def test_earth_relief_60():
     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(), -8596)
+    npt.assert_allclose(data.min(), -8592)
     npt.assert_allclose(data.max(), 5559)
 
 
@@ -72,5 +72,5 @@ def test_earth_relief_30():
     assert data.shape == (361, 721)
     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(), -9458)
+    npt.assert_allclose(data.min(), -9460)
     npt.assert_allclose(data.max(), 5888)

pygmt/tests/baseline/test_grdcontour_file.png

@@ -59,13 +59,13 @@ def test_grdinfo():
     "Make sure grd info works as expected"
     grid = load_earth_relief()
     result = grdinfo(grid, L=0, C="n")
-    assert result.strip() == "-180 180 -90 90 -8596 5559 1 1 361 181"
+    assert result.strip() == "-180 180 -90 90 -8592 5559 1 1 361 181"
 
 
 def test_grdinfo_file():
     "Test grdinfo with file input"
     result = grdinfo("@earth_relief_60m", L=0, C="n")
-    assert result.strip() == "-180 180 -90 90 -8596 5559 1 1 361 181"
+    assert result.strip() == "-180 180 -90 90 -8592 5559 1 1 361 181"
 
 
 def test_grdinfo_fails():