Using skyfield, the astronomical python module, I want to find the latitude and longitude for the sun and the moon subpoints, that is, the location on Earth where the body is directly overhead. I need to be able to compute them for any time of day and any day of the year. I've looked through the examples the documentation provides but I didn't find these specific calculations nor any mention of them in the text.
Is there a way of doing such a computation?