changed ryan's skypos_uni_vol() method to the more general pareto_hrss() method which draws from a pareto distribution in hrss uniformly over the sky, keeping only those events for which network_snr>min_snr