Extend ancil_create_shapefile.py to support rotating coordinates to a domain

[arjclark]

✨ Feature Request

For regional usage, it is desirable to specify the coordinates of a shapefile once in the json but be able to use it for differing rotated pole domains. Rotated pole coordinates can be provided so the points in the json just need rotating the polygon to the new coords before generating the shapefile from them (shapefiles as we have them don't know about coordinate systems so if you tell it it's at X,Y that's what's used when it's applied).

Requirements

• Add a --target-lsm option to ancil_create_shapefile.py to take in a land sea mask for the rotated pole domain

• Extract the rotated pole coordinate details from there

• Take the points defined in the json file and use cartopy to rotate them to the new rotated pole coordinate system

• Generate the shapefile as before but using the rotated points

• Need to make sure it's clear in documentation that the routine assumes the source json is on a standard (unrotated) lat-lon grid

• Logging/outputs should make it clear when a rotation occurs

• Don't do anything to the points should a non-rotated pole lsm get fed in (realistic possibility given this will get plumbed into workflows) - some sort of warning/error should be raised in this case

Context

This functionality will be used in the Regional Ancillary Suite.