theta = linspace(0,60,7);
A = matrix(length(theta),200);

for (i = 1:length(theta)){

switchtolayout;
setnamed('source1','plane wave type','BFAST'); # enable BFAST
setnamed('source1','angle theta',theta(i));
run;
T = -transmission('T');
R = transmission('R');
A(i,1:200) = 1-R-T;
?i;
}

# revert to original
switchtolayout;
setnamed('source1','plane wave type','Bloch/periodic'); # disable BFAST
setnamed('source1','angle theta',0);

f = getdata('R','f');
image(c*1e6/f,theta,transpose(A),'wavelength (um)','theta','A');

f = getdata('R','f');
image(theta,c*1e6/f,(A),'wavelength (um)','theta','A');