lat=37.2;
G0dm=c(2766, 3491, 4494, 5912, 6989, 7742, 7919, 7027, 5369, 3562, 2814,
2179)
Ta=c(10, 14.1, 15.6, 17.2, 19.3, 21.2, 28.4, 29.9, 24.3, 18.2, 17.2, 15.2)
prom=list(G0dm=G0dm, Ta=Ta)
###Different tracker methods
ProdFixed<-prodGCPV(lat=lat,dataRad=prom, keep.night=FALSE)
Prod2x<-prodGCPV(lat=lat, dataRad=prom, modeTrk='two', keep.night=FALSE)
ProdHoriz<-prodGCPV(lat=lat,dataRad=prom, modeTrk='horiz', keep.night=FALSE)
prod <- mergesolaR(ProdFixed, Prod2x, ProdHoriz)
head(prod)
## Due to changes in SIAR webpage this code no longer works
EstMadrid <- subset(RedEstaciones, NomProv=='Madrid')
nEstMadrid <- nrow(EstMadrid)
namesMadrid <- EstMadrid$NomEst
prodMadrid <- lapply(1:nEstMadrid,
function(x){try(prodGCPV(lat=41, modeRad='siar',
dataRad=list(prov=28, est=x,
start='01/01/2009', end='31/12/2010'))
)})
names(prodMadrid) <- namesMadrid
okMadrid <- lapply(prodMadrid, class)!='try-error'
prodMadrid <- prodMadrid[okMadrid]
YfMadrid <- do.call(mergesolaR, prodMadrid)
horizonplot(YfMadrid-rowMeans(YfMadrid),
origin=0,
scales=list(y=list(relation='same')),
colorkey=TRUE)
TargetDiagram(YfMadrid, end=as.POSIXct('2010-12-31'), ndays=c(10, 20,
30, 40, 50, 60), cex=0.6)
Run the code above in your browser using DataLab