Add an inline example for load_earth_mask

pygmt/datasets/earth_mask.py

@@ -57,6 +57,19 @@ def load_earth_mask(resolution="01d", region=None, registration=None):
         - 3: Islands in lakes in the land areas
         - 4: Smaller lakes in islands that are found within lakes
           inside the land area
+
+    Examples
+    --------
+
+    >>> from pygmt.datasets import load_earth_mask
+    >>> # load the default grid (gridline-registered 1 arc-degree grid)
+    >>> grid = load_earth_mask()
+    >>> # Location (120E, 50N) is in land area (1)
+    >>> grid.sel(lon=120, lat=50).values
+    array(1, dtype=int8)
+    >>> # Location (170E, 50N) is in oceanic area (0)
+    >>> grid.sel(lon=170, lat=50).values
+    array(0, dtype=int8)
     """
     grid = _load_remote_dataset(
         dataset_name="earth_mask",