% Make eigenfunction intensity map
decim=1;
let_len=20; % letter length in grid spacings in numerics
let_len_mkm=1; %1.35; % letter length in mkm
res=0.01/1;
NN=99;
 [map,ens] = eigmode_map(ff(1:NN,[1:decim:end]),idos(1,2:2:2*NN),res);
 figure
 xxxx=size(map,2);
 imagesc([1:xxxx]*let_len_mkm/(let_len/1/decim),ens,map); axis xy
 xlabel('x [\mum]','Interpreter','Tex')
 xlabel('x [letters]','Interpreter','Tex')
 ylabel('$E [meV]$','Interpreter','LaTex')
 clear xxxx