Day-light terminators and other sunlight parameters
coast(limits=:global, frame=(annot=:a,ticks=:a,grid=:a), proj=:EckertVI,
area=5000, shore=0.5, borders=(type=1, pen=(0.5,:gray)),
water=[175 210 255], par=(:MAP_FRAME_TYPE, :plain), figsize=16)
solar!(terminators=(term=:d, date="2016-02-09T16:00:00"), fill="navy@95")
solar!(terminators=(term=:c, date="2016-02-09T16:00:00"), fill="navy@85")
solar!(terminators=(term=:n, date="2016-02-09T16:00:00"), fill="navy@80")
solar!(terminators=(term=:a, date="2016-02-09T16:00:00"), fill="navy@80")
# Get Sun position at at that date and time an plot it.
t = solar(sun=(date="2016-02-09T16:00:00",), formated=true)
plot!(t, symbol="k@sunglasses/1.5", fill=:yellow, show=true)
See also GMT ex46