<div dir="ltr">Boa tarde, Éder!!<div><br></div><div>Mais uma vez, lhe agradeço pelas contribuições. Foram de grande valia.</div><div>Como havia comentado, segue o script que programei para fazer o calculo diário da ETP através dos métodos de Penman-Monteith (padronizado por Allen, 1998) e através do método proposto por Priestly-Taylor.</div><div><br></div><div>Abraços!</div><div><br></div><div>rm(list = ls())</div><div><br></div><div>setwd("DIRETÓRIO")</div><div>a = read.table("ARQUIVO_CLIMATICO'); head(a)</div><div>a$Data = as.Date(a$Data)</div><div>a$DiaJuliano = as.numeric(format(a$Data, trim = T, '%j')); head(a)</div><div><br></div><div>altitude = 259.38</div><div>Patm = 101.3*((293-0.0065*altitude)/293)**5.26</div><div>gama = 0.665*(10**-3)*Patm</div><div>latitude = -7.53</div><div>corr = pi/180</div><div>decl = 23.45*sin(corr*((a$DiaJuliano-80)*360/365)); decl</div><div>hn = 1/corr*acos(-tan(corr*lat)*tan(corr*decl))</div><div>a$N = 2*hn/15</div><div>s = (1+0.033*cos(a$DiaJuliano*360/365))</div><div>t = corr*hn*sin(corr*lat)*sin(corr*decl)</div><div>u = cos(corr*lat)*cos(corr*decl)*sin(corr*hn)</div><div>a$Qo = 37.6*s*(t+u)</div><div>head(a)</div><div>x = 0.25</div><div>b = 0.50</div><div>n.est = (a$RS/a$Qo-x)*a$N/b; n.est</div><div>a$n<- ifelse((n.est)<0,1, n.est); head(a)</div><div>a$n.N = a$n/a$N; head(a)</div><div>r = 0.23</div><div><br></div><div>a$Tmed = (a$Tmax+a$Tmin)/2; head(a)</div><div>a$esMax = 0.6100*exp((17.3*a$Tmax)/(237.3+a$Tmax))</div><div>a$esMin = 0.6100*exp((17.3*a$Tmin)/(237.3+a$Tmin))</div><div>a$esMed = (a$esMax+a$esMin)/2</div><div>a$ea = (a$esMed*a$UR)/100</div><div>a$DPV = a$esMed-a$ea</div><div>a$Delta = (4098*a$esMed)/(a$Tmed+237.3)**2</div><div>a$BOC = a$RS*(1-r)</div><div>a$BOL = -(0.903*(10**-9)*((a$Tmed+273)**4)*(0.34-0.14*(sqrt(a$ea)))*(0.1+(0.9*a$n.N)))</div><div>a$Rn = (a$BOC+a$BOL)</div><div>a$ETP.Penman = ((0.408*a$Delta*a$Rn)+(gama*(900/(a$Tmed+273))*a$Vel*a$DPV))/(a$Delta+(gama*(1+0.34*a$u2)))<br></div><div>#View(a)</div><div><br></div><div>a$W = ifelse(a$Tmed<16,(0.407+0.0145*a$Tmed),(0.483+0.01*a$Tmed))</div><div>a$ETP.Priestley = (1.26*a$W*a$Rn)/2.45</div><div><br></div><div>plot(a$ETP.Penman,a$ETP.Priestley)</div><div><br></div><div>write.table(a, "ETP.txt") </div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div>Yury Duarte<br></div>Engenheiro Agrônomo - ESALQ/USP<br></div></div></div>
<br><div class="gmail_quote">Em 28 de abril de 2016 14:51, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Yuri, boa tarde!<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Você deve editar as constantes conforme os dados disponíveis e fórmula utilizada.<br></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br><b>Penman-Monteith diário utilizando "Rs":</b><br><span style="font-family:monospace,monospace"># ET.PenmanMonteith(data, constants, ts="daily", solar="data", wind="yes", crop="short", ...)<br>labels <- c("Elev", "lambda", "lat_rad", "Gsc", "z", "sigma", "G")</span><br><br><b>Penman-Monteith diário utilizando "n":</b><br><span style="font-family:monospace,monospace"># ET.PenmanMonteith(data, constants, ts="daily", solar="sunshine hours", wind="yes", crop="short", ...)<br>labels <- c("Elev", "lambda", "lat_rad", "Gsc", "z", "sigma", "G", "as", "bs")</span><br><br><b>Priestley-Taylor diário utilizando "Rs":</b><br><span style="font-family:monospace,monospace"># ET.PriestleyTaylor(data, constants, ts="daily", solar="data", alpha=0.23, ...)<br>labels <- c("Elev", "lambda", "lat_rad", "Gsc", "sigma", "G", "alphaPT")</span><br><br><b>Priestley-Taylor diário utilizando "n":</b><br><span style="font-family:monospace,monospace"># ET.PriestleyTaylor(data, constants, ts="daily", solar="sunshine hours", alpha=0.23, ...)<br>labels <- c("Elev", "lambda", "lat_rad", "Gsc", "sigma", "G", "alphaPT", "as", "bs")</span><br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">No caso de comparar resultados atente para diferenças de padrão das fórmulas para crop="short" (FAO-56) ou crop="tall" (ASCE-EWRI).<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Você pode avaliar diretamente as fórmulas empregadas no pacote! Se você digitar a função sem estar seguida de parênteses, terá acesso ao código utilizado e poderá comparar com teu script:<br><br><span style="font-family:monospace,monospace">> ET.PenmanMonteith<br># function (data, constants, ts = "daily", solar = "sunshine hours", <br>#           wind = "yes", crop = "short", ...) <br># {<br>#   if (is.null(data$Tmax) | is.null(data$Tmin)) {<br>#     stop("Required data missing for 'Tmax.daily' and 'Tmin.daily', or 'Temp.subdaily'")<br>#   }<br>#   if (is.null(data$RHmax) | is.null(data$RHmin)) {<br>#     stop("Required data missing for 'RHmax.daily' and 'RHmin.daily', or 'RH.subdaily'")<br>#   }<br># ...</span><br></div><div class="gmail_extra"><span class=""><br clear="all"><div><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font style="font-size:small" face="arial, helvetica, sans-serif"><div style="font-family:'trebuchet ms',sans-serif;display:inline">​</div></font><font style="font-size:12.8px" face="arial, helvetica, sans-serif"><span style="font-size:small">================================================</span><br>Éder Comunello</font><div style="font-size:12.8px"><span style="font-family:arial,helvetica,sans-serif">Agronomist (UEM), MSc in Environ. Sciences (UEM)</span><font face="arial, helvetica, sans-serif"><br></font></div><div style="font-size:12.8px"><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">DSc in Agricultural Systems Engineering (USP/Esalq)</span><br></div><div><span style="font-size:small">Brazilian Agricultural Research Corporation (</span><font style="font-size:small" face="arial, helvetica, sans-serif">Embrapa)</font></div><div><div><font style="font-size:small" face="arial, helvetica, sans-serif">Dourados, MS, Brazil </font><font face="arial, helvetica, sans-serif" size="2">|<O>|</font></div><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">================================================</span><br></div></div></div><div style="font-size:small"><font face="arial, helvetica, sans-serif">GEO, </font>-22.2752, -54.8182, 408m<span style="font-family:arial,helvetica,sans-serif"><br></span></div><div style="font-size:12.8px"><span style="font-size:small;font-family:arial,helvetica,sans-serif">UTC-04:00 / DST: UTC-03:00</span></div><div><div><div><br></div><div><br></div></div><div style="font-size:small"><br></div></div></div></div></div></div></div></div></div></div></div>
<br></span><div><div class="h5"><div class="gmail_quote">Em 28 de abril de 2016 13:19, Yury Duarte <span dir="ltr"><<a href="mailto:yurynepomuceno@gmail.com" target="_blank">yurynepomuceno@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Obrigado, Éder!!<div><br></div><div>De fato, não consegui entender a necessidade de informar todas essas constantes que vc abordou no código quando li o "help" do pacote.<div>Irei rodar seus comandos e acertarei as constantes segundo minhas situações.</div><div>Estou finalizando a programação das fórmulas de PriestleyTaylor e PenmanMonteith para comparar com as saídas das evapo fornecidas pelo pacote. Assim que tiver os scripts prontos os enviarei nessa mesma conversa!</div></div><div><br></div><div>Abs</div></div><div class="gmail_extra"><span><br clear="all"><div><div><div dir="ltr"><div>Yury Duarte<br></div>Engenheiro Agrônomo - ESALQ/USP<br></div></div></div>
<br></span><div><div><div class="gmail_quote">Em 28 de abril de 2016 12:18, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div style="font-family:verdana,sans-serif">Yuri, bom dia!<br><br></div><div style="font-family:verdana,sans-serif">Esse pacote é um pouco melindroso quanto ao formato de entrada dos dados e muita coisa não está documentada.<br><br></div><div style="font-family:verdana,sans-serif">O objeto de entrada deve ser uma lista, iniciando com Date e J (data e dia juliano), sendo seguida pelas variáveis climáticas (na classe {zoo}).<br><br></div><div style="font-family:verdana,sans-serif">O script tá rodando a partir dos seus dados de exemplo, mas você tem que utilizar as constantes adequadas pro seu caso.<br></div><div style="font-family:verdana,sans-serif"><br><span style="font-family:monospace,monospace">### <code r><br>rm(list=ls())<br>require(Evapotranspiration)<br><br># URL <- "<a href="http://r-br.2285057.n4.nabble.com/attachment/4666053/0/a.txt" target="_blank">http://r-br.2285057.n4.nabble.com/attachment/4666053/0/a.txt</a>"<br><br>setwd("~/LAB/LEARN") ### Alterar!<br>df <- read.table("a.txt", sep=";", head=T, <a href="http://as.is" target="_blank">as.is</a>=T)<br>df$Data <- as.Date(df$Data)<br><br>climate <- lapply(as.list(df)[2:7], zoo, df$Data)<br>J       <- as.numeric(format(df$Data, "%j"))<br>data    <- c(list(Date.daily=df$Data, J=J), climate)<br><br>names(data)<br># [1] "Date.daily" "J"          "Rs"         "Tmax"       "Tmin"       "RHmax"      "u2"         "RHmin"     <br>str(data)<br><br># Editar adequadamente as constantes utilizadas!!!<br># As constantes a serem utilizadas variam com a formulação escolhida... Ver ajuda...<br># ?ET.PenmanMonteith<br># ?ET.PristleyTaylor<br># pi/180*-23.45 # [1] -0.4092797<br>myConst <- list(lambda = 2.45, sigma = 4.903e-09, Gsc = 0.082, <br>                          lat = -23.45, lat_rad = -0.40928, as = 0.25, <br>                          bs = 0.55, Elev = 480, z = 2, Roua = 1.2, Ca = 0.001013, G = 0, <br>                          alphaA = 0.14, alphaPT = 1.26, ap = 2.4, fz = 28, b0 = 1, <br>                          a_0 = 11.9, b_0 = -0.15, c_0 = -0.25, d_0 = -0.0107, e0 = 0.81917, <br>                          e1 = -0.0040922, e2 = 1.0705, e3 = 0.065649, e4 = -0.0059684, <br>                          e5 = -0.0005967, gammaps = 0.66, epsilonMo = 0.92, PA = 285.8, <br>                          alphaMo = 17.27, betaMo = 237.3, sigmaMo = 5.67e-08, lambdaMo = 28.5, <br>                          b1 = 14, b2 = 1.2)<br><br>res1 <- ET.PenmanMonteith(data,  myConst, ts="daily", solar="data", wind="yes", crop = "short")<br>res2 <- ET.PriestleyTaylor(data, myConst, ts="daily", solar="data", alpha=.23)<br><br>res <- cbind(PM=res1$ET.Daily, PT=res2$ET.Daily)<br>head(res)<br>fit <- lm(res$PM~res$PT-1); summary(fit)<br>plot(res$PM~res$PT); abline(fit, col=2); abline(a=0, b=1, col=3, lty=2) # 1:1<br>### </code></span><br></div></div><div class="gmail_extra"><br clear="all"><div><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font style="font-size:small" face="arial, helvetica, sans-serif"><div style="font-family:'trebuchet ms',sans-serif;display:inline">​</div></font><font style="font-size:12.8px" face="arial, helvetica, sans-serif"><span style="font-size:small">================================================</span><span><font color="#888888"><br>Éder Comunello</font></span></font><span><font color="#888888"><div style="font-size:12.8px"><span style="font-family:arial,helvetica,sans-serif">Agronomist (UEM), MSc in Environ. Sciences (UEM)</span><font face="arial, helvetica, sans-serif"><br></font></div><div style="font-size:12.8px"><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">DSc in Agricultural Systems Engineering (USP/Esalq)</span><br></div><div><span style="font-size:small">Brazilian Agricultural Research Corporation (</span><font style="font-size:small" face="arial, helvetica, sans-serif">Embrapa)</font></div><div><div><font style="font-size:small" face="arial, helvetica, sans-serif">Dourados, MS, Brazil </font><font face="arial, helvetica, sans-serif" size="2">|<O>|</font></div><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">================================================</span><br></div></div></div><div style="font-size:small"><font face="arial, helvetica, sans-serif">GEO, </font>-22.2752, -54.8182, 408m<span style="font-family:arial,helvetica,sans-serif"><br></span></div><div style="font-size:12.8px"><span style="font-size:small;font-family:arial,helvetica,sans-serif">UTC-04:00 / DST: UTC-03:00</span></div><div><div><div><br></div><div><br></div></div><div style="font-size:small"><br></div></div></font></span></div></div></div></div></div></div></div></div></div><div><div>
<br><div class="gmail_quote">Em 27 de abril de 2016 18:09, Jônatan <span dir="ltr"><<a href="mailto:jdtatsch@gmail.com" target="_blank">jdtatsch@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Qual erro? Qual os parâmetros usados no comando? <div>Exponha um exemplo reproduzível.</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>2016-04-27 15:02 GMT-03:00 Yury Duarte <span dir="ltr"><<a href="mailto:yurynepomuceno@gmail.com" target="_blank">yurynepomuceno@gmail.com</a>></span>:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><span style="font-size:12.8px">Boa tarde colegas programadores!</span><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Estou com algumas dificuldades em rodar o comando "ET.PenmanMonteith" do pacote "Evapotranspiration". Minha intenção é de calcular a evapotranspiração potencial para a escala diária. Segundo o pacote, as variáveis exigidas são:</div><div style="font-size:12.8px">Tmax; Tmin; RHmax; RHmin; Rs e u2.</div><div style="font-size:12.8px">Todas as variáveis exigidas pelo pacote estão compreendidas no meu arquivo de entrada (que segue em anexo), entretanto o comando segue dando erro.</div><div style="font-size:12.8px">Caso alguém seja familiar com o pacote e/ou com o comando e puder me dar uma luz, seria de muita ajuda!</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Desde já agradeço pela colaboração de todos!</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Abraços</div><span><font color="#888888"><div><div><div dir="ltr"><div>Yury Duarte<br></div>Engenheiro Agrônomo - ESALQ/USP<br></div></div></div>
</font></span></div>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<span><font color="#888888"><br></font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><div><div dir="ltr"><div><font style="font-family:arial" size="2">###############################################################</font><br><span style="font-family:arial;font-size:small">##  </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Jônatan Dupont Tatsch</span></div><div><span style="font-family:arial;font-size:small">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Professor do Departamento de Física</span><font style="font-size:small" face="arial"><br></font></div><div><span style="font-family:arial;font-size:small">## </span><span style="font-family:arial;font-size:small"> </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.9090909957886px;line-height:12.9971590042114px;text-align:center;background-color:rgb(239,239,239)">Coordenador Substituto do Programa de Pós-Graduação em Meteorologia (PPGMET)</span><br></div><div><span style="font-size:small;font-family:arial">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Centro de Ciências Exatas e Naturais (CCNE)</span><br></div><div><span style="font-size:small;font-family:arial">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Universidade Federal de Santa Maria - UFSM</span></div><div><span style="font-family:arial;font-size:small">##  </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Faixa de Camobi, Prédio 13 - Campus UFSM - </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Santa Maria, RS, Brasil - 97105-900</span><br></div><div style="font-family:arial;font-size:small">##  <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Telefone: <a href="tel:%2B55%2855%2933012083" value="+555533012083" target="_blank">+55(55)33012083</a></span><br></div><div style="font-family:arial;font-size:small">##  <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px"><a href="http://www.ufsm.br/meteorologia" target="_blank">www.ufsm.br/meteorologia</a></span></div><div style="font-family:arial;font-size:small">###############################################################</div></div></div></div></div>
</font></span></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div></div>
</blockquote></div><br></div></div></div></div>
</blockquote></div><br></div>