Daily Path of the Sun through the Year 

Plot the position of the Sun for your location between sunrise and sunset for the solstices and equinox.

show complete Wolfram Language inputhide input

In[1]:=

X summersolstice = DateObject[{2013, 6, 20}]; autumnalequinox = DateObject[{2013, 9, 22}]; wintersolstice = DateObject[{2013, 12, 21}]; coordsToCartesian = Evaluate[CoordinateTransformData["Spherical" -> "Cartesian", "Mapping", {1, \[Pi]/2 - (#2 Degree), 2 Pi - (#1 Degree)}]] &; getValues = Function[{series}, Map[QuantityMagnitude, series["Values"], {2}]]; getSunriseToSunset = DateRange[Sunrise[#], Sunset[#], {30, "Minute"}] &; latitude = QuantityMagnitude[Latitude[$GeoLocation], "Radians"]; Graphics3D[{ Black, Table[ Text[{"W", "S", "E", "N"}[[j]], 1.1 {Cos[Pi/2 j], Sin[Pi/2 j], 0}], {j, 4}], Rotate[{ Arrow[{{0, 0, 0}, {0, 0, 1.4}}], Arrowheads[.02], Arrow@Table[{.1 Cos[t], .1 Sin[t], 1}, {t, 0, -3 Pi/2, -.2}]}, latitude, {0, 1, 0}], , Thickness[.007], Transpose[{{RGBColor[.49, 0, 0], RGBColor[.88, .63, .23], RGBColor[.27, .2, .75]}, (Line@(coordsToCartesian @@@ (getValues@ SunPosition[ getSunriseToSunset[#]]))) & /@ {summersolstice, autumnalequinox, wintersolstice}}], RGBColor[.6, .73, .36], Cylinder[{{0, 0, -.005}, {0, 0, 0}}], Opacity[.1], Sphere[]}, PlotRange -> {0, 1.4}, Boxed -> False, SphericalRegion -> True, ViewAngle -> \[Pi]/12, ViewCenter -> {.5, .5, .2}]

Out[1]=