clear all
load f_correct.mat

left_figure = subplot(131);
left_pos = get(left_figure, 'Position');
plot(x,[phi_0,phi,phi_1],'linewidth',3,'color',[9,77,140]/256);
title('phase field \phi');
axis([R-0.5,R+0.5,-0.05,1]);
set(gca,'xtick',[2.4:0.2:3.5],'fontsize',15);
set(gca,'ytick',[0.0:0.2:1],'fontsize',15);
grid on;
hold on;

load f_wrong.mat
plot(x,[phi_0,phi,phi_1],'--','linewidth',3,'color',[20,139,67]/256);
plot([3.15,3.42,3.42,3.15,3.15],[-0.025,-0.025,0.025,0.025,-0.025],'-','linewidth',1,'color',[90,90,90]/256);
quiver(3.26,0.05,0,0.1,'MaxHeadSize',1.0,'linewidth',2,'color',[90,90,90]/256);
h_legend = legend('\phi in new model','\phi in old model');
set(h_legend,'Location','NorthWest','FontSize',17);


hold off;

ax1 = axes('Position',left_pos);
plot(ax1,x,[phi_0,phi,phi_1],'--','linewidth',3,'color',[20,139,67]/256);hold on;
load f_correct.mat
plot(ax1,x,[phi_0,phi,phi_1],'-','linewidth',3,'color',[9,77,140]/256);
set(ax1, 'Units', 'normalized', 'Position', [0.27 0.30 0.07 0.1]);
set(ax1, 'XLim',[R+0.3,R+0.5], 'YLim',[-1e-3 1e-3],'xtick',[3.3,3.4],'ytick',[0,1e-3],'fontsize',12);
hold off

mid_figure = subplot(132);
mid_pos = get(mid_figure, 'Position');
plot(r,force_surf,':','linewidth',2,'color',[9,77,140]/256);hold on;
plot(r,force_LJ,'-.','linewidth',2,'color',[20,139,67]/256);
plot(r,force_CFL,'--','linewidth',2,'color',[204,81,39]/256);
plot(r,force_ttl,'-','linewidth',2,'color',[90,90,90]/256);
plot([2.89,3.03,3.03,2.89,2.89],[-0.025,-0.025,0.025,0.025,-0.025],'-','linewidth',1,'color',[90,90,90]/256);
quiver(2.88,-0.03,-0.05,-0.05,'MaxHeadSize',1.9,'linewidth',2,'color',[90,90,90]/256)
title('three forces in new model')
axis([R-0.5,R+0.5,-0.2,0.4]);
set(gca,'xtick',[2.4:0.2:3.5],'fontsize',15);
set(gca,'ytick',[-0.4:0.1:0.4],'fontsize',15);
h2_legend = legend('tension','LJ','elec','total');
set(h2_legend,'Location','NorthWest','FontSize',17);
grid on;
hold off;

ax2 = axes('Position',mid_pos);
plot(ax2,r,force_ttl,'-','linewidth',2,'color',[90,90,90]/256);
set(ax2, 'Units', 'normalized', 'Position', [0.43 0.17 0.07 0.1]);
set(ax2, 'XLim',[R-0.1,R+0.1], 'YLim',[-2e-4 1e-3],'xtick',[2.9,3],'ytick',[0,1e-3],'fontsize',12);
hold off


load f_wrong.mat

right_figure = subplot(133);
right_pos = get(right_figure, 'Position');
plot(r,force_surf,':','linewidth',2,'color',[9,77,140]/256);hold on;
plot(r,force_LJ,'-.','linewidth',2,'color',[20,139,67]/256);
plot(r,force_CFL,'--','linewidth',2,'color',[204,81,39]/256);
plot(r,force_ttl,'-','linewidth',2,'color',[90,90,90]/256);
plot([2.89,3.03,3.03,2.89,2.89],[-0.025,-0.025,0.025,0.025,-0.025],'-','linewidth',1,'color',[90,90,90]/256);
quiver(2.88,-0.03,-0.05,-0.05,'MaxHeadSize',1.9,'linewidth',2,'color',[90,90,90]/256)
title('three forces in old model')
axis([R-0.5,R+0.5,-0.2,0.4]);
set(gca,'xtick',[2.4:0.2:3.5],'fontsize',15);
set(gca,'ytick',[-0.4:0.1:0.4],'fontsize',15);
h2_legend = legend('tension','LJ','elec','total');
set(h2_legend,'Location','NorthWest','FontSize',17);
grid on;
hold off;

ax3 = axes('Position',right_pos);
plot(ax3,r,force_ttl,'-','linewidth',2,'color',[90,90,90]/256);
set(ax3, 'Units', 'normalized', 'Position', [0.71 0.17 0.07 0.1]);
set(ax3, 'XLim',[R-0.1,R+0.1], 'YLim',[-5e-3 1e-3],'xtick',[2.9,3],'ytick',[-5e-3,0],'fontsize',12);
hold off

Published with MATLAB® R2016b