function RangeDopplerRadar_plot(fid, titlestr, dumpstr, ...
                         m, n, s, SNR, mu, GERR, FERR, ...
                         gammatil_S,abs_theta_S, gammatil, absthetatil, ...
                         tau, alpha, abstheta, tautil_S,alphatil_S)
%%
% figure(21); plot(gammatil(S),abs(theta(S)), 'ro', gamma+zeta_til, abs_theta_til,'bx'); axis('tight');
% save2pdf('fig_RangeDopplerRadar_fixpt_gridfixON_OMPsc1.pdf', gcf); %saveas2('fig_RangeDopplerRadar_fixpt_gridfixON1.png');
% 
% figure(23); imagesc(tau, alpha, reshape(abs_theta_til,n,n)'); 
% xlabel('Time delay'); ylabel('Doppler shift');
% caxis([0, ceil(max([abs_theta_til;abs(theta)]))]); colormap('summer');
% hold on; plot(tautil(S_xi),alphatil(S_eta),'ro', 'LineWidth', 2); hold off;
% title({'Range-Doppler Image, off-grid fixing on, OMPsc. '; infostr(m,n*n,s,SNR,coherence(A)) });   
% colorbar; axis('square'); 
% save2pdf('fig_RangeDopplerRadar_fixpt_gridfixON_OMPsc2.pdf', gcf); %saveas2('fig_RangeDopplerRadar_fixpt_gridfixON2.png');
%   

%% Misc
%infostr = @(m,n,s,SNR,mu) strcat('m=', int2str(m), ...
%                                 ', n=', int2str(n), ...
%                                 ', s=', int2str(s), ...
%                                 ', SNR=', num2str(SNR), ...
%                                 ', \mu=', num2str(mu) ); 
infostr = @(m,n,s,SNR,mu) strcat('m=', int2str(m), ...
                                 ', n=', int2str(n), ...
                                 ', s=', int2str(s), ...
                                 ', \mu=', num2str(mu), ...
                                 ', GERR=', num2str(GERR), ...
                                 ', FERR=', num2str(FERR)  ); 
%%

figure(fid+1); %plot(gammatil_S,abs_theta_S, 'ro', gammatil, absthetatil,'bx'); axis('tight');
stem(gammatil_S,abs_theta_S, 'ro');
hold on
stem(gammatil, absthetatil, 'bx'); 
xlabel('\gamma'); ylabel('|\theta|');
axis('tight');
hold off
save2pdf( sprintf('%s_fig%d.pdf',dumpstr,1) , gcf); %saveas2('fig_RangeDopplerRadar_fixpt_gridfixON1.png');


figure(fid+2); 
imagesc(tau, alpha, reshape(absthetatil,n,n)'); 
xlabel('Time delay'); ylabel('Doppler shift');
caxis([0, ceil(max([absthetatil;abstheta]))]); colormap('summer');
hold on; plot(tautil_S, alphatil_S,'ro', 'LineWidth', 2); hold off;
title({titlestr; infostr(m,n*n,s,SNR,mu) });   
colorbar; axis('square'); 
save2pdf( sprintf('%s_fig%d.pdf',dumpstr,2), gcf); %saveas2('fig_RangeDopplerRadar_fixpt_gridfixON2.png');

%%
