clear all; close all; clc

load('D50_a.mat');
load('D50_b.mat');

load('D100_a.mat');
load('D100_b.mat');
load('D200_a.mat');
load('D200_b.mat');
load('D300_a.mat');
load('D300_a.mat');

Tb1_monthly_avg=[Tb1_monthly_avg_a Tb1_monthly_avg_b Tb1_monthly_avg_c Tb1_monthly_avg_d];
Tb1_std=[Tb1_std_a Tb1_std_a Tb1_std_a Tb1_std_a]; 
% errorbar(1:12, Tb1_monthly_avg_c, Tb1_std_c, 'o-', 'LineWidth', 1.5, 'DisplayName', 'Borehole Wall Temperature');
% 
% legend('show','location','best');
% xlabel('Time [month]');
% ylabel('Temperature [°C]');
% set(gca,'fontsize',30)
% 
