function figure2()
  model_config;

  [~, out, f0, f1, f2, f3] = analyze('2019-05-01-14-18-13', true, '');
  
  figure;
  ax0 = subplot(3,1,1);
  plot(out.T, out.C, '.');
  hold on
  plot(out.t, out.r1, '.');
  idx = find(f1);
  plot(out.t(idx), out.r1(idx), '*');
  ylabel('A or mm/s')
  title('Nominal');
  xlim([80 220]);
  ylim([-0.4 0.4])
  legend('current (A)', 'r1 (mm/s)', 'r1 faults')
  grid on
  
  [~, out, f0, f1, f2, f3] = analyze('2019-04-30-08-07-19', false, '');
  
  ax0 = subplot(3,1,2);
  plot(out.T, out.C, '.');
  hold on
  plot(out.t, out.r1, '.');
  idx = find(f1);
  plot(out.t(idx), out.r1(idx), '*');
  ylabel('A or mm/s')
  title('Shorter fault against aft limit')
  xlim([105 117]);
  ylim([-0.4 0.4])
  legend('current (A)', 'r1 (mm/s)', 'r1 faults')
  grid on

  [~, out, f0, f1, f2, f3] = analyze('2019-05-07-16-25-48', true, '');
  
  ax0 = subplot(3,1,3);
  plot(out.T, out.C, '.');
  hold on
  plot(out.t, out.r1, '.');
  idx = find(f1);
  plot(out.t(idx), out.r1(idx), '*');
  ylabel('A or mm/s')
  title('Longer fault against aft limit')
  xlim([175 192]);
  ylim([-0.4 0.4])
  legend('current (A)', 'r1 (mm/s)', 'r1 faults')
  grid on

  set(gcf,'position',[10,10,600,1800])
  print('2.png', '-dpng', '-r300')


end
