close all;
fSize = 16;

data_mu = mean(data,3);
data_filt = medfilt2(data_mu,5*[1 1]);
clevels = 0.25:0.25:5;
contour(Thetas,Sigs*V1*eps*1e3,data_filt,clevels); hold on;
contour(Thetas,Sigs*V1*eps*1e3,data_filt,[1 1],'b','Linewidth',3);
% subplot(121);
% surfc(Thetas,Sigs*100,data_mu); 
colorbar;
set(gca,'FontSize',fSize);
xlabel('Phase offset (degrees)');
ylabel('Error in V_2 (mV)');
xlim([0 3]);
ylim([0 30]);
zlim([0 10]);
title('Peak-to-peak amplitude (um)');

% subplot(122);
% surfc(Thetas,Sigs*100,data_filt); 
% set(gca,'FontSize',fSize);
% xlabel('Phase offset (degrees)');
% ylabel('Ratio jitter (%)');
% zlim([0 10]);
% title('Smoothed');