library('pacman')
p_load('ggplot2')
p_load('solaR')
p_load('stringr')
p_load('lubridate')

lat=52.23

sol = readRDS(file="sol.Rda")

sol$timestamp = as.Date(sol$timestamp)

class(sol$timestamp)
x = dfI2Meteo(sol, time.col = 'timestamp',lat=lat,format='%Y-%m-%d %H:%M:%S')
head(x)
g = calcGef(lat,
            modeTrk='fixed',
            modeRad='bdI',
            dfI2Meteo(sol, time.col = 'timestamp',lat=lat,format='%Y-%m-%d %H:%M:%S'),
            keep.night=TRUE,
            sunGeometry='michalsky',
            beta=90,
            alfa=36.29,
            corr='none',
            iS=2,
            alb=0.0)
   
head(g)

#g = as.data.frameI(g)[c('Gef')]

#head(g$timestamp)
#rownames(g)