clear;

load data/Rita_ent_dat.mat;
h1=plot(dth_e,dq,'b*');
hold on;
load data/Harvey_ent_dat.mat;
h2=plot(dth_e,dq,'ro');
load data/Michael_ent_dat.mat;
h3=plot(dth_e,dq,'gs');
load data/Patricia_ent_dat.mat;
h4=plot(dth_e,dq,'k>');
plot([-12 0],[-4 0],'--','color',[0.6 0.6 0.6]);
h=legend([h1,h2,h3,h4],'Rita 50 legs','Harvey 22 legs','Michael 34 legs','Patricia 7 legs');
set(h,'location','southeast','box','off','fontsize',13);
axis([-12,0,-4,0]);
xlabel('$\Delta \theta_{e}$ (K)','Interpreter','latex');
ylabel('$\overline{\theta}\Delta q / \alpha$ (K)','Interpreter','latex');
h=text(-10,-0.7,'Unstable');set(h,'fontsize',13);
h=text(-9,-3.5,'Stable');set(h,'fontsize',13);
h=title('Entrainment instability criterion $\Delta \theta_{e} < \overline{\theta}\Delta q / \alpha$','Interpreter','latex');

hold off;

print -dpng obs_entrain_plot
