load 'uh82janT'
load 'uh83janT'
load 'uh84janT'
load 'uh85janT'
load 'uh86janT'
load 'uh87janT'
load 'uh88janT'
load 'uh89janT'
load 'uh90janT'
load 'uh91janT'
load 'uh92janT'
load 'uh93janT'
load 'uh94janT'
load 'uh95janT'
load 'uh96janT'
load 'uh97janT'
load 'uh98janT'
load 'uh99janT'
load 'uh00janT'
load 'uh01janT'
load 'uh02janT'
load 'uh03janT'
load 'uh07janT'
load 'uh08janT'
load 'uh09janT'
load 'uh10janT'


maxjan = [uh82janT;uh83janT;uh84janT;uh85janT;uh86janT;uh87janT;uh88janT;uh89janT;uh90janT;uh91janT;uh92janT;uh93janT;uh94janT;uh95janT;uh96janT;uh97janT;uh98janT;uh99janT;uh00janT;uh01janT;uh02janT;uh03janT;uh07janT;uh08janT;uh09janT;uh10janT];

[jan0 jan1 jan2 jan3 jan4 jan5 jan6 jan7 jan8 jan9 jan10 jan11 jan12 jan13 jan14 jan15 jan16 jan17 jan18 jan19 jan20 jan21 jan22 jan23] = extracthours (maxjan);

mjan0 = nanmean (jan0);
mjan1 = nanmean (jan1);
mjan2 = nanmean (jan2);
mjan3 = nanmean (jan3);
mjan4 = nanmean (jan4);
mjan5 = nanmean (jan5);
mjan6 = nanmean (jan6);
mjan7 = nanmean (jan7);
mjan8 = nanmean (jan8);
mjan9 = nanmean (jan9);
mjan10 = nanmean (jan10);
mjan11 = nanmean (jan11);
mjan12 = nanmean (jan12);
mjan13 = nanmean (jan13);
mjan14 = nanmean (jan14);
mjan15 = nanmean (jan15);
mjan16 = nanmean (jan16);
mjan17 = nanmean (jan17);
mjan18 = nanmean (jan18);
mjan19 = nanmean (jan19);
mjan20 = nanmean (jan20);
mjan21 = nanmean (jan21);
mjan22 = nanmean (jan22);
mjan23 = nanmean (jan23);

maxjan0 = nanmax (jan0);
maxjan1 = nanmax (jan1);
maxjan2 = nanmax (jan2);
maxjan3 = nanmax (jan3);
maxjan4 = nanmax (jan4);
maxjan5 = nanmax (jan5);
maxjan6 = nanmax (jan6);
maxjan7 = nanmax (jan7);
maxjan8 = nanmax (jan8);
maxjan9 = nanmax (jan9);
maxjan10 = nanmax (jan10);
maxjan11 = nanmax (jan11);
maxjan12 = nanmax (jan12);
maxjan13 = nanmax (jan13);
maxjan14 = nanmax (jan14);
maxjan15 = nanmax (jan15);
maxjan16 = nanmax (jan16);
maxjan17 = nanmax (jan17);
maxjan18 = nanmax (jan18);
maxjan19 = nanmax (jan19);
maxjan20 = nanmax (jan20);
maxjan21 = nanmax (jan21);
maxjan22 = nanmax (jan22);
maxjan23 = nanmax (jan23);

minjan0 = nanmin (jan0);
minjan1 = nanmin (jan1);
minjan2 = nanmin (jan2);
minjan3 = nanmin (jan3);
minjan4 = nanmin (jan4);
minjan5 = nanmin (jan5);
minjan6 = nanmin (jan6);
minjan7 = nanmin (jan7);
minjan8 = nanmin (jan8);
minjan9 = nanmin (jan9);
minjan10 = nanmin (jan10);
minjan11 = nanmin (jan11);
minjan12 = nanmin (jan12);
minjan13 = nanmin (jan13);
minjan14 = nanmin (jan14);
minjan15 = nanmin (jan15);
minjan16 = nanmin (jan16);
minjan17 = nanmin (jan17);
minjan18 = nanmin (jan18);
minjan19 = nanmin (jan19);
minjan20 = nanmin (jan20);
minjan21 = nanmin (jan21);
minjan22 = nanmin (jan22);
minjan23 = nanmin (jan23);


tjan0 = [mjan0,maxjan0(1,4),minjan0(1,4)];
tjan1 = [mjan1,maxjan1(1,4),minjan1(1,4)];
tjan2 = [mjan2,maxjan2(1,4),minjan2(1,4)];
tjan3 = [mjan3,maxjan3(1,4),minjan3(1,4)];
tjan4 = [mjan4,maxjan4(1,4),minjan4(1,4)];
tjan5 = [mjan5,maxjan5(1,4),minjan5(1,4)];
tjan6 = [mjan6,maxjan6(1,4),minjan6(1,4)];
tjan7 = [mjan7,maxjan7(1,4),minjan7(1,4)];
tjan8 = [mjan8,maxjan8(1,4),minjan8(1,4)];
tjan9 = [mjan9,maxjan9(1,4),minjan9(1,4)];
tjan10 = [mjan10,maxjan10(1,4),minjan10(1,4)];
tjan11 = [mjan11,maxjan11(1,4),minjan11(1,4)];
tjan12 = [mjan12,maxjan12(1,4),minjan12(1,4)];
tjan13 = [mjan13,maxjan13(1,4),minjan13(1,4)];
tjan14 = [mjan14,maxjan14(1,4),minjan14(1,4)];
tjan15 = [mjan15,maxjan15(1,4),minjan15(1,4)];
tjan16 = [mjan16,maxjan16(1,4),minjan16(1,4)];
tjan17 = [mjan17,maxjan17(1,4),minjan17(1,4)];
tjan18 = [mjan18,maxjan18(1,4),minjan18(1,4)];
tjan19 = [mjan19,maxjan19(1,4),minjan19(1,4)];
tjan20 = [mjan20,maxjan20(1,4),minjan20(1,4)];
tjan21 = [mjan21,maxjan21(1,4),minjan21(1,4)];
tjan22 = [mjan22,maxjan22(1,4),minjan22(1,4)];
tjan23 = [mjan23,maxjan23(1,4),minjan23(1,4)];

tjan = [tjan0;tjan1;tjan2;tjan3;tjan4;tjan5;tjan6;tjan7;tjan8;tjan9;tjan10;tjan11;tjan12;tjan13;tjan14;tjan15;tjan16;tjan17;tjan18;tjan19;tjan20;tjan21;tjan22;tjan23];


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

load 'uh82febT'
load 'uh83febT'
load 'uh84febT'
load 'uh85febT'
load 'uh86febT'
load 'uh87febT'
load 'uh88febT'
load 'uh89febT'
load 'uh90febT'
load 'uh91febT'
load 'uh92febT'
load 'uh93febT'
load 'uh94febT'
load 'uh95febT'
load 'uh96febT'
load 'uh97febT'
load 'uh98febT'
load 'uh99febT'
load 'uh00febT'
load 'uh01febT'
load 'uh02febT'
load 'uh03febT'
load 'uh07febT'
load 'uh08febT'
load 'uh09febT'
load 'uh10febT'


maxfeb = [uh82febT;uh83febT;uh84febT;uh85febT;uh86febT;uh87febT;uh88febT;uh89febT;uh90febT;uh91febT;uh92febT;uh93febT;uh94febT;uh95febT;uh96febT;uh97febT;uh98febT;uh99febT;uh00febT;uh01febT;uh02febT;uh03febT;uh07febT;uh08febT;uh09febT;uh10febT];

[feb0 feb1 feb2 feb3 feb4 feb5 feb6 feb7 feb8 feb9 feb10 feb11 feb12 feb13 feb14 feb15 feb16 feb17 feb18 feb19 feb20 feb21 feb22 feb23] = extracthours (maxfeb);

mfeb0 = nanmean (feb0);
mfeb1 = nanmean (feb1);
mfeb2 = nanmean (feb2);
mfeb3 = nanmean (feb3);
mfeb4 = nanmean (feb4);
mfeb5 = nanmean (feb5);
mfeb6 = nanmean (feb6);
mfeb7 = nanmean (feb7);
mfeb8 = nanmean (feb8);
mfeb9 = nanmean (feb9);
mfeb10 = nanmean (feb10);
mfeb11 = nanmean (feb11);
mfeb12 = nanmean (feb12);
mfeb13 = nanmean (feb13);
mfeb14 = nanmean (feb14);
mfeb15 = nanmean (feb15);
mfeb16 = nanmean (feb16);
mfeb17 = nanmean (feb17);
mfeb18 = nanmean (feb18);
mfeb19 = nanmean (feb19);
mfeb20 = nanmean (feb20);
mfeb21 = nanmean (feb21);
mfeb22 = nanmean (feb22);
mfeb23 = nanmean (feb23);

maxfeb0 = nanmax (feb0);
maxfeb1 = nanmax (feb1);
maxfeb2 = nanmax (feb2);
maxfeb3 = nanmax (feb3);
maxfeb4 = nanmax (feb4);
maxfeb5 = nanmax (feb5);
maxfeb6 = nanmax (feb6);
maxfeb7 = nanmax (feb7);
maxfeb8 = nanmax (feb8);
maxfeb9 = nanmax (feb9);
maxfeb10 = nanmax (feb10);
maxfeb11 = nanmax (feb11);
maxfeb12 = nanmax (feb12);
maxfeb13 = nanmax (feb13);
maxfeb14 = nanmax (feb14);
maxfeb15 = nanmax (feb15);
maxfeb16 = nanmax (feb16);
maxfeb17 = nanmax (feb17);
maxfeb18 = nanmax (feb18);
maxfeb19 = nanmax (feb19);
maxfeb20 = nanmax (feb20);
maxfeb21 = nanmax (feb21);
maxfeb22 = nanmax (feb22);
maxfeb23 = nanmax (feb23);

minfeb0 = nanmin (feb0);
minfeb1 = nanmin (feb1);
minfeb2 = nanmin (feb2);
minfeb3 = nanmin (feb3);
minfeb4 = nanmin (feb4);
minfeb5 = nanmin (feb5);
minfeb6 = nanmin (feb6);
minfeb7 = nanmin (feb7);
minfeb8 = nanmin (feb8);
minfeb9 = nanmin (feb9);
minfeb10 = nanmin (feb10);
minfeb11 = nanmin (feb11);
minfeb12 = nanmin (feb12);
minfeb13 = nanmin (feb13);
minfeb14 = nanmin (feb14);
minfeb15 = nanmin (feb15);
minfeb16 = nanmin (feb16);
minfeb17 = nanmin (feb17);
minfeb18 = nanmin (feb18);
minfeb19 = nanmin (feb19);
minfeb20 = nanmin (feb20);
minfeb21 = nanmin (feb21);
minfeb22 = nanmin (feb22);
minfeb23 = nanmin (feb23);


tfeb0 = [mfeb0,maxfeb0(1,4),minfeb0(1,4)];
tfeb1 = [mfeb1,maxfeb1(1,4),minfeb1(1,4)];
tfeb2 = [mfeb2,maxfeb2(1,4),minfeb2(1,4)];
tfeb3 = [mfeb3,maxfeb3(1,4),minfeb3(1,4)];
tfeb4 = [mfeb4,maxfeb4(1,4),minfeb4(1,4)];
tfeb5 = [mfeb5,maxfeb5(1,4),minfeb5(1,4)];
tfeb6 = [mfeb6,maxfeb6(1,4),minfeb6(1,4)];
tfeb7 = [mfeb7,maxfeb7(1,4),minfeb7(1,4)];
tfeb8 = [mfeb8,maxfeb8(1,4),minfeb8(1,4)];
tfeb9 = [mfeb9,maxfeb9(1,4),minfeb9(1,4)];
tfeb10 = [mfeb10,maxfeb10(1,4),minfeb10(1,4)];
tfeb11 = [mfeb11,maxfeb11(1,4),minfeb11(1,4)];
tfeb12 = [mfeb12,maxfeb12(1,4),minfeb12(1,4)];
tfeb13 = [mfeb13,maxfeb13(1,4),minfeb13(1,4)];
tfeb14 = [mfeb14,maxfeb14(1,4),minfeb14(1,4)];
tfeb15 = [mfeb15,maxfeb15(1,4),minfeb15(1,4)];
tfeb16 = [mfeb16,maxfeb16(1,4),minfeb16(1,4)];
tfeb17 = [mfeb17,maxfeb17(1,4),minfeb17(1,4)];
tfeb18 = [mfeb18,maxfeb18(1,4),minfeb18(1,4)];
tfeb19 = [mfeb19,maxfeb19(1,4),minfeb19(1,4)];
tfeb20 = [mfeb20,maxfeb20(1,4),minfeb20(1,4)];
tfeb21 = [mfeb21,maxfeb21(1,4),minfeb21(1,4)];
tfeb22 = [mfeb22,maxfeb22(1,4),minfeb22(1,4)];
tfeb23 = [mfeb23,maxfeb23(1,4),minfeb23(1,4)];

tfeb = [tfeb0;tfeb1;tfeb2;tfeb3;tfeb4;tfeb5;tfeb6;tfeb7;tfeb8;tfeb9;tfeb10;tfeb11;tfeb12;tfeb13;tfeb14;tfeb15;tfeb16;tfeb17;tfeb18;tfeb19;tfeb20;tfeb21;tfeb22;tfeb23];


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



load 'uh82marT'
load 'uh83marT'
load 'uh84marT'
load 'uh85marT'
load 'uh86marT'
load 'uh87marT'
load 'uh88marT'
load 'uh89marT'
load 'uh90marT'
load 'uh91marT'
load 'uh92marT'
load 'uh93marT'
load 'uh94marT'
load 'uh95marT'
load 'uh96marT'
load 'uh97marT'
load 'uh98marT'
load 'uh99marT'
load 'uh00marT'
load 'uh01marT'
load 'uh02marT'
load 'uh03marT'
load 'uh07marT'
load 'uh08marT'
load 'uh09marT'
load 'uh10marT'


maxmar = [uh82marT;uh83marT;uh84marT;uh85marT;uh86marT;uh87marT;uh88marT;uh89marT;uh90marT;uh91marT;uh92marT;uh93marT;uh94marT;uh95marT;uh96marT;uh97marT;uh98marT;uh99marT;uh00marT;uh01marT;uh02marT;uh03marT;uh07marT;uh08marT;uh09marT;uh10marT];

[mar0 mar1 mar2 mar3 mar4 mar5 mar6 mar7 mar8 mar9 mar10 mar11 mar12 mar13 mar14 mar15 mar16 mar17 mar18 mar19 mar20 mar21 mar22 mar23] = extracthours (maxmar);

mmar0 = nanmean (mar0);
mmar1 = nanmean (mar1);
mmar2 = nanmean (mar2);
mmar3 = nanmean (mar3);
mmar4 = nanmean (mar4);
mmar5 = nanmean (mar5);
mmar6 = nanmean (mar6);
mmar7 = nanmean (mar7);
mmar8 = nanmean (mar8);
mmar9 = nanmean (mar9);
mmar10 = nanmean (mar10);
mmar11 = nanmean (mar11);
mmar12 = nanmean (mar12);
mmar13 = nanmean (mar13);
mmar14 = nanmean (mar14);
mmar15 = nanmean (mar15);
mmar16 = nanmean (mar16);
mmar17 = nanmean (mar17);
mmar18 = nanmean (mar18);
mmar19 = nanmean (mar19);
mmar20 = nanmean (mar20);
mmar21 = nanmean (mar21);
mmar22 = nanmean (mar22);
mmar23 = nanmean (mar23);

maxmar0 = nanmax (mar0);
maxmar1 = nanmax (mar1);
maxmar2 = nanmax (mar2);
maxmar3 = nanmax (mar3);
maxmar4 = nanmax (mar4);
maxmar5 = nanmax (mar5);
maxmar6 = nanmax (mar6);
maxmar7 = nanmax (mar7);
maxmar8 = nanmax (mar8);
maxmar9 = nanmax (mar9);
maxmar10 = nanmax (mar10);
maxmar11 = nanmax (mar11);
maxmar12 = nanmax (mar12);
maxmar13 = nanmax (mar13);
maxmar14 = nanmax (mar14);
maxmar15 = nanmax (mar15);
maxmar16 = nanmax (mar16);
maxmar17 = nanmax (mar17);
maxmar18 = nanmax (mar18);
maxmar19 = nanmax (mar19);
maxmar20 = nanmax (mar20);
maxmar21 = nanmax (mar21);
maxmar22 = nanmax (mar22);
maxmar23 = nanmax (mar23);

minmar0 = nanmin (mar0);
minmar1 = nanmin (mar1);
minmar2 = nanmin (mar2);
minmar3 = nanmin (mar3);
minmar4 = nanmin (mar4);
minmar5 = nanmin (mar5);
minmar6 = nanmin (mar6);
minmar7 = nanmin (mar7);
minmar8 = nanmin (mar8);
minmar9 = nanmin (mar9);
minmar10 = nanmin (mar10);
minmar11 = nanmin (mar11);
minmar12 = nanmin (mar12);
minmar13 = nanmin (mar13);
minmar14 = nanmin (mar14);
minmar15 = nanmin (mar15);
minmar16 = nanmin (mar16);
minmar17 = nanmin (mar17);
minmar18 = nanmin (mar18);
minmar19 = nanmin (mar19);
minmar20 = nanmin (mar20);
minmar21 = nanmin (mar21);
minmar22 = nanmin (mar22);
minmar23 = nanmin (mar23);


tmar0 = [mmar0,maxmar0(1,4),minmar0(1,4)];
tmar1 = [mmar1,maxmar1(1,4),minmar1(1,4)];
tmar2 = [mmar2,maxmar2(1,4),minmar2(1,4)];
tmar3 = [mmar3,maxmar3(1,4),minmar3(1,4)];
tmar4 = [mmar4,maxmar4(1,4),minmar4(1,4)];
tmar5 = [mmar5,maxmar5(1,4),minmar5(1,4)];
tmar6 = [mmar6,maxmar6(1,4),minmar6(1,4)];
tmar7 = [mmar7,maxmar7(1,4),minmar7(1,4)];
tmar8 = [mmar8,maxmar8(1,4),minmar8(1,4)];
tmar9 = [mmar9,maxmar9(1,4),minmar9(1,4)];
tmar10 = [mmar10,maxmar10(1,4),minmar10(1,4)];
tmar11 = [mmar11,maxmar11(1,4),minmar11(1,4)];
tmar12 = [mmar12,maxmar12(1,4),minmar12(1,4)];
tmar13 = [mmar13,maxmar13(1,4),minmar13(1,4)];
tmar14 = [mmar14,maxmar14(1,4),minmar14(1,4)];
tmar15 = [mmar15,maxmar15(1,4),minmar15(1,4)];
tmar16 = [mmar16,maxmar16(1,4),minmar16(1,4)];
tmar17 = [mmar17,maxmar17(1,4),minmar17(1,4)];
tmar18 = [mmar18,maxmar18(1,4),minmar18(1,4)];
tmar19 = [mmar19,maxmar19(1,4),minmar19(1,4)];
tmar20 = [mmar20,maxmar20(1,4),minmar20(1,4)];
tmar21 = [mmar21,maxmar21(1,4),minmar21(1,4)];
tmar22 = [mmar22,maxmar22(1,4),minmar22(1,4)];
tmar23 = [mmar23,maxmar23(1,4),minmar23(1,4)];

tmar = [tmar0;tmar1;tmar2;tmar3;tmar4;tmar5;tmar6;tmar7;tmar8;tmar9;tmar10;tmar11;tmar12;tmar13;tmar14;tmar15;tmar16;tmar17;tmar18;tmar19;tmar20;tmar21;tmar22;tmar23];


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

load 'uh82aprT'
load 'uh83aprT'
load 'uh84aprT'
load 'uh85aprT'
load 'uh86aprT'
load 'uh87aprT'
load 'uh88aprT'
load 'uh89aprT'
load 'uh90aprT'
load 'uh91aprT'
load 'uh92aprT'
load 'uh93aprT'
load 'uh94aprT'
load 'uh95aprT'
load 'uh96aprT'
load 'uh97aprT'
load 'uh98aprT'
load 'uh99aprT'
load 'uh00aprT'
load 'uh01aprT'
load 'uh02aprT'
load 'uh03aprT'
load 'uh07aprT'
load 'uh08aprT'
load 'uh09aprT'
load 'uh10aprT'


maxapr = [uh82aprT;uh83aprT;uh84aprT;uh85aprT;uh86aprT;uh87aprT;uh88aprT;uh89aprT;uh90aprT;uh91aprT;uh92aprT;uh93aprT;uh94aprT;uh95aprT;uh96aprT;uh97aprT;uh98aprT;uh99aprT;uh00aprT;uh01aprT;uh02aprT;uh03aprT;uh07aprT;uh08aprT;uh09aprT;uh10aprT];

[apr0 apr1 apr2 apr3 apr4 apr5 apr6 apr7 apr8 apr9 apr10 apr11 apr12 apr13 apr14 apr15 apr16 apr17 apr18 apr19 apr20 apr21 apr22 apr23] = extracthours (maxapr);

mapr0 = nanmean (apr0);
mapr1 = nanmean (apr1);
mapr2 = nanmean (apr2);
mapr3 = nanmean (apr3);
mapr4 = nanmean (apr4);
mapr5 = nanmean (apr5);
mapr6 = nanmean (apr6);
mapr7 = nanmean (apr7);
mapr8 = nanmean (apr8);
mapr9 = nanmean (apr9);
mapr10 = nanmean (apr10);
mapr11 = nanmean (apr11);
mapr12 = nanmean (apr12);
mapr13 = nanmean (apr13);
mapr14 = nanmean (apr14);
mapr15 = nanmean (apr15);
mapr16 = nanmean (apr16);
mapr17 = nanmean (apr17);
mapr18 = nanmean (apr18);
mapr19 = nanmean (apr19);
mapr20 = nanmean (apr20);
mapr21 = nanmean (apr21);
mapr22 = nanmean (apr22);
mapr23 = nanmean (apr23);

maxapr0 = nanmax (apr0);
maxapr1 = nanmax (apr1);
maxapr2 = nanmax (apr2);
maxapr3 = nanmax (apr3);
maxapr4 = nanmax (apr4);
maxapr5 = nanmax (apr5);
maxapr6 = nanmax (apr6);
maxapr7 = nanmax (apr7);
maxapr8 = nanmax (apr8);
maxapr9 = nanmax (apr9);
maxapr10 = nanmax (apr10);
maxapr11 = nanmax (apr11);
maxapr12 = nanmax (apr12);
maxapr13 = nanmax (apr13);
maxapr14 = nanmax (apr14);
maxapr15 = nanmax (apr15);
maxapr16 = nanmax (apr16);
maxapr17 = nanmax (apr17);
maxapr18 = nanmax (apr18);
maxapr19 = nanmax (apr19);
maxapr20 = nanmax (apr20);
maxapr21 = nanmax (apr21);
maxapr22 = nanmax (apr22);
maxapr23 = nanmax (apr23);

minapr0 = nanmin (apr0);
minapr1 = nanmin (apr1);
minapr2 = nanmin (apr2);
minapr3 = nanmin (apr3);
minapr4 = nanmin (apr4);
minapr5 = nanmin (apr5);
minapr6 = nanmin (apr6);
minapr7 = nanmin (apr7);
minapr8 = nanmin (apr8);
minapr9 = nanmin (apr9);
minapr10 = nanmin (apr10);
minapr11 = nanmin (apr11);
minapr12 = nanmin (apr12);
minapr13 = nanmin (apr13);
minapr14 = nanmin (apr14);
minapr15 = nanmin (apr15);
minapr16 = nanmin (apr16);
minapr17 = nanmin (apr17);
minapr18 = nanmin (apr18);
minapr19 = nanmin (apr19);
minapr20 = nanmin (apr20);
minapr21 = nanmin (apr21);
minapr22 = nanmin (apr22);
minapr23 = nanmin (apr23);


tapr0 = [mapr0,maxapr0(1,4),minapr0(1,4)];
tapr1 = [mapr1,maxapr1(1,4),minapr1(1,4)];
tapr2 = [mapr2,maxapr2(1,4),minapr2(1,4)];
tapr3 = [mapr3,maxapr3(1,4),minapr3(1,4)];
tapr4 = [mapr4,maxapr4(1,4),minapr4(1,4)];
tapr5 = [mapr5,maxapr5(1,4),minapr5(1,4)];
tapr6 = [mapr6,maxapr6(1,4),minapr6(1,4)];
tapr7 = [mapr7,maxapr7(1,4),minapr7(1,4)];
tapr8 = [mapr8,maxapr8(1,4),minapr8(1,4)];
tapr9 = [mapr9,maxapr9(1,4),minapr9(1,4)];
tapr10 = [mapr10,maxapr10(1,4),minapr10(1,4)];
tapr11 = [mapr11,maxapr11(1,4),minapr11(1,4)];
tapr12 = [mapr12,maxapr12(1,4),minapr12(1,4)];
tapr13 = [mapr13,maxapr13(1,4),minapr13(1,4)];
tapr14 = [mapr14,maxapr14(1,4),minapr14(1,4)];
tapr15 = [mapr15,maxapr15(1,4),minapr15(1,4)];
tapr16 = [mapr16,maxapr16(1,4),minapr16(1,4)];
tapr17 = [mapr17,maxapr17(1,4),minapr17(1,4)];
tapr18 = [mapr18,maxapr18(1,4),minapr18(1,4)];
tapr19 = [mapr19,maxapr19(1,4),minapr19(1,4)];
tapr20 = [mapr20,maxapr20(1,4),minapr20(1,4)];
tapr21 = [mapr21,maxapr21(1,4),minapr21(1,4)];
tapr22 = [mapr22,maxapr22(1,4),minapr22(1,4)];
tapr23 = [mapr23,maxapr23(1,4),minapr23(1,4)];

tapr = [tapr0;tapr1;tapr2;tapr3;tapr4;tapr5;tapr6;tapr7;tapr8;tapr9;tapr10;tapr11;tapr12;tapr13;tapr14;tapr15;tapr16;tapr17;tapr18;tapr19;tapr20;tapr21;tapr22;tapr23];


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


load 'uh82mayT'
load 'uh83mayT'
load 'uh84mayT'
load 'uh85mayT'
load 'uh86mayT'
load 'uh87mayT'
load 'uh88mayT'
load 'uh89mayT'
load 'uh90mayT'
load 'uh91mayT'
load 'uh92mayT'
load 'uh93mayT'
load 'uh94mayT'
load 'uh95mayT'
load 'uh96mayT'
load 'uh97mayT'
load 'uh98mayT'
load 'uh99mayT'
load 'uh00mayT'
load 'uh01mayT'
load 'uh02mayT'
load 'uh03mayT'
load 'uh07mayT'
load 'uh08mayT'
load 'uh09mayT'
load 'uh10mayT'


maxmay = [uh82mayT;uh83mayT;uh84mayT;uh85mayT;uh86mayT;uh87mayT;uh88mayT;uh89mayT;uh90mayT;uh91mayT;uh92mayT;uh93mayT;uh94mayT;uh95mayT;uh96mayT;uh97mayT;uh98mayT;uh99mayT;uh00mayT;uh01mayT;uh02mayT;uh03mayT;uh07mayT;uh08mayT;uh09mayT;uh10mayT];

[may0 may1 may2 may3 may4 may5 may6 may7 may8 may9 may10 may11 may12 may13 may14 may15 may16 may17 may18 may19 may20 may21 may22 may23] = extracthours (maxmay);

mmay0 = nanmean (may0);
mmay1 = nanmean (may1);
mmay2 = nanmean (may2);
mmay3 = nanmean (may3);
mmay4 = nanmean (may4);
mmay5 = nanmean (may5);
mmay6 = nanmean (may6);
mmay7 = nanmean (may7);
mmay8 = nanmean (may8);
mmay9 = nanmean (may9);
mmay10 = nanmean (may10);
mmay11 = nanmean (may11);
mmay12 = nanmean (may12);
mmay13 = nanmean (may13);
mmay14 = nanmean (may14);
mmay15 = nanmean (may15);
mmay16 = nanmean (may16);
mmay17 = nanmean (may17);
mmay18 = nanmean (may18);
mmay19 = nanmean (may19);
mmay20 = nanmean (may20);
mmay21 = nanmean (may21);
mmay22 = nanmean (may22);
mmay23 = nanmean (may23);

maxmay0 = nanmax (may0);
maxmay1 = nanmax (may1);
maxmay2 = nanmax (may2);
maxmay3 = nanmax (may3);
maxmay4 = nanmax (may4);
maxmay5 = nanmax (may5);
maxmay6 = nanmax (may6);
maxmay7 = nanmax (may7);
maxmay8 = nanmax (may8);
maxmay9 = nanmax (may9);
maxmay10 = nanmax (may10);
maxmay11 = nanmax (may11);
maxmay12 = nanmax (may12);
maxmay13 = nanmax (may13);
maxmay14 = nanmax (may14);
maxmay15 = nanmax (may15);
maxmay16 = nanmax (may16);
maxmay17 = nanmax (may17);
maxmay18 = nanmax (may18);
maxmay19 = nanmax (may19);
maxmay20 = nanmax (may20);
maxmay21 = nanmax (may21);
maxmay22 = nanmax (may22);
maxmay23 = nanmax (may23);

minmay0 = nanmin (may0);
minmay1 = nanmin (may1);
minmay2 = nanmin (may2);
minmay3 = nanmin (may3);
minmay4 = nanmin (may4);
minmay5 = nanmin (may5);
minmay6 = nanmin (may6);
minmay7 = nanmin (may7);
minmay8 = nanmin (may8);
minmay9 = nanmin (may9);
minmay10 = nanmin (may10);
minmay11 = nanmin (may11);
minmay12 = nanmin (may12);
minmay13 = nanmin (may13);
minmay14 = nanmin (may14);
minmay15 = nanmin (may15);
minmay16 = nanmin (may16);
minmay17 = nanmin (may17);
minmay18 = nanmin (may18);
minmay19 = nanmin (may19);
minmay20 = nanmin (may20);
minmay21 = nanmin (may21);
minmay22 = nanmin (may22);
minmay23 = nanmin (may23);


tmay0 = [mmay0,maxmay0(1,4),minmay0(1,4)];
tmay1 = [mmay1,maxmay1(1,4),minmay1(1,4)];
tmay2 = [mmay2,maxmay2(1,4),minmay2(1,4)];
tmay3 = [mmay3,maxmay3(1,4),minmay3(1,4)];
tmay4 = [mmay4,maxmay4(1,4),minmay4(1,4)];
tmay5 = [mmay5,maxmay5(1,4),minmay5(1,4)];
tmay6 = [mmay6,maxmay6(1,4),minmay6(1,4)];
tmay7 = [mmay7,maxmay7(1,4),minmay7(1,4)];
tmay8 = [mmay8,maxmay8(1,4),minmay8(1,4)];
tmay9 = [mmay9,maxmay9(1,4),minmay9(1,4)];
tmay10 = [mmay10,maxmay10(1,4),minmay10(1,4)];
tmay11 = [mmay11,maxmay11(1,4),minmay11(1,4)];
tmay12 = [mmay12,maxmay12(1,4),minmay12(1,4)];
tmay13 = [mmay13,maxmay13(1,4),minmay13(1,4)];
tmay14 = [mmay14,maxmay14(1,4),minmay14(1,4)];
tmay15 = [mmay15,maxmay15(1,4),minmay15(1,4)];
tmay16 = [mmay16,maxmay16(1,4),minmay16(1,4)];
tmay17 = [mmay17,maxmay17(1,4),minmay17(1,4)];
tmay18 = [mmay18,maxmay18(1,4),minmay18(1,4)];
tmay19 = [mmay19,maxmay19(1,4),minmay19(1,4)];
tmay20 = [mmay20,maxmay20(1,4),minmay20(1,4)];
tmay21 = [mmay21,maxmay21(1,4),minmay21(1,4)];
tmay22 = [mmay22,maxmay22(1,4),minmay22(1,4)];
tmay23 = [mmay23,maxmay23(1,4),minmay23(1,4)];

tmay = [tmay0;tmay1;tmay2;tmay3;tmay4;tmay5;tmay6;tmay7;tmay8;tmay9;tmay10;tmay11;tmay12;tmay13;tmay14;tmay15;tmay16;tmay17;tmay18;tmay19;tmay20;tmay21;tmay22;tmay23];


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


load 'uh82junT'
load 'uh83junT'
load 'uh84junT'
load 'uh85junT'
load 'uh86junT'
load 'uh87junT'
load 'uh88junT'
load 'uh89junT'
load 'uh90junT'
load 'uh91junT'
load 'uh92junT'
load 'uh93junT'
load 'uh94junT'
load 'uh95junT'
load 'uh96junT'
load 'uh97junT'
load 'uh98junT'
load 'uh99junT'
load 'uh00junT'
load 'uh01junT'
load 'uh02junT'
load 'uh03junT'
load 'uh07junT'
load 'uh08junT'
load 'uh09junT'
load 'uh10junT'


maxjun = [uh82junT;uh83junT;uh84junT;uh85junT;uh86junT;uh87junT;uh88junT;uh89junT;uh90junT;uh91junT;uh92junT;uh93junT;uh94junT;uh95junT;uh96junT;uh97junT;uh98junT;uh99junT;uh00junT;uh01junT;uh02junT;uh03junT;uh07junT;uh08junT;uh09junT;uh10junT];

[jun0 jun1 jun2 jun3 jun4 jun5 jun6 jun7 jun8 jun9 jun10 jun11 jun12 jun13 jun14 jun15 jun16 jun17 jun18 jun19 jun20 jun21 jun22 jun23] = extracthours (maxjun);

mjun0 = nanmean (jun0);
mjun1 = nanmean (jun1);
mjun2 = nanmean (jun2);
mjun3 = nanmean (jun3);
mjun4 = nanmean (jun4);
mjun5 = nanmean (jun5);
mjun6 = nanmean (jun6);
mjun7 = nanmean (jun7);
mjun8 = nanmean (jun8);
mjun9 = nanmean (jun9);
mjun10 = nanmean (jun10);
mjun11 = nanmean (jun11);
mjun12 = nanmean (jun12);
mjun13 = nanmean (jun13);
mjun14 = nanmean (jun14);
mjun15 = nanmean (jun15);
mjun16 = nanmean (jun16);
mjun17 = nanmean (jun17);
mjun18 = nanmean (jun18);
mjun19 = nanmean (jun19);
mjun20 = nanmean (jun20);
mjun21 = nanmean (jun21);
mjun22 = nanmean (jun22);
mjun23 = nanmean (jun23);

maxjun0 = nanmax (jun0);
maxjun1 = nanmax (jun1);
maxjun2 = nanmax (jun2);
maxjun3 = nanmax (jun3);
maxjun4 = nanmax (jun4);
maxjun5 = nanmax (jun5);
maxjun6 = nanmax (jun6);
maxjun7 = nanmax (jun7);
maxjun8 = nanmax (jun8);
maxjun9 = nanmax (jun9);
maxjun10 = nanmax (jun10);
maxjun11 = nanmax (jun11);
maxjun12 = nanmax (jun12);
maxjun13 = nanmax (jun13);
maxjun14 = nanmax (jun14);
maxjun15 = nanmax (jun15);
maxjun16 = nanmax (jun16);
maxjun17 = nanmax (jun17);
maxjun18 = nanmax (jun18);
maxjun19 = nanmax (jun19);
maxjun20 = nanmax (jun20);
maxjun21 = nanmax (jun21);
maxjun22 = nanmax (jun22);
maxjun23 = nanmax (jun23);

minjun0 = nanmin (jun0);
minjun1 = nanmin (jun1);
minjun2 = nanmin (jun2);
minjun3 = nanmin (jun3);
minjun4 = nanmin (jun4);
minjun5 = nanmin (jun5);
minjun6 = nanmin (jun6);
minjun7 = nanmin (jun7);
minjun8 = nanmin (jun8);
minjun9 = nanmin (jun9);
minjun10 = nanmin (jun10);
minjun11 = nanmin (jun11);
minjun12 = nanmin (jun12);
minjun13 = nanmin (jun13);
minjun14 = nanmin (jun14);
minjun15 = nanmin (jun15);
minjun16 = nanmin (jun16);
minjun17 = nanmin (jun17);
minjun18 = nanmin (jun18);
minjun19 = nanmin (jun19);
minjun20 = nanmin (jun20);
minjun21 = nanmin (jun21);
minjun22 = nanmin (jun22);
minjun23 = nanmin (jun23);


tjun0 = [mjun0,maxjun0(1,4),minjun0(1,4)];
tjun1 = [mjun1,maxjun1(1,4),minjun1(1,4)];
tjun2 = [mjun2,maxjun2(1,4),minjun2(1,4)];
tjun3 = [mjun3,maxjun3(1,4),minjun3(1,4)];
tjun4 = [mjun4,maxjun4(1,4),minjun4(1,4)];
tjun5 = [mjun5,maxjun5(1,4),minjun5(1,4)];
tjun6 = [mjun6,maxjun6(1,4),minjun6(1,4)];
tjun7 = [mjun7,maxjun7(1,4),minjun7(1,4)];
tjun8 = [mjun8,maxjun8(1,4),minjun8(1,4)];
tjun9 = [mjun9,maxjun9(1,4),minjun9(1,4)];
tjun10 = [mjun10,maxjun10(1,4),minjun10(1,4)];
tjun11 = [mjun11,maxjun11(1,4),minjun11(1,4)];
tjun12 = [mjun12,maxjun12(1,4),minjun12(1,4)];
tjun13 = [mjun13,maxjun13(1,4),minjun13(1,4)];
tjun14 = [mjun14,maxjun14(1,4),minjun14(1,4)];
tjun15 = [mjun15,maxjun15(1,4),minjun15(1,4)];
tjun16 = [mjun16,maxjun16(1,4),minjun16(1,4)];
tjun17 = [mjun17,maxjun17(1,4),minjun17(1,4)];
tjun18 = [mjun18,maxjun18(1,4),minjun18(1,4)];
tjun19 = [mjun19,maxjun19(1,4),minjun19(1,4)];
tjun20 = [mjun20,maxjun20(1,4),minjun20(1,4)];
tjun21 = [mjun21,maxjun21(1,4),minjun21(1,4)];
tjun22 = [mjun22,maxjun22(1,4),minjun22(1,4)];
tjun23 = [mjun23,maxjun23(1,4),minjun23(1,4)];

tjun = [tjun0;tjun1;tjun2;tjun3;tjun4;tjun5;tjun6;tjun7;tjun8;tjun9;tjun10;tjun11;tjun12;tjun13;tjun14;tjun15;tjun16;tjun17;tjun18;tjun19;tjun20;tjun21;tjun22;tjun23];


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




load 'uh82julT'
load 'uh83julT'
load 'uh84julT'
load 'uh85julT'
load 'uh86julT'
load 'uh87julT'
load 'uh88julT'
load 'uh89julT'
load 'uh90julT'
load 'uh91julT'
load 'uh92julT'
load 'uh93julT'
load 'uh94julT'
load 'uh95julT'
load 'uh96julT'
load 'uh97julT'
load 'uh98julT'
load 'uh99julT'
load 'uh00julT'
load 'uh01julT'
load 'uh02julT'
load 'uh03julT'
load 'uh07julT'
load 'uh08julT'
load 'uh09julT'
load 'uh10julT'


maxjul = [uh82julT;uh83julT;uh84julT;uh85julT;uh86julT;uh87julT;uh88julT;uh89julT;uh90julT;uh91julT;uh92julT;uh93julT;uh94julT;uh95julT;uh96julT;uh97julT;uh98julT;uh99julT;uh00julT;uh01julT;uh02julT;uh03julT;uh07julT;uh08julT;uh09julT;uh10julT];

[jul0 jul1 jul2 jul3 jul4 jul5 jul6 jul7 jul8 jul9 jul10 jul11 jul12 jul13 jul14 jul15 jul16 jul17 jul18 jul19 jul20 jul21 jul22 jul23] = extracthours (maxjul);

mjul0 = nanmean (jul0);
mjul1 = nanmean (jul1);
mjul2 = nanmean (jul2);
mjul3 = nanmean (jul3);
mjul4 = nanmean (jul4);
mjul5 = nanmean (jul5);
mjul6 = nanmean (jul6);
mjul7 = nanmean (jul7);
mjul8 = nanmean (jul8);
mjul9 = nanmean (jul9);
mjul10 = nanmean (jul10);
mjul11 = nanmean (jul11);
mjul12 = nanmean (jul12);
mjul13 = nanmean (jul13);
mjul14 = nanmean (jul14);
mjul15 = nanmean (jul15);
mjul16 = nanmean (jul16);
mjul17 = nanmean (jul17);
mjul18 = nanmean (jul18);
mjul19 = nanmean (jul19);
mjul20 = nanmean (jul20);
mjul21 = nanmean (jul21);
mjul22 = nanmean (jul22);
mjul23 = nanmean (jul23);

maxjul0 = nanmax (jul0);
maxjul1 = nanmax (jul1);
maxjul2 = nanmax (jul2);
maxjul3 = nanmax (jul3);
maxjul4 = nanmax (jul4);
maxjul5 = nanmax (jul5);
maxjul6 = nanmax (jul6);
maxjul7 = nanmax (jul7);
maxjul8 = nanmax (jul8);
maxjul9 = nanmax (jul9);
maxjul10 = nanmax (jul10);
maxjul11 = nanmax (jul11);
maxjul12 = nanmax (jul12);
maxjul13 = nanmax (jul13);
maxjul14 = nanmax (jul14);
maxjul15 = nanmax (jul15);
maxjul16 = nanmax (jul16);
maxjul17 = nanmax (jul17);
maxjul18 = nanmax (jul18);
maxjul19 = nanmax (jul19);
maxjul20 = nanmax (jul20);
maxjul21 = nanmax (jul21);
maxjul22 = nanmax (jul22);
maxjul23 = nanmax (jul23);

minjul0 = nanmin (jul0);
minjul1 = nanmin (jul1);
minjul2 = nanmin (jul2);
minjul3 = nanmin (jul3);
minjul4 = nanmin (jul4);
minjul5 = nanmin (jul5);
minjul6 = nanmin (jul6);
minjul7 = nanmin (jul7);
minjul8 = nanmin (jul8);
minjul9 = nanmin (jul9);
minjul10 = nanmin (jul10);
minjul11 = nanmin (jul11);
minjul12 = nanmin (jul12);
minjul13 = nanmin (jul13);
minjul14 = nanmin (jul14);
minjul15 = nanmin (jul15);
minjul16 = nanmin (jul16);
minjul17 = nanmin (jul17);
minjul18 = nanmin (jul18);
minjul19 = nanmin (jul19);
minjul20 = nanmin (jul20);
minjul21 = nanmin (jul21);
minjul22 = nanmin (jul22);
minjul23 = nanmin (jul23);


tjul0 = [mjul0,maxjul0(1,4),minjul0(1,4)];
tjul1 = [mjul1,maxjul1(1,4),minjul1(1,4)];
tjul2 = [mjul2,maxjul2(1,4),minjul2(1,4)];
tjul3 = [mjul3,maxjul3(1,4),minjul3(1,4)];
tjul4 = [mjul4,maxjul4(1,4),minjul4(1,4)];
tjul5 = [mjul5,maxjul5(1,4),minjul5(1,4)];
tjul6 = [mjul6,maxjul6(1,4),minjul6(1,4)];
tjul7 = [mjul7,maxjul7(1,4),minjul7(1,4)];
tjul8 = [mjul8,maxjul8(1,4),minjul8(1,4)];
tjul9 = [mjul9,maxjul9(1,4),minjul9(1,4)];
tjul10 = [mjul10,maxjul10(1,4),minjul10(1,4)];
tjul11 = [mjul11,maxjul11(1,4),minjul11(1,4)];
tjul12 = [mjul12,maxjul12(1,4),minjul12(1,4)];
tjul13 = [mjul13,maxjul13(1,4),minjul13(1,4)];
tjul14 = [mjul14,maxjul14(1,4),minjul14(1,4)];
tjul15 = [mjul15,maxjul15(1,4),minjul15(1,4)];
tjul16 = [mjul16,maxjul16(1,4),minjul16(1,4)];
tjul17 = [mjul17,maxjul17(1,4),minjul17(1,4)];
tjul18 = [mjul18,maxjul18(1,4),minjul18(1,4)];
tjul19 = [mjul19,maxjul19(1,4),minjul19(1,4)];
tjul20 = [mjul20,maxjul20(1,4),minjul20(1,4)];
tjul21 = [mjul21,maxjul21(1,4),minjul21(1,4)];
tjul22 = [mjul22,maxjul22(1,4),minjul22(1,4)];
tjul23 = [mjul23,maxjul23(1,4),minjul23(1,4)];

tjul = [tjul0;tjul1;tjul2;tjul3;tjul4;tjul5;tjul6;tjul7;tjul8;tjul9;tjul10;tjul11;tjul12;tjul13;tjul14;tjul15;tjul16;tjul17;tjul18;tjul19;tjul20;tjul21;tjul22;tjul23];


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



load 'uh82augT'
load 'uh83augT'
load 'uh84augT'
load 'uh85augT'
load 'uh86augT'
load 'uh87augT'
load 'uh88augT'
load 'uh89augT'
load 'uh90augT'
load 'uh91augT'
load 'uh92augT'
load 'uh93augT'
load 'uh94augT'
load 'uh95augT'
load 'uh96augT'
load 'uh97augT'
load 'uh98augT'
load 'uh99augT'
load 'uh00augT'
load 'uh01augT'
load 'uh02augT'
load 'uh03augT'
load 'uh07augT'
load 'uh08augT'
load 'uh09augT'
load 'uh10augT'


maxaug = [uh82augT;uh83augT;uh84augT;uh85augT;uh86augT;uh87augT;uh88augT;uh89augT;uh90augT;uh91augT;uh92augT;uh93augT;uh94augT;uh95augT;uh96augT;uh97augT;uh98augT;uh99augT;uh00augT;uh01augT;uh02augT;uh03augT;uh07augT;uh08augT;uh09augT;uh10augT];

[aug0 aug1 aug2 aug3 aug4 aug5 aug6 aug7 aug8 aug9 aug10 aug11 aug12 aug13 aug14 aug15 aug16 aug17 aug18 aug19 aug20 aug21 aug22 aug23] = extracthours (maxaug);

maug0 = nanmean (aug0);
maug1 = nanmean (aug1);
maug2 = nanmean (aug2);
maug3 = nanmean (aug3);
maug4 = nanmean (aug4);
maug5 = nanmean (aug5);
maug6 = nanmean (aug6);
maug7 = nanmean (aug7);
maug8 = nanmean (aug8);
maug9 = nanmean (aug9);
maug10 = nanmean (aug10);
maug11 = nanmean (aug11);
maug12 = nanmean (aug12);
maug13 = nanmean (aug13);
maug14 = nanmean (aug14);
maug15 = nanmean (aug15);
maug16 = nanmean (aug16);
maug17 = nanmean (aug17);
maug18 = nanmean (aug18);
maug19 = nanmean (aug19);
maug20 = nanmean (aug20);
maug21 = nanmean (aug21);
maug22 = nanmean (aug22);
maug23 = nanmean (aug23);

maxaug0 = nanmax (aug0);
maxaug1 = nanmax (aug1);
maxaug2 = nanmax (aug2);
maxaug3 = nanmax (aug3);
maxaug4 = nanmax (aug4);
maxaug5 = nanmax (aug5);
maxaug6 = nanmax (aug6);
maxaug7 = nanmax (aug7);
maxaug8 = nanmax (aug8);
maxaug9 = nanmax (aug9);
maxaug10 = nanmax (aug10);
maxaug11 = nanmax (aug11);
maxaug12 = nanmax (aug12);
maxaug13 = nanmax (aug13);
maxaug14 = nanmax (aug14);
maxaug15 = nanmax (aug15);
maxaug16 = nanmax (aug16);
maxaug17 = nanmax (aug17);
maxaug18 = nanmax (aug18);
maxaug19 = nanmax (aug19);
maxaug20 = nanmax (aug20);
maxaug21 = nanmax (aug21);
maxaug22 = nanmax (aug22);
maxaug23 = nanmax (aug23);

minaug0 = nanmin (aug0);
minaug1 = nanmin (aug1);
minaug2 = nanmin (aug2);
minaug3 = nanmin (aug3);
minaug4 = nanmin (aug4);
minaug5 = nanmin (aug5);
minaug6 = nanmin (aug6);
minaug7 = nanmin (aug7);
minaug8 = nanmin (aug8);
minaug9 = nanmin (aug9);
minaug10 = nanmin (aug10);
minaug11 = nanmin (aug11);
minaug12 = nanmin (aug12);
minaug13 = nanmin (aug13);
minaug14 = nanmin (aug14);
minaug15 = nanmin (aug15);
minaug16 = nanmin (aug16);
minaug17 = nanmin (aug17);
minaug18 = nanmin (aug18);
minaug19 = nanmin (aug19);
minaug20 = nanmin (aug20);
minaug21 = nanmin (aug21);
minaug22 = nanmin (aug22);
minaug23 = nanmin (aug23);


taug0 = [maug0,maxaug0(1,4),minaug0(1,4)];
taug1 = [maug1,maxaug1(1,4),minaug1(1,4)];
taug2 = [maug2,maxaug2(1,4),minaug2(1,4)];
taug3 = [maug3,maxaug3(1,4),minaug3(1,4)];
taug4 = [maug4,maxaug4(1,4),minaug4(1,4)];
taug5 = [maug5,maxaug5(1,4),minaug5(1,4)];
taug6 = [maug6,maxaug6(1,4),minaug6(1,4)];
taug7 = [maug7,maxaug7(1,4),minaug7(1,4)];
taug8 = [maug8,maxaug8(1,4),minaug8(1,4)];
taug9 = [maug9,maxaug9(1,4),minaug9(1,4)];
taug10 = [maug10,maxaug10(1,4),minaug10(1,4)];
taug11 = [maug11,maxaug11(1,4),minaug11(1,4)];
taug12 = [maug12,maxaug12(1,4),minaug12(1,4)];
taug13 = [maug13,maxaug13(1,4),minaug13(1,4)];
taug14 = [maug14,maxaug14(1,4),minaug14(1,4)];
taug15 = [maug15,maxaug15(1,4),minaug15(1,4)];
taug16 = [maug16,maxaug16(1,4),minaug16(1,4)];
taug17 = [maug17,maxaug17(1,4),minaug17(1,4)];
taug18 = [maug18,maxaug18(1,4),minaug18(1,4)];
taug19 = [maug19,maxaug19(1,4),minaug19(1,4)];
taug20 = [maug20,maxaug20(1,4),minaug20(1,4)];
taug21 = [maug21,maxaug21(1,4),minaug21(1,4)];
taug22 = [maug22,maxaug22(1,4),minaug22(1,4)];
taug23 = [maug23,maxaug23(1,4),minaug23(1,4)];

taug = [taug0;taug1;taug2;taug3;taug4;taug5;taug6;taug7;taug8;taug9;taug10;taug11;taug12;taug13;taug14;taug15;taug16;taug17;taug18;taug19;taug20;taug21;taug22;taug23];


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


load 'uh82sepT'
load 'uh83sepT'
load 'uh84sepT'
load 'uh85sepT'
load 'uh86sepT'
load 'uh87sepT'
load 'uh88sepT'
load 'uh89sepT'
load 'uh90sepT'
load 'uh91sepT'
load 'uh92sepT'
load 'uh93sepT'
load 'uh94sepT'
load 'uh95sepT'
load 'uh96sepT'
load 'uh97sepT'
load 'uh98sepT'
load 'uh99sepT'
load 'uh00sepT'
load 'uh01sepT'
load 'uh02sepT'
load 'uh03sepT'
load 'uh07sepT'
load 'uh08sepT'
load 'uh09sepT'
load 'uh10sepT'


maxsep = [uh82sepT;uh83sepT;uh84sepT;uh85sepT;uh86sepT;uh87sepT;uh88sepT;uh89sepT;uh90sepT;uh91sepT;uh92sepT;uh93sepT;uh94sepT;uh95sepT;uh96sepT;uh97sepT;uh98sepT;uh99sepT;uh00sepT;uh01sepT;uh02sepT;uh03sepT;uh07sepT;uh08sepT;uh09sepT;uh10sepT];

[sep0 sep1 sep2 sep3 sep4 sep5 sep6 sep7 sep8 sep9 sep10 sep11 sep12 sep13 sep14 sep15 sep16 sep17 sep18 sep19 sep20 sep21 sep22 sep23] = extracthours (maxsep);

msep0 = nanmean (sep0);
msep1 = nanmean (sep1);
msep2 = nanmean (sep2);
msep3 = nanmean (sep3);
msep4 = nanmean (sep4);
msep5 = nanmean (sep5);
msep6 = nanmean (sep6);
msep7 = nanmean (sep7);
msep8 = nanmean (sep8);
msep9 = nanmean (sep9);
msep10 = nanmean (sep10);
msep11 = nanmean (sep11);
msep12 = nanmean (sep12);
msep13 = nanmean (sep13);
msep14 = nanmean (sep14);
msep15 = nanmean (sep15);
msep16 = nanmean (sep16);
msep17 = nanmean (sep17);
msep18 = nanmean (sep18);
msep19 = nanmean (sep19);
msep20 = nanmean (sep20);
msep21 = nanmean (sep21);
msep22 = nanmean (sep22);
msep23 = nanmean (sep23);

maxsep0 = nanmax (sep0);
maxsep1 = nanmax (sep1);
maxsep2 = nanmax (sep2);
maxsep3 = nanmax (sep3);
maxsep4 = nanmax (sep4);
maxsep5 = nanmax (sep5);
maxsep6 = nanmax (sep6);
maxsep7 = nanmax (sep7);
maxsep8 = nanmax (sep8);
maxsep9 = nanmax (sep9);
maxsep10 = nanmax (sep10);
maxsep11 = nanmax (sep11);
maxsep12 = nanmax (sep12);
maxsep13 = nanmax (sep13);
maxsep14 = nanmax (sep14);
maxsep15 = nanmax (sep15);
maxsep16 = nanmax (sep16);
maxsep17 = nanmax (sep17);
maxsep18 = nanmax (sep18);
maxsep19 = nanmax (sep19);
maxsep20 = nanmax (sep20);
maxsep21 = nanmax (sep21);
maxsep22 = nanmax (sep22);
maxsep23 = nanmax (sep23);

minsep0 = nanmin (sep0);
minsep1 = nanmin (sep1);
minsep2 = nanmin (sep2);
minsep3 = nanmin (sep3);
minsep4 = nanmin (sep4);
minsep5 = nanmin (sep5);
minsep6 = nanmin (sep6);
minsep7 = nanmin (sep7);
minsep8 = nanmin (sep8);
minsep9 = nanmin (sep9);
minsep10 = nanmin (sep10);
minsep11 = nanmin (sep11);
minsep12 = nanmin (sep12);
minsep13 = nanmin (sep13);
minsep14 = nanmin (sep14);
minsep15 = nanmin (sep15);
minsep16 = nanmin (sep16);
minsep17 = nanmin (sep17);
minsep18 = nanmin (sep18);
minsep19 = nanmin (sep19);
minsep20 = nanmin (sep20);
minsep21 = nanmin (sep21);
minsep22 = nanmin (sep22);
minsep23 = nanmin (sep23);


tsep0 = [msep0,maxsep0(1,4),minsep0(1,4)];
tsep1 = [msep1,maxsep1(1,4),minsep1(1,4)];
tsep2 = [msep2,maxsep2(1,4),minsep2(1,4)];
tsep3 = [msep3,maxsep3(1,4),minsep3(1,4)];
tsep4 = [msep4,maxsep4(1,4),minsep4(1,4)];
tsep5 = [msep5,maxsep5(1,4),minsep5(1,4)];
tsep6 = [msep6,maxsep6(1,4),minsep6(1,4)];
tsep7 = [msep7,maxsep7(1,4),minsep7(1,4)];
tsep8 = [msep8,maxsep8(1,4),minsep8(1,4)];
tsep9 = [msep9,maxsep9(1,4),minsep9(1,4)];
tsep10 = [msep10,maxsep10(1,4),minsep10(1,4)];
tsep11 = [msep11,maxsep11(1,4),minsep11(1,4)];
tsep12 = [msep12,maxsep12(1,4),minsep12(1,4)];
tsep13 = [msep13,maxsep13(1,4),minsep13(1,4)];
tsep14 = [msep14,maxsep14(1,4),minsep14(1,4)];
tsep15 = [msep15,maxsep15(1,4),minsep15(1,4)];
tsep16 = [msep16,maxsep16(1,4),minsep16(1,4)];
tsep17 = [msep17,maxsep17(1,4),minsep17(1,4)];
tsep18 = [msep18,maxsep18(1,4),minsep18(1,4)];
tsep19 = [msep19,maxsep19(1,4),minsep19(1,4)];
tsep20 = [msep20,maxsep20(1,4),minsep20(1,4)];
tsep21 = [msep21,maxsep21(1,4),minsep21(1,4)];
tsep22 = [msep22,maxsep22(1,4),minsep22(1,4)];
tsep23 = [msep23,maxsep23(1,4),minsep23(1,4)];

tsep = [tsep0;tsep1;tsep2;tsep3;tsep4;tsep5;tsep6;tsep7;tsep8;tsep9;tsep10;tsep11;tsep12;tsep13;tsep14;tsep15;tsep16;tsep17;tsep18;tsep19;tsep20;tsep21;tsep22;tsep23];


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

load 'uh82octT'
load 'uh83octT'
load 'uh84octT'
load 'uh85octT'
load 'uh86octT'
load 'uh87octT'
load 'uh88octT'
load 'uh89octT'
load 'uh90octT'
load 'uh91octT'
load 'uh92octT'
load 'uh93octT'
load 'uh94octT'
load 'uh95octT'
load 'uh96octT'
load 'uh97octT'
load 'uh98octT'
load 'uh99octT'
load 'uh00octT'
load 'uh01octT'
load 'uh02octT'
load 'uh03octT'
load 'uh07octT'
load 'uh08octT'
load 'uh09octT'
load 'uh10octT'


maxoct = [uh82octT;uh83octT;uh84octT;uh85octT;uh86octT;uh87octT;uh88octT;uh89octT;uh90octT;uh91octT;uh92octT;uh93octT;uh94octT;uh95octT;uh96octT;uh97octT;uh98octT;uh99octT;uh00octT;uh01octT;uh02octT;uh03octT;uh07octT;uh08octT;uh09octT;uh10octT];

[oct0 oct1 oct2 oct3 oct4 oct5 oct6 oct7 oct8 oct9 oct10 oct11 oct12 oct13 oct14 oct15 oct16 oct17 oct18 oct19 oct20 oct21 oct22 oct23] = extracthours (maxoct);

moct0 = nanmean (oct0);
moct1 = nanmean (oct1);
moct2 = nanmean (oct2);
moct3 = nanmean (oct3);
moct4 = nanmean (oct4);
moct5 = nanmean (oct5);
moct6 = nanmean (oct6);
moct7 = nanmean (oct7);
moct8 = nanmean (oct8);
moct9 = nanmean (oct9);
moct10 = nanmean (oct10);
moct11 = nanmean (oct11);
moct12 = nanmean (oct12);
moct13 = nanmean (oct13);
moct14 = nanmean (oct14);
moct15 = nanmean (oct15);
moct16 = nanmean (oct16);
moct17 = nanmean (oct17);
moct18 = nanmean (oct18);
moct19 = nanmean (oct19);
moct20 = nanmean (oct20);
moct21 = nanmean (oct21);
moct22 = nanmean (oct22);
moct23 = nanmean (oct23);

maxoct0 = nanmax (oct0);
maxoct1 = nanmax (oct1);
maxoct2 = nanmax (oct2);
maxoct3 = nanmax (oct3);
maxoct4 = nanmax (oct4);
maxoct5 = nanmax (oct5);
maxoct6 = nanmax (oct6);
maxoct7 = nanmax (oct7);
maxoct8 = nanmax (oct8);
maxoct9 = nanmax (oct9);
maxoct10 = nanmax (oct10);
maxoct11 = nanmax (oct11);
maxoct12 = nanmax (oct12);
maxoct13 = nanmax (oct13);
maxoct14 = nanmax (oct14);
maxoct15 = nanmax (oct15);
maxoct16 = nanmax (oct16);
maxoct17 = nanmax (oct17);
maxoct18 = nanmax (oct18);
maxoct19 = nanmax (oct19);
maxoct20 = nanmax (oct20);
maxoct21 = nanmax (oct21);
maxoct22 = nanmax (oct22);
maxoct23 = nanmax (oct23);

minoct0 = nanmin (oct0);
minoct1 = nanmin (oct1);
minoct2 = nanmin (oct2);
minoct3 = nanmin (oct3);
minoct4 = nanmin (oct4);
minoct5 = nanmin (oct5);
minoct6 = nanmin (oct6);
minoct7 = nanmin (oct7);
minoct8 = nanmin (oct8);
minoct9 = nanmin (oct9);
minoct10 = nanmin (oct10);
minoct11 = nanmin (oct11);
minoct12 = nanmin (oct12);
minoct13 = nanmin (oct13);
minoct14 = nanmin (oct14);
minoct15 = nanmin (oct15);
minoct16 = nanmin (oct16);
minoct17 = nanmin (oct17);
minoct18 = nanmin (oct18);
minoct19 = nanmin (oct19);
minoct20 = nanmin (oct20);
minoct21 = nanmin (oct21);
minoct22 = nanmin (oct22);
minoct23 = nanmin (oct23);


toct0 = [moct0,maxoct0(1,4),minoct0(1,4)];
toct1 = [moct1,maxoct1(1,4),minoct1(1,4)];
toct2 = [moct2,maxoct2(1,4),minoct2(1,4)];
toct3 = [moct3,maxoct3(1,4),minoct3(1,4)];
toct4 = [moct4,maxoct4(1,4),minoct4(1,4)];
toct5 = [moct5,maxoct5(1,4),minoct5(1,4)];
toct6 = [moct6,maxoct6(1,4),minoct6(1,4)];
toct7 = [moct7,maxoct7(1,4),minoct7(1,4)];
toct8 = [moct8,maxoct8(1,4),minoct8(1,4)];
toct9 = [moct9,maxoct9(1,4),minoct9(1,4)];
toct10 = [moct10,maxoct10(1,4),minoct10(1,4)];
toct11 = [moct11,maxoct11(1,4),minoct11(1,4)];
toct12 = [moct12,maxoct12(1,4),minoct12(1,4)];
toct13 = [moct13,maxoct13(1,4),minoct13(1,4)];
toct14 = [moct14,maxoct14(1,4),minoct14(1,4)];
toct15 = [moct15,maxoct15(1,4),minoct15(1,4)];
toct16 = [moct16,maxoct16(1,4),minoct16(1,4)];
toct17 = [moct17,maxoct17(1,4),minoct17(1,4)];
toct18 = [moct18,maxoct18(1,4),minoct18(1,4)];
toct19 = [moct19,maxoct19(1,4),minoct19(1,4)];
toct20 = [moct20,maxoct20(1,4),minoct20(1,4)];
toct21 = [moct21,maxoct21(1,4),minoct21(1,4)];
toct22 = [moct22,maxoct22(1,4),minoct22(1,4)];
toct23 = [moct23,maxoct23(1,4),minoct23(1,4)];

toct = [toct0;toct1;toct2;toct3;toct4;toct5;toct6;toct7;toct8;toct9;toct10;toct11;toct12;toct13;toct14;toct15;toct16;toct17;toct18;toct19;toct20;toct21;toct22;toct23];


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

load 'uh82novT'
load 'uh83novT'
load 'uh84novT'
load 'uh85novT'
load 'uh86novT'
load 'uh87novT'
load 'uh88novT'
load 'uh89novT'
load 'uh90novT'
load 'uh91novT'
load 'uh92novT'
load 'uh93novT'
load 'uh94novT'
load 'uh95novT'
load 'uh96novT'
load 'uh97novT'
load 'uh98novT'
load 'uh99novT'
load 'uh00novT'
load 'uh01novT'
load 'uh02novT'
load 'uh03novT'
load 'uh07novT'
load 'uh08novT'
load 'uh09novT'
load 'uh10novT'


maxnov = [uh82novT;uh83novT;uh84novT;uh85novT;uh86novT;uh87novT;uh88novT;uh89novT;uh90novT;uh91novT;uh92novT;uh93novT;uh94novT;uh95novT;uh96novT;uh97novT;uh98novT;uh99novT;uh00novT;uh01novT;uh02novT;uh03novT;uh07novT;uh08novT;uh09novT;uh10novT];

[nov0 nov1 nov2 nov3 nov4 nov5 nov6 nov7 nov8 nov9 nov10 nov11 nov12 nov13 nov14 nov15 nov16 nov17 nov18 nov19 nov20 nov21 nov22 nov23] = extracthours (maxnov);

mnov0 = nanmean (nov0);
mnov1 = nanmean (nov1);
mnov2 = nanmean (nov2);
mnov3 = nanmean (nov3);
mnov4 = nanmean (nov4);
mnov5 = nanmean (nov5);
mnov6 = nanmean (nov6);
mnov7 = nanmean (nov7);
mnov8 = nanmean (nov8);
mnov9 = nanmean (nov9);
mnov10 = nanmean (nov10);
mnov11 = nanmean (nov11);
mnov12 = nanmean (nov12);
mnov13 = nanmean (nov13);
mnov14 = nanmean (nov14);
mnov15 = nanmean (nov15);
mnov16 = nanmean (nov16);
mnov17 = nanmean (nov17);
mnov18 = nanmean (nov18);
mnov19 = nanmean (nov19);
mnov20 = nanmean (nov20);
mnov21 = nanmean (nov21);
mnov22 = nanmean (nov22);
mnov23 = nanmean (nov23);

maxnov0 = nanmax (nov0);
maxnov1 = nanmax (nov1);
maxnov2 = nanmax (nov2);
maxnov3 = nanmax (nov3);
maxnov4 = nanmax (nov4);
maxnov5 = nanmax (nov5);
maxnov6 = nanmax (nov6);
maxnov7 = nanmax (nov7);
maxnov8 = nanmax (nov8);
maxnov9 = nanmax (nov9);
maxnov10 = nanmax (nov10);
maxnov11 = nanmax (nov11);
maxnov12 = nanmax (nov12);
maxnov13 = nanmax (nov13);
maxnov14 = nanmax (nov14);
maxnov15 = nanmax (nov15);
maxnov16 = nanmax (nov16);
maxnov17 = nanmax (nov17);
maxnov18 = nanmax (nov18);
maxnov19 = nanmax (nov19);
maxnov20 = nanmax (nov20);
maxnov21 = nanmax (nov21);
maxnov22 = nanmax (nov22);
maxnov23 = nanmax (nov23);

minnov0 = nanmin (nov0);
minnov1 = nanmin (nov1);
minnov2 = nanmin (nov2);
minnov3 = nanmin (nov3);
minnov4 = nanmin (nov4);
minnov5 = nanmin (nov5);
minnov6 = nanmin (nov6);
minnov7 = nanmin (nov7);
minnov8 = nanmin (nov8);
minnov9 = nanmin (nov9);
minnov10 = nanmin (nov10);
minnov11 = nanmin (nov11);
minnov12 = nanmin (nov12);
minnov13 = nanmin (nov13);
minnov14 = nanmin (nov14);
minnov15 = nanmin (nov15);
minnov16 = nanmin (nov16);
minnov17 = nanmin (nov17);
minnov18 = nanmin (nov18);
minnov19 = nanmin (nov19);
minnov20 = nanmin (nov20);
minnov21 = nanmin (nov21);
minnov22 = nanmin (nov22);
minnov23 = nanmin (nov23);


tnov0 = [mnov0,maxnov0(1,4),minnov0(1,4)];
tnov1 = [mnov1,maxnov1(1,4),minnov1(1,4)];
tnov2 = [mnov2,maxnov2(1,4),minnov2(1,4)];
tnov3 = [mnov3,maxnov3(1,4),minnov3(1,4)];
tnov4 = [mnov4,maxnov4(1,4),minnov4(1,4)];
tnov5 = [mnov5,maxnov5(1,4),minnov5(1,4)];
tnov6 = [mnov6,maxnov6(1,4),minnov6(1,4)];
tnov7 = [mnov7,maxnov7(1,4),minnov7(1,4)];
tnov8 = [mnov8,maxnov8(1,4),minnov8(1,4)];
tnov9 = [mnov9,maxnov9(1,4),minnov9(1,4)];
tnov10 = [mnov10,maxnov10(1,4),minnov10(1,4)];
tnov11 = [mnov11,maxnov11(1,4),minnov11(1,4)];
tnov12 = [mnov12,maxnov12(1,4),minnov12(1,4)];
tnov13 = [mnov13,maxnov13(1,4),minnov13(1,4)];
tnov14 = [mnov14,maxnov14(1,4),minnov14(1,4)];
tnov15 = [mnov15,maxnov15(1,4),minnov15(1,4)];
tnov16 = [mnov16,maxnov16(1,4),minnov16(1,4)];
tnov17 = [mnov17,maxnov17(1,4),minnov17(1,4)];
tnov18 = [mnov18,maxnov18(1,4),minnov18(1,4)];
tnov19 = [mnov19,maxnov19(1,4),minnov19(1,4)];
tnov20 = [mnov20,maxnov20(1,4),minnov20(1,4)];
tnov21 = [mnov21,maxnov21(1,4),minnov21(1,4)];
tnov22 = [mnov22,maxnov22(1,4),minnov22(1,4)];
tnov23 = [mnov23,maxnov23(1,4),minnov23(1,4)];

tnov = [tnov0;tnov1;tnov2;tnov3;tnov4;tnov5;tnov6;tnov7;tnov8;tnov9;tnov10;tnov11;tnov12;tnov13;tnov14;tnov15;tnov16;tnov17;tnov18;tnov19;tnov20;tnov21;tnov22;tnov23];


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



load 'uh82decT'
load 'uh83decT'
load 'uh84decT'
load 'uh85decT'
load 'uh86decT'
load 'uh87decT'
load 'uh88decT'
load 'uh89decT'
load 'uh90decT'
load 'uh91decT'
load 'uh92decT'
load 'uh93decT'
load 'uh94decT'
load 'uh95decT'
load 'uh96decT'
load 'uh97decT'
load 'uh98decT'
load 'uh99decT'
load 'uh00decT'
load 'uh01decT'
load 'uh02decT'
load 'uh03decT'
load 'uh07decT'
load 'uh08decT'
load 'uh09decT'
load 'uh10decT'


maxdec = [uh82decT;uh83decT;uh84decT;uh85decT;uh86decT;uh87decT;uh88decT;uh89decT;uh90decT;uh91decT;uh92decT;uh93decT;uh94decT;uh95decT;uh96decT;uh97decT;uh98decT;uh99decT;uh00decT;uh01decT;uh02decT;uh03decT;uh07decT;uh08decT;uh09decT;uh10decT];

[dec0 dec1 dec2 dec3 dec4 dec5 dec6 dec7 dec8 dec9 dec10 dec11 dec12 dec13 dec14 dec15 dec16 dec17 dec18 dec19 dec20 dec21 dec22 dec23] = extracthours (maxdec);

mdec0 = nanmean (dec0);
mdec1 = nanmean (dec1);
mdec2 = nanmean (dec2);
mdec3 = nanmean (dec3);
mdec4 = nanmean (dec4);
mdec5 = nanmean (dec5);
mdec6 = nanmean (dec6);
mdec7 = nanmean (dec7);
mdec8 = nanmean (dec8);
mdec9 = nanmean (dec9);
mdec10 = nanmean (dec10);
mdec11 = nanmean (dec11);
mdec12 = nanmean (dec12);
mdec13 = nanmean (dec13);
mdec14 = nanmean (dec14);
mdec15 = nanmean (dec15);
mdec16 = nanmean (dec16);
mdec17 = nanmean (dec17);
mdec18 = nanmean (dec18);
mdec19 = nanmean (dec19);
mdec20 = nanmean (dec20);
mdec21 = nanmean (dec21);
mdec22 = nanmean (dec22);
mdec23 = nanmean (dec23);

maxdec0 = nanmax (dec0);
maxdec1 = nanmax (dec1);
maxdec2 = nanmax (dec2);
maxdec3 = nanmax (dec3);
maxdec4 = nanmax (dec4);
maxdec5 = nanmax (dec5);
maxdec6 = nanmax (dec6);
maxdec7 = nanmax (dec7);
maxdec8 = nanmax (dec8);
maxdec9 = nanmax (dec9);
maxdec10 = nanmax (dec10);
maxdec11 = nanmax (dec11);
maxdec12 = nanmax (dec12);
maxdec13 = nanmax (dec13);
maxdec14 = nanmax (dec14);
maxdec15 = nanmax (dec15);
maxdec16 = nanmax (dec16);
maxdec17 = nanmax (dec17);
maxdec18 = nanmax (dec18);
maxdec19 = nanmax (dec19);
maxdec20 = nanmax (dec20);
maxdec21 = nanmax (dec21);
maxdec22 = nanmax (dec22);
maxdec23 = nanmax (dec23);

mindec0 = nanmin (dec0);
mindec1 = nanmin (dec1);
mindec2 = nanmin (dec2);
mindec3 = nanmin (dec3);
mindec4 = nanmin (dec4);
mindec5 = nanmin (dec5);
mindec6 = nanmin (dec6);
mindec7 = nanmin (dec7);
mindec8 = nanmin (dec8);
mindec9 = nanmin (dec9);
mindec10 = nanmin (dec10);
mindec11 = nanmin (dec11);
mindec12 = nanmin (dec12);
mindec13 = nanmin (dec13);
mindec14 = nanmin (dec14);
mindec15 = nanmin (dec15);
mindec16 = nanmin (dec16);
mindec17 = nanmin (dec17);
mindec18 = nanmin (dec18);
mindec19 = nanmin (dec19);
mindec20 = nanmin (dec20);
mindec21 = nanmin (dec21);
mindec22 = nanmin (dec22);
mindec23 = nanmin (dec23);


tdec0 = [mdec0,maxdec0(1,4),mindec0(1,4)];
tdec1 = [mdec1,maxdec1(1,4),mindec1(1,4)];
tdec2 = [mdec2,maxdec2(1,4),mindec2(1,4)];
tdec3 = [mdec3,maxdec3(1,4),mindec3(1,4)];
tdec4 = [mdec4,maxdec4(1,4),mindec4(1,4)];
tdec5 = [mdec5,maxdec5(1,4),mindec5(1,4)];
tdec6 = [mdec6,maxdec6(1,4),mindec6(1,4)];
tdec7 = [mdec7,maxdec7(1,4),mindec7(1,4)];
tdec8 = [mdec8,maxdec8(1,4),mindec8(1,4)];
tdec9 = [mdec9,maxdec9(1,4),mindec9(1,4)];
tdec10 = [mdec10,maxdec10(1,4),mindec10(1,4)];
tdec11 = [mdec11,maxdec11(1,4),mindec11(1,4)];
tdec12 = [mdec12,maxdec12(1,4),mindec12(1,4)];
tdec13 = [mdec13,maxdec13(1,4),mindec13(1,4)];
tdec14 = [mdec14,maxdec14(1,4),mindec14(1,4)];
tdec15 = [mdec15,maxdec15(1,4),mindec15(1,4)];
tdec16 = [mdec16,maxdec16(1,4),mindec16(1,4)];
tdec17 = [mdec17,maxdec17(1,4),mindec17(1,4)];
tdec18 = [mdec18,maxdec18(1,4),mindec18(1,4)];
tdec19 = [mdec19,maxdec19(1,4),mindec19(1,4)];
tdec20 = [mdec20,maxdec20(1,4),mindec20(1,4)];
tdec21 = [mdec21,maxdec21(1,4),mindec21(1,4)];
tdec22 = [mdec22,maxdec22(1,4),mindec22(1,4)];
tdec23 = [mdec23,maxdec23(1,4),mindec23(1,4)];

tdec = [tdec0;tdec1;tdec2;tdec3;tdec4;tdec5;tdec6;tdec7;tdec8;tdec9;tdec10;tdec11;tdec12;tdec13;tdec14;tdec15;tdec16;tdec17;tdec18;tdec19;tdec20;tdec21;tdec22;tdec23];


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