%dfname = 'dm_max_data.dat';
dfname = 'dp_max_data.dat';
%dfname = 'fm_max_data.dat';
%dfname = 'fp_max_data.dat';

DD = load(dfname);

rhovec = union(DD(:,4),DD(:,4));
disp(sprintf('There are %d possible rho values',length(rhovec)));
count = input('Which entries do you want to plot? input vector of indices :  ');

for i = 1:length(count);
  figure;
  idx = find(DD(:,4) == rhovec(count(i)));
  scatter3(DD(idx,1),DD(idx,2),DD(idx,3),'.');
  titlestr = sprintf('%s, rho = %d',dfname(1:2),rhovec(count(i)));
  xlabel('k','FontSize',18);
  ylabel('B','FontSize',18);
  zlabel('Max(\Re(\sigma))','FontSize',18);
  title(titlestr,'FontSize',18);
end

