code: def main():
      # URL of NASA GIBS
      url = (
          "https://gibs.earthdata.nasa.gov/wmts/epsg4326/best/wmts.cgi"
      )
      wmts = WebMapTileService(url)
  
      # Layers for MODIS true color and snow RGB
      layers = [
          "MODIS_Terra_SurfaceReflectance_Bands143",
          "MODIS_Terra_CorrectedReflectance_Bands367",
      ]
  
      date_str = "2016-02-05"
  
      # Plot setup
      plot_crs = ccrs.Mercator()
      geodetic_crs = ccrs.Geodetic()
      x0, y0 = plot_crs.transform_point(4.6, 43.1, geodetic_crs)
      x1, y1 = plot_crs.transform_point(11.0, 47.4, geodetic_crs)
      ysize = 8
      xsize = 2 * ysize * (x1 - x0) / (y1 - y0)
      fig = plt.figure(figsize=(xsize, ysize), dpi=100)
  
      for layer, offset in zip(layers, [0, 0.5]):
          ax = fig.add_axes([offset, 0, 0.5, 1], projection=plot_crs)
          ax.set_xlim((x0, x1))
          ax.set_ylim((y0, y1))
          ax.add_wmts(wmts, layer, wmts_kwargs={"time": date_str})
          txt = ax.text(
              4.7,
              43.2,
              wmts[layer].title,
              fontsize=18,
              color="wheat",
              transform=geodetic_crs,
          )
          txt.set_path_effects(
              [patheffects.withStroke(linewidth=5, foreground="black")]
          )
      plt.show()
  
imports: from matplotlib import patheffects
  import matplotlib.pyplot
  from owslib.wmts import WebMapTileService
  import cartopy.crs
label_desc: Use image node_id 1 for calling a global function with example usage: main() without return types