Building a numpy array of geometries had an undesired behaviour when passing on only one (shapely) MultiPolygon. This happened for the first time when running Gibraltar, which has one data unit composed of one MultiPolygon with two parts.
This MR fixes this by using lists instead of numpy arrays.