#Hohenpeissenberg, Temperatur, Monatlich, mit 100-Jahres Monatsmittel 
#R: text( ) nur innerhalb des Plot-Rahmens 

hohve<- scan("hohtpm.dat")

postscript(file="HPmo.ps",height=18,width=10,horizontal=F)
par(mfrow=c(2,1))

hohve<- hohve/10
hohv<- hohve[1021:1200]

length(hohv)
hohme<- mean(hohve); hohme
#Monatsmittel
hohma<- matrix(scan("hohtpm.dat"),ncol=12,byrow=T)
hohma<- hohma/10 
 
mo<- 1:12
for (j in 1:12)
{mo[j]<- mean(hohma[,j])}
mo 
mose<- rep(mo,times=4)

seq<- 1:180
plot(seq,hohv,type="l",lty=1,xlim=c(1,180),ylim=c(-5,20),
   xlab="Monat Jan 1905 - Dez 2004",ylab="Temperatur [Celsius]")
title(main="Hohenpeissenberg, Mittlere Monatstemperaturen",cex=0.85)
abline(h=hohme,lty=3)
text(150,-4.5,"- - - - - - -   100-Jahres Mittel",cex=0.65)
abline(h=c(0,5,10,15),lty=2)
text(7,-7.5,"1990",cex=0.7); text(187,-7.5,"2005",cex=0.7)
text(67,-7.5,"1995",cex=0.7); text(127,-7.5,"2000",cex=0.7)

#die letzten 4 Jahre
seq<- 1:48
su<- seq; yt<- seq

for (t in 1:48)
{ yt[t]<-  hohve[t+1152]}
 
for (t in 7:42)
{su[t]<- 0
{ for (k in -5:5)
 su[t]<-  su[t]+ hohve[t+1152+k] }}
 
plot(seq,yt,type="l",lty=1,xlim=c(1,48),ylim=c(-5,20),
  xlab="Monat Jan 2001 - Dez 2004",ylab="Temperatur [Celsius]")
points(seq,yt,pch=4,cex=0.6)
abline(h=c(0,5,10,15),lty=2)
lines(seq,mose,lty=2)
text(38,-4.5,".........   100-Jahre Monatsmittel",cex=0.65)
text(7,-1,"2001",cex=0.7); text(19,-1,"2002",cex=0.7)
text(31,-1,"2003",cex=0.7); text(43,-1,"2004",cex=0.7)

dev.off()
rm(hohve,hohv,hohme)
rm(yt,su,seq)
rm(hohma,mo,mose)