#read Agrawal
I=sqrt(-1);
mydata =readdata("AgrawalFig4_21a.csv");
sz = size(mydata);
numsamples = sz(1);
xmatlab = mydata(1:numsamples,1);
ymatlab = mydata(1:numsamples,2); #+I*mydata(1:numsamples,3);

timeWindow=getnamed("","time window");
width=getnamed("GAUSS_1","width");
fc=getnamed("GAUSS_1","frequency");
chirp=getnamed("GAUSS_1","chirp");
P0in=getnamed("GAUSS_1","power");

Tfwhm=width*timeWindow;
T0g=Tfwhm/(2*sqrt(log(2)));

#read IC
myr=getresult("OOSC_4","mode 1/signal");
x=myr.getparameter(myr.getparameter);
y=myr.getattribute(myr.getattribute);
tc=( x(1) + x(length(x)) ) /2;
P0out = max(abs(y)^2);

#plot IC Intensity
plot((x-tc)/T0g,abs(y)^2/P0in);
setplot("title","Gaussian, Chirp="+num2str(chirp));
setplot("x min",-5);
setplot("x max",5);
setplot("y max",0.4);
setplot("x label","(t-t_c)/T_0");
setplot("y label", "Normalized Pulse Power");
holdon;

#plot Agrawal matlab Intensity
plot(xmatlab,abs(ymatlab));
setplot("title","Gaussian, Chirp="+num2str(chirp));
setplot("x min",-5);
setplot("x max",5);
setplot("y max",0.4);
setplot("x label","(t-t_c)/T_0");
setplot("y label", "Normalized Pulse Power");
legend("a) INTERCONNECT","b) Ref.[1] Fig 4.21a");
holdoff;

#plot IC Phase
#plot((x-tc)/T0g,angle(y));
#setplot("title","Gaussian, Chirp="+num2str(chirp));
#setplot("x min",-5);
#setplot("x max",5);
#setplot("y max",1.05);
#setplot("x label","(t-t_c)/T_0");
#setplot("y label", "Phase");
#holdon;

#plot Agrawal matlab Phase
#plot(xmatlab,angle(ymatlab));
#setplot("title","Gaussian, Chirp="+num2str(chirp));
#setplot("x min",-5);
#setplot("x max",5);
#setplot("y max",1.05);
#setplot("x label","(t-t_c)/T_0");
#setplot("y label", "Phase");
#holdoff;