Added option to calculate Gehrels (1986) upper confidence limit. Small improvements.
readx1d() now makes use of astropy.table.Table objects.