Sensibilidade à altura da torre de medida
Introdução
As alturas zwind e zmet, são definidas como a altura de referência da medida da velocidade do vento (reference height for wind measurement) e altura de referência das medidas de umidade e temperatura (reference height for temperature, humidity measurement), respectivamente.
Com o presente se planeja apresentar a influência de zwind e zmet na saídas dos fluxos de energia.
Conserto no código:
Na subrotina VARCAL, as seguentes linhas preescreviam sempre zwind e zmet a 45 metros.
if (z2 + 10.0 * dd + 0.1 < 2.0) then
zwind = 2.0 ! (data from Chinese Weather Bureau)
zmet = 2.0 ! (usually it is 2.00 m)
else
!c! added by tatsch sep/2010, if vtype=1 this give unrealistic values
!c! zwind = z2+10.*dd+0.1 ! (data from Chinese Weather Bureau)
!c! zmet = z2+10.*dd+0.1 ! (usually it is 2.00 m)
zwind = 45.00 ! (data from Chinese Weather Bureau)
zmet = 45.00 ! (usually it is 2.00 m)
endifSolução proposta
if (z2 + 10.0 * dd + 0.1 < 2.0) then
zwind = 2.0 ! (data from Chinese Weather Bureau)
zmet = 2.0 ! (usually it is 2.00 m)
! else
!c! added by tatsch sep/2010, if vtype=1 this give unrealistic values
!c! zwind = z2+10.*dd+0.1 ! (data from Chinese Weather Bureau)
!c! zmet = z2+10.*dd+0.1 ! (usually it is 2.00 m)
! zwind = 45.00 ! (data from Chinese Weather Bureau)
! zmet = 45.00 ! (usually it is 2.00 m)
endifQue fazer?
Se fizeram 5 rodadas do SiB2:
- zmet 45 e zwind 45 , a que seria até o momento como referência.
zmet45_zwind45 <- do.call(sib2,
as.list(modifyList(formals(sib2),
list(date1=date1,date2=date2,
infile2 = infile,
zmet = 45,
zwind = 45,
vars_out=c("LE","H","rn","ustar"),
dir_out = dir_out)))
) %>% selectByDate(start="2009-07-02",end="2012-06-30")- zmet 21 e zwind 21 , altura de medida de temperatura e umidade, além de precipitação Cabral et al., 2015 e velocidade do vento da Rocha et al., 2002
zmet21_zwind21 <- do.call(sib2,
as.list(modifyList(formals(sib2),
list(date1=date1,date2=date2,
infile2 = infile,
zmet = 21,
zwind = 21,
vars_out=c("LE","H","rn","ustar"),
dir_out = dir_out)))
) %>% selectByDate(start="2009-07-02",end="2012-06-30")- zmet 11 e zwind 11 , altura da torre de medida dos fluxos (H e LE) Cabral et al., 2015
zmet11_zwind11 <- do.call(sib2,
as.list(modifyList(formals(sib2),
list(date1=date1,date2=date2,
infile2 = infile,
zmet = 11,
zwind = 11,
vars_out=c("LE","H","rn","ustar"),
dir_out = dir_out)))
) %>% selectByDate(start="2009-07-02",end="2012-06-30")- zmet 11 e zwind 21 , para testar sensibilidade à diferença da altura dos instrumentos no SiB2
zmet11_zwind21 <- do.call(sib2,
as.list(modifyList(formals(sib2),
list(date1=date1,date2=date2,
infile2 = infile,
zmet = 11,
zwind = 21,
vars_out=c("LE","H","rn","ustar"),
dir_out = dir_out)))
) %>% selectByDate(start="2009-07-02",end="2012-06-30")- zmet 21 e zwind 11 , para testar sensibilidade à diferença da altura dos instrumentos no SiB2
zmet21_zwind11 <- do.call(sib2,
as.list(modifyList(formals(sib2),
list(date1=date1,date2=date2,
infile2 = infile,
zmet = 21,
zwind = 11,
vars_out=c("LE","H","rn","ustar"),
dir_out = dir_out)))
) %>% selectByDate(start="2009-07-02",end="2012-06-30")Comparando com observação
RN
LE
H
Ustar
Sensibilidade com LHOAT
Inicialização do LHOAT
### Goodness-of-fit function, either customized or pre-defined from hydroGOF
gof.Ini="2009-07-02"
gof.Fin="2012-06-30"
gof.FUN <- "NSE"
gof.FUN.args <- list()
dir_out = paste0("/home/",system("echo $USER",intern = TRUE), "/ADBHM/trash")
model.FUN.args <- as.list(modifyList(formals(sib2),
list(date1=date1,
date2=date2,
infile2 = infile,
vars_out= "rn",
dir_out = dir_out)))
param.ranges <- "./PSO.in/ParamRanges-Zs.txt"Execução
lhoat4sib2(
fn="sib2"
,
model.FUN="sib2"
,
model.FUN.args=model.FUN.args
,
obs = list(
gof.Ini = gof.Ini,
gof.Fin = gof.Fin,
gof.obs = obs)
,
control=list(
N=10,
f=0.25,
drty.out="LH_OAT_rn",
param.ranges=param.ranges,
gof.name=gof.FUN,
REPORT = 10,
do.plots=FALSE,
write2disk=TRUE,
parallel = "multicore",
verbose= TRUE) ###END control options
) -> out.rn
# file.remove("./LH_OAT_rn/LH_OAT-out.txt")
saveRDS(out,"./LH_OAT_rn/out.rds")Intervalo definido
| MinValue | MaxValue | |
|---|---|---|
| zwind | 11 | 21 |
| zmet | 11 | 21 |
LHOAT para o Saldo de Radiação
| RankingNmbr | ParameterName | RelativeImportance | RelativeImportance.Norm |
|---|---|---|---|
| 1 | zwind | 0.1260441 | 0.7477528 |
| 2 | zmet | 0.0425198 | 0.2522472 |
LHOAT para Evapotranspiração
| RankingNmbr | ParameterName | RelativeImportance | RelativeImportance.Norm |
|---|---|---|---|
| 1 | zmet | 3.132970 | 0.6414654 |
| 2 | zwind | 1.751113 | 0.3585346 |
LHOAT para Calor Sensível
| RankingNmbr | ParameterName | RelativeImportance | RelativeImportance.Norm |
|---|---|---|---|
| 1 | zmet | 7.842108 | 0.5963432 |
| 2 | zwind | 5.308220 | 0.4036568 |