title: Geopandas plot shapefile on xarray with same legend
body: <p>I'm trying to create some maps of precipitation data (xarray) with a shapefile of the region of interest on top. However, when Python plots the figures, I get two seperate figures:</p>
  <p><a href="https://i.stack.imgur.com/HoWeE.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/HoWeE.png" alt="The precipitation data" /></a></p>
  <p><a href="https://i.stack.imgur.com/aInww.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/aInww.png" alt="The region of interest" /></a></p>
  <p>When I open the data in QGIS they do appear on top of each other, so the coordinate systems do check out. Then I have an additional bonus question: I have to create multiple precipitation maps,  on for a visual analysis it would be ideal if I could have the same legend (thus the same min/max for the colorbar) for each map. Anyone an idea how to proceed further?</p>
  <p>My code so far:</p>
  <pre class="lang-py prettyprint-override"><code>def chirps_to_map(input1, input2, title):
  
      projection = input1 + input2
  
      plt.figure(figsize=(9, 9))
  
      projection['pr'].plot()
  
      watershed.plot()
  
      plt.title(title)
  
      plt.show()
  
      plt.close()
  
      projection.to_netcdf(str(path)+str(title)+&quot;.nc&quot;)
  
      return projection
  
  </code></pre>