new functions: "zsquare" and "z2fold"

similar to "epfold", these new functions allow to calculate the Z^2_m statistics for a test period to search for pulsations over a range of periods
minor bugfix for get_sron_colors: allow arrays as input, if "cont" qualifier is given