R/calcShd.R

In solaR: Solar Photovoltaic Systems

 # Copyright (C) 2010 Oscar Perpiñán Lamigueiro
 #
 # This program is free software; you can redistribute it and/or
 # modify it under the terms of the GNU General Public License
 # as published by the Free Software Foundation; either version 2
 # of the License, or (at your option) any later version.
 #
 # This program is distributed in the hope that it will be useful,
 # but WITHOUT ANY WARRANTY; without even the implied warranty of
 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 # GNU General Public License for more details.
 #
 # You should have received a copy of the GNU General Public License
 # along with this program; if not, write to the Free Software
 # Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
 #/
calcShd<-function(radEf,##class='Gef'
                  modeTrk='fixed',     #c('two','horiz','fixed')
                  modeShd='prom',      #modeShd=c('area','bt','prom')
                  struct=list(), #list(W=23.11, L=9.8, Nrow=2, Ncol=8), 
                  distances=data.frame() #data.frame(Lew=40, Lns=30, H=0)){
                  ){
  stopifnot(is.list(struct), is.data.frame(distances))

  ##Por ahora sólo utilizo modeShd='area'
  ##Con diferentes modeShd (por definir) podré calcular Gef de diferente forma
  ##Ver tesis de Macagnan
  prom=("prom"  %in%  modeShd)
  prev <- as.zooI(radEf, complete=TRUE)
###Cálculo de sombras
  sol <- prev[,c('AzS', 'AlS')]
  theta <- prev[,c('Beta', 'Alfa', 'cosTheta')]
  AngGen <- CBIND(sol, theta, index=indexI(radEf))
  Shd<-fSombra(AngGen, distances, struct, modeTrk, prom)
  FS=coredata(Shd)
###Cálculo de irradiancia
  ## gef0 <- prev[,c('G', 'D', 'B', 'R',
  ##              'FTb', 'FTd', 'FTr',
  ##              'Dief', 'Dcef',
  ##              'Gef', 'Def', 'Bef', 'Ref')]
  gef0 <- radEf@GefI
  Bef0=gef0$Bef
  Dcef0=gef0$Dcef
  Gef0=gef0$Gef
  Dief0=gef0$Dief
  Ref0=gef0$Ref
  ##Cálculos
  Bef=Bef0*(1-FS)
  Dcef=Dcef0*(1-FS)
  Def=Dief0+Dcef
  Gef=Dief0+Ref0+Bef+Dcef               #Incluyendo sombras
  ##Cambio nombres
  nms=c('Gef', 'Def', 'Dcef', 'Bef')
  nmsIndex <- which(names(gef0) %in% nms)
  names(gef0)[nmsIndex]<- paste(names(gef0)[nmsIndex], '0', sep='')
  ##Entrego zoo con resultados, incluyendo previos sin sombras
  GefShd=CBIND(gef0, data.frame(Gef, Def, Bef, Dcef, FS), index=indexI(radEf))


###Valores diarios, mensuales y anuales
  DayOfMonth=c(31,28,31,30,31,30,31,31,30,31,30,31) ###OJO
  
  if (radEf@type=='prom') {
    Gefdm=aggregate(GefShd[,c('Gef0', 'Def0', 'Bef0', 'G', 'D', 'B', 'Gef', 'Def', 'Bef')]/1000,
      by=as.yearmon, FUN=P2E, radEf@sample)       #kWh
    names(Gefdm)=paste(names(Gefdm), 'd', sep='')

    GefD=Gefdm*1000                  #Wh
    index(GefD) <- indexD(radEf)     ##para que sea compatible con G0D
    
    Gefy=zoo(t(colSums(Gefdm*DayOfMonth)),
      unique(year(index(Gefdm))))
  } else {
    GefD=aggregate(GefShd[,c('Gef0', 'Def0', 'Bef0', 'G', 'D', 'B', 'Gef', 'Def', 'Bef')],
      by=truncDay, FUN=P2E, radEf@sample)         #Wh
    names(GefD)=paste(names(GefD), 'd', sep='')

    Gefdm=aggregate(GefD/1000, by=as.yearmon, mean, na.rm=1)
    Gefy=aggregate(GefD/1000, by=year, sum, na.rm=1)
  }

###Entrego un objecto de clase Gef
###modificando los slots 'modeShd', 'GefI', 'GefD', 'Gefdm', y 'Gefy'
###del objecto original radEf
  radEf@modeShd=modeShd
  radEf@GefI=GefShd
  radEf@GefD=GefD
  radEf@Gefdm=Gefdm
  radEf@Gefy=Gefy
  return(radEf)
}