-
Notifications
You must be signed in to change notification settings - Fork 0
/
_2a_long_timeseries-Hofstetten.R
148 lines (114 loc) · 4.23 KB
/
_2a_long_timeseries-Hofstetten.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
library(tidyverse)
library(lubridate)
################### runoff
path <- "/home/christoph/Dokumente/BOKU/Masterarbeit/Daten/Stationsdaten"
if( .Platform$OS.type == "windows" )
path <- "C:/Users/Russ/Desktop/master/daten/Stationsdaten/"
setwd(path)
fileq <- "Q207852-TM-Hofstetten(Bad).dat"
rff <- read_table(fileq, col_names = F, skip = 27, cols(
X1 = col_date(format = "%d.%m.%Y"),
X2 = col_time(format = ""),
X3 = col_double()
))
#subsetting
rff <- rff[as_date(rff$X1)<as_date("2016-01-01"), ]
rff <- rff[as_date(rff$X1)>as_date("1990-12-31"), ]
mean(na.omit(rff$X3))
# plot
plot_rff<- ggplot(rff,aes(x =X1 , y = X3)) +
xlab("Time [d]") + ylab("Runoff [m³/s]") +
ylim(0,60)+
ggtitle("Runoff (Gauge Hofstetten)")+
geom_line(stat="identity", size=0.05) +
geom_hline(yintercept=mean(rff$X3), linetype="dashed", color="red",size=1.2) +
geom_hline(yintercept=6.47, color="blue",size=1)
plot_rff
################## air temperature
path <- "/home/christoph/Dokumente/BOKU/Masterarbeit/Daten/Stationsdaten"
if( .Platform$OS.type == "windows" )
path <- "C:/Users/Russ/Desktop/master/daten/Stationsdaten/"
setwd(path)
fileat <- "LT107300.dat"
at <- read_table(fileat, col_names = F, skip = 20, cols(
X1 = col_date(format = "%d.%m.%Y"),
X2 = col_time(format = ""),
X3 = col_double()
))
#subsetting
at <- at[as_date(at$X1)<as_date("2016-01-01"), ]
plot_at<- ggplot(at,aes(x =X1 , y = X3)) +
xlab("time [d]") + ylab("air temperature") +
expand_limits(y=c(-10,20))+
ggtitle("Station 107300")+
geom_line(stat="identity", size=0.05)
plot_at
################## Water temperature
path <- "/home/christoph/Dokumente/BOKU/Masterarbeit/Daten/Stationsdaten"
if( .Platform$OS.type == "windows" )
path <- "C:/Users/Russ/Desktop/master/daten/Stationsdaten/"
setwd(path)
filewt <- "WT-Tagesmitte-Hofstetten(Bad).dat"
wt <- read_table(filewt, col_names = F, skip = 31, cols(
X1 = col_date(format = "%d.%m.%Y"),
X2 = col_time(format = ""),
X3 = col_double()
))
wt <- wt[as_date(wt$X1)<as_date("2016-01-01"), ]
wt <- wt[as_date(wt$X1)>as_date("1990-12-31"), ]
# #subsetting
# wt <- wt[as_date(wt$X1)>as_date("2008-01-01"), ]
# plot
plot_wt<- ggplot(wt,aes(x =X1 , y = X3)) +
xlab("Time [d]") + ylab("Water temperature [°C]") +
ggtitle("Water Temperature (Gauge Hofstetten)")+
geom_line(stat="identity", size=0.05)
plot_wt # 864 missing observations
####################################################
# 2 years of data
path <- "/home/christoph/Dokumente/BOKU/Masterarbeit/Daten/Stationsdaten"
if( .Platform$OS.type == "windows" )
path <- "C:/Users/Russ/Desktop/master/daten/Stationsdaten/"
setwd(path)
file <- "LT115642.dat"
at <- read_table(file, col_names = F, skip = 24, cols(
X1 = col_date(format = "%d.%m.%Y"),
X2 = col_time(format = ""),
X3 = col_double()
))
#subsetting
at2 <- at[as_date(at$X1)<as_date("2015-01-01"), ]
at2 <- at2[as_date(at2$X1)>as_date("2013-01-01"), ]
# plot
plot_at<- ggplot(at2,aes(x =X1 , y = X3)) +
xlab("time") + ylab("air temperature") +
expand_limits(y=c(-10,30))+
ggtitle("air temperature curve")+
geom_line(stat="identity")
plot_at
####
path <- "/home/christoph/Dokumente/BOKU/Masterarbeit/Daten/Stationsdaten"
if( .Platform$OS.type == "windows" )
path <- "C:/Users/Russ/Desktop/master/daten/Stationsdaten/"
setwd(path)
file <- "WT-Tagesmitte-Hofstetten(Bad).dat"
wt <- read_table(file, col_names = F, skip = 31, cols(
X1 = col_date(format = "%d.%m.%Y"),
X2 = col_time(format = ""),
X3 = col_double()
))
#subsetting
wt <- wt[as_date(wt$X1)<as_date("2015-01-01"), ]
wt <- wt[as_date(wt$X1)>as_date("2013-01-01"), ]
# plot
plot_wt<- ggplot(wt,aes(x =X1 , y = X3)) +
xlab("time") + ylab("water temperature") +
expand_limits(y=c(-10,30))+
ggtitle("water temperature curve")+
geom_line(stat="identity")
plot_wt
#### gemeinsam in einem plot
plot_at + geom_line(data = wt, aes(x=X1, y=X3), stat="identity", colour="red")+
xlab("Time [d]") + ylab("Temperature [°C]") +
expand_limits(y=c(-10,30))+
ggtitle("Air and Water Temperature (Gauge Hofstetten)")