I’m trying to create a python application to automatically take pictures during the next solar eclipse. One of the things I still need to implement is to find the times of first contact (C1), start of totality (C2), maximun eclipse (M), end of totality (C3) and end of eclipse (C4).
Is there an easy way in sunpy to calculate these times accurately?
Thanks a lot!
I guess the example at Obtaining solar-eclipse information — sunpy 5.1.1 documentation
shows everything I need. I guess the maximum eclipse is exactly in the middle of the start and end of the total phase.
Glad you found that example!
I’ll just warn you that I have checked sunpy’s calculations against what’s computed by various eclipse websites, and the contact times can be discrepant by up to a few seconds. I don’t know who is more accurate, and of course the true contact times can be tricky anyway due to the irregularity of the lunar limb profile. Please do account for some uncertainty when using sunpy’s output to set up your automatic picture taking.
I was indeed interested in the accuracies of the timings. I also found them to be a few seconds off, and I don’t know why that is the case.