clear all
clc

% subplot(5,2,9);
% West1z=[
%     0
%     25
%     25.01
%     50
%     50.01
%     75
%     75.01
%     100
%     ];
% West1c=[
%     3.4
%     3.4
%     4.7
%     4.7
%     2.16
%     2.16
%     2.73
%     2.73
%     ];
% 
% Xerror9=[
%     12.5
%     37.5
%     62.5
%     87.5];
% Yerror9=[3.4
%     4.7
%     2.16
%     2.73];
% Error9=[1.02
% 0.811
% 0.875
% 1.02
% ];
%     
% plot(West1z,West1c,'LineWidth',2);
% grid on
% xlim([0 100]);
% ylim([0 7]);
% hold on
% 
% for i=1:1:4  
%     errorbar(Xerror9(i,1),Yerror9(i,1),Error9(i,1),'r','LineWidth',2);
%     hold on
% end
% 
% title('9 West 1','fontsize',8,'FontWeight','Bold');
% orient tall
% view(90,90)
% ylabel('Cesium (Bq/kg DM)','fontsize',8);
% 
% % ------------------------------------------------
% 
% subplot(5,2,10);
% West2z=[
%     0
%     23
%     23.01
%     46
%     46.01
%     69
%     69.01
%     95
%     ];
% West2c=[
%     3.45
%     3.45
%     5.19
%     5.19
%     4.61
%     4.61
%     4.53
%     4.53
%     ];
% 
% Xerror10=[
%     11.5
%     34.5
%     57.5
%     80.5];
% Yerror10=[
%     3.45
%     5.19
%     4.61
%     4.53];
% Error10=[
%   0.937
% 1.03
% 0.85
% 1.01
% ];
% 
% plot(West2z,West2c,'LineWidth',2);
% grid on
% xlim([0 100]);
% ylim([0 7]);
% hold on
% for i=1:1:4  
%     errorbar(Xerror10(i,1),Yerror10(i,1),Error10(i,1),'r','LineWidth',2);
%     hold on
% end
% ylabel('Cesium (Bq/kg DM)','fontsize',8);
% view(90,90)
% title('10 West 2','fontsize',8,'FontWeight','Bold');
% orient tall

% -------------------------------------------------------

subplot(1,2,2);
PBoutz=[
    0
    21
    21.01
    42
    42.01
    63
     ];
PBoutc=[
    4.42
    4.42
    2.7
    2.7
    2.09
    2.09
    ];

Xerror8=[
    10.5
    31.5
    53.5
    ];
Yerror8=[
    4.42
    2.7
    2.09];
Error8=[
0.909
0.428
0.741
];

plot(PBoutz,PBoutc,'LineWidth',4);
grid on
xlim([0 60]);
ylim([0 6]);
hold on
for i=1:1:3  
    errorbar(Xerror8(i,1),Yerror8(i,1),Error8(i,1),'r','LineWidth',3);
    hold on
end
ylabel('Cesium content (Bq/kg DM)','fontsize',12);
xlabel('Core depth (cm)','fontsize',12);
title('8 Proefput B out','fontsize',12,'FontWeight','Bold');
%orient tall
view(90,90)

%-------------------------------------------------------------

subplot(1,2,1);
PBinz=[
    0
    18
    18.01
    36
    36.01
    54
     ];
PBinc=[
    17.3
    17.3
    20.9
    20.9
    20.5
    20.5
    ];
Xerror7=[
    9
    27
    45
    ];
Yerror7=[
    17.3
    20.9
    20.5];
Error7=[
1.27
0.187
1.7
];


plot(PBinz,PBinc,'LineWidth',4);
grid on
xlim([0 60]);
ylim([0 25]);
hold on
for i=1:1:3
    errorbar(Xerror7(i,1),Yerror7(i,1),Error7(i,1),'r','LineWidth',4);
    hold on
end

ylabel('Cesium content (Bq/kg DM)','fontsize',12);
title('7 Proefput B in','fontsize',12,'FontWeight','Bold');
xlabel('Core depth (cm)','fontsize',12);
%orient tall
view(90,90)


% % -----------------------------------------------------------------
% 
% subplot(5,2,6);
% P1z=[
%     0
%     15
%     15.01
%     30
%      ];
% P1c=[
%     4.45
%     4.45
%     2.1
%     2.1
%     ];
% Xerror6=[
%     7.5
%     22.5
%     ];
% Yerror6=[
%     4.45
%     2.1
%     ];
% Error6=[
% 0.157
% 0.712
% ];
% 
% plot(P1z,P1c,'LineWidth',2);
% grid on
% xlim([0 40]);
% ylim([0 7]);
% hold on
% for i=1:1:2  
%     errorbar(Xerror6(i,1),Yerror6(i,1),Error6(i,1),'r','LineWidth',2);
%     hold on
% end
% %xlabel('Cesium (Bq)','fontsize',8);
% title('6 Pilot Pit 1','fontsize',8,'FontWeight','Bold');
% orient tall
% view(90,90)
% %-----------------------------------------------------------------
% 
% subplot(5,2,5);
% P2z=[
%     0
%     26
%     26.01
%     52
%     52.01
%     78
%     ];
% P2c=[
%     1.97
%     1.97
%     4.12
%     4.12
%     1.8
%     1.8
%     ];
% Xerror5=[
%     13
%     39
%     65
%     ];
% Yerror5=[
%     1.97
%     4.12
%     1.8];
% Error5=[
% 0.489
% 0.934
% 0.83
% ];
% 
% plot(P2z,P2c,'LineWidth',2);
% grid on
% xlim([0 100]);
% ylim([0 7]);
% hold on
% for i=1:1:3  
%     errorbar(Xerror5(i,1),Yerror5(i,1),Error5(i,1),'r','LineWidth',2);
%     hold on
% end
% %xlabel('Cesium (Bq)','fontsize',8);
% xlabel('Core Depth (cm)','fontsize',8);
% title('5 Pilot Pit 2','fontsize',8,'FontWeight','Bold');
% orient tall
% view(90,90)
% 
% %----------------------------------------------------------------------
% 
% subplot(5,2,4);
% P3z=[
%     0
%     23
%     23.01
%     46
%     46.01
%     69
%      ];
% P3c=[
%     3.06
%     3.06
%     3.38
%     3.38
%     3.09
%     3.09
%     ];
% Xerror4=[
%     11.5
%     34.5
%     57.5
%     ];
% Yerror4=[
%     3.06
%     3.38
%     3.09];
% Error4=[
% 0.189
% 0.207
% 0.548
% ];
% 
% plot(P3z,P3c,'LineWidth',2);
% grid on
% xlim([0 80]);
% ylim([0 7]);
% hold on
% for i=1:1:3  
%     errorbar(Xerror4(i,1),Yerror4(i,1),Error4(i,1),'r','LineWidth',2);
%     hold on
% end
% %xlabel('Cesium (Bq)','fontsize',8);
% title('4 Pilot Pit 3','fontsize',8,'FontWeight','Bold');
% orient tall
% view(90,90)
% %-----------------------------------------------------------------
% 
% subplot(5,2,3);
% PAoutz=[
%     0
%     25
%     25.01
%     50
%       ];
% PAoutc=[
%     5.5
%     5.5
%     3.14
%     3.14
%     ];
% Xerror3=[
%     12.5
%     37.5
%     ];
% Yerror3=[
%     5.5
%     3.14
%     ];
% Error3=[
% 0.811
% 0.684
% ];
% plot(PAoutz,PAoutc,'LineWidth',2);
% grid on
% xlim([0 60]);
% ylim([0 7]);
% hold on
% for i=1:1:2  
%     errorbar(Xerror3(i,1),Yerror3(i,1),Error3(i,1),'r','LineWidth',2);
%     hold on
% end
% %xlabel('Cesium (Bq)','fontsize',8);
% title('3 Proefput A out','fontsize',8,'FontWeight','Bold');
% orient tall
% view(90,90)
% %-----------------------------------------------------------------
% 
% subplot(5,2,2);
% PAinz=[
%     0
%     25
%       ];
% PAinc=[
%     14.1
%     14.1
%     ];
% Xerror2=[
%     12.5
%     ];
% Yerror2=[
%     14.1
%     ];
% Error2=[
% 1.37
% ];
% plot(PAinz,PAinc,'LineWidth',2);
% grid on
% xlim([0 30]);
% ylim([0 25]);
% hold on
% for i=1:1:1
%     errorbar(Xerror2(i,1),Yerror2(i,1),Error2(i,1),'r','LineWidth',2);
%     hold on
% end
% %xlabel('Cesium (Bq)','fontsize',8);
% title('2 Proefput A in','fontsize',8,'FontWeight','Bold');
% orient tall
% view(90,90)
% %------------------------------------------------------------------
% 
% subplot(5,2,1);
% FL42z=[
%     0
%     23
%     23.01
%     46
%     46.01
%     69
%       ];
% FL42c=[
%     2.01
%     2.01
%     2.35
%     2.35
%     1.58
%     1.58
%     ];
% Xerror1=[
%     11.5
%     34.5
%     57.5
%     ];
% Yerror1=[
%     2.01
%     2.35
%     1.58
%     ];
% Error1=[
% 0.764
% 0.831
% 0.775
% ];
% 
% plot(FL42z,FL42c,'LineWidth',2);
% grid on
% xlim([0 80]);
% ylim([0 7]);
% hold on
% for i=1:1:3
%     errorbar(Xerror1(i,1),Yerror1(i,1),Error1(i,1),'r','LineWidth',2);
%     hold on
% end
% %xlabel('Cesium (Bq)','fontsize',8);
% title('1 FL42','fontsize',8,'FontWeight','Bold');
% orient tall
% view(90,90)