require(openair)
## Loading required package: openair
require(knitr)
## Loading required package: knitr
setwd("~/ADBHM/pauliceia/simctl_roilan_marcelo/SOURCE/R_SCRIPTS/")
date <- seq.POSIXt(from = as.POSIXct("2009-01-01 00:00:00"),
to = as.POSIXct("2014-12-31 23:00:00"),
by = "hour", tz = "GMT")
dados <- readRDS("~/Dissertação/DADOS/nee_hh_r_qc_flags_pdg.rds")
vars <- dados[,c("date","swc1m", "swc2m" )]
vars$sum.int <- vars$swc2m + vars$swc1m
timePlot(vars, c( "swc1m", "swc2m", "sum.int"), ylab = "Conteúdo de água [mm]", xlab = "Data",
group = T, lty = 1, lwd = 2,
name.pol = c("Camada 1 m", "Camada 2 m", "Suma integrada"),
main = "Conteúdo de água no solo")
umid1m <- aggregate(vars$swc1m,
by = list(format(vars$date, "%Y-%m-%d %H")),
FUN = mean,
na.rm = F)
range(umid1m$Group.1)
## [1] "2009-07-08 00" "2014-07-01 00"
#Preenchendo com NA as datas sem observações
umid1m <- c(rep(NA, times = 182 * 24), umid1m$x)
umid1m <- c(umid1m, rep(NA, times = length(date) - length(umid1m)))
length(umid1m) == length(date)
## [1] TRUE
#Lendo o nome das pastas das rodadas e organizando por o numero no arquivo
# "rodadas"
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"
# Lendo os arquivos de umidade do solo.
umidade <- sapply(1:length(list.Simulacoes),
function(i){ # i = 12
read.table(file = paste("../../run/",list.Simulacoes[i],
"/OutStn_www.txt", sep = ""),
header = F,skip = 1,
na.strings = "-9999.0000000000000000")[,6] * 10})
umidade <- as.data.frame(umidade)
umidade$date <- date
umidade$obv <- umid1m
# str(umidade)
names(umidade) <- gsub("V", "",names(umidade))
rodadas <- read.table(file = "~/Dropbox/rodadas", header = F)
rodadas
## V1 V2 V3 V4 V5
## 1 1 0 1 0 1
## 2 2 1 1 0 1
## 3 3 2 1 0 1
## 4 4 3 1 0 1
## 5 5 0 6 0 1
## 6 6 1 6 0 1
## 7 7 2 6 0 1
## 8 8 3 6 0 1
## 9 9 0 11 0 1
## 10 10 1 11 0 1
## 11 11 2 11 0 1
## 12 12 3 11 0 1
## 13 13 0 16 0 1
## 14 14 1 16 0 1
## 15 15 2 16 0 1
## 16 16 3 16 0 1
## 17 17 0 21 0 1
## 18 18 1 21 0 1
## 19 19 2 21 0 1
## 20 20 3 21 0 1
## 21 21 0 26 0 1
## 22 22 1 26 0 1
## 23 23 2 26 0 1
## 24 24 3 26 0 1
# Simulações com distribuição Tang09
hour.ind <- rodadas[rodadas[,2] == 0,1]
timePlot(umidade, c("obv",as.character(hour.ind)), group = T, lty = 1,
avg.time = "day", stack = F ,
name.pol = c("Observada",paste("Ks",as.character(rodadas[hour.ind,3] * 3.36))),
ylab = "Conteúdo de agua no solo", xlab = "Data",
main = "Conteúdo de água na camada 1 metro, distribuição Tang09", auto.text = T
)
# Simulações com as distribuição 1 hora
hour.ind <- rodadas[rodadas[,2] == 1,1]
timePlot(umidade, c("obv",as.character(hour.ind)), group = T, lty = 1,
avg.time = "day", stack = F ,
name.pol = c("Observada",paste("Ks",as.character(rodadas[hour.ind,3] * 3.36))),
ylab = "Conteúdo de agua no solo", xlab = "Data",
main = "Conteúdo de água na camada 1 metro, distribuição 1 hora", auto.text = T
)
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
# Simulações com as distribuição media em 24 horas
hour.ind <- rodadas[rodadas[,2] == 2,1]
timePlot(umidade, c("obv",as.character(hour.ind)), group = T, lty = 1,
avg.time = "day", stack = F ,
name.pol = c("Observada",paste("Ks",as.character(rodadas[hour.ind,3] * 3.36))),
ylab = "Conteúdo de agua no solo", xlab = "Data",
main = "Conteúdo de água na camada 1 metro, distribuição 24 hora", auto.text = T
)
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
# Simulações com a condutividade hidraulica inicial
hour.ind <- rodadas[rodadas[,3] == 1,1]
timePlot(umidade, c("obv",as.character(hour.ind)), group = T, lty = 1,
avg.time = "day", stack = F ,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[hour.ind,2]))),
ylab = "Conteúdo de agua no solo", xlab = "Data",
main = paste("Conteúdo de água na camada 1 metro, ks=", rodadas[hour.ind[1],3] * 3.3), auto.text = T
)
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
# Simulações com a condutividade hidraulica inicial
hour.ind <- rodadas[rodadas[,3] == 16,1]
timePlot(umidade, c("obv",as.character(hour.ind)), group = T, lty = 1,
avg.time = "day", stack = F ,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[hour.ind,2]))),
ylab = "Conteúdo de agua no solo", xlab = "Data",
main = paste("Conteúdo de água na camada 1 metro, ks=", rodadas[hour.ind[1],3] * 3.3), auto.text = T
)
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
# Simulações com a condutividade hidraulica inicial
hour.ind <- rodadas[rodadas[,3] == 26,1]
timePlot(umidade, c("obv",as.character(hour.ind)), group = T, lty = 1,
avg.time = "day", stack = F ,
name.pol = c("Observada",paste("Desagregação",as.character(rodadas[hour.ind,2]))),
ylab = "Conteúdo de agua no solo", xlab = "Data",
main = paste("Conteúdo de água na camada 1 metro, ks=", rodadas[hour.ind[1],3] * 3.3), auto.text = T
)
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
