# NOT RUN {
#Calculating solar hour angle for two consecutive days
DOY <- seq(0, 2, .05)
h <- HourAngle(DOY, Lon=0, SLon=0, DS=60)
#Note: only the difference between Lon and SLon matters not each value
plot(DOY, h)
# }
Run the code above in your browser using DataLab