function gradtagklassisch
%Berechnung des klassischen Gradtagfaktors mit Lufftemperatur und Schmelze
%zuerst einlesen der Pegeldaten und Wetterstationsdaten von Matritsch
%(benoetigt nur das metinput mit korrigierten Werten)
[year,jday,month,day,hour,minute,T,RH,Swin]=madritsch_met
[yyyy04,mm04,dd04,pegel07_04,pegel08_04,pegel09_04,pegel10_04,pegel11_04,pegel12_04]=pegeldaten2004
[yyyy05,mm05,dd05,pegel07_05,pegel0705_05,pegel08_05,pegel0805_05,pegel09_05,pegel10_05,pegel11_05,pegel12_05]=pegeldaten2005
[yyyy06,mm06,dd06,pegel07_06,pegel08_06,pegel09_06,pegel10_06,pegel11_06,pegel12_06]=pegeldaten2006
[yyyy07,mm07,dd07,pegel07_07,pegel08_07,pegel09_07,pegel10_07,pegel11_07,pegel12_07]=pegeldaten2007
[yyyy08,mm08,dd08,pegel07_08,pegel08_08,pegel09_08,pegel10_08,pegel11_08,pegel12_08]=pegeldaten2008
[yyyy09,mm09,dd09,pegel07_09,pegel08_09,pegel09_09,pegel10_09,pegel11_09,pegel12_09]=pegeldaten2009
[yyyy10,mm10,dd10,pegel07_10,pegel08_10]=pegeldaten2010
[yyyy11,mm11,dd11,pegel07_11,pegel08_11,pegel09_11,pegel10_11,pegel11_11]=pegeldaten2011
[yyyy12,mm12,dd12,pegel07_12,pegel08_12,pegel09_12,pegel10_12,pegel11_12,pegel12_12]=pegeldaten2012
[yyyy13,mm13,dd13,pegel07_13,pegel08_13,pegel09_13]=pegeldaten2013
[Tmeanh]=berechnung_temperaturmean

function [year,jday,month,day,hour,minute,T,RH,Swin]=madritsch_met
load('metinput_mad_10.mat')
%sdate=datenum(metinput_mad_10(:,[1,3,4,5,6]));
year=metinput_mad_10(:,1);
jday=metinput_mad_10(:,2);
month=metinput_mad_10(:,3);
day=metinput_mad_10(:,4);
hour=metinput_mad_10(:,5);
minute=metinput_mad_10(:,6);
T=metinput_mad_10(:,7);
RH=metinput_mad_10(:,8);
Swin=metinput_mad_10(:,9);




function [yyyy04,mm04,dd04,pegel07_04,pegel08_04,pegel09_04,pegel10_04,pegel11_04,pegel12_04]=pegeldaten2004
[yyyy04,mm04,dd04,pegel07_04,pegel08_04,pegel09_04,pegel10_04,pegel11_04,pegel12_04]=textread('pegeldaten_LGF2004_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f %f',-1,'headerlines',1);
%einlesen des files zufritt
%berechnung datumsformat
date=[yyyy04 mm04 dd04];
% sdate=datenum(date);
% sdatelim=[datenum(1973,01,01) datenum(2013,12,31)];
% sdatelim1=[datenum(2004,01,01) datenum(2013,12,31)];

function [yyyy05,mm05,dd05,pegel07_05,pegel0705_05,pegel08_05,pegel0805_05,pegel09_05,pegel10_05,pegel11_05,pegel12_05]=pegeldaten2005
[yyyy05,mm05,dd05,pegel07_05,pegel0705_05,pegel08_05,pegel0805_05,pegel09_05,pegel10_05,pegel11_05,pegel12_05]=textread('pegeldaten_LGF2005_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f %f %f %f',-1,'headerlines',1);
date=[yyyy05,mm05,dd05];

function [yyyy06,mm06,dd06,pegel07_06,pegel08_06,pegel09_06,pegel10_06,pegel11_06,pegel12_06]=pegeldaten2006
[yyyy06,mm06,dd06,pegel07_06,pegel08_06,pegel09_06,pegel10_06,pegel11_06,pegel12_06]=textread('pegeldaten_LGF2006_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f %f',-1,'headerlines',1);

function [yyyy07,mm07,dd07,pegel07_07,pegel08_07,pegel09_07,pegel10_07,pegel11_07,pegel12_07]=pegeldaten2007
[yyyy07,mm07,dd07,pegel07_07,pegel08_07,pegel09_07,pegel10_07,pegel11_07,pegel12_07]=textread('pegeldaten_LGF2007_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f %f',-1,'headerlines',1);

function [yyyy08,mm08,dd08,pegel07_08,pegel08_08,pegel09_08,pegel10_08,pegel11_08,pegel12_08]=pegeldaten2008
[yyyy08,mm08,dd08,pegel07_08,pegel08_08,pegel09_08,pegel10_08,pegel11_08,pegel12_08]=textread('pegeldaten_LGF2008_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f %f',-1,'headerlines',1);

function [yyyy09,mm09,dd09,pegel07_09,pegel08_09,pegel09_09,pegel10_09,pegel11_09,pegel12_09]=pegeldaten2009
[yyyy09,mm09,dd09,pegel07_09,pegel08_09,pegel09_09,pegel10_09,pegel11_09,pegel12_09]=textread('pegeldaten_LGF2009_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f %f',-1,'headerlines',1);

function [yyyy10,mm10,dd10,pegel07_10,pegel08_10]=pegeldaten2010
[yyyy10,mm10,dd10,pegel07_10,pegel08_10]=textread('pegeldaten_LGF2010_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f',-1,'headerlines',1);

function [yyyy11,mm11,dd11,pegel07_11,pegel08_11,pegel09_11,pegel10_11,pegel11_11]=pegeldaten2011
[yyyy11,mm11,dd11,pegel07_11,pegel08_11,pegel09_11,pegel10_11,pegel11_11]=textread('pegeldaten_LGF2011_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f',-1,'headerlines',1);

function [yyyy12,mm12,dd12,pegel07_12,pegel08_12,pegel09_12,pegel10_12,pegel11_12,pegel12_12]=pegeldaten2012
[yyyy12,mm12,dd12,pegel07_12,pegel08_12,pegel09_12,pegel10_12,pegel11_12,pegel12_12]=textread('pegeldaten_LGF2012_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f %f %f %f',-1,'headerlines',1);

function [yyyy13,mm13,dd13,pegel07_13,pegel08_13,pegel09_13]=pegeldaten2013
[yyyy13,mm13,dd13,pegel07_13,pegel08_13,pegel09_13]=textread('pegeldaten_LGF2013_pegel7bis12.txt','%f %*1c %f %*1c %f %f %f %f',-1,'headerlines',1);

function [Tmeanh]=berechnung_temperaturmean(year,jday,month,day,hour,minute,T,RH,Swin)
%date=[year month day]
%sdate=datenum(sdate);

ss = 0 
sdate=datenum(year,month,day,hour,minute,ss); 
sdatelim=[datenum(2004,07,15,0,0,0) datenum(2004,08,04,23,59,0)]
dt=1/24;
sdate1=sdatelim(1):dt:sdatelim(2);
for i = 1:length(sdate1)
    j=find(sdate >= sdate1(i) & sdate < sdate1(i)+dt);
    Tmeanh(i)=mean(T(j));
end
%http://www.gomatlab.de/viewtopic,p,141157.html#141157