Need function which returns day/night longitude for a given latitude. May be it is possible to get even 2 values: night and shadow border?