LeafWeb/Core/CharterOld/CurveParamSet.cs

using System.Collections.Generic;
using System.IO;
using LeafWeb.Core.Utility;

namespace LeafWeb.Core.CharterOld
{
	public class CurveParamSet
	{
		public List<XyPoint> AnetMeasChloro1Data { get; } = new List<XyPoint>(); // y=AnetMeas column, x=PCO2c, for PointLimitType=1
		public List<XyPoint> AnetMeasChloro2Data { get; } = new List<XyPoint>(); // y=AnetMeas column, x=PCO2c, for PointLimitType=2
		public List<XyPoint> AnetMeasChloro3Data { get; } = new List<XyPoint>(); // y=AnetMeas column, x=PCO2c, for PointLimitType=3
		public List<XyPoint> AnetMeasInter1Data { get; } = new List<XyPoint>(); // y=AnetMeas column, x=PCO2i, for PointLimitType=1
		public List<XyPoint> AnetMeasInter2Data { get; } = new List<XyPoint>(); // y=AnetMeas column, x=PCO2i, for PointLimitType=2
		public List<XyPoint> AnetMeasInter3Data { get; } = new List<XyPoint>(); // y=AnetMeas column, x=PCO2i, for PointLimitType=3

		public List<XyPoint> AcChloroData { get; } = new List<XyPoint>();
		public List<XyPoint> AjChloroData { get; } = new List<XyPoint>();
		public List<XyPoint> AtChloroData { get; } = new List<XyPoint>();
		public List<XyPoint> AcInterData { get; } = new List<XyPoint>();
		public List<XyPoint> AjInterData { get; } = new List<XyPoint>();
		public List<XyPoint> AtInterData { get; } = new List<XyPoint>();
		
		public CurveParamSet(TextReader sr, ref int lineNbr, ref string curveId)
		{
			bool curveIdSet = false, doneWithAnet = false;
			string line;
			List<string> phrases;

			while (!doneWithAnet)
			{
				lineNbr++;
				line = sr.ReadLine();
				if (line == null)
				{
					throw new ParseException("Unexpected end-of-file at line " + lineNbr);
				}

				phrases = SplitCsvLine(line);

				var firstField = phrases[0];
				if (firstField.Equals("CO2i"))
				{
					doneWithAnet = true;
				}
				else
				{
					// The fields on the line:
					// Column  Name
					//   0     CurveID
					//   1     ChlFlUse
					//   2     FitGi
					//   3     FitGamma
					//   4     FitKco
					//   5     FitRd
					//   6     FitAlpha
					//   7     LimitCombina
					//   8     PCO2i
					//   9     PCO2c
					//  10     AnetMeas
					//  11     AnetCal
					//  12     weitedrms
					//  13     PointLimitType

					if (!curveIdSet)
					{
						curveId = firstField;
						curveIdSet = true;
					}
					var xyPoint1 = new XyPoint(phrases[9], phrases[10]); // AnetMeas(y), PCO2c(x)
					var xyPoint2 = new XyPoint(phrases[8], phrases[10]);
					var pointLimitType = int.Parse(phrases[13]);
					switch (pointLimitType)
					{
						case 1:
							AnetMeasChloro1Data.Add(xyPoint1);
							AnetMeasInter1Data.Add(xyPoint2);
							break;
						case 2:
							AnetMeasChloro2Data.Add(xyPoint1);
							AnetMeasInter2Data.Add(xyPoint2);
							break;
						case 3:
							AnetMeasChloro3Data.Add(xyPoint1);
							AnetMeasInter3Data.Add(xyPoint2);
							break;
					}
				}
			}

			// The next set of lines will have three pairs of x,y-coordinates to save.
			// A blank line signals the end of the data.

			var moreData = true;
			while (moreData)
			{
				// The fields on the line:
				// Column  Name
				//   0     CO2i
				//   1     CO2cc
				//   2     Ac
				//   3     CO2cj
				//   4     Aj
				//   5     CO2ct
				//   6     At

				lineNbr++;
				line = sr.ReadLine();
				if (line == null)
				{
					throw new ParseException("Unexpected end-of-file at line " + lineNbr);
				}

				if (line.Length == 0)
					moreData = false;
				else
				{
					phrases = SplitCsvLine(line);
					var xyPoint1 = new XyPoint(phrases[1],phrases[2]); // Ac(y),CO2cc(x)
					if (xyPoint1.YIsInRange(-20.0, 50.0))
						AcChloroData.Add(xyPoint1);
					xyPoint1 = new XyPoint(phrases[3], phrases[4]); // Aj(y),CO2cj(x)
					if (xyPoint1.YIsInRange(-20.0, 50.0))
						AjChloroData.Add(xyPoint1);
					xyPoint1 = new XyPoint(phrases[5], phrases[6]); // At(y),CO2ct(x)
					if (xyPoint1.YIsInRange(-20.0, 50.0))
						AtChloroData.Add(xyPoint1);

					xyPoint1 = new XyPoint(phrases[0], phrases[2]); // Ac(y),CO2i(x)
					if (xyPoint1.YIsInRange(-20.0, 50.0))
						AcInterData.Add(xyPoint1);
					xyPoint1 = new XyPoint(phrases[0], phrases[4]); // Aj(y),CO2i(x)
					if (xyPoint1.YIsInRange(-20.0, 50.0))
						AjInterData.Add(xyPoint1);
					xyPoint1 = new XyPoint(phrases[0], phrases[6]); // At(y),CO2i(x)
					if (xyPoint1.YIsInRange(-20.0, 50.0))
						AtInterData.Add(xyPoint1);
				}
			}
		}

		/// <summary>
		///     This method assumes that the argument is a comma-, blank-, or
		///     tab-separated set of strings. It returns an ArrayList of those
		///     quantities. Consecutive commas will be returned as an empty string,
		///     but all empty strings  at the end of the line will be thrown away.
		/// </summary>
		/// <param name="line"></param>
		/// <returns></returns>
		private static List<string> SplitCsvLine(string line)
		{
			int i;
			var separator = new [] {',', ' ', '\t'};
			var phrases = line.Split(separator);

			var retPhrases = new List<string>();
			for (i = 0; i < phrases.Length; i++)
			{
				var phrase = phrases[i].Trim();
				if (phrase.Length > 0)
					retPhrases.Add(phrase);
			}

			return retPhrases;
		}
	}
}