露点温度 海平面气压 场面气压 计算公式 类

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace LuDianOperation
{
    public class Ludian
    {
        //基于马哥拉斯改进公式计算水汽压与露点温度

        public double T { get; set; }
        public double F { get; set; }
        private double Pa { get; set; }
        private double a { get; set; }
        private double b { get; set; }

        public static readonly double E0 = 6.1078;

        //海平面气压
        // P0=Ps×10^[h/18400(1+Tm/273)]

        //P0为海平面气压，单位hPa；
        //Ps为本站气压，单位hPa；
        //h为气压传感器拔海高度，单位为m；
        //Tm为气柱平均温度，单位为℃，Tm=（t+t12)/2+h/400；
        //t为干球温度，单位℃；
        //t12为观测前12小时气温，单位℃。

        /// <summary>
        /// 获得海平面气压
        /// </summary>
        /// <param name="CurrentT">当前时间的温度</param>
        /// <param name="T12">12个小时之前的温度</param>
        /// <param name="CurrentStationPS">当前站点的气压</param>
        /// <param name="HaibaH">气压传感器拔海高度（海拔高度）</param>
        /// <returns>海平面气压（double）</returns>
        public double GetHaiPingMianQiYa(double CurrentT, double T12, double CurrentStationPS, double HaibaH)
        {
            try
            {
                //Tm=（t+t12)/2+h/400；
                double Tm = (CurrentT + T12) / 2 + HaibaH / 400;
                //[h/18400(1+Tm/273)]
                double centerDataOne = HaibaH / (18400 * (1 + Tm / 273));
                //10^[h/18400(1+Tm/273)]
                double centerDataTwo = Math.Pow(10.0, centerDataOne);
                //Ps×10^[h/18400(1+Tm/273)]
                return Math.Round(CurrentStationPS * centerDataTwo, 1);
            }
            catch (Exception ex)
            {
                Console.WriteLine(ex.Message);
                return 0.0;
            }
        }

        /// <summary>
        /// 无参数构造函数
        /// </summary>
        public Ludian()
        { 
            
        }

        /// <summary>
        /// 构造函数
        /// </summary>
        /// <param name="t">温度</param>
        /// <param name="f">湿度</param>
        public Ludian(double t, double f)
        {
            this.T = t;
            this.F = f;
            if (this.T > 0)
            {
                this.a = 7.69;
                this.b = 243.92;
            }
            else
            {
                this.a = 9.5;
                this.b = 265.5;
            }
        }

        #region 露点温度以及饱和水气压
        /// <summary>
        /// 获得饱和水蒸气压
        /// </summary>
        /// <returns></returns>
        public double GetES()
        {
            double Center = (a * this.T) / (b + this.T);
            return Math.Round(E0 * Math.Pow(10.0, Center), 1);
        }

        /// <summary>
        /// 计算中间变量
        /// </summary>
        /// <returns></returns>
        private double E()
        {
            return GetES() * F * 0.01;
        }

        /// <summary>
        /// 计算出露点温度
        /// </summary>
        /// <returns>double</returns>
        public double GetTd()
        {
            double ESS = E();
            double AAA = b * Math.Log10(ESS / E0);
            double BBB = a - Math.Log10(ESS / E0);
            return Math.Round(AAA / BBB, 1);
        }

        #endregion



        #region 水气压
        /// <summary>
        /// 根据水面温度获得水汽压
        /// </summary>
        /// <returns></returns>
        public double GetEShuiQIYA(double Ts)
        {
            double Two = 17.269 * ((Ts - 273.16) / (Ts - 35.86));
            return Math.Round(610.78 * Math.Exp(Two), 1);
        }

        /// <summary>
        /// 干湿球算法求水汽压
        /// </summary>
        /// <param name="FSS">相对湿度</param> 
        /// <returns></returns>
        public double GetEShuiQIYA()
        {
            return Math.Round((this.F * GetES() / 100), 2);
        }

        #endregion

    }
}