%% load data
% clear all
fn_data = 'Fig3BDF.mat';
data = load(fn_data);

freq_sel = data.freq_sel;
x = data.x;
y = data.y;
prs = data.prs;

%% 
obs.radius = 0.094;
obs.x_center = 0;
obs.y_center = 0.095;
obs.phi_trace = linspace(0, 2*pi, 100).';
obs.rho_trace = obs.radius;
obs.x_trace = obs.rho_trace .* cos(obs.phi_trace) + obs.x_center;
obs.y_trace = obs.rho_trace .* sin(obs.phi_trace) + obs.y_center;

% prs((y-obs.y_center).^2 + (x - obs.x_center).^2 <= obs.radius^2) = 0;


%% plot 
for i = 1:length(freq_sel)
    figure;
    pp = pcolor(x{i}, y{i}, abs(prs{i}));
    
    colormap(MyColor('vik').rgb);
    shading interp
    clb = colorbar;
    xlabel('x (m)')
    ylabel('y (m)')
    
    xlim([-360, 360]/1e3)
    ylim([0, 600]/1e3)
    pbaspect([720,600,1])
    
    if freq_sel(i) > 10e3
        clim([0, 400]);
    else
        clim([0, .1]);
    end
    hold on
    fill(obs.x_trace, obs.y_trace, [1,1,1], 'edgecolor', 'none')
    
    set(gca, 'xtick', [-360, -200, 0, 200, 360]/1e3)
    fontsize(22, 'points');
        xtickangle(0)
end