示例#1
0
        private static Series HGHTarget(FloodControlPoint pt, double forecastValue, DateTime t1, DateTime t2)
        {
            HydrometRuleCurve m_ruleCurve = RuleCurveFactory.Create(pt, 7100);
            Series            target      = new Series("Target");

            string   flag = "";
            DateTime t    = t1;

            while (t < t2.AddDays(-1))
            {
                t = t.AddDays(1).Date;
                double val = -m_ruleCurve.LookupRequiredSpace(t, forecastValue, out flag) * pt.PercentSpace / 100 + pt.TotalUpstreamActiveSpace;
                target.Add(t, val);
            }

            return(target);
        }
示例#2
0
        /// <summary>
        /// Create a list of targts
        /// </summary>
        /// <param name="pt"></param>
        /// <param name="waterYear"></param>
        /// <param name="start"></param>
        /// <param name="optionalPercents"></param>
        /// <param name="m_ruleCurve"></param>
        /// <param name="t2"></param>
        /// <param name="forecastMonth"></param>
        /// <param name="forecastValue"></param>
        /// <returns></returns>
        private static SeriesList GetTargets(FloodControlPoint pt, int waterYear, Point start, int[] optionalPercents,
                                             int forecastMonth, double forecastValue, bool dashed)
        {
            SeriesList        rval        = new SeriesList();
            HydrometRuleCurve m_ruleCurve = RuleCurveFactory.Create(pt, 7100, dashed);
            var t1 = new DateTime(waterYear, pt.ForecastMonthStart, 1);
            var t2 = new DateTime(waterYear, pt.ForecastMonthEnd, 1).EndOfMonth();

            var avg30yrQU = Get30YearAverageSeries(pt.DailyStationQU, "qu", forecastMonth);

            // sum volume for the forecast period
            var    t = new DateTime(start.DateTime.Year, start.DateTime.Month, 1);
            double historicalAverageResidual = SumResidual(avg30yrQU, t, t2);

            // capture if running the process before forecast season starts (before Jan)
            if (DateTime.Now.Month > 9)
            {
                historicalAverageResidual = SumResidual(avg30yrQU, t.AddYears(1), t2.AddYears(1));
            }
            double percent = forecastValue / historicalAverageResidual;

            //get thirty year average QU from daily
            //avg30yrQU = Get30YearAverageSeries(pt.DailyStationQU, "qu", forecastMonth);
            Series targetx = CalculateTarget(pt, percent, waterYear, start, m_ruleCurve, avg30yrQU, t2, forecastValue);

            targetx.Name = "Forecast (" + (100 * percent).ToString("F0") + "%)" + (forecastValue / 1000.0).ToString("F0");
            //targetx.Add(start);
            rval.Add(targetx);

            for (int i = 0; i < optionalPercents.Length; i++)
            {
                var fc = historicalAverageResidual * optionalPercents[i] / 100.0;
                targetx = CalculateTarget(pt, optionalPercents[i] / 100.0, waterYear, start,
                                          m_ruleCurve, avg30yrQU
                                          , t2, fc);
                targetx.Name = "Target (" + optionalPercents[i].ToString("F0") + "%) ";
                //targetx.Add(start);
                rval.Add(targetx);
            }

            return(rval);
        }