# data
s11 = getresult("Optical Network Analyzer","input 1/mode 1/transmission");
S11 = getattribute(s11, "TE transmission");
s21 = getresult("Optical Network Analyzer","input 2/mode 1/transmission");
S21 = getattribute(s21, "TE transmission");
f   = getparameter(s11, "frequency");

# power
plot(c/f*1e6,abs(S11)^2,abs(S21)^2,"wavelength (microns)","Power");
legend("R","T");

# amplitude
plot(c/f*1e6,abs(S11),abs(S21),"wavelength (microns)","|S|");
legend("|S11|","|S21|");

# phase
S11 = pinch(S11);
S21 = pinch(S21);
phi0_21 = angle(S21(1))*(f-f(1))*2*pi/c*540e-6;
phi0_11 = angle(S11(1))*(f-f(1))*2*pi/c*500e-6;
plot(c/f*1e6,unwrap(angle(S11)-phi0_11)+phi0_11,unwrap(angle(S21)-phi0_21)+phi0_21,"wavelength (microns)","angle(S)");
legend("angle(S11)","angle(S21)");