public bool CalculateWindSpeed(double shearht)
{
try
{
//make sure there are axis vals set
if (_xaxis == null || _yaxis == null)
{
return(false);
}
Console.WriteLine(" x axis " + _xaxis.AxisValues.Length + " y axis " + _yaxis.AxisValues.Length);
IAxis hourAxis = null;
IAxis monthAxis = null;
IAxis wdAxis = null;
//add axes to a list for easier parsing
List <IAxis> AssignedAxis = new List <IAxis>();
AssignedAxis.Add(_xaxis);
AssignedAxis.Add(_yaxis);
bool useWD = false;
foreach (IAxis a in AssignedAxis)
{
switch (a.AxisType)
{
case AxisType.Hour:
hourAxis = a;
break;
case AxisType.Month:
monthAxis = a;
break;
case AxisType.WD:
wdAxis = a;
useWD = true;
throw new NotImplementedException("shear by WD not yet developed");
default:
break;
}
}
Console.WriteLine("assigned axis types");
_data.Sort = Settings.Default.TimeStampName;
DateTime t = DateTime.Now;
double missing = Settings.Default.MissingValue;
//get a date at the beginning of the dataset provided
DateTime start = DateTime.Parse(_data[2][_collection.DateIndex].ToString());
//get a list of the columns for this config sorted ascending by height
SortedDictionary <double, ISessionColumn> orderedWsCols
= _collection.GetColumnsByType(SessionColumnType.WSAvg, start, true);
//sort columns descending by height
orderedWsCols.Reverse();
//put the column indexes into a list for coincident values constructor
SortedDictionary <double, int> TopCols = new SortedDictionary <double, int>();
foreach (KeyValuePair <double, ISessionColumn> kv in orderedWsCols)
{
TopCols.Add(kv.Key, kv.Value.ColIndex);
}
List <double> SensorHts = TopCols.Keys.ToList();
SensorHts.Sort();
//Console.WriteLine("first row after sort " + _data[1][0] + " " + _[SensorHts [1]] + " " + _data[1][SensorHts[0]]);
//pass the column list and the dataview to get only coincident rows for the columns passed
XbyYShearCoincidentRowCollection coincident
= new XbyYShearCoincidentRowCollection(_collection.DateIndex, TopCols, _data);
int wdColIndex = _collection["WDAvgComposite"].ColIndex;
Console.WriteLine("wd index=" + wdColIndex);
List <XbyYCoincidentRow> coincidentValues = coincident.GetRows(wdColIndex);
Console.WriteLine("coincident rows count " + coincidentValues.Count);
//collection to hold averages at each x y point
List <double[, ]> GridCollection = new List <double[, ]>();
//get number of elements in each axis
int hour = 0;
int month = 0;
int x = _xaxis.AxisValues.Length;
int y = _yaxis.AxisValues.Length;
Console.WriteLine("assigned axis lengths " + x + ", " + y);
//arrays to hold upper and lower avgs at each x,y point
double[,] UpperAvgGrid = new double[x, y];
double[,] LowerAvgGrid = new double[x, y];
//list to hold avg values for each sensor grid at each x and y bin
List <double> upperWsVals = new List <double>();
List <double> lowerWSVals = new List <double>();
for (int a = 0; a < x; a++)
{
for (int j = 0; j < y; j++)
{
int TestCount = 0;
if (useWD)
{
}
else
{
if (_xaxis.AxisType == AxisType.Hour)
{
hour = a; month = j;
}
else
{
hour = j; month = a;
}
hour = (int)hourAxis.AxisValues[hour];
month = (int)monthAxis.AxisValues[month];
//Console.WriteLine("hour " + hour + ", Month " + month);
//find the row collection for this x,y point
var r = from s in coincidentValues.AsEnumerable()
where s.Date.Month == month && s.Date.Hour == hour &&
s.WD >= 0
select s;
//Console.WriteLine(month + " " + hour + " value: " + r.Count());
foreach (var row in r)
{
//Console.WriteLine("upper ws " + (double)row.UpperWS );
upperWsVals.Add((double)row.UpperWS);
lowerWSVals.Add((double)row.LowerWS);
TestCount++;
}
}
//Console.WriteLine("upper count at " + month + ", " + hour + ": " +
//TestCount);
//add each collection avg to the avg grids
//Console.WriteLine("lower average at " + month + ", " + hour + ": " + lowerWSVals.Average() + "/n");
UpperAvgGrid[a, j] = upperWsVals.Average();
LowerAvgGrid[a, j] = lowerWSVals.Average();
upperWsVals.Clear();
lowerWSVals.Clear();
//Console.WriteLine("Upper average at " + a + ", " + hour + " "+ UpperAvgGrid[a, j]);
//Console.WriteLine("Lower average at " + month + ", " + hour + " " + LowerAvgGrid[a, j]);
}
}
GridCollection.Add(LowerAvgGrid);
GridCollection.Add(UpperAvgGrid);
Console.WriteLine("Grid Collection now contains " + GridCollection.Count + " grids");
//calculate each alpha value and assign to a single grid
double[,] alphaGrid = new double[x, y];
//fill alpha table
for (int a = 0; a < x; a++)
{
for (int j = 0; j < y; j++)
{
//Console.WriteLine(" x=" + a.ToString() + " y=" + j.ToString());
alphaGrid[a, j] = Math.Log(GridCollection[1][a, j] / GridCollection[0][a, j])
/ Math.Log(SensorHts[1] / SensorHts[0]);
Console.WriteLine("Alpha in grid=" + alphaGrid[a, j].ToString() +
" upper ws =" + GridCollection[1][a, j].ToString() +
" Lower ws =" + GridCollection[0][a, j].ToString() +
" upper ht " + SensorHts[1] + " lower ht " + SensorHts[0] +
" month=" + (a + 1) + " hour=" + j);
}
}
GridCollection.Add(alphaGrid);
//create a single multi-d array from the original dictionary of values
//values must stay in order
List <int> CombineCols = TopCols.Values.ToList();
//height and value of the valid ws comp at each row ordered by datetime
double[,] result = CombineSensorVals(CombineCols);
Console.WriteLine("Rows in full data " + _data.Count + " rows in combined dataset " + result.GetLength(1));
//array to store sheared ws
List <double> ShearResult = new List <double>(result.GetLength(1));
//add colummn to dataset and capture colindex
string localcolname = shearht.ToString();
localcolname += Enum.GetName(typeof(SessionColumnType), SessionColumnType.WSAvgShear);
localcolname += x + "by" + y;
Console.WriteLine("X by Y shear column will be named " + localcolname);
int shearIndex = AddShearColumnToTable(localcolname, shearht);
double thisAlpha;
int xVal;
int yVal;
DateTime thisDate;
double thisWD;
//loop combined values
for (int i = 0; i < result.GetLength(1); i++)
{
thisDate = (DateTime)_data[i][_collection.DateIndex];
if (useWD)
{
thisWD = (double)_data[i][wdColIndex];
}
else
{
xVal = _xaxis.ReturnAxisValue(thisDate);
yVal = _yaxis.ReturnAxisValue(thisDate);
thisAlpha = GridCollection[2][xVal, yVal];
if (result[1, i] >= 0)
{
_data[i][shearIndex] = (result[1, i] * Math.Pow((shearht / result[0, i]), thisAlpha));
}
else
{
_data[i][shearIndex] = missing;
}
}
}
_data.Table.AcceptChanges();
DateTime t_end = DateTime.Now;
TimeSpan dur = t_end - t;
Console.WriteLine("Calculating and adding " + x + " by " + y + " shear ws takes " + dur);
return(true);
}
catch (Exception e)
{
throw e;
}
}
