hydrology/discharge.R at main · dmkahler/hydrology · GitHub

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
# analyze river flow in Kruger
library(readr)
library(dplyr)
library(tidyr)
library(lubridate)
library(ggplot2)
library(latex2exp)

mamba <- read_csv("Olifants_Mamba_B7H015.csv", skip = 9, n_max = 12174) # upstream
mamba <- mamba %>%
      mutate(dn = as.numeric(date)) %>%
      select(-date) %>%
      mutate(year = floor(dn/10000)) %>%
      mutate(mon = floor((dn-year*10000)/100)) %>%
      mutate(day = dn-(10000*year+100*mon)) %>%
      mutate(dt = ymd(dn)) %>%
      select(-dn)

balule <- read_csv("Olifants_Balule_B7H026.csv", skip = 9) # downstream, near Olifants Camp
balule <- balule %>%
      mutate(dn = as.numeric(date)) %>%
      select(-date) %>%
      mutate(year = floor(dn/10000)) %>%
      mutate(mon = floor((dn-year*10000)/100)) %>%
      mutate(day = dn-(10000*year+100*mon)) %>%
      mutate(dt = ymd(dn)) %>%
      select(-dn)

# Start comparison table
start <- as.numeric(min(c(min(mamba$dt), min(balule$dt))))
stop <- as.numeric(max(c(max(mamba$dt), max(mamba$dt))))
dur <- as.numeric(stop) - as.numeric(start) + 1
dn <- c(start:stop)
mamba.discharge <- array(NA, dim = dur)
balule.discharge <- array(NA, dim = dur)
for (i in 1:nrow(mamba)) {
      dateposition <- as.numeric(mamba$dt[i]) - start + 1
      mamba.discharge[dateposition] <- mamba$discharge[i]/100
}
for (i in 1:nrow(balule)) {
      dateposition <- as.numeric(balule$dt[i]) - start + 1
      balule.discharge[dateposition] <- balule$discharge[i]/100
}
comp <- data.frame(dn, mamba.discharge, balule.discharge)
comp.l <- rename(comp, Mamba = mamba.discharge, Balule = balule.discharge)
comp.l <- pivot_longer(comp.l, cols = c(Mamba, Balule), names_to = "Site", values_to = "Discharge")

# plot
ggplot(comp, aes(x = as_date(dn))) +
      geom_point(aes(y = mamba.discharge)) +
      geom_point(aes(y = balule.discharge))

ggplot(comp, aes(x = as_date(dn), y = Discharge, color = Site)) +
      geom_line(linetype = "dotted") +
      xlab("Date") +
      ylab(TeX('Discharge $(m^3/s)$')) +
      theme(panel.background = element_rect(fill = "white", colour = "black")) +
      theme(legend.position="right") +
      theme(aspect.ratio = 1) +
      theme(axis.text = element_text(face = "plain", size = 12))