Open ClaudiaGEscribano opened 5 years ago
He calculado productividad anual a patir de los datos diarios obtenidos con el código anterior y creo que los resultados no están bien, me parece que son valores bajos.
Por actualizar y dar algún dato más,
Haciendo los cálculos para una celda, la 200:
cell <- 200
x <- as.vector(SISS[200])
xx <- as.vector(Tas[200])
lat <- y[200]
data <- c(lat, x, xx)
Usando datos intradiarios obtengo una productividad anual de ~ 805 kWh/kWp para el (año 2000) y de 1500 si en lugar de datos cada 3 horas uso datos de diarios.
Este es el código para hacer un ejemplo del cálculo de productividad diaria a partir de datos cada 3 horas. Creo que los resultados tienen sentido. Quizá lo único un poco raro es si el índice temporal cada tres horas se interpreta bien por ProdCGPV en la forma en la que está.
Hay que ejecutar ProdparaleloI.R que llama a fooIprod.R y usa los datos en data.