计算给定经纬度的真太阳时怎么做呢？ about astro HOT 2 CLOSED

package main

import (
	"fmt"
	"github.com/starainrt/astro/basic"
	"github.com/starainrt/astro/sun"
	"math"
	"time"
)

func main() {
	cst := time.FixedZone("CST", 8*3600)
	date := time.Date(2022, 1, 1, 0, 0, 0, 0, cst)
	fmt.Println(date, "\n",
		ApparentSolarTime(date, 115), "\n",
		ApparentSolarTime2(date, 115))
}

// 方法一：真太阳时=当地平太阳时+均时差
func ApparentSolarTime(date time.Time, lon float64) time.Time {
	//取世界时，用于计算均时差
	jde := basic.Date2JDE(date.UTC())
	LocalTimezone := int(lon * 3600.00 / 15.0)
	//pTime:获取当地经度时区的平太阳时
	pTime := date.In(time.FixedZone("LTZ", LocalTimezone))
	//真太阳时=平太阳时+均时差
	return pTime.Add(time.Duration(basic.SunTime(jde) * 3600 * 1000 * 1000 * 1000))
}


// 方法二 真太阳时=太阳时角+12小时
func ApparentSolarTime2(date time.Time, lon float64) time.Time {
	//真太阳时=太阳时角+12小时
	trueTime := (sun.HourAngle(date, lon, 0) + 180) / 15
	if trueTime > 24 {
		trueTime -= 24
	}
	//真太阳时的分
	minute := (trueTime - math.Floor(trueTime)) * 60
	//真太阳时的秒
	second := (minute - math.Floor(minute)) * 60
	//当地经度下的本地时区
	trueSunTime := date.In(time.FixedZone("LTZ", int(lon*3600.00/15.0)))
	if trueSunTime.Hour()-int(trueTime) > 12 {
		trueSunTime = trueSunTime.Add(time.Hour * 24)
	} else if int(trueTime)-trueSunTime.Hour() > 12 {
		trueSunTime = trueSunTime.Add(-time.Hour * 24)
	}
	return time.Date(trueSunTime.Year(), trueSunTime.Month(), trueSunTime.Day(),
		int(trueTime), int(minute), int(second), int((second-math.Floor(second))*1000000000),
		time.FixedZone("LTZ", int(lon*3600.00/15.0)))
}