Untitled
Roilan Hernández Valdés
July 26, 2015
pacotes <- c("knitr","raster", "dygraphs", "openair","EcoHydRology")
invisible(sapply(pacotes, FUN = require, character.only = TRUE, quietly = TRUE))##
## Attaching package: 'operators'
##
## The following object is masked from 'package:dygraphs':
##
## %>%
##
## The following object is masked from 'package:base':
##
## options
##
##
## DEoptim package
## Differential Evolution algorithm in R
## Authors: D. Ardia, K. Mullen, B. Peterson and J. Ulrichsim.dir <- "~/ADBHM/pauliceia/simctl_roilan_marcelo/SOURCE/R_SCRIPTS"
setwd(sim.dir )
start.date <- as.POSIXct("2009-01-01", tz = "GMT")
end.date <- as.POSIXct("2014-12-31", tz = "GMT")
date <- seq.POSIXt(start.date, end.date, by = "day")Condições iniciais
Rasters da bacia
# Só precisamos uma rodada para os raster de condições inicias
list.Simulacoes <- basename(list.dirs(path = "../../run", full.names = F, recursive = F))
sims <- as.numeric(gsub(pattern = "_",replacement = "",x = substr(list.Simulacoes, start = 1,stop = 3)))
list.Simulacoes <- list.Simulacoes[order(sims)]
list.Simulacoes## [1] "1_____0____1____0____1" "2_____1____1____0____1"
## [3] "3_____2____1____0____1" "4_____3____1____0____1"
## [5] "5_____0____6____0____1" "6_____1____6____0____1"
## [7] "7_____2____6____0____1" "8_____3____6____0____1"
## [9] "9_____0___11____0____1" "10____1___11____0____1"
## [11] "11____2___11____0____1" "12____3___11____0____1"
## [13] "13____0___16____0____1" "14____1___16____0____1"
## [15] "15____2___16____0____1" "16____3___16____0____1"
## [17] "17____0___21____0____1" "18____1___21____0____1"
## [19] "19____2___21____0____1" "20____3___21____0____1"
## [21] "21____0___26____0____1" "22____1___26____0____1"
## [23] "23____2___26____0____1" "24____3___26____0____1" basin <- raster(paste("../../run/",list.Simulacoes[1],"/basin1.asc", sep = ""))
plot(basin)
wsdem <- raster(paste("../../run/",list.Simulacoes[1],"/wsdem_pauliceia.asc", sep = ""))
plot(wsdem)
Rasters criados dos parametros
# list.derive <- list.files(path = "../../Derive_D", pattern = ".asc$" , full.names = T)
# del <- grep(pattern = "stn", x = list.derive)
# list.derive <- list.derive[-del]
#
# derive <- stack(list.derive)
# plot(derive[[c(2,5,6,8:10)]])
#
# all.basin <- stack(bacia,derive)Vazão no exutório
vazao <- read.table(file = paste("../../run/",list.Simulacoes[1],"/Exutorio_d.txt", sep = ""),
header = T,
na.strings = "-9999.0000000000000000")
vazao <- data.frame(date, Qd_obv = vazao$Qd_obv * 24 * 3.6e6 / 11.94e6)
timePlot(vazao, "Qd_obv")
vaza0 <- sapply(1:length(list.Simulacoes), function(i){
read.table(file = paste("../../run/",list.Simulacoes[i],
"/Exutorio_d.txt", sep = ""),
header = T,
na.strings = "-9999.0000000000000000")[,4] * 24 * 3.6e6 / 11.94e6
})
vaza0 <- as.data.frame(vaza0)
names(vaza0) <- list.Simulacoes
vazao <- cbind(vaza0,vazao)
str(vazao)## 'data.frame': 2191 obs. of 26 variables:
## $ 1_____0____1____0____1: num 27.9 7.5 6.18 6.38 3.79 ...
## $ 2_____1____1____0____1: num 37.8 26.4 27.3 41.4 28.2 ...
## $ 3_____2____1____0____1: num 23.48 5.59 2.72 1.97 2.47 ...
## $ 4_____3____1____0____1: num 28.9 17.4 14 25.4 19 ...
## $ 5_____0____6____0____1: num 24.46 6.54 4.36 4.38 4.24 ...
## $ 6_____1____6____0____1: num 28.08 14.38 19.65 14.95 9.79 ...
## $ 7_____2____6____0____1: num 25.05 8.99 6.32 6.03 6.98 ...
## $ 8_____3____6____0____1: num 24.93 9.03 6.99 8.12 7.62 ...
## $ 9_____0___11____0____1: num 25.39 8.32 5.77 5.55 4.95 ...
## $ 10____1___11____0____1: num 26.11 12.35 8.82 8.27 8.55 ...
## $ 11____2___11____0____1: num 26.05 11.76 8.76 7.69 8.34 ...
## $ 12____3___11____0____1: num 26.63 7.85 5.59 10.52 7.7 ...
## $ 13____0___16____0____1: num 38.36 9.23 6.75 6.16 5.33 ...
## $ 14____1___16____0____1: num 27.2 14.8 10 56.4 16.3 ...
## $ 15____2___16____0____1: num 27.26 14.09 9.69 8.85 8.66 ...
## $ 16____3___16____0____1: num 27.94 14.37 11.26 10.64 7.26 ...
## $ 17____0___21____0____1: num 27.4 10.88 6.97 6.37 5.92 ...
## $ 18____1___21____0____1: num 101.95 21.33 6.89 6.5 6.05 ...
## $ 19____2___21____0____1: num 28.32 15.83 10.82 9.15 9.09 ...
## $ 20____3___21____0____1: num 28.24 16.06 39.1 11.65 7.25 ...
## $ 21____0___26____0____1: num 28.43 11.3 11.44 6.46 6.27 ...
## $ 22____1___26____0____1: num 30.02 17.77 12.32 9.93 10.99 ...
## $ 23____2___26____0____1: num 29.67 17.26 10.99 9.47 9.64 ...
## $ 24____3___26____0____1: num 29.65 17.12 81.81 19.92 8.59 ...
## $ date : POSIXct, format: "2009-01-01" "2009-01-02" ...
## $ Qd_obv : num NA NA NA NA NA NA NA NA NA NA ... # Arquivo que contem a descripção das configurações
rodadas <- read.table("../../rodadas", header = F)
names(rodadas) <- c("ID", "Tprec", "K_sFac", "Decay", "GHcFac")Vazão acumulada
cumvazao <- cumsum(vazao[!is.na(vazao$Qd_obv),names(vazao) != "date"])
cumvazao$date <- date[!is.na(vazao$Qd_obv)] vars <- which(rodadas$Tprec == 0)
timePlot(vazao, c("Qd_obv",names(vazao[,vars])),
group = T, stack = F, lty = 1,
name.pol =c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão", xlab = "Date",
main = "Vazão [mm], distribuição Tang09", log = T
)
timePlot(cumvazao, c("Qd_obv", names(cumvazao[,vars])),
group = T, stack = F, lty = 1, lwd = 2,
name.pol = c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão acumulada" , xlab = "Date",
main = "Vazão acumulada [mm], distribuição Tang09"
)
vars <- which(rodadas$Tprec == 1)
timePlot(vazao, c("Qd_obv",names(vazao[,vars])),
group = T, stack = F, lty = 1,
name.pol =c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão", xlab = "Date",
main = "Vazão [mm], distribuição 1 hora", log = T
)
timePlot(cumvazao, c("Qd_obv", names(cumvazao[,vars])),
group = T, stack = F, lty = 1, lwd = 2,
name.pol = c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão acumulada" , xlab = "Date",
main = "Vazão acumulada [mm], distribuição 1 hora"
) 
vars <- which(rodadas$Tprec == 2)
timePlot(vazao, c("Qd_obv",names(vazao[,vars])),
group = T, stack = F, lty = 1,
name.pol =c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão", xlab = "Date",
main = "Vazão [mm], distribuição 24 horas" , log = T
)
timePlot(cumvazao, c("Qd_obv", names(cumvazao[,vars])),
group = T, stack = F, lty = 1, lwd = 2,
name.pol = c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão acumulada" , xlab = "Date",
main = "Vazão acumulada [mm], distribuição 24 horas"
) 
vars <- which(rodadas$Tprec == 3)
timePlot(vazao, c("Qd_obv",names(vazao[,vars])),
group = T, stack = F, lty = 1,
name.pol = c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão", xlab = "Date",
main = "Vazão [mm], distribuição PRatio" , log = T
)
timePlot(cumvazao, c("Qd_obv", names(cumvazao[,vars])),
group = T, stack = F, lty = 1, lwd = 2,
name.pol = c("Observada",paste("Ks",as.character(rodadas[vars,3]*3.3))),
ylab = "Vazão acumulada" , xlab = "Date",
main = "Vazão acumulada [mm], distribuição PRatio"
) 
vars <- which(rodadas$K_sFac == 1)
timePlot(vazao, c("Qd_obv",names(vazao[,vars])),
group = T, stack = F, lty = 1,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[vars,2]))),
ylab = "Vazão", xlab = "Date",
main = paste("Vazão [mm], Condutividade Hidraúlica",unique(rodadas[vars,3])*3.3) , log = T
)
timePlot(cumvazao, c("Qd_obv", names(cumvazao[,vars])),
group = T, stack = F, lty = 1, lwd = 2,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[vars,2]))),
ylab = "Vazão acumulada" , xlab = "Date",
main = paste("Vazão acumulada [mm], Condutividade Hidraúlica",unique(rodadas[vars,3])*3.3)
) 
vars <- which(rodadas$K_sFac == 16)
timePlot(vazao, c("Qd_obv",names(vazao[,vars])),
group = T, stack = F, lty = 1,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[vars,2]))),
ylab = "Vazão", xlab = "Date",
main = paste("Vazão [mm], Condutividade Hidraúlica",unique(rodadas[vars,3])*3.3) , log = T
)
timePlot(cumvazao, c("Qd_obv", names(cumvazao[,vars])),
group = T, stack = F, lty = 1, lwd = 2,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[vars,2]))),
ylab = "Vazão acumulada" , xlab = "Date",
main = paste("Vazão acumulada [mm], Condutividade Hidraúlica",unique(rodadas[vars,3])*3.3)
)
vars <- which(rodadas$K_sFac == 26)
timePlot(vazao, c("Qd_obv",names(vazao[,vars])),
group = T, stack = F, lty = 1,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[vars,2]))),
ylab = "Vazão", xlab = "Date",
main = paste("Vazão [mm], Condutividade Hidraúlica",unique(rodadas[vars,3])*3.3) , log = T
)
timePlot(cumvazao, c("Qd_obv", names(cumvazao[,vars])),
group = T, stack = F, lty = 1, lwd = 2,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[vars,2]))),
ylab = "Vazão acumulada" , xlab = "Date",
main = paste("Vazão acumulada [mm], Condutividade Hidraúlica",unique(rodadas[vars,3])*3.3)
) 
# table(is.na(vazao.exutorio$Qd))
#
# escoamentos <- BaseflowSeparation(vazao.exutorio$Qd)
# vazao.exutorio$bt_s <- escoamentos$bt
# vazao.exutorio$qft_s <- escoamentos$qft
#
# # Qd_obv <- esc
# # escoamentos <- BaseflowSeparation(Qd_obv)
#
#
# timePlot(vazao.exutorio, c("bt", "Qd_obv"))