modisaqua.py

import xarray as xr
# https://oceandata.sci.gsfc.nasa.gov/ob/getfile/AQUA_MODIS.20230606T215500.L2.SST.NRT.nc

with xr.open_dataset('/Users/jklymak/Downloads/AQUA_MODIS.20230606T215500.L2.OC.NRT.nc',group='geophysical_data') as ds, xr.open_dataset('/Users/jklymak/Downloads/AQUA_MODIS.20230606T215500.L2.OC.NRT.nc',group='navigation_data') as nav:
    print(nav)
    ds['lon'] = (('number_of_lines', 'pixels_per_line'), nav['longitude'].values)
    ds['lat'] = (('number_of_lines', 'pixels_per_line'), nav['latitude'].values)
    for j in range(0, 2030):
        if np.any(~np.isfinite(ds.lat[j, :])):
            print(j)
    ds =  ds.sel(number_of_lines=slice(0, 2019))