function fig = plot_cal_diff(enh,cal_diff,flags,mark)
    u = repmat(1:16,16,1);
    mask = zeros(size(u));
    mask(u == mark) = 1;
    mask = mask(:);
    idxMask = find(mask);
    idxMask_not = find(~mask);
    fig = figure;
    scatter(enh(idxMask,1),enh(idxMask,2),120,cal_diff(idxMask),'filled','LineWidth',0.1,'Marker','s');
    hold on
    scatter(enh(idxMask_not,1),enh(idxMask_not,2),80,cal_diff(idxMask_not),'filled','LineWidth',0.1);
    scatter(enh(find(flags & mask),1),enh(find(flags & mask),2),120,[0,0,0],'filled','Marker','s');
    scatter(enh(find(flags & ~mask),1),enh(find(flags & ~mask),2),80,[0,0,0],'filled');
    caxis([-60,60])
    axis([-20,20,-20,20]);
    axis equal
    cb = colorbar;colormap('jet')
    title('Difference to Reference Phase Calibration, in Degrees')
    xlabel('East/m');
    ylabel('North/m');
end

