def METHOD_NAME():
    region = db.get_timestep("test_tipsy/tiny.000640").load_region(pynbody.filt.Sphere(2000,[1000,1000,1000]))
    assert (region['iord']==[ 9980437, 10570437, 10630437, 10640437, 10770437, 10890437,
       10900437, 10960437, 11030437, 11090437, 11480437, 11490437,
       11550437, 11740437, 12270437, 12590437, 12600437, 12920437,
       13380437, 13710437]).all()