Hi folks, I saw that ECMWF have developed this: https://github.com/ecmwf/polytope  It might be a useful package to watch in terms of managing complex subsetting to complex shapes.