AbsorptionEnable=1;
setnamed("absorption","enabled",1);
switchtolayout;
setnamed("source","polarization angle",0);# x polarization
#if "absorption analysis group is enbaled:
if (AbsorptionEnable==1){
setnamed("absorption","Q export filename","Q_heat_x");}
run;
Rx=transmission("R");
t=getresult("T","T");
Tx=-t.T;
f=t.f;
if (AbsorptionEnable==1){
runanalysis;}

switchtolayout;
setnamed("source","polarization angle",90);# y polarization
#if "absorption analysis group is enbaled:
if (AbsorptionEnable==1){
setnamed("absorption","Q export filename","Q_heat_y");}
run;
Ry=transmission("R");
t=getresult("T","T");
Ty=-t.T;
if (AbsorptionEnable==1){
runanalysis;}
##########################
plot(c/f*1e6,Rx,Ry,Tx,Ty,"w um","R T","Reflection/transmission");
legend("R polarization 0","R polarization 90","T polarization","T polarization 90");
setplot("x min",0.4);
setplot("x max",0.8);

plot(c/f*1e6,1-Rx-Tx,1-Ry-Ty,"w um","A","Absorption");
legend("polarization 0","polarization 90");
setplot("x min",0.4);
setplot("x max",0.8);