function phi = cycs(z)
% Phase angle in periods: phi = cyc(unwrap(angle(z)))

  phi = cyc(unwrap(angle(z)));
  
  return