OPS C1442068492-GES_DISC (S5P_L2__CLOUD_)

[podaac-cicd]

New association between l2ss-py and C1442068492-GES_DISC found in OPS.
Beginning verification of collection.

[github-actions]

Test Results for C1442068492-GES_DISC

2 tests   1 ✅  1m 48s ⏱️
1 suites  0 💤
1 files    1 ❌

For more details on these failures, see this check.

Results for commit 7349a04.

♻️ This comment has been updated with latest results.

[github-actions]

Action Needed
Tests have failed or encountered an unknown error. Open the status checks to view the logs and review the failure.

[ocsmit]

Leaving a comment here for reference.

The failure here is caused by the autotester is subsetting a portion of the granule which just simply does not have any data.

I.e. it grabs a science variable (cloud_base_height), masks the lat and lon variables where the science variable is NaN, and then checks if there is any coordinate data for the non NaN pixels. Problem is is that there is no science data in the subsetted portion being used so when mask is applied to the lat/lon vars they end up with completely empty arrays. (see image)

I verified with the autotester locally that when running with a larger subset BBox that it passes fully.

I also noticed that there is a break statement in the get_bounding_box method that stops the method from getting the full bounding box of the entire granule in the case that there are multiple polygons (such as with this collection). If this is intended behavior to get only a portion of the bounding box even before making a smaller bounding box with create_smaller_bounding_box, then it should be made clearer. Removing this break statement (L218) and getting full spatial extent bbox for the granule before calling create_smaller_bounding_box also results in a passing spatial subset test.

l2ss-py-autotest/tests/verify_collection.py

Lines 212 to 218 in 4d065af

	if polygons:
	for polygon in polygons:
	points = polygon['Boundary']['Points']
	for point in points:
	longitude_list.append(point.get('Longitude'))
	latitude_list.append(point.get('Latitude'))
	break