figure;
set(gcf,'Color','w');
plot(disp/r0/2,Force/r0/sigmay,'-.');
hold;
plot(xydat(:,1)*0.5/r0/2,2*xydat(:,2)/r0/sigmay,'o');

xlabel('\Delta/D','FontSize',16);
ylabel('P/D/S_u','FontSize',16);
set(gca,'XLim',[0,0.12]);
set(gca,'YLim',[0,16]);
legend('Currrent analysis','FE analysis');

dlmwrite('piledisp.dat',disp,'\t');
dlmwrite('pileload.dat',Force,'\t');
dlmwrite('r_points.dat',ri,'\t');
dlmwrite('theta_points.dat',ti,'\t');

% dlmwrite('q.dat',q,'\t');
dlmwrite('U.dat',U,'\t');
dlmwrite('V.dat',V,'\t');

dlmwrite('Uinterp.dat',Uinterp,'\t');
dlmwrite('Vinterp.dat',Vinterp,'\t');


dlmwrite('erE.dat',erE,'\t');
dlmwrite('etE.dat',etE,'\t');
dlmwrite('ertE.dat',ertE,'\t');
dlmwrite('ezE.dat',ezE,'\t');

dlmwrite('erP.dat',erPc,'\t');
dlmwrite('etP.dat',etPc,'\t');
dlmwrite('ertP.dat',ertPc,'\t');
dlmwrite('ezP.dat',ezPc,'\t');

dlmwrite('Fr1.dat',Fr1,'\t');
dlmwrite('Fr2.dat',Fr2,'\t');
dlmwrite('Ft1.dat',Ft1,'\t');
dlmwrite('Ft2.dat',Ft2,'\t');

dlmwrite('derP.dat',derP,'\t');
dlmwrite('etP.dat',detP,'\t');
dlmwrite('dertP.dat',dertP,'\t');
dlmwrite('dezP.dat',dezP,'\t');

    
% % ZZ2=[zc-R:R/100:zc+R];
% % XX2=(R^2-(zc-ZZ2).^2).^0.5;
% % 
% % 
% % figure;
% % [C,h] = contour(x,z,gama*100);
% % hold;
% % [C,h] = contour(-x,z,gama*100);
% % set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2);
% % colormap cool;
% % plot(XX2,ZZ2,'Color','black','Linewidth',2);
% % plot(-XX2,ZZ2,'Color','black','Linewidth',2);
% % Create polar data
[r,t] = meshgrid(ri,tf);
% Convert to Cartesian
x = r.*cos(t);
y = r.*sin(t);
figure;
% h = polar(x,y);
[C,h] =contour(x,y,sqrt(Uinterp.^2+Vinterp.^2));
set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2);
%  set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2);
% % Hide the POLAR function data and leave annotations
% set(h,'Visible','off')
% % Turn off axes and set square aspect ratio
% axis off
% axis image