How complicated is it to calculate a surface normal on the spherical approximation of the earths surface pointing towards the sun at a given point in time?
What I try do is to highlight a small area on a world map where the radiation from the sun tangential hits a imaginary sphere around the earth.
This is just to get an idea where this area was at the given time, so this information has not to be very accurate. As a mater of fact I will just incorporated this information in the map if it is easy enough to calculate.