I found a detailed page with examples of how to compute solar positions and twilight times.
Bob Miller (Michigan State) gave me a computer code in the early 2000s that tabulated moon poisition and phase. I used it to compute these data until ??? For now, here is an example of how I use it:
[sco@mcs ~]$ utdatemoon.sh Usage: utdatemoon.sh 20190515 arg1 - UT date (YYYYMMDD) arg2 - run in debug/verbose mode [sco@mcs ~]$ utdatemoon.sh 20190515 Y Debug flag is ON in utdatemoon.sh utdate = 20190515 Enter any key to continue: 81.500 12:19:00.90 +03:16:39.0 20190515 (illum,ra,dec,UTdate)