[year,jan,feb,mar,apr,may,jun,jul,aug,sep,oct,nov,dec,annmean]=...
   textread('Tempelhofmonthtemp.txt','%d %f %f  %f %f %f %f %f %f %f %f %f %f %f','headerlines',7);

badval=year==-999.9;
year(badval)=NaN;
badval=jan==-999.9;
jan(badval)=NaN;
badval=feb==-999.9;
feb(badval)=NaN;
badval=mar==-999.9;
mar(badval)=NaN;
badval=apr==-999.9;
apr(badval)=NaN;
badval=may==-999.9;
may(badval)=NaN;
badval=jun==-999.9;
jun(badval)=NaN;
badval=jul==-999.9;
jul(badval)=NaN;
badval=aug==-999.9;
aug(badval)=NaN;
badval=sep==-999.9;
sep(badval)=NaN;
badval=oct==-999.9;
oct(badval)=NaN;
badval=nov==-999.9;
nov(badval)=NaN;
badval=dec==-999.9;
dec(badval)=NaN;
badval=annmean==-999.9;
annmean(badval)=NaN;

%Festlegung: kleine Eiszeit von 1675-1715, Industrialisierungsphase in
%Deutschland 1830-1873, 1990er Jahre 1990-1993

mittelwertallejahre=nanmean(annmean);
data(data==-999.9)=NaN;
klez=annmean(1:15);
mnklez=nanmean(klez); %durchschnittliche Jahresmitteltemperatur während 
%der kleinen Eiszeit
ind=annmean(130:173);
mnind=nanmean(ind); %durchschnittliche Jahresmitteltemperatur während der 
%Industrialisierungsphase in Deutschland
neun=annmean(290:293);
mnneun=nanmean(neun); %durchschnittliche Jahresmitteltemperatur während der
%1990er Jahre

[r,p]=corr(year,annmean);
[r1,p1]=corr(year,jan);
[r2,p2]=corr(year,feb);
[r3,p3]=corr(year,mar);
[r4,p4]=corr(year,apr);
[r5,p5]=corr(year,may);
[r6,p6]=corr(year,jun);
[r7,p7]=corr(year,jul);
[r8,p8]=corr(year,aug);





