pk.bf.stats function

Calculate baseflow peak statistics