function zpkTF = SR_TF_parametric(SRphi, zero_F)

% circle of poles parameters
circleEdge = -zero_F;
circleCenter = -zero_F/2-4;
circleRadius = circleCenter - circleEdge;

% defin position of pole along the circle
pole_circle = circleCenter + circleRadius*exp(1j*(pi-SRphi));

% convert into pole frequency and quality factor
pole_F = abs(pole_circle)
pole_Q = pole_F/abs(real(pole_circle))/2

% convert into zpk model used in matlab
zpkTF = virgo2zpk(zero_F, 0, pole_F, pole_Q, 1, 1)

