This routine uses Forsythe et al. 1995.
daylength(doy, latitude)nested list with daylength (daylength) and solar elevation (solar_elev) elements
if (FALSE) { # \dontrun{
# calcualte the hours of sunlight and solar elevation on day of year 1
# and latitude 51
ephem <- daylength(1, 51)
print(ephem)
} # }