myr=getresult("ONA_1","input 1/mode 1/group delay");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y,"frequency (THz)","Group Delay (s)","");
holdon;

myr=getresult("ONA_2","input 1/mode 1/group delay");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y,"frequency (THz)","Group Delay (s)","");


myr=getresult("ONA_3","input 1/mode 1/group delay");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y,"frequency (THz)","Group Delay (s)","");
legend("ONA_1 scattering analysis","ONA_2 impulse reponse, no frac. delay,", "ONA_3 impulse reponse, with frac. delay");  
holdoff;

myr=getresult("ONA_4","input 1/mode 1/peak/free spectral range");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y/1e12,"frequency (THz)","free spectral range (THz)","");
holdon;

myr=getresult("ONA_5","input 1/mode 1/peak/free spectral range");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y/1e12,"frequency (THz)","free spectral range (THz)","");


myr=getresult("ONA_6","input 1/mode 1/peak/free spectral range");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y/1e12,"frequency (THz)","free spectral range (THz)","");
legend("ONA_4 scattering analysis","ONA_5 impulse reponse, no frac. delay,", "ONA_6: impulse reponse, with frac. delay");  
holdoff;