path="/home/marc/dev/scilab/Tempelhof/";
f=mopen(path+"Tempelhofmonthtemp_2_.txt");

//---------------------------------------------------------
//-- Datenaufnahme aus TXT Datei
//---------------------------------------------------------
printf("\nDaten auslesen:\n");
[n,year,jan,feb,mar,apr,may,jun,jul,aug,sep,oct,nov,dec,annmean]...
       =mfscanf(-1,f,'%d %f %f  %f %f %f %f %f %f %f %f %f %f %f');
printf("Fertig mit Daten einlesen.\n");
mclose(f);

//---------------------------------------------------------
//-- Gesamtverlauf erfassen 
//-- ??? Keine Information ???
//---------------------------------------------------------
verlauf=[]; cc=0;
for m=1:size(jan,1) do
  cc=cc+1;  verlauf(cc)=jan(m);
  cc=cc+1;  verlauf(cc)=feb(m);
  cc=cc+1;  verlauf(cc)=mar(m);
  cc=cc+1;  verlauf(cc)=apr(m);
  cc=cc+1;  verlauf(cc)=may(m);
  cc=cc+1;  verlauf(cc)=jun(m);
  cc=cc+1;  verlauf(cc)=jul(m);
  cc=cc+1;  verlauf(cc)=aug(m);
  cc=cc+1;  verlauf(cc)=sep(m);
  cc=cc+1;  verlauf(cc)=oct(m);
  cc=cc+1;  verlauf(cc)=nov(m);
  cc=cc+1;  verlauf(cc)=dec(m);
end

//---------------------------------------------------------
//-- NAN Bereich ersetzen mit 10
//-- und Fehlerfeld erzeugen
//---------------------------------------------------------
fehler=zeros(size(annmean,1),1);
nan=find(annmean<-998);
annmean(nan)=10;
fehler(nan)=6;

//---------------------------------------------------------
//-- Exponentielle Glaettung
//---------------------------------------------------------
ex=[];
ex(1)=annmean(1);
alpha=0.1;
for i=2:size(annmean,1) do
  ex(i)=alpha*annmean(i)+(1-alpha)*ex(i-1);
end

//---------------------------------------------------------
//-- Graphik
//---------------------------------------------------------
clf();
plot2d(year,[annmean,ex,fehler],[3,2,5],rect=[year(1),5,year(size(year,1)),12]);
legend(['annmean','expo.','Fehler']);
    
