Xarray with zarr and kerchunk

import earthaccess
import fsspec
import xarray as xr

def load_data():
    src = "earthaccess_data/20020601090000-JPL-L4_GHRSST-SSTfnd-MUR-GLOB-v02.0-fv04.1.json"
    variable = ["analysed_sst"]
    earthaccess.login()
    s3_fs = earthaccess.get_s3fs_session(daac="PODAAC")
    storage_options = s3_fs.storage_options.copy()
    fsspec_caching = {
        "cache_type": "none",
    }
    fs = fsspec.filesystem("reference", fo=src, **fsspec_caching)
    m = fs.get_mapper("")
    da = xr.open_dataset(
        m,
        engine="kerchunk",
        chunks={},
        storage_options=storage_options,
    )[variable]
    return da.load()

if __name__ == "__main__":
    da = load_data()