daylength function

Compute sunrise and sunset times, and daylength