prompt: We are currently solving the following issue within our repository. Here is the issue text:
  --- BEGIN ISSUE ---
  Remove "cfeature.LAND" from "Feature Creation" example ?
  In the nice example provided on the [Feature Creation page](https://scitools.org.uk/cartopy/docs/latest/gallery/feature_creation.html), you could probably remove the useless line:
  B) `ax.add_feature(cfeature.LAND)`
  because you have already called
  A) `ax.stock_img()`
  
  As a matter of fact, B) being called after A), it should theoretically be plotted above A). I see you are using a default `zorder=-1` for LAND, so it is probably below everything else (did not know you could use negative zorder values)
  ```
  >>> cfeature.LAND.kwargs
  {'edgecolor': 'face', 'zorder': -1, 'facecolor': array([0.9375  , 0.9375  , 0.859375])}
  ```
  If I use `ax.add_feature(cfeature.LAND, zorder=10)`, I do get land with a uniform color above the stock image on land (and the stock ocean elsewhere)
  
  --- END ISSUE ---
  
  Below are some code segments, each from a relevant file. One or more of these files may contain bugs.
  --- BEGIN FILES ---
  Path: `examples/lines_and_polygons/feature_creation.py`
  Content:
  ```
  1 """
  2 Feature Creation
  3 ----------------
  4 
  5 This example manually instantiates a
  6 :class:`cartopy.feature.NaturalEarthFeature` to access administrative
  7 boundaries (states and provinces).
  8 
  9 Note that this example is intended to illustrate the ability to construct
  10 Natural Earth features that cartopy does not necessarily know about
  11 *a priori*.
  12 In this instance however, it would be possible to make use of the
  13 pre-defined :data:`cartopy.feature.STATES` constant.
  14 
  15 """
  16 from matplotlib.offsetbox import AnchoredText
  17 import matplotlib.pyplot as plt
  18 
  19 import cartopy.crs as ccrs
  20 import cartopy.feature as cfeature
  21 
  22 
  23 def main():
  24     fig = plt.figure()
  25     ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
  26     ax.set_extent([80, 170, -45, 30], crs=ccrs.PlateCarree())
  27 
  28     # Put a background image on for nice sea rendering.
  29     ax.stock_img()
  30 
  31     # Create a feature for States/Admin 1 regions at 1:50m from Natural Earth
  32     states_provinces = cfeature.NaturalEarthFeature(
  33         category='cultural',
  34         name='admin_1_states_provinces_lines',
  35         scale='50m',
  36         facecolor='none')
  37 
  38     SOURCE = 'Natural Earth'
  39     LICENSE = 'public domain'
  40 
  41     ax.add_feature(cfeature.LAND)
  42     ax.add_feature(cfeature.COASTLINE)
  43     ax.add_feature(states_provinces, edgecolor='gray')
  44 
  45     # Add a text annotation for the license information to the
  46     # the bottom right corner.
  47     text = AnchoredText('\u00A9 {}; license: {}'
  48                         ''.format(SOURCE, LICENSE),
  49                         loc=4, prop={'size': 12}, frameon=True)
  50     ax.add_artist(text)
  51 
  52     plt.show()
  53 
  54 
  55 if __name__ == '__main__':
  56     main()
  57
  ```
  --- END FILES ---
  
  Please first localize the bug based on the issue statement, and then generate a patch according to the `git diff` format fenced by three backticks.
  
  Here is an example:
  
  ```diff
  diff --git a/examples/server_async.py b/examples/server_async.py
  --- a/examples/server_async.py
  +++ b/examples/server_async.py
  @@ -313,4 +313,4 @@
   
   
   if __name__ == "__main__":
  -    asyncio.run(run_async_server("."), debug=True)
  +    asyncio.run(run_async_server(), debug=True)
  diff --git a/examples/server_sync.py b/examples/server_sync.py
  --- a/examples/server_sync.py
  +++ b/examples/server_sync.py
  @@ -313,5 +313,5 @@
   
   
   if __name__ == "__main__":
  -    server = run_sync_server(".")
  +    server = run_sync_server()
       server.shutdown()
  
  ```
  
golden_diff: diff --git a/examples/lines_and_polygons/feature_creation.py b/examples/lines_and_polygons/feature_creation.py
  --- a/examples/lines_and_polygons/feature_creation.py
  +++ b/examples/lines_and_polygons/feature_creation.py
  @@ -28,7 +28,7 @@
       # Put a background image on for nice sea rendering.
       ax.stock_img()
   
  -    # Create a feature for States/Admin 1 regions at 1:50m from Natural Earth
  +    # Create a feature for States/Admin 1 regions at 1:50m from Natural Earth.
       states_provinces = cfeature.NaturalEarthFeature(
           category='cultural',
           name='admin_1_states_provinces_lines',
  @@ -38,9 +38,11 @@
       SOURCE = 'Natural Earth'
       LICENSE = 'public domain'
   
  -    ax.add_feature(cfeature.LAND)
  -    ax.add_feature(cfeature.COASTLINE)
  +    # Add our states feature.
       ax.add_feature(states_provinces, edgecolor='gray')
  +    # Add land feature, overriding the default negative zorder so it shows
  +    # above the background image.
  +    ax.add_feature(cfeature.LAND, zorder=1, edgecolor='k')
   
       # Add a text annotation for the license information to the
       # the bottom right corner.
  
verification_info: {"golden_diff": "diff --git a/examples/lines_and_polygons/feature_creation.py b/examples/lines_and_polygons/feature_creation.py\n--- a/examples/lines_and_polygons/feature_creation.py\n+++ b/examples/lines_and_polygons/feature_creation.py\n@@ -28,7 +28,7 @@\n     # Put a background image on for nice sea rendering.\n     ax.stock_img()\n \n-    # Create a feature for States/Admin 1 regions at 1:50m from Natural Earth\n+    # Create a feature for States/Admin 1 regions at 1:50m from Natural Earth.\n     states_provinces = cfeature.NaturalEarthFeature(\n         category='cultural',\n         name='admin_1_states_provinces_lines',\n@@ -38,9 +38,11 @@\n     SOURCE = 'Natural Earth'\n     LICENSE = 'public domain'\n \n-    ax.add_feature(cfeature.LAND)\n-    ax.add_feature(cfeature.COASTLINE)\n+    # Add our states feature.\n     ax.add_feature(states_provinces, edgecolor='gray')\n+    # Add land feature, overriding the default negative zorder so it shows\n+    # above the background image.\n+    ax.add_feature(cfeature.LAND, zorder=1, edgecolor='k')\n \n     # Add a text annotation for the license information to the\n     # the bottom right corner.\n", "issue": "Remove \"cfeature.LAND\" from \"Feature Creation\" example ?\nIn the nice example provided on the [Feature Creation page](https://scitools.org.uk/cartopy/docs/latest/gallery/feature_creation.html), you could probably remove the useless line:\r\nB) `ax.add_feature(cfeature.LAND)`\r\nbecause you have already called\r\nA) `ax.stock_img()`\r\n\r\nAs a matter of fact, B) being called after A), it should theoretically be plotted above A). I see you are using a default `zorder=-1` for LAND, so it is probably below everything else (did not know you could use negative zorder values)\r\n```\r\n>>> cfeature.LAND.kwargs\r\n{'edgecolor': 'face', 'zorder': -1, 'facecolor': array([0.9375  , 0.9375  , 0.859375])}\r\n```\r\nIf I use `ax.add_feature(cfeature.LAND, zorder=10)`, I do get land with a uniform color above the stock image on land (and the stock ocean elsewhere)\n", "before_files": [{"content": "\"\"\"\nFeature Creation\n----------------\n\nThis example manually instantiates a\n:class:`cartopy.feature.NaturalEarthFeature` to access administrative\nboundaries (states and provinces).\n\nNote that this example is intended to illustrate the ability to construct\nNatural Earth features that cartopy does not necessarily know about\n*a priori*.\nIn this instance however, it would be possible to make use of the\npre-defined :data:`cartopy.feature.STATES` constant.\n\n\"\"\"\nfrom matplotlib.offsetbox import AnchoredText\nimport matplotlib.pyplot as plt\n\nimport cartopy.crs as ccrs\nimport cartopy.feature as cfeature\n\n\ndef main():\n    fig = plt.figure()\n    ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())\n    ax.set_extent([80, 170, -45, 30], crs=ccrs.PlateCarree())\n\n    # Put a background image on for nice sea rendering.\n    ax.stock_img()\n\n    # Create a feature for States/Admin 1 regions at 1:50m from Natural Earth\n    states_provinces = cfeature.NaturalEarthFeature(\n        category='cultural',\n        name='admin_1_states_provinces_lines',\n        scale='50m',\n        facecolor='none')\n\n    SOURCE = 'Natural Earth'\n    LICENSE = 'public domain'\n\n    ax.add_feature(cfeature.LAND)\n    ax.add_feature(cfeature.COASTLINE)\n    ax.add_feature(states_provinces, edgecolor='gray')\n\n    # Add a text annotation for the license information to the\n    # the bottom right corner.\n    text = AnchoredText('\\u00A9 {}; license: {}'\n                        ''.format(SOURCE, LICENSE),\n                        loc=4, prop={'size': 12}, frameon=True)\n    ax.add_artist(text)\n\n    plt.show()\n\n\nif __name__ == '__main__':\n    main()\n", "path": "examples/lines_and_polygons/feature_creation.py"}], "after_files": [{"content": "\"\"\"\nFeature Creation\n----------------\n\nThis example manually instantiates a\n:class:`cartopy.feature.NaturalEarthFeature` to access administrative\nboundaries (states and provinces).\n\nNote that this example is intended to illustrate the ability to construct\nNatural Earth features that cartopy does not necessarily know about\n*a priori*.\nIn this instance however, it would be possible to make use of the\npre-defined :data:`cartopy.feature.STATES` constant.\n\n\"\"\"\nfrom matplotlib.offsetbox import AnchoredText\nimport matplotlib.pyplot as plt\n\nimport cartopy.crs as ccrs\nimport cartopy.feature as cfeature\n\n\ndef main():\n    fig = plt.figure()\n    ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())\n    ax.set_extent([80, 170, -45, 30], crs=ccrs.PlateCarree())\n\n    # Put a background image on for nice sea rendering.\n    ax.stock_img()\n\n    # Create a feature for States/Admin 1 regions at 1:50m from Natural Earth.\n    states_provinces = cfeature.NaturalEarthFeature(\n        category='cultural',\n        name='admin_1_states_provinces_lines',\n        scale='50m',\n        facecolor='none')\n\n    SOURCE = 'Natural Earth'\n    LICENSE = 'public domain'\n\n    # Add our states feature.\n    ax.add_feature(states_provinces, edgecolor='gray')\n    # Add land feature, overriding the default negative zorder so it shows\n    # above the background image.\n    ax.add_feature(cfeature.LAND, zorder=1, edgecolor='k')\n\n    # Add a text annotation for the license information to the\n    # the bottom right corner.\n    text = AnchoredText('\\u00A9 {}; license: {}'\n                        ''.format(SOURCE, LICENSE),\n                        loc=4, prop={'size': 12}, frameon=True)\n    ax.add_artist(text)\n\n    plt.show()\n\n\nif __name__ == '__main__':\n    main()\n", "path": "examples/lines_and_polygons/feature_creation.py"}]}