% load data to matlab from txt files
% txt files produced in Ferret


%clearn out extreme inventories
filename = 'data_2010/lev_inv.txt';
data3 = importdata(filename);
INV3 = data3(:,4);
INV3(INV3>=9E14)=NaN;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Ave depth of max conc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
filename = 'data_2010/filtering.txt';
data1 = importdata(filename);
pinp1 = data1(:,3);
prise1 = data1(:,4);
pbc1 = data1(:,5);
pagg1 = data1(:,6);
kp1 = data1(:,7);
zpref1 = data1(:,8);
ping1 = data1(:,9);
MAXD1 = data1(:,2);
RUN1 = data1(:,1);

% Filter for input rates greater than 40% tot production
pinp1(pinp1>0.4)=NaN;
% Filter for inventories greater than integrated input at 2010
pinp1(isnan(INV3))=NaN;
MAXD1(isnan(pinp1))=NaN;
% Filter for aggregation fraction less than 10%
pagg1(pagg1>0.1)=NaN;
MAXD1(isnan(pagg1))=NaN;

RUN1(isnan(MAXD1))=NaN;
RUN2 = RUN1(~isnan(RUN1))
pbc1(isnan(MAXD1))=NaN;
pagg1(isnan(MAXD1))=NaN;
kp1(isnan(MAXD1))=NaN;
zpref1(isnan(MAXD1))=NaN;
ping1(isnan(MAXD1))=NaN;
prise1(isnan(MAXD1))=NaN;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Now make the plot
figure
positionVector1 = [0.1 0.9 0.07 0.05] ;
positionVector2 = [0.2 0.9 0.07 0.05] ;
positionVector3 = [0.3 0.9 0.07 0.05] ;
positionVector4 = [0.4 0.9 0.07 0.05] ;
positionVector5 = [0.5 0.9 0.07 0.05] ;
positionVector6 = [0.6 0.9 0.07 0.05] ;
positionVector7 = [0.7 0.9 0.07 0.05];

positionVector8 = [0.1 0.8 0.07 0.05] ;
positionVector9 = [0.2 0.8 0.07 0.05] ;
positionVector10 = [0.3 0.8 0.07 0.05];
positionVector11 = [0.4 0.8 0.07 0.05];
positionVector12 = [0.5 0.8 0.07 0.05];
positionVector13 = [0.6 0.8 0.07 0.05];
positionVector14 = [0.7 0.8 0.07 0.05];

positionVector15 = [0.1 0.7 0.07 0.05];
positionVector16 = [0.2 0.7 0.07 0.05];
positionVector17 = [0.3 0.7 0.07 0.05];
positionVector18 = [0.4 0.7 0.07 0.05];
positionVector19 = [0.5 0.7 0.07 0.05] ;
positionVector20 = [0.6 0.7 0.07 0.05] ;
positionVector21 = [0.7 0.7 0.07 0.05] ;

positionVector22 = [0.1 0.6 0.07 0.05] ;
positionVector23 = [0.2 0.6 0.07 0.05] ;
positionVector24 = [0.3 0.6 0.07 0.05] ;
positionVector25 = [0.4 0.6 0.07 0.05];
positionVector26 = [0.5 0.6 0.07 0.05] ;
positionVector27 = [0.6 0.6 0.07 0.05] ;
positionVector28 = [0.7 0.6 0.07 0.05];

positionVector29 = [0.1 0.5 0.07 0.05];
positionVector30 = [0.2 0.5 0.07 0.05];
positionVector31 = [0.3 0.5 0.07 0.05];
positionVector32 = [0.4 0.5 0.07 0.05];
positionVector33 = [0.5 0.5 0.07 0.05];
positionVector34 = [0.6 0.5 0.07 0.05];
positionVector35 = [0.7 0.5 0.07 0.05];

positionVector36 = [0.1 0.4 0.07 0.05];
positionVector37 = [0.2 0.4 0.07 0.05];
positionVector38 = [0.3 0.4 0.07 0.05];
positionVector39 = [0.4 0.4 0.07 0.05];
positionVector40 = [0.5 0.4 0.07 0.05];
positionVector41 = [0.6 0.4 0.07 0.05];
positionVector42 = [0.7 0.4 0.07 0.05];

positionVector43 = [0.1 0.3 0.07 0.05];
positionVector44 = [0.2 0.3 0.07 0.05];
positionVector45 = [0.3 0.3 0.07 0.05];
positionVector46 = [0.4 0.3 0.07 0.05];
positionVector47 = [0.5 0.3 0.07 0.05];
positionVector48 = [0.6 0.3 0.07 0.05];
positionVector49 = [0.7 0.3 0.07 0.05];

positionVector50 = [0.1 0.2 0.07 0.05];
positionVector51 = [0.2 0.2 0.07 0.05];
positionVector52 = [0.3 0.2 0.07 0.05];
positionVector53 = [0.4 0.2 0.07 0.05];
positionVector54 = [0.5 0.2 0.07 0.05];
positionVector55 = [0.6 0.2 0.07 0.05];
positionVector56 = [0.7 0.2 0.07 0.05];


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot('Position',positionVector1);

scatter(prise1,MAXD1,30,'blue','filled'); 
%xlim( [0 0.1])
%ylim([200 600])
ylabel({'Ave Max'; 'Conc'; 'Depth (m)'},'FontSize',10)
hold on
scatter(prise1(71),MAXD1(71),150,'cyan','filled','p');
scatter(prise1(173),MAXD1(173),150,'magenta','filled','p');
scatter(prise1(226),MAXD1(226),150,'green','filled','p');

subplot('Position',positionVector2);
scatter(pinp1,MAXD1,30,'blue','filled');
%ylim([200 600])
hold on
scatter(pinp1(71),MAXD1(71),150,'cyan','filled','p');
scatter(pinp1(173),MAXD1(173),150,'magenta','filled','p');
scatter(pinp1(226),MAXD1(226),150,'green','filled','p');

subplot('Position',positionVector3);
scatter(pbc1,MAXD1,30,'blue','filled'); 
%ylim([200 600])
hold on
scatter(pbc1(71),MAXD1(71),150,'cyan','filled','p');
scatter(pbc1(173),MAXD1(173),150,'magenta','filled','p');
scatter(pbc1(226),MAXD1(226),150,'green','filled','p');

subplot('Position',positionVector4);
scatter(pagg1,MAXD1,30,'blue','filled');
%ylim([200 600])
hold on
scatter(pagg1(71),MAXD1(71),150,'cyan','filled','p');
scatter(pagg1(173),MAXD1(173),150,'magenta','filled','p');
scatter(pagg1(226),MAXD1(226),150,'green','filled','p');

subplot('Position',positionVector5);
scatter(kp1,MAXD1,30,'blue','filled');
%ylim([200 600])
hold on
scatter(kp1(71),MAXD1(71),150,'cyan','filled','p');
scatter(kp1(173),MAXD1(173),150,'magenta','filled','p');
scatter(kp1(226),MAXD1(226),150,'green','filled','p');

subplot('Position',positionVector6);
scatter(zpref1,MAXD1,30,'blue','filled');
%ylim([200 600])
hold on
scatter(zpref1(71),MAXD1(71),150,'cyan','filled','p');
scatter(zpref1(173),MAXD1(173),150,'magenta','filled','p');
scatter(zpref1(226),MAXD1(226),150,'green','filled','p');

subplot('Position',positionVector7);
scatter(ping1,MAXD1,30,'blue','filled');
%ylim([200 600])
hold on
scatter(ping1(71),MAXD1(71),150,'cyan','filled','p');
scatter(ping1(173),MAXD1(173),150,'magenta','filled','p');
scatter(ping1(226),MAXD1(226),150,'green','filled','p');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% Percent total MP inventory by depth
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
filename = 'data_2010/lev_inv.txt';
data1 = importdata(filename);
SURF1 = data1(:,1);
MID1 = data1(:,2);
DEEP1 = data1(:,3);
INV1 = data1(:,4);

SURF1(isnan(MAXD1))=NaN;
MID1(isnan(MAXD1))=NaN;
DEEP1(isnan(MAXD1))=NaN;
INV1(isnan(MAXD1))=NaN;

subplot('Position',positionVector8);
scatter(prise1,INV1,30,'blue','filled');
%xlim( [0 0.1])
%ylim([0 5E17])
ylabel({'Global';'MP inv';'(part)'},'FontSize',10)
%xlabel('F_R','FontSize',10) 
hold on
scatter(prise1(71),INV1(71),150,'cyan','filled','p');
scatter(prise1(173),INV1(173),150,'magenta','filled','p');
scatter(prise1(226),INV1(226),150,'green','filled','p');

subplot('Position',positionVector9);
scatter(pinp1,INV1,30,'blue','filled');
%xlabel('F_T','FontSize',10)
%ylim([0 5E17])
hold on
scatter(pinp1(71),INV1(71),150,'cyan','filled','p');
scatter(pinp1(173),INV1(173),150,'magenta','filled','p');
scatter(pinp1(226),INV1(226),150,'green','filled','p');

subplot('Position',positionVector10);
scatter(pbc1,INV1,30,'blue','filled');
%ylim([0 5E17]) 
%xlabel('F_B','FontSize',10)
hold on
scatter(pbc1(71),INV1(71),150,'cyan','filled','p');
scatter(pbc1(173),INV1(173),150,'magenta','filled','p');
scatter(pbc1(226),INV1(226),150,'green','filled','p');


subplot('Position',positionVector11);
scatter(pagg1,INV1,30,'blue','filled');
%ylim([0 5E17])
%xlabel('F_A','FontSize',10)
hold on
scatter(pagg1(71),INV1(71),150,'cyan','filled','p');
scatter(pagg1(173),INV1(173),150,'magenta','filled','p');
scatter(pagg1(226),INV1(226),150,'green','filled','p');

subplot('Position',positionVector12);
scatter(kp1,INV1,30,'blue','filled');
%ylim([0 5E17])
%xlabel('k_P','FontSize',10)
hold on
scatter(kp1(71),INV1(71),150,'cyan','filled','p');
scatter(kp1(173),INV1(173),150,'magenta','filled','p');
scatter(kp1(226),INV1(226),150,'green','filled','p');

subplot('Position',positionVector13);
scatter(zpref1,INV1,30,'blue','filled');
%ylim([0 5E17])
%xlabel('\psi_{MP}','FontSize',10)
hold on
scatter(zpref1(71),INV1(71),150,'cyan','filled','p');
scatter(zpref1(173),INV1(173),150,'magenta','filled','p');
scatter(zpref1(226),INV1(226),150,'green','filled','p');

subplot('Position',positionVector14);
scatter(ping1,INV1,30,'blue','filled');
%ylim([0 5E17])
%xlabel('R_{F:MP}','FontSize',10)
hold on
scatter(ping1(71),INV1(71),150,'cyan','filled','p');
scatter(ping1(173),INV1(173),150,'magenta','filled','p');
scatter(ping1(226),INV1(226),150,'green','filled','p');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot('Position',positionVector15);
scatter(prise1,SURF1,30,'blue','filled');
%xlim( [0 0.1])
ylim([0 1])
ylabel({'Frac';'< 100 m'},'FontSize',10)
hold on
scatter(prise1(71),SURF1(71),150,'cyan','filled','p');
scatter(prise1(173),SURF1(173),150,'magenta','filled','p');
scatter(prise1(226),SURF1(226),150,'green','filled','p');

subplot('Position',positionVector16);
scatter(pinp1,SURF1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pinp1(71),SURF1(71),150,'cyan','filled','p');
scatter(pinp1(173),SURF1(173),150,'magenta','filled','p');
scatter(pinp1(226),SURF1(226),150,'green','filled','p');

subplot('Position',positionVector17);
scatter(pbc1,SURF1,30,'blue','filled'); 
ylim([0 1])
hold on
scatter(pbc1(71),SURF1(71),150,'cyan','filled','p');
scatter(pbc1(173),SURF1(173),150,'magenta','filled','p');
scatter(pbc1(226),SURF1(226),150,'green','filled','p');

subplot('Position',positionVector18);
scatter(pagg1,SURF1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pagg1(71),SURF1(71),150,'cyan','filled','p');
scatter(pagg1(173),SURF1(173),150,'magenta','filled','p');
scatter(pagg1(226),SURF1(226),150,'green','filled','p');

subplot('Position',positionVector19);
scatter(kp1,SURF1,30,'blue','filled');
ylim([0 1])
hold on
scatter(kp1(71),SURF1(71),150,'cyan','filled','p');
scatter(kp1(173),SURF1(173),150,'magenta','filled','p');
scatter(kp1(226),SURF1(226),150,'green','filled','p');

subplot('Position',positionVector20);
scatter(zpref1,SURF1,30,'blue','filled');
ylim([0 1])
hold on
scatter(zpref1(71),SURF1(71),150,'cyan','filled','p');
scatter(zpref1(173),SURF1(173),150,'magenta','filled','p');
scatter(zpref1(226),SURF1(226),150,'green','filled','p');

subplot('Position',positionVector21);
scatter(ping1,SURF1,30,'blue','filled');
ylim([0 1])
hold on
scatter(ping1(71),SURF1(71),150,'cyan','filled','p');
scatter(ping1(173),SURF1(173),150,'magenta','filled','p');
scatter(ping1(226),SURF1(226),150,'green','filled','p');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

subplot('Position',positionVector22);
scatter(prise1,MID1,30,'blue','filled');
ylim([0 1])
xlim( [0 0.1])
ylabel({'Frac';'100-500';'m'},'FontSize',10)
hold on
scatter(prise1(71),MID1(71),150,'cyan','filled','p');
scatter(prise1(173),MID1(173),150,'magenta','filled','p');
scatter(prise1(226),MID1(226),150,'green','filled','p');

subplot('Position',positionVector23);
scatter(pinp1,MID1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pinp1(71),MID1(71),150,'cyan','filled','p');
scatter(pinp1(173),MID1(173),150,'magenta','filled','p');
scatter(pinp1(226),MID1(226),150,'green','filled','p');

subplot('Position',positionVector24);
scatter(pbc1,MID1,30,'blue','filled'); 
ylim([0 1])
hold on
scatter(pbc1(71),MID1(71),150,'cyan','filled','p');
scatter(pbc1(173),MID1(173),150,'magenta','filled','p');
scatter(pbc1(226),MID1(226),150,'green','filled','p');

subplot('Position',positionVector25);
scatter(pagg1,MID1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pagg1(71),MID1(71),150,'cyan','filled','p');
scatter(pagg1(173),MID1(173),150,'magenta','filled','p');
scatter(pagg1(226),MID1(226),150,'green','filled','p');

subplot('Position',positionVector26);
scatter(kp1,MID1,30,'blue','filled');
ylim([0 1])
hold on
scatter(kp1(71),MID1(71),150,'cyan','filled','p');
scatter(kp1(173),MID1(173),150,'magenta','filled','p');
scatter(kp1(226),MID1(226),150,'green','filled','p');

subplot('Position',positionVector27);
scatter(zpref1,MID1,30,'blue','filled');
ylim([0 1])
hold on
scatter(zpref1(71),MID1(71),150,'cyan','filled','p');
scatter(zpref1(173),MID1(173),150,'magenta','filled','p');
scatter(zpref1(226),MID1(226),150,'green','filled','p');

subplot('Position',positionVector28);
scatter(ping1,MID1,30,'blue','filled');
ylim([0 1])
hold on
scatter(ping1(71),MID1(71),150,'cyan','filled','p');
scatter(ping1(173),MID1(173),150,'magenta','filled','p');
scatter(ping1(226),MID1(226),150,'green','filled','p');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

subplot('Position',positionVector29);
scatter(prise1,DEEP1,30,'blue','filled'); 
ylim([0 1])
xlim( [0 0.1])
ylabel({'Frac';'> 500 m'},'FontSize',10)
hold on
scatter(prise1(71),DEEP1(71),150,'cyan','filled','p');
scatter(prise1(173),DEEP1(173),150,'magenta','filled','p');
scatter(prise1(226),DEEP1(226),150,'green','filled','p');

subplot('Position',positionVector30);
scatter(pinp1,DEEP1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pinp1(71),DEEP1(71),150,'cyan','filled','p');
scatter(pinp1(173),DEEP1(173),150,'magenta','filled','p');
scatter(pinp1(226),DEEP1(226),150,'green','filled','p');

subplot('Position',positionVector31);
scatter(pbc1,DEEP1,30,'blue','filled'); 
ylim([0 1])
hold on
scatter(pbc1(71),DEEP1(71),150,'cyan','filled','p');
scatter(pbc1(173),DEEP1(173),150,'magenta','filled','p');
scatter(pbc1(226),DEEP1(226),150,'green','filled','p');

subplot('Position',positionVector32);
scatter(pagg1,DEEP1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pagg1(71),DEEP1(71),150,'cyan','filled','p');
scatter(pagg1(173),DEEP1(173),150,'magenta','filled','p');
scatter(pagg1(226),DEEP1(226),150,'green','filled','p');

subplot('Position',positionVector33);
scatter(kp1,DEEP1,30,'blue','filled');
ylim([0 1])
hold on
scatter(kp1(71),DEEP1(71),150,'cyan','filled','p');
scatter(kp1(173),DEEP1(173),150,'magenta','filled','p');
scatter(kp1(226),DEEP1(226),150,'green','filled','p');


subplot('Position',positionVector34);
scatter(zpref1,DEEP1,30,'blue','filled');
ylim([0 1])
hold on
scatter(zpref1(71),DEEP1(71),150,'cyan','filled','p');
scatter(zpref1(173),DEEP1(173),150,'magenta','filled','p');
scatter(zpref1(226),DEEP1(226),150,'green','filled','p');

subplot('Position',positionVector35);
scatter(ping1,DEEP1,30,'blue','filled');
ylim([0 1])
hold on
scatter(ping1(71),DEEP1(71),150,'cyan','filled','p');
scatter(ping1(173),DEEP1(173),150,'magenta','filled','p');
scatter(ping1(226),DEEP1(226),150,'green','filled','p');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Percent total MP inventory by latitude
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
filename = 'data_2010/lo_lat.txt';
data1 = importdata(filename);

LL1 = data1(:,1);
LL1(isnan(MAXD1))=NaN;


subplot('Position',positionVector36);
scatter(prise1,LL1,30,'blue','filled');
xlim( [0 0.1])
ylim([0 1])
ylabel({'Frac'; '40S-40N'},'FontSize',10)
hold on
scatter(prise1(71),LL1(71),150,'cyan','filled','p');
scatter(prise1(173),LL1(173),150,'magenta','filled','p');
scatter(prise1(226),LL1(226),150,'green','filled','p');

subplot('Position',positionVector37);
scatter(pinp1,LL1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pinp1(71),LL1(71),150,'cyan','filled','p');
scatter(pinp1(173),LL1(173),150,'magenta','filled','p');
scatter(pinp1(226),LL1(226),150,'green','filled','p');

subplot('Position',positionVector38);
scatter(pbc1,LL1,30,'blue','filled'); 
ylim([0 1])
hold on
scatter(pbc1(71),LL1(71),150,'cyan','filled','p');
scatter(pbc1(173),LL1(173),150,'magenta','filled','p');
scatter(pbc1(226),LL1(226),150,'green','filled','p');

subplot('Position',positionVector39);
scatter(pagg1,LL1,30,'blue','filled');
ylim([0 1])
hold on
scatter(pagg1(71),LL1(71),150,'cyan','filled','p');
scatter(pagg1(173),LL1(173),150,'magenta','filled','p');
scatter(pagg1(226),LL1(226),150,'green','filled','p');

subplot('Position',positionVector40);
scatter(kp1,LL1,30,'blue','filled');
ylim([0 1])
hold on
scatter(kp1(71),LL1(71),150,'cyan','filled','p');
scatter(kp1(173),LL1(173),150,'magenta','filled','p');
scatter(kp1(226),LL1(226),150,'green','filled','p');


subplot('Position',positionVector41);
scatter(zpref1,LL1,30,'blue','filled');
ylim([0 1])
hold on
scatter(zpref1(71),LL1(71),150,'cyan','filled','p');
scatter(zpref1(173),LL1(173),150,'magenta','filled','p');
scatter(zpref1(226),LL1(226),150,'green','filled','p');

subplot('Position',positionVector42);
scatter(ping1,LL1,30,'blue','filled');
ylim([0 1])
hold on
scatter(ping1(71),LL1(71),150,'cyan','filled','p');
scatter(ping1(173),LL1(173),150,'magenta','filled','p');
scatter(ping1(226),LL1(226),150,'green','filled','p');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
filename = 'data_2010/hi_lat.txt';
data1 = importdata(filename);

SO1 = data1(:,2);
AR1 = data1(:,3);

SO1(isnan(MAXD1))=NaN;
AR1(isnan(MAXD1))=NaN;

subplot('Position',positionVector43);
scatter(prise1,SO1,30,'blue','filled');
xlim( [0 0.1])
ylim([0 0.02])
ylabel({'Frac';'> 60S'},'FontSize',10)
hold on
scatter(prise1(71),SO1(71),150,'cyan','filled','p');
scatter(prise1(173),SO1(173),150,'magenta','filled','p');
scatter(prise1(226),SO1(226),150,'green','filled','p');

subplot('Position',positionVector44);
scatter(pinp1,SO1,30,'blue','filled');
ylim([0 0.02])
hold on
scatter(pinp1(71),SO1(71),150,'cyan','filled','p');
scatter(pinp1(173),SO1(173),150,'magenta','filled','p');
scatter(pinp1(226),SO1(226),150,'green','filled','p');

subplot('Position',positionVector45);
scatter(pbc1,SO1,30,'blue','filled'); 
ylim([0 0.02])
hold on
scatter(pbc1(71),SO1(71),150,'cyan','filled','p');
scatter(pbc1(173),SO1(173),150,'magenta','filled','p');
scatter(pbc1(226),SO1(226),150,'green','filled','p');

subplot('Position',positionVector46);
scatter(pagg1,SO1,30,'blue','filled');
ylim([0 0.02])
hold on
scatter(pagg1(71),SO1(71),150,'cyan','filled','p');
scatter(pagg1(173),SO1(173),150,'magenta','filled','p');
scatter(pagg1(226),SO1(226),150,'green','filled','p');

subplot('Position',positionVector47);
scatter(kp1,SO1,30,'blue','filled');
ylim([0 0.02])
hold on
scatter(kp1(71),SO1(71),150,'cyan','filled','p');
scatter(kp1(173),SO1(173),150,'magenta','filled','p');
scatter(kp1(226),SO1(226),150,'green','filled','p');

subplot('Position',positionVector48);
scatter(zpref1,SO1,30,'blue','filled');
ylim([0 0.02])
hold on
scatter(zpref1(71),SO1(71),150,'cyan','filled','p');
scatter(zpref1(173),SO1(173),150,'magenta','filled','p');
scatter(zpref1(226),SO1(226),150,'green','filled','p');

subplot('Position',positionVector49);
scatter(ping1,SO1,30,'blue','filled');
ylim([0 0.02])
hold on
scatter(ping1(71),SO1(71),150,'cyan','filled','p');
scatter(ping1(173),SO1(173),150,'magenta','filled','p');
scatter(ping1(226),SO1(226),150,'green','filled','p');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot('Position',positionVector50);
scatter(prise1,AR1,30,'blue','filled');
xlim( [0 0.1])
ylim([0 0.05])
xlabel('F_R','FontSize',10) 
ylabel({'Frac';'> 60N'},'FontSize',10)
hold on
scatter(prise1(71),AR1(71),150,'cyan','filled','p');
scatter(prise1(173),AR1(173),150,'magenta','filled','p');
scatter(prise1(226),AR1(226),150,'green','filled','p');

subplot('Position',positionVector51);
scatter(pinp1,AR1,30,'blue','filled');
ylim([0 0.05])
xlabel('F_T','FontSize',10) 
hold on
scatter(pinp1(71),AR1(71),150,'cyan','filled','p');
scatter(pinp1(173),AR1(173),150,'magenta','filled','p');
scatter(pinp1(226),AR1(226),150,'green','filled','p');

subplot('Position',positionVector52);
scatter(pbc1,AR1,30,'blue','filled'); 
ylim([0 0.05])
xlabel('F_B','FontSize',10)
hold on
scatter(pbc1(71),AR1(71),150,'cyan','filled','p');
scatter(pbc1(173),AR1(173),150,'magenta','filled','p');
scatter(pbc1(226),AR1(226),150,'green','filled','p');

subplot('Position',positionVector53);
scatter(pagg1,AR1,30,'blue','filled');
ylim([0 0.05])
xlabel('F_A','FontSize',10)
hold on
scatter(pagg1(71),AR1(71),150,'cyan','filled','p');
scatter(pagg1(173),AR1(173),150,'magenta','filled','p');
scatter(pagg1(226),AR1(226),150,'green','filled','p');

subplot('Position',positionVector54);
scatter(kp1,AR1,30,'blue','filled');
ylim([0 0.05])
xlabel('k_P','FontSize',10)
hold on
scatter(kp1(71),AR1(71),150,'cyan','filled','p');
scatter(kp1(173),AR1(173),150,'magenta','filled','p');
scatter(kp1(226),AR1(226),150,'green','filled','p');

subplot('Position',positionVector55);
scatter(zpref1,AR1,30,'blue','filled');
ylim([0 0.05])
xlabel('\psi_{MP}','FontSize',10)
hold on
scatter(zpref1(71),AR1(71),150,'cyan','filled','p');
scatter(zpref1(173),AR1(173),150,'magenta','filled','p');
scatter(zpref1(226),AR1(226),150,'green','filled','p');

subplot('Position',positionVector56);
scatter(ping1,AR1,30,'blue','filled');
ylim([0 0.05])
xlabel('R_{F:MP}','FontSize',10)
hold on
scatter(ping1(71),AR1(71),150,'cyan','filled','p');
scatter(ping1(173),AR1(173),150,'magenta','filled','p');
scatter(ping1(226),AR1(226),150,'green','filled','p');
