def METHOD_NAME(self, domain, organization, desired_nodes):
    icesat2.init(domain, organization=organization, desired_nodes=desired_nodes, bypass_dns=True)
    resource = "ATL03_20181017222812_02950102_005_01.h5"
    region = sliderule.toregion(os.path.join(TESTDIR, "data/grandmesa.geojson"))
    parms = {
        "poly": region["poly"],
        "srt": icesat2.SRT_LAND,
        "cnf": icesat2.CNF_SURFACE_HIGH,
        "ats": 10.0,
        "cnt": 10,
        "len": 40.0,
        "res": 20.0,
        "output": { "path": "testfile4.parquet", "format": "parquet", "open_on_complete": True } }
    gdf = icesat2.atl03sp(parms, resources=[resource])
    os.remove("testfile4.parquet")
    assert len(gdf) == 20642
    assert gdf.index.values.min() == numpy.datetime64('2018-10-17T22:31:17.349347328')
    assert gdf.index.values.max() == numpy.datetime64('2018-10-17T22:31:19.582347520')