Merge pull request #1372 from knutfrode/dev

Including now also Antarctica coastline (as well as ice sheet-water-b…
OpenDrift · Jul 26, 2024 · 8022e94 · 8022e94
2 parents 84786ea + 5f20104
commit 8022e94
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/opendrift/readers/basereader/__init__.py b/opendrift/readers/basereader/__init__.py
@@ -268,7 +268,7 @@ def plot(self, variable=None, vmin=None, vmax=None, time=None,
             lscale = s.scale_from_extent([lonmin, lonmax, latmin, latmax])
 
         # GSHHS coastlines
-        f = cfeature.GSHHSFeature(scale=lscale, levels=[1])
+        f = cfeature.GSHHSFeature(scale=lscale, levels=[1,5,6])
         f._geometries_cache = {}
         ax.add_geometries(
             #f.intersecting_geometries([lonmin, lonmax, latmin, latmax]),

diff --git a/opendrift/readers/reader_global_landmask.py b/opendrift/readers/reader_global_landmask.py
@@ -108,7 +108,7 @@ def show_landmask_roaring(roaring):
     if fast:
         show_landmask_roaring(get_mask())
     else:
-        land = LandmaskFeature(scale=lscale, facecolor=land_color, globe=globe)
+        land = LandmaskFeature(scale=lscale, facecolor=land_color, globe=globe, levels=[1,5,6])
 
         ax.add_feature(land, zorder=2,
                        facecolor=land_color,

diff --git a/opendrift/readers/reader_schism_native.py b/opendrift/readers/reader_schism_native.py
@@ -770,7 +770,7 @@ def plot_mesh(self, variable=None, vmin=None, vmax=None,
             #              resolution='c', projection='cyl')
 
         # GSHHS coastlines
-        f = cfeature.GSHHSFeature(scale=lscale, levels=[1],
+        f = cfeature.GSHHSFeature(scale=lscale, levels=[1,5,6],
                                   facecolor=cfeature.COLORS['land'])
         ax.add_geometries(
             f.intersecting_geometries([lonmin, lonmax, latmin, latmax]),