powered by
This routine uses Forsythe et al. 1995.
daylength(doy, latitude)
nested list with daylength (daylength) and solar elevation (solar_elev) elements
a vector with doy values 1 - 365(6)
a given latitude
if (FALSE) { # calcualte the hours of sunlight and solar elevation on day of year 1 # and latitude 51 ephem <- daylength(1, 51) print(ephem) }
Run the code above in your browser using DataLab