# install.packages("suntools")
# install.packages("lubridate")

library(suntools)
library(lubridate)


lat <- 44.7912
lon <- -0.6078
tz  <- "Europe/Paris"
d_deb <- "2026-01-05"
d_fin <- "2026-01-10"
date_debut <- as.POSIXct(d_deb, tz = tz)
date_fin   <- as.POSIXct(d_fin, tz = tz)

sequence_temps <- seq(from = date_debut, to = date_fin, by = "15 min")

unix_time <- as.numeric(sequence_temps)

coords <- matrix(c(lon, lat), nrow = 1)

positions <- solarpos(coords, sequence_temps)

df_soleil <- data.frame(
  timestamp = sequence_temps,
  utime = unix_time,
  azimut = positions[, 1],
  elevation = positions[, 2]
)

csv_name <-c("solar_data_", d_deb ,"_to_", d_fin ,".csv")
csv_name <-paste(csv_name, collapse = "")
print(csv_name)

store_path  <- "../data/solar_pos_data/"
csv_path_name <- paste(c(store_path, csv_name), collapse = "")
print(csv_path_name)

write.csv(df_soleil, csv_path_name, row.names = FALSE)

print(paste(csv_name,"succesfully generated"))