#run;  runanalysis;

m="solar_generation";

# total absorption
P=getresult(m,"Pabs_total");
plot(c/P.f*1e6,P.Pabs_total,"wavelength (um)","fraction of sourcepower","Absorption");

# spatial absorption
P=getresult(m,"Pabs");
data = P.Pabs;
data = pinch(data,4,find(P.f,c/1.1e-6));
data = pinch(data,2,find(P.y,0));
image(P.x*1e6,P.z*1e6,data,"x (um)","z (um)","Absorption at 1.1um");
data = P.Pabs;
data = pinch(data,4,find(P.f,c/0.4e-6));
data = pinch(data,2,find(P.y,0));
image(P.x*1e6,P.z*1e6,data,"x (um)","z (um)","Absorption at 0.4um");

# generation rate from simulation
G=getresult(m,"G");
image(G.x*1e6,G.z*1e6,pinch(G.G,2,find(G.y,0)),"x (um)","z (um)","Generation rate");

G=getresult(m,"G_export");
image(G.x*1e6,G.z*1e6,pinch(pinch(G.G,2,12)),"x (um)","z (um)","Generation rate exported");

Jsc=getresult(m,"Jsc");
?num2str(Jsc/10)+" mA/cm^2";