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);

format short;
setplot("title","Time Window = "+num2str(timeWindow/1e-12)+" ps");
format long;
#setplot("x min",0);
#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) );
holdoff;