0001 lambda=stM.lambda;
0002
0003 figure('Name','Q-coefficients and average LFIEF on surface');
0004 plot(lambda,stM.Qext,'b',lambda,stM.Qsca,'k',lambda,stM.Qabs,'k:');
0005 ax1 = gca;
0006 set(ax1,'XColor','k','YColor','b')
0007 set(get(ax1,'XLabel'),'String','Wavelength [nm]');
0008 set(get(ax1,'YLabel'),'String','Q_{ext},Q_{sca},Q_{abs}');
0009 legend({'Q_{ext}','Q_{sca}','Q_{abs}'});
0010 title(['Mie for PWE: nNmax=' int2str(stM.nNmax) ', a=' num2str(stM.a) 'nm']);
0011 ax2=axes('Position',get(ax1,'Position'),...
0012 'XAxisLocation','top','XColor','k', ...
0013 'YAxisLocation','right','Ycolor','r', ...
0014 'Color','none', ...
0015 'YScale','Log', ...
0016 'XTickLabel',{});
0017 set(get(ax2,'Ylabel'),'String','Average LFIEF on surface');
0018
0019 line(lambda,stM.MLocAve,'Color','r');