mdsumner read_geoid()