load 'tjan'
load 'tfeb' 
load 'tmar'
load 'tapr'
load 'tmay'
load 'tjun'
load 'tjul'
load 'taug'
load 'tsep'
load 'toct'
load 'tnov'
load 'tdec'

hour0 = [tjan(1,:); tfeb(1,:); tmar(1,:); tapr(1,:); tmay(1,:); tjun(1,:); tjul(1,:); taug(1,:); tsep(1,:); toct(1,:); tnov(1,:); tdec(1,:)];

hour1 = [tjan(2,:); tfeb(2,:); tmar(2,:); tapr(2,:); tmay(2,:); tjun(2,:); tjul(2,:); taug(2,:); tsep(2,:); toct(2,:); tnov(2,:); tdec(2,:)];

hour2 = [tjan(3,:); tfeb(3,:); tmar(3,:); tapr(3,:); tmay(3,:); tjun(3,:); tjul(3,:); taug(3,:); tsep(3,:); toct(3,:); tnov(3,:); tdec(3,:)];

hour3 = [tjan(4,:); tfeb(4,:); tmar(4,:); tapr(4,:); tmay(4,:); tjun(4,:); tjul(4,:); taug(4,:); tsep(4,:); toct(4,:); tnov(4,:); tdec(4,:)];

hour4 = [tjan(5,:); tfeb(5,:); tmar(5,:); tapr(5,:); tmay(5,:); tjun(5,:); tjul(5,:); taug(5,:); tsep(5,:); toct(5,:); tnov(5,:); tdec(5,:)];

hour5 = [tjan(6,:); tfeb(6,:); tmar(6,:); tapr(6,:); tmay(6,:); tjun(6,:); tjul(6,:); taug(6,:); tsep(6,:); toct(6,:); tnov(6,:); tdec(6,:)];

hour6 = [tjan(7,:); tfeb(7,:); tmar(7,:); tapr(7,:); tmay(7,:); tjun(7,:); tjul(7,:); taug(7,:); tsep(7,:); toct(7,:); tnov(7,:); tdec(7,:)];

hour7 = [tjan(8,:); tfeb(8,:); tmar(8,:); tapr(8,:); tmay(8,:); tjun(8,:); tjul(8,:); taug(8,:); tsep(8,:); toct(8,:); tnov(8,:); tdec(8,:)];

hour8 = [tjan(9,:); tfeb(9,:); tmar(9,:); tapr(9,:); tmay(9,:); tjun(9,:); tjul(9,:); taug(9,:); tsep(9,:); toct(9,:); tnov(9,:); tdec(9,:)];

hour9 = [tjan(10,:); tfeb(10,:); tmar(10,:); tapr(10,:); tmay(10,:); tjun(10,:); tjul(10,:); taug(10,:); tsep(10,:); toct(10,:); tnov(10,:); tdec(10,:)];

hour10 = [tjan(11,:); tfeb(11,:); tmar(11,:); tapr(11,:); tmay(11,:); tjun(11,:); tjul(11,:); taug(11,:); tsep(11,:); toct(11,:); tnov(11,:); tdec(11,:)];

hour11 = [tjan(12,:); tfeb(12,:); tmar(12,:); tapr(12,:); tmay(12,:); tjun(12,:); tjul(12,:); taug(12,:); tsep(12,:); toct(12,:); tnov(12,:); tdec(12,:)];

hour12 = [tjan(13,:); tfeb(13,:); tmar(13,:); tapr(13,:); tmay(13,:); tjun(13,:); tjul(13,:); taug(13,:); tsep(13,:); toct(13,:); tnov(13,:); tdec(13,:)];

hour13 = [tjan(14,:); tfeb(14,:); tmar(14,:); tapr(14,:); tmay(14,:); tjun(14,:); tjul(14,:); taug(14,:); tsep(14,:); toct(14,:); tnov(14,:); tdec(14,:)];

hour14 = [tjan(15,:); tfeb(15,:); tmar(15,:); tapr(15,:); tmay(15,:); tjun(15,:); tjul(15,:); taug(15,:); tsep(15,:); toct(15,:); tnov(15,:); tdec(15,:)];

hour15 = [tjan(16,:); tfeb(16,:); tmar(16,:); tapr(16,:); tmay(16,:); tjun(16,:); tjul(16,:); taug(16,:); tsep(16,:); toct(16,:); tnov(16,:); tdec(16,:)];

hour16 = [tjan(17,:); tfeb(17,:); tmar(17,:); tapr(17,:); tmay(17,:); tjun(17,:); tjul(17,:); taug(17,:); tsep(17,:); toct(17,:); tnov(17,:); tdec(17,:)];

hour17 = [tjan(18,:); tfeb(18,:); tmar(18,:); tapr(18,:); tmay(18,:); tjun(18,:); tjul(18,:); taug(18,:); tsep(18,:); toct(18,:); tnov(18,:); tdec(18,:)];

hour18 = [tjan(19,:); tfeb(19,:); tmar(19,:); tapr(19,:); tmay(19,:); tjun(19,:); tjul(19,:); taug(19,:); tsep(19,:); toct(19,:); tnov(19,:); tdec(19,:)];

hour19 = [tjan(20,:); tfeb(20,:); tmar(20,:); tapr(20,:); tmay(20,:); tjun(20,:); tjul(20,:); taug(20,:); tsep(20,:); toct(20,:); tnov(20,:); tdec(20,:)];

hour20 = [tjan(21,:); tfeb(21,:); tmar(21,:); tapr(21,:); tmay(21,:); tjun(21,:); tjul(21,:); taug(21,:); tsep(21,:); toct(21,:); tnov(21,:); tdec(21,:)];

hour21 = [tjan(22,:); tfeb(22,:); tmar(22,:); tapr(22,:); tmay(22,:); tjun(22,:); tjul(22,:); taug(22,:); tsep(22,:); toct(22,:); tnov(22,:); tdec(22,:)];

hour22 = [tjan(23,:); tfeb(23,:); tmar(23,:); tapr(23,:); tmay(23,:); tjun(23,:); tjul(23,:); taug(23,:); tsep(23,:); toct(23,:); tnov(23,:); tdec(23,:)];

hour23 = [tjan(24,:); tfeb(24,:); tmar(24,:); tapr(24,:); tmay(24,:); tjun(24,:); tjul(24,:); taug(24,:); tsep(24,:); toct(24,:); tnov(24,:); tdec(24,:)];



hour0m = nanmean (hour0);
hour1m = nanmean (hour1);
hour2m = nanmean (hour2);
hour3m = nanmean (hour3);
hour4m = nanmean (hour4);
hour5m = nanmean (hour5);
hour6m = nanmean (hour6);
hour7m = nanmean (hour7);
hour8m = nanmean (hour8);
hour9m = nanmean (hour9);
hour10m = nanmean (hour10);
hour11m = nanmean (hour11);
hour12m = nanmean (hour12);
hour13m = nanmean (hour13);
hour14m = nanmean (hour14);
hour15m = nanmean (hour15);
hour16m = nanmean (hour16);
hour17m = nanmean (hour17);
hour18m = nanmean (hour18);
hour19m = nanmean (hour19);
hour20m = nanmean (hour20);
hour21m = nanmean (hour21);
hour22m = nanmean (hour22);
hour23m = nanmean (hour23);


gen = [hour0m;hour1m;hour2m;hour3m;hour4m;hour5m;hour6m;hour7m;hour8m;hour9m;hour10m;hour11m;hour12m;hour13m;hour14m;hour15m;hour16m;hour17m;hour18m;hour19m;hour20m;hour21m;hour22m;hour23m];

general = gen (:,3:end);

save ('/Users/sarahstuart/Desktop/MK/validated/temperature/diurnal cycle/general','general');