#' @export
Geracao_de_Energia_Estocastico5<-function(long,long.st,lat,A,
DadosAmbientais,Pdc0,Beta,gama,
TipoDePainel,pot_inv,n_inv,n_non,perdasTot,G.sc=1367,TempRef=25,IrrRef=1000,n_ref=0.9637,
tempo.inicial="01/01/2015 00:00", tempo.final="31/12/2015 23:00", passo=60, Hora_do_passo="Mais_tarde",
AlbedoFixo=0.2,TipoAlbedo="fixo",VVFixo=1.7,TipoVV="fixo",TaFixo=25,TipoTA="fixo",modo_radiacao="ceu_limpo", ModeloSintese="Estocastico",ModeloComposicao="Orgill",
TipoDeClima,TipoModulo_Montagem,DadosMensais, Precos){
g<-geometrico(lat, long, long.st, Beta, gama, passo, Hora_do_passo, tempo.inicial, tempo.final) #saida ok
r<-Rad.Extraterrestre(g$tempos, g$theta.z, g$HaSol, G.sc) #saida ok
c<-CeuLimpo(A,TipoDeClima,g$theta,g$theta.z,g$HaSol, r$G.on, r$G.o)
s<-Sintetico(g$omega.sunset,g$omega.t,g$omega.b_delta.t,
DadosMensais, r$G.on, g$delta, g$omega.1, g$omega.2,
lat, g$HaSol, g$tempos, ModeloComposicao, ModeloSintese, g$omega.t_delta.t, g$theta.z) # Orgill Estocastico semi-validado
T.Sintetico<-TemperaturaSintetico(g$tempos,DadosMensais)
# Possivelmente alterar
a<-Ambientais(modo_radiacao, c$G.cnb, c$G.cd, g$theta.z, r$G.o, #K.T tem NaNs. To Do: Corrigir.
DadosAmbientais, g$HaSol,
AlbedoFixo, TipoAlbedo, VVFixo, TipoVV, TaFixo, TipoTA,
s$Ibn, s$Id, T.Sintetico, tempos, DadosMensais)
pl<-Pl.Inc_DirRef(g$theta,g$theta.z,g$HaSol,a$ro.g,Beta, a$G.b, a$G)
hdkr<-HDKR(g$theta,g$theta.z,a$G.nb,a$G.b,a$G.d,a$G, Beta, pl$R.b,r$G.on,r$G.o) # to do #' @export
liu<-Isotropico(a$G.d,Beta)
i<-Perez(a$G.d, a$G.nb, g$theta, g$theta.z, r$G.on, Beta)
G.t<-i$G.Td+pl$G.T.b+pl$G.T.refl
#write.table(G.t, file="perez_R.csv", append="FALSE",sep=";", dec=",")
p<-CaracteristicPainel(TipoDePainel)
e<-RadiacaoEfetiva(p$eta.Ar,p$eta.VidroAR, p$eta.Vidro,g$theta,pl$G.T.b,i$G.Td,pl$G.T.refl)
m<-ModelagemTermica(TipoModulo_Montagem)
###Pac!!?!
ce<-ConversaoEnergetica(perdasTot,passo,TempRef,p$TempCoef,n_inv*pot_inv,Pdc0,IrrRef,e$G.ef.T,a$TA, m$a, m$b, m$DeltaT, a$VV ,n_non/n_ref, n_non, g$Delta.T)
Valor<-ValorEnergetico(ce$E, Precos)
Precos<-as.matrix(Precos)
# saida<-list("G.t"=G.t,"E"=E, "s"=s, "a"=a, "pl"=pl, "i"=i,"e"=e, "T.Sin"=T.Sintetico)
saida2<-data.frame("tempos"=g$tempos,
"HaSol"=g$HaSol,
"theta"=g$theta,
"G.on"=r$G.on,
"G.nb"=a$G.nb,
"G.d"=a$G.d,
"G"=a$G,
"G.t"=G.t,
"Ta"=a$TA,
"WS"=a$VV,
"G.ef"=e$G.ef.T,
"Tm"=ce$Tm,
"E"=ce$E,
"liu"=liu,
"hay"=hdkr,
"perez"=i,
"E.dc"=ce$E_dc)
detach(f)
return(saida2)
}
#' @export
head.list <- function(obj, n = 6L, ...)
{
stopifnot(length(n) == 1L)
origN <- n
n <- if (n < 0L)
max(length(obj) + n, 0L)
else min(n, length(obj))
lapply(obj[seq_len(n)], function(x)
{
tryCatch({
head(x, origN, ...)
}, error = function(e) {
x
})
})
}
environment(head.list) <- asNamespace('utils')
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.