C# 计算太阳倾角、太阳高度角

	/// 计算太阳高度角ho
        /// </summary>
        /// <param name="φ">当地纬度，度(如34.64)</param>
        /// <param name="t">北京时间，时(如9)</param>
        /// <param name="δ">太阳倾角，度(如-23.1)</param>
        /// <param name="λ">当地经度，度(如34.64)</param>
        /// <returns>返回太阳高度角，度</returns>
        public double ho(double φ, int t, double δ, double λ)
        {
            //公式：ho=arcsin[sinφsinδ+cosφcosδcos(15t+λ-300)]
            double a = 15 * t + λ - 300;//时角
            double b = Math.Cos(φ * Math.PI / 180d) * Math.Cos(δ * Math.PI / 180d) * Math.Cos(a * Math.PI / 180d);
            double c = Math.Sin(φ * Math.PI / 180) * Math.Sin(δ * Math.PI / 180);
            double d = c + b;
            double ho = Math.Asin(d) / (Math.PI /