Added threshold and reldiff_cube to cubist_numpy