public static void ageReclass(map8 m)
{
m.dim((uint)snr, (uint)snc);
m.rename("Age class representation");
for (uint j = 1; j < map8.MapmaxValue; j++)
{
m.assignLeg(j, "");
}
string str;
//J.Yang hard coding changing itr*sites.TimeStep to itr
//J.Yang maxLeg is defined as 256 in map8.h, therefore, maximum age cohorts it can output is 254
for (uint i = 1; i < map8.MaxValueforLegend - 4; i++)
{
str = string.Format("{0: } - {1: } yr", (i - 1) * time_step + 1, i * time_step);
m.assignLeg(i, str);
}
m.assignLeg(0, "NoSpecies");
m.assignLeg(map8.MaxValueforLegend - 1, "N/A");
m.assignLeg(map8.MaxValueforLegend - 2, "Water");
m.assignLeg(map8.MaxValueforLegend - 3, "NonForest");
str = string.Format(" >{0: } yr", (map8.MaxValueforLegend - 4 - 1) * time_step);
m.assignLeg(map8.MaxValueforLegend - 4, str);
for (int i = snr; i >= 1; i--)
{
for (int j = 1; j <= snc; j++)
{
if (Program.sites.locateLanduPt(i, j).active())
{
m[(uint)i, (uint)j] = 0;
uint myage = 0;
Site local_site = Program.sites[i, j];
Specie s = local_site.first();
while (s != null)
{
uint temp = (uint)s.oldest();
if (temp > myage)
{
myage = temp;
}
s = local_site.next();
}
m[(uint)i, (uint)j] = (ushort)(myage / time_step);
}
else if (Program.sites.locateLanduPt(i, j).lowland())
{
m[(uint)i, (uint)j] = (ushort)(map8.MaxValueforLegend - 3);
}
else if (Program.sites.locateLanduPt(i, j).water())
{
m[(uint)i, (uint)j] = (ushort)(map8.MaxValueforLegend - 2);
}
else
{
m[(uint)i, (uint)j] = (ushort)(map8.MaxValueforLegend - 1);
}
}
}
}
//This will copy the contents of another map into the current map.
//This includes legends, titles, and current display dimensions.
public void copy(map8 n)
{
largeCell = n.largeCell;
dim(n.numRows, n.numCols);
title = n.title;
for (int i = 0; i < maxLeg; i++) //the original part might be problematic : for (int i = 0; i < 22; i++)
{
legend[i] = n.legend[i];
}
Array.Copy(n.data, data, n.data.Length);
}
private static void putHarvestOutput(int itr, double[] wAdfGeoTransform)
{
map8 m = new map8(sites.Header);
string str = null;
string str1 = null;
string str_htyp = null;
string str_htyp1 = null;
string str_dec = null;
string str_dec1 = null;
Harvest.GlobalFunctions.writeStandReport();
if (itr * freq[3] <= parameters.numberOfIterations)
{
str_htyp1 = string.Format("{0}/{1}/htyp{2:D}.img", parameters.outputDir, "Harvest", itr * freq[3] * sites.TimeStep_Harvest); //*
str_dec1 = string.Format("{0}/{1}/hdec{2:D}.img", parameters.outputDir, "Harvest", itr * freq[3] * sites.TimeStep_Harvest); //*
Harvest.GlobalFunctions.output_harvest_Dec_Type(itr, str_htyp, str_htyp1, str_dec, str_dec1, wAdfGeoTransform); //*
for (int i = 0; i < sites.specNum; i++)
{
str1 = string.Format("{0}/{1}/{2}_BACut{3:D}.img", parameters.outputDir, "Harvest", speciesAttrs[i + 1].name, itr * freq[3] * sites.TimeStep_Harvest); //*
Harvest.GlobalFunctions.PutOutput_harvestBACut_spec(str, str1, i, wAdfGeoTransform); //*
}
str1 = string.Format("{0}/{1}/BACut{2:D}.img", parameters.outputDir, "Harvest", itr * freq[3] * sites.TimeStep_Harvest); //*
Harvest.GlobalFunctions.PutOutput_harvestBACut(str, str1, wAdfGeoTransform);
}
if (parameters.numberOfIterations % freq[3] != 0)
{
str_htyp1 = string.Format("{0}/{1}/htyp{2:D}.img", parameters.outputDir, "Harvest", parameters.numberOfIterations * sites.TimeStep_Harvest); //*
str_dec1 = string.Format("{0}/{1}/hdec{2:D}.img", parameters.outputDir, "Harvest", parameters.numberOfIterations * sites.TimeStep_Harvest); //*
Harvest.GlobalFunctions.output_harvest_Dec_Type(itr, str_htyp, str_htyp1, str_dec, str_dec1, wAdfGeoTransform); //*
for (int i = 0; i < sites.specNum; i++)
{
str1 = string.Format("{0}/{1}/{2}_BACut{3:D}.img", parameters.outputDir, "Harvest", speciesAttrs[i + 1].name, parameters.numberOfIterations * sites.TimeStep_Harvest); //*
Harvest.GlobalFunctions.PutOutput_harvestBACut_spec(str, str1, i, wAdfGeoTransform); //*
}
str1 = string.Format("{0}/{1}/BACut{2:D}.img", parameters.outputDir, "Harvest", parameters.numberOfIterations * sites.TimeStep_Harvest); //*
Harvest.GlobalFunctions.PutOutput_harvestBACut(str, str1, wAdfGeoTransform);
}
}
public static void speciesAgeMap(map8 m, string ageFile)
{
int curSp = Program.speciesAttrs.current(ageFile);
m.dim((uint)snr, (uint)snc);
m.rename(ageFile);
string str;
for (uint i = 1; i < map8.maxLeg - 4; i++)
{
str = string.Format("{0: } - {1: } yr", (i - 1) * time_step + 1, i * time_step);
m.assignLeg(i, str);
}
m.assignLeg(0, "NotPresent");
m.assignLeg(map8.MaxValueforLegend - 1, "N/A");
m.assignLeg(map8.MaxValueforLegend - 2, "Water");
m.assignLeg(map8.MaxValueforLegend - 3, "NonForest");
str = string.Format(" >{0} yr", (map8.maxLeg - 4 - 1) * time_step);
m.assignLeg(map8.MaxValueforLegend - 4, str);
for (int i = snr; i >= 1; i--)
{
for (int j = 1; j <= snc; j++)
{
if (Program.sites.locateLanduPt(i, j) == null)
{
throw new Exception("Invalid landunit error\n");
}
if (Program.sites.locateLanduPt(i, j).active())
{
m[(uint)i, (uint)j] = 0; //where species not presents
if (Program.sites[i, j] == null)
{
throw new Exception("No site\n");
}
Specie s = Program.sites[i, j].current(curSp);
if (s == null)
{
Console.WriteLine("{0}\n", curSp);
throw new Exception("No Species\n");
}
if (s.query())
{
m[(uint)i, (uint)j] = (ushort)(s.oldest() / time_step); //compare ageReclass which uses +3 there???
if (m[(uint)i, (uint)j] > map8.MaxValueforLegend - 4) //maximum longevity is 640 years// Notice 66 means 640 years
{
m[(uint)i, (uint)j] = (ushort)(map8.MaxValueforLegend - 4);
}
}
}
else if (Program.sites.locateLanduPt(i, j).water())
{
m[(uint)i, (uint)j] = (ushort)(map8.MaxValueforLegend - 2);
}
else if (Program.sites.locateLanduPt(i, j).lowland())
{
m[(uint)i, (uint)j] = (ushort)(map8.MaxValueforLegend - 3);
}
else
{
m[(uint)i, (uint)j] = (ushort)(map8.MaxValueforLegend - 1);
}
}
}
}