I want to plot contourf inside a country shape. So I use the following code:
def get_geometries(*readers): geoms = [] for reader in readers: for r