$h_0=\arcsin[\sin\phi\sin\sigma+\cos\phi\cos\sigma\cos(15t+\lambda-300+E_Q/4)]$
参数说明:
$\theta=\frac{2 \pi \times (N+△N-N_0)}{365.2422}$
$\sigma=0.3723+23.2567 \sin \theta+0.1149 \sin 2\theta−0.1712 \sin 3\theta-0.7580 \cos \theta+0.3656 \cos2\theta+0.0201 \cos3\theta$