Hello Marie,
If found the cause for this problem. OLCI RR Orbits can contain the north pole. And this wasn’t properly handled for all cases in our code (SNAP-759). This is fixed and will be released with the next update. Not sure when will be published, maybe in two week or so.
As a temporary workaround, you can do a pixel based subset first and afterwards a geo-region subset.
The pixel based subset might look like this:
from snappy import Rectangle
op1 = SubsetOp()
op1.setSourceProduct(p)
size = p.getSceneRasterSize()
op1.setRegion(Rectangle(0, 500, size.width, size,height))
pixelSubset = op1.getTargetProduct()
and in your second op you use pixelSubset
instead of p
op.setSourceProduct(pixelSubset)