Example #1
0
func GetSunEclipticLatitudeForEarth(jd float64) float64 {
	t := calendarutil.GetJulianThousandYears(jd)
	B0 := GetEarthB0(t)
	B1 := GetEarthB1(t)
	B2 := GetEarthB2(t)
	B3 := GetEarthB3(t)
	B4 := GetEarthB4(t)
	B := ((((B4*t)+B3)*t+B2)*t+B1)*t + B0
	return B
}
Example #2
0
/**
 * 按照儒略日计算地球和太阳的距离
 *
 * @param jd
 *            儒略日
 * @return 地球和太阳的距离,单位是天文单位(au)
 */
func GetSunRadiusForEarth(jd float64) float64 {
	t := calendarutil.GetJulianThousandYears(jd)
	R0 := GetEarthR0(t)
	R1 := GetEarthR1(t)
	R2 := GetEarthR2(t)
	R3 := GetEarthR3(t)
	R4 := GetEarthR4(t)
	R5 := GetEarthR5(t)
	R := ((((R5*t+R4)*t+R3)*t+R2)*t+R1)*t + R0
	return R
}
Example #3
0
func GetSunEclipticLongitudeForEarth(jd float64) float64 {
	t := calendarutil.GetJulianThousandYears(jd)
	L0 := GetEarthL0(t)
	L1 := GetEarthL1(t)
	L2 := GetEarthL2(t)
	L3 := GetEarthL3(t)
	L4 := GetEarthL4(t)
	L5 := GetEarthL5(t)
	L := ((((L5*t+L4)*t+L3)*t+L2)*t+L1)*t + L0
	return mathutil.Mod2Pi(L)
}