import cartopy.crs as ccrs

gmp_2010 = ccrs.RotatedPole(pole_longitude=72.21, pole_latitude=80.08)

import matplotlib.pyplot as plt

plt.figure(figsize=(14, 14))
ax = plt.axes(projection=gmp_2010)
ax.coastlines()
ax.gridlines()
plt.show()

import matplotlib.pyplot as plt

plt.figure(figsize=(14, 14))

ax1 = plt.subplot(1, 2, 1, projection=ccrs.NorthPolarStereo())
ax1.coastlines()
ax1.gridlines(gmp_2010, color='red')


ax1 = plt.subplot(1, 2, 2, projection=ccrs.NorthPolarStereo())
ax1.coastlines()
ax1.gridlines()


plt.show()