I=sqrt(-1);
timeWindow=getnamed("","time window");
width=getnamed("GAUSS_1","width");
inputPower = getnamed("GAUSS_1","power");
fc=getnamed("GAUSS_1","frequency");
chirp=getnamed("GAUSS_1","chirp");


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

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


#plot Input Intensity	

myr=getresult("OOSC_1","mode 1/signal");
x=myr.getparameter(myr.getparameter);
y=myr.getattribute(myr.getattribute);

tc= ( x(1) + x(length(x)) ) /2.0;
plot((x-tc)/T0g,abs(y)^2/inputPower);
holdon;

myr=getresult("OOSC_2","mode 1/signal");
x=myr.getparameter(myr.getparameter);
y=myr.getattribute(myr.getattribute);

# plot Intensity IC 
tc= ( x(1) + x(length(x)) ) /2.0;
plot((x-tc)/T0g,abs(y)^2/inputPower);

#plot Intensity Agrawal
#plot(xmatlab,abs(ymatlab)^2);

setplot("title","Gaussian");
setplot("x min",-10);
setplot("x max",10);
setplot("y max",1.05);
setplot("x label","(t-t_c)/T_0");
setplot("y label", "Pulse Power Shape");
legend("a) Input Pulse, z/L_d=0","b) INTERCONNECT, z/L_d="+num2str(normLength),"c) Agrawal Fig 3.1, z/L_d="+num2str(normLength) );
holdoff;

#plot Phase
#IC
#tc= ( x(1) + x(length(x)) ) /2.0;
#plot((x-tc)/T0g,angle(y));
#holdon;
##Agrawal
#plot(xmatlab,angle(ymatlab));

#setplot("title","Gaussian, Chirp="+num2str(chirp));
#setplot("x min",-10);
#setplot("x max",10);
#setplot("y max",1.05);
#setplot("x label","(t-t_c)/T_0");
#setplot("y label", "Pulse Phasee");
#holdoff;