Dear Sander,

I am sure that the reply to my question will be something very obvious, which right now eludes me.

For reasons I do not wish to bother you with, I would like to output the scanline from the operational S5P product. Harp merges scanline and ground_pixel as

**scan_subindex** *description* the scanline and pixel dimensions are collapsed into a temporal dimension; the index of the pixel within the scanline is computed as the index on the temporal dimension modulo the number of scanlines

and indeed, I can turn the ingested product index into the ground pixel with modulo, for e.g.

In [16]:product.index.data % 450

Out[16]: array([14, 15, 13, 14, 15, 16, 12, 13, 14, 15, 16, 12, 13, 14, 15, 16, 12,

13, 14, 15, 16, 12, 13, 14, 15, 16, 13, 14, 15], dtype=int32)

In [19]: product.scan_subindex.data

Out[19]: array([14, 15, 13, 14, 15, 16, 12, 13, 14, 15, 16, 12, 13, 14, 15, 16, 12,

13, 14, 15, 16, 12, 13, 14, 15, 16, 13, 14, 15], dtype=int16)

Is there a similar way to find out the original scanline? these should be just a few consecutive scanlines since they refer to a very small part of the orbit [some 20km around a location]

Many thanks,

MariLiza