//---------------------------------------------------------------------
protected override Dictionary <int, IDynamicInputRecord[, ]> Parse()
{
//InputVar<string> landisData = new InputVar<string>("LandisData");
//ReadVar(landisData);
//if (landisData.Value.Actual != FileName)
// throw new InputValueException(landisData.Value.String, "The value is not \"{0}\"", PlugIn.ExtensionName);
ReadLandisDataVar();
Dictionary <int, IDynamicInputRecord[, ]> allData = new Dictionary <int, IDynamicInputRecord[, ]>();
//---------------------------------------------------------------------
//Read in establishment probability data:
InputVar <int> year = new InputVar <int>("Time step for updating values");
InputVar <string> ecoregionName = new InputVar <string>("Ecoregion Name");
InputVar <string> speciesName = new InputVar <string>("Species Name");
InputVar <double> pest = new InputVar <double>("Probability of Establishment");
while (!AtEndOfInput)
{
StringReader currentLine = new StringReader(CurrentLine);
ReadValue(year, currentLine);
int yr = year.Value.Actual;
if (!allData.ContainsKey(yr))
{
IDynamicInputRecord[,] inputTable = new IDynamicInputRecord[EcoregionData.ModelCore.Species.Count, EcoregionData.ModelCore.Ecoregions.Count];
allData.Add(yr, inputTable);
EcoregionData.ModelCore.UI.WriteLine(" Dynamic Input Parser: Add new year = {0}.", yr);
}
ReadValue(ecoregionName, currentLine);
IEcoregion ecoregion = GetEcoregion(ecoregionName.Value);
ReadValue(speciesName, currentLine);
ISpecies species = GetSpecies(speciesName.Value);
IDynamicInputRecord dynamicInputRecord = new DynamicInputRecord();
ReadValue(pest, currentLine);
dynamicInputRecord.ProbEst = pest.Value;
allData[yr][species.Index, ecoregion.Index] = dynamicInputRecord;
CheckNoDataAfter("the " + pest.Name + " column",
currentLine);
GetNextLine();
}
return(allData);
}
public static IDynamicInputRecord Clone(IDynamicInputRecord oldRecord)
{
IDynamicInputRecord newRecord = new DynamicInputRecord();
newRecord.EcoIgnitionNum = oldRecord.EcoIgnitionNum;
newRecord.FallFMCHi = oldRecord.FallFMCHi;
newRecord.FallFMCHiProp = oldRecord.FallFMCHiProp;
newRecord.FallFMCLo = oldRecord.FallFMCLo;
newRecord.FireRegionSites = oldRecord.FireRegionSites;
newRecord.Index = oldRecord.Index;
newRecord.MapCode = oldRecord.MapCode;
newRecord.MaxSize = oldRecord.MaxSize;
newRecord.MeanSize = oldRecord.MeanSize;
newRecord.MinSize = oldRecord.MinSize;
newRecord.Name = oldRecord.Name;
newRecord.OpenFuelType = oldRecord.OpenFuelType;
newRecord.SpringFMCHi = oldRecord.SpringFMCHi;
newRecord.SpringFMCHiProp = oldRecord.SpringFMCHiProp;
newRecord.SpringFMCLo = oldRecord.SpringFMCLo;
newRecord.StandardDeviation = oldRecord.StandardDeviation;
newRecord.SummerFMCHi = oldRecord.SummerFMCHi;
newRecord.SummerFMCHiProp = oldRecord.SummerFMCHiProp;
newRecord.SummerFMCLo = oldRecord.SummerFMCLo;
return newRecord;
}
//---------------------------------------------------------------------
public static DataTable ReadWeatherFile(string path, IDynamicInputRecord[] regionRecords, ISeasonParameters[] seasonParms)
{
PlugIn.ModelCore.UI.WriteLine(" Loading Weather Data...");
CSVParser weatherParser = new CSVParser();
DataTable weatherTable = weatherParser.ParseToDataTable(path);
int recordCount = 0;
for (int i = 0; i <= 2; i++)
{
string seasName = seasonParms[i].NameOfSeason.ToString();
foreach (IDynamicInputRecord fire_region in regionRecords)
{
string ecoName = fire_region.Name;
//PlugIn.ModelCore.Log.WriteLine("Read Weather File: Season={0}, FireRegion={1}.", seasName, ecoName);
string selectText = ("FireRegion = '" + ecoName + "' AND Season = '" + seasName + "'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
DataRow[] foundRows = weatherTable.Select(selectText);
if (foundRows.Length == 0)
{
selectText = ("FireRegion = 'All' AND Season = '" + seasName + "'");
foundRows = weatherTable.Select(selectText);
}
if ((foundRows.Length == 0) && (seasonParms[i].FireProbability > 0) && (fire_region.EcoIgnitionNum > 0))
{
PlugIn.ModelCore.UI.WriteLine("WARNING (Ln 430): FireRegion " + ecoName + ", Season " + seasName + " has fire probability > 0, but 0 weather records. Using alternate season.");
if (seasName == "Spring")
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Summer'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
if (foundRows.Length == 0)
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Fall'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
}
if (foundRows.Length == 0)
{
PlugIn.ModelCore.UI.WriteLine("WARNING (Ln 445): FireRegion " + ecoName + " has fire probability > 0, but 0 weather records. No fires will occur in this fire region.");
}
}
if (seasName == "Fall")
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Summer'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
if (foundRows.Length == 0)
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Spring'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
}
if (foundRows.Length == 0)
{
PlugIn.ModelCore.UI.WriteLine("WARNING (Ln 464): FireRegion " + ecoName + " has fire probability > 0, but 0 weather records. No fires will occur in this fire region.");
}
}
if (seasName == "Summer")
{
if (PlugIn.ModelCore.GenerateUniform() >= 0.5)
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Spring'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
if (foundRows.Length == 0)
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Fall'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
}
if (foundRows.Length == 0)
{
PlugIn.ModelCore.UI.WriteLine("WARNING (Ln 485): FireRegion " + ecoName + " has fire probability > 0, but 0 weather records. No fires will occur in this fire region.");
}
}
else
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Fall'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
if (foundRows.Length == 0)
{
selectText = ("FireRegion = '" + ecoName + "' AND Season = 'Spring'");
//PlugIn.ModelCore.Log.WriteLine("Read Weather File SelectText = {0}.", selectText);
foundRows = weatherTable.Select(selectText);
}
if (foundRows.Length == 0)
{
PlugIn.ModelCore.UI.WriteLine("WARNING (Ln 503): FireRegion " + ecoName + " has fire probability > 0, but 0 weather records. No fires will occur in this fire region.");
}
}
}
}
if(foundRows.Length > 0)
{
//Input validation
double WSV, FFMC, BUI;
int WINDDIR, FWIBIN;
for(int j = 0; j < foundRows.Length; j ++) //weatherDataSet.Tables["Table"].Rows)
{
DataRow myDataRow = foundRows[j];
//WSV = (double)myDataRow["WSV"];
WSV = Convert.ToDouble(myDataRow["WSV"]);
//Console.WriteLine("WSV: {0}", WSV);
if (WSV < 0.0)
{
throw new System.ApplicationException("Error: Wind Speed < 0: FireRegion = " + ecoName + "; Season = " + seasName);
}
FFMC = Convert.ToDouble(myDataRow["FFMC"]);
//Console.WriteLine("FFMC: {0}", FFMC);
if (FFMC < 0.0)
{
throw new System.ApplicationException("Error: FFMC < 0: FireRegion = " + ecoName + "; Season = " + seasName);
}
else if (FFMC > 100.0)
{
throw new System.ApplicationException("Error: FFMC > 100: FireRegion = " + ecoName + "; Season = " + seasName);
}
BUI = Convert.ToDouble(myDataRow["BUI"]);
//Console.WriteLine("BUI: {0}", BUI);
if (BUI < 0.0)
{
throw new System.ApplicationException("Error: BUI < 0: FireRegion = " + ecoName + "; Season = " + seasName);
}
WINDDIR = (int) myDataRow["WindDir"];
//Console.WriteLine("WindDir: {0}", WINDDIR);
if (WINDDIR < 0)
{
throw new System.ApplicationException("Error: WINDDIR < 0: FireRegion = " + ecoName + "; Season = " + seasName);
}
else if (WINDDIR > 360)
{
throw new System.ApplicationException("Error: WINDDIR > 360: FireRegion = " + ecoName + "; Season = " + seasName);
}
FWIBIN = (int) myDataRow["FWIBin"];
//Console.WriteLine("FWIBIN: {0}", FWIBIN);
if (FWIBIN < 1)
{
throw new System.ApplicationException("Error: FWIBIN < 1: FireRegion = " + ecoName + "; Season = " + seasName);
}
else if (FWIBIN > 5)
{
throw new System.ApplicationException("Error: FWIBIN > 5: FireRegion = " + ecoName + "; Season = " + seasName);
}
}
}
if(seasName == "Fall")
{
fire_region.FallRecords = recordCount;
}
else if (seasName == "Spring")
{
fire_region.SpringRecords = recordCount;
}
else if (seasName == "Summer")
{
fire_region.SummerRecords = recordCount;
}
}
}
return weatherTable;
}
//---------------------------------------------------------------------
public static int GenerateFMC(ISeasonParameters season, IDynamicInputRecord fire_region)
{
int FMC = 0;
if (season.NameOfSeason == SeasonName.Spring)
{
if (PlugIn.ModelCore.GenerateUniform() < fire_region.SpringFMCHiProp)
FMC = fire_region.SpringFMCHi;
else
FMC = fire_region.SpringFMCLo;
}
if (season.NameOfSeason == SeasonName.Summer)
{
if (PlugIn.ModelCore.GenerateUniform() < fire_region.SummerFMCHiProp)
FMC = fire_region.SummerFMCHi;
else
FMC = fire_region.SummerFMCLo;
}
if (season.NameOfSeason == SeasonName.Fall)
{
if (PlugIn.ModelCore.GenerateUniform() < fire_region.FallFMCHiProp)
FMC = fire_region.FallFMCHi;
else
FMC = fire_region.FallFMCLo;
}
return FMC;
}
//---------------------------------------------------------------------
public static DataRow GenerateDataRow(ISeasonParameters season, IDynamicInputRecord fire_region, int sizeBin)
{
int weatherRandomizer = PlugIn.WeatherRandomizer;
string seasonName = season.NameOfSeason.ToString();
string ecoName = fire_region.Name;
int weatherBin = 0;
if(weatherRandomizer == 0)
weatherBin = sizeBin;
else
{
if (weatherRandomizer == 4)
{
weatherBin = -1;
}
else
{
// First, tally the available bins and assign a probability based on their size
// (number of records).
// Bins can only be 1 - 5.
int minBin = (int)System.Math.Max(1, sizeBin - weatherRandomizer);
int maxBin = (int)System.Math.Min(5, sizeBin + weatherRandomizer);
double[] binProbability = new double[maxBin - minBin + 1];
int[] binCount = new int[maxBin - minBin + 1];
int binSum = 0;
for (int bin = minBin; bin <= maxBin; bin++)
{
string selectString = "FWIBin = '" + bin + "' AND Season = '" + seasonName + "' AND FireRegion = '" + ecoName + "'";
DataRow[] rows = PlugIn.WeatherDataTable.Select(selectString);
if (rows.Length == 0)
{
selectString = "FWIBin = '" + bin + "' AND Season = '" + seasonName + "' AND FireRegion = 'All'";
rows = PlugIn.WeatherDataTable.Select(selectString);
}
int numRecs = rows.Length;
binCount[bin - minBin] = numRecs;
binSum += numRecs;
}
for (int bin = minBin; bin <= maxBin; bin++)
binProbability[bin - minBin] = (double)binCount[bin - minBin] / (double)binSum;
// Now randomly select from the available bins:
double randomBinNum = PlugIn.ModelCore.GenerateUniform();
double minProb = 0.0;
double maxProb = 0.0;
for (int bin = minBin; bin <= maxBin; bin++)
{
maxProb += binProbability[bin - minBin];
if (randomBinNum >= minProb && randomBinNum < maxProb)
{
weatherBin = bin;
break;
}
else
{
minProb = binProbability[bin - minBin];
}
}
}
}
if(weatherBin == 0)
weatherBin = sizeBin;
//throw new System.ApplicationException("No Weather Bin randomly selected. FireRegion = "+ecoName+", Season = "+seasonName+", sizeBin = "+sizeBin);
int rowCount = 0;
int loopCount = 0;
int firstDir = 0;
DataRow[] foundRows = null;
string seasonOrig = seasonName;
if (PlugIn.ModelCore.GenerateUniform() >= 0.5)
firstDir = 1; // Direction (+ or -) to check first if target bin is not available (+1 or -1)
else
firstDir = -1;
while (rowCount == 0)
{
string selectString = "FWIBin = '" + weatherBin + "' AND Season = '" + seasonName + "' AND FireRegion = '" + ecoName + "'";
if (weatherBin == -1)
{
selectString = "Season = '" + seasonName + "' AND FireRegion = '" + ecoName + "'";
}
foundRows = PlugIn.WeatherDataTable.Select(selectString);
rowCount = foundRows.Length;
if (rowCount == 0)
{
selectString = "FWIBin = '" + weatherBin + "' AND Season = '" + seasonName + "' AND FireRegion = 'All'";
{
selectString = "Season = '" + seasonName + "' AND FireRegion = 'All'";
}
foundRows = PlugIn.WeatherDataTable.Select(selectString);
rowCount = foundRows.Length;
}
if (rowCount == 0)
{
//PlugIn.ModelCore.Log.WriteLine(" weatherBin "+weatherBin+" Not Found. Using alternate weatherBin.");
if (sizeBin == 5)
weatherBin = weatherBin - 1;
else if (sizeBin == 1)
weatherBin = weatherBin + 1;
else
{
if (loopCount == 0)
weatherBin = sizeBin + firstDir;
else if (loopCount == 1)
weatherBin = sizeBin - firstDir;
else if (loopCount == 2)
weatherBin = sizeBin + (firstDir * 2);
else if (loopCount == 3)
weatherBin = sizeBin - (firstDir * 2);
else if (loopCount == 4)
weatherBin = sizeBin + (firstDir * 3);
else if (loopCount == 5)
weatherBin = sizeBin - (firstDir * 3);
else if (loopCount == 6)
weatherBin = sizeBin + (firstDir * 4);
else
weatherBin = sizeBin - (firstDir * 4);
}
loopCount ++;
if (loopCount > 100)
{
PlugIn.ModelCore.UI.WriteLine(" No Weather Rows Selected. Using alternate season.");
if (seasonOrig == "Spring" || seasonOrig == "Fall")
seasonName = "Summer";
else // Summer
{
if (firstDir == 1)
{
seasonName = "Spring";
}
else
{
seasonName = "Fall";
}
}
}
if (loopCount > 200)
{
PlugIn.ModelCore.UI.WriteLine(" No Weather Rows Selected. Using alternate season 2.");
if (seasonOrig == "Spring")
seasonName = "Fall";
else
if (seasonOrig == "Fall")
seasonName = "Spring";
else // Summer
{
if (firstDir == 1)
{
seasonName = "Fall";
}
else
{
seasonName = "Spring";
}
}
}
if (loopCount > 300)
{
PlugIn.ModelCore.UI.WriteLine("WARNING (Ln 350): FireRegion " + ecoName + " has fire probability > 0, but 0 weather records. No fires will occur in this fire region.");
break;
}
}
}
if (rowCount > 0)
{
int newRandNum = (int)(Math.Round(PlugIn.ModelCore.GenerateUniform() * (rowCount - 1)));
DataRow weatherRow = foundRows[newRandNum];
return weatherRow;
}
else
{
return null;
}
}
