Hi all,
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?
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.