clear;

#Run the multifrequency simulation and collect Ey at the source 
#and after propagation through the wire
switchtolayout;
select("source");
set("multifrequency mode calculation",1);
set("frequency points",15);

run;

Ey_multi_injected=getresult("Injected","Ey");
Ey_multi_injected=pinch(Ey_multi_injected);
Ey_multi_propagated=getresult("Propagated","Ey");
Ey_multi_propagated=pinch(Ey_multi_propagated);
t=getresult("Injected","t")*10e12;

plot(t,Ey_multi_injected,Ey_multi_propagated,"t [ps]","Ey","Multi frequency mode source(15 frequency points)");
legend("Ey time profile of the injected pulse","Ey time profile after the impulse propagated through the coated wire");

#Run single frequency simulation and collect Ey at the source 
#and after propagation through the wire
switchtolayout;
select("source");
set("multifrequency mode calculation",0);

run;

Ey_single_injected=getresult("Injected","Ey");
Ey_single_injected=pinch(Ey_single_injected);
Ey_single_propagated=getresult("Propagated","Ey");
Ey_single_propagated=pinch(Ey_single_propagated);
t=getresult("Injected","t")*10e12;

plot(t,Ey_single_injected,Ey_single_propagated,"t [ps]","Ey","Single frequency mode source(1 frequency point)");
legend("Ey time profile of the injected pulse","Ey time profile after the impulse propagated through the coated wire");