/*
* public static void updateLandtypeMap8(BinaryReader ltMapFile)
* {
* uint[] dest = new uint[32];
*
* //LDfread((char*)dest, 4, 32, ltMapFile);
* for (int i = 0; i < 32; i++)
* dest[i] = ltMapFile.ReadUInt32();
*
*
* int b16or8;
* if ((dest[1] & 0xff0000) == 0x020000)
* b16or8 = 16;
* else if ((dest[1] & 0xff0000) == 0)
* b16or8 = 8;
* else
* {
* b16or8 = -1;
* throw new Exception("Error: IO: Landtype map is neither 16 bit or 8 bit.");
* }
*
* uint nCols = dest[4];
* uint nRows = dest[5];
*
*
* uint xDim = gl_sites.numColumns;
* uint yDim = gl_sites.numRows;
*
* if ((nCols != xDim) || (nRows != yDim))
* throw new Exception("landtype map and species map do not match.");
*
*
*
* if (b16or8 == 8) //8 bit
* {
* for (uint i = yDim; i > 0; i--)
* {
* for (uint j = 1; j <= xDim; j++)
* {
* int coverType = ltMapFile.Read();
*
* if (coverType >= 0)
* gl_sites.fillinLanduPt(i, j, gl_landUnits[coverType]);
* else
* throw new Exception("illegel landtype class found.");
* }
* }
* }
* else if (b16or8 == 16) //16 bit
* {
* for (uint i = yDim; i > 0; i--)
* {
* for (uint j = 1; j <= xDim; j++)
* {
* int coverType = ltMapFile.ReadUInt16();
*
* if (coverType >= 0)
* gl_sites.fillinLanduPt(i, j, gl_landUnits[coverType]);
* else
* throw new Exception("illegel landtype class found.");
* }
* }
* }
*
* }
*/
//Main program. This contains start and shut down procedures as well as the main iteration loop.
public override void Run()
{
int i = modelCore.TimeSinceStart;
int i_d_timestep = i / gl_sites.SuccessionTimeStep;
for (int r = 0; r < gl_landUnits.Num_Landunits; ++r)
{
gl_landUnits[r].MinShade = Land_type_Attributes.get_min_shade(r);
float[] rd = new float[4];
rd[0] = Land_type_Attributes.get_gso(0, r);
rd[1] = Land_type_Attributes.get_gso(1, r);
rd[2] = Land_type_Attributes.get_gso(2, r);
rd[3] = gl_land_attrs.get_maxgso(i, r);
gl_landUnits[r].MaxRDArrayItem = rd;
gl_landUnits[r].MaxRD = rd[3];
}
//Simulation loops////////////////////////////////////////////////
//for (int i = 1; i <= numbOfIter * gl_sites.TimeStep; i++)
//{
if (i % gl_sites.SuccessionTimeStep == 0)
{
if (gl_param.FlagforSECFile == 3)
{
int index = i_d_timestep - 1;
if (index == 0)
{
SEC_landtypefiles.Clear();
SEC_gisfiles.Clear();
/*
* using (StreamReader FpSECfile = new StreamReader(gl_param.VarianceSECFile))
* {
* int num_of_files = system1.read_int(FpSECfile);
*
* for (int ii_count_num = 0; ii_count_num < num_of_files; ii_count_num++)
* {
* string SECfileMapGIS = system1.read_string(FpSECfile);
* string SECfileNametemp = system1.read_string(FpSECfile);
*
* SEC_landtypefiles.Add(SECfileNametemp);
* SEC_gisfiles.Add(SECfileMapGIS);
* }
* }
*/
if (index < gl_land_attrs.year_arr.Count)
{
/*string SECfileNametemp = SEC_landtypefiles[index];
*
* string SECfileMapGIS = SEC_gisfiles[index];
*
* gl_param.LandUnitFile = SECfileNametemp;
*
* gl_landUnits.read(gl_param.LandUnitFile);
*/
Console.Write("\nEnvironment parameter Updated.\n");
string SECfileMapGIS = gl_land_attrs.Get_new_landtype_map(index);
//strcpy(parameters.landUnitMapFile, SECfileMapGIS);
gl_param.LandImgMapFile = SECfileMapGIS;
/*Dataset simgFile = null;
*
* if ((simgFile = (Dataset)Gdal.Open(gl_param.LandImgMapFile, Access.GA_ReadOnly)) == null) //* landtype.img
* {
* Console.Write("Land Map Img file %s not found.\n", gl_param.LandImgMapFile);
*
* throw new Exception(gl_param.LandImgMapFile);
* }
* else
* {
* double[] adfGeoTransform = new double[6];
*
* simgFile.GetGeoTransform(adfGeoTransform);
* {
* for (int ii = 0; ii < 6; ii++)
* {
* wAdfGeoTransform[ii] = adfGeoTransform[ii];
* }
* }
*
* updateLandtypeImg8(simgFile);
*
* simgFile = null;
* }
*/
Console.WriteLine("\nEnvironment map Updated.");
landunit SECLog_use = gl_landUnits.first();
int ii_count = 0;
using (StreamWriter fpLogFileSEC = new StreamWriter(fpLogFileSEC_name))
{
fpLogFileSEC.Write("Year: {0}\n", i);
for (; ii_count < gl_landUnits.Num_Landunits; ii_count++)
{
fpLogFileSEC.Write("Landtype{0}:\n", ii_count);
for (int jj_count = 1; jj_count <= specAtNum; jj_count++)
{
fpLogFileSEC.Write("spec{0}: {1:N6}, ", jj_count, SECLog_use.probRepro(jj_count));
}
SECLog_use = gl_landUnits.next();
fpLogFileSEC.Write("\n");
}
}
}
}
if (index > 0)
{
if (index < SEC_landtypefiles.Count)
{
/*gl_param.LandUnitFile = SEC_landtypefiles[index];
* //gl_param.landUnitMapFile = SEC_gisfiles[index];
* gl_param.LandImgMapFile = SEC_gisfiles[index];
*
* gl_landUnits.read(gl_param.LandUnitFile);
*/
gl_param.LandImgMapFile = gl_land_attrs.Get_new_landtype_map(index);
Console.WriteLine("\nEnvironment parameter Updated.");
/*using(BinaryReader GISmap = new BinaryReader(File.Open(gl_param.landUnitMapFile, FileMode.Open)))
* {
* updateLandtypeMap8(GISmap);
* }*/
/*Dataset ltimgFile = null;
*
* if ((ltimgFile = (Dataset)Gdal.Open(gl_param.LandImgMapFile, Access.GA_ReadOnly)) == null) //* landtype.img
* {
* Console.WriteLine("Land Map Img file {0} not found.", gl_param.LandImgMapFile);
*
* throw new Exception(gl_param.LandImgMapFile);
* }
* else
* {
* double[] adfGeoTransform = new double[6];
*
* ltimgFile.GetGeoTransform(adfGeoTransform);
* {
* for (int ii = 0; ii < 6; ii++)
* {
* wAdfGeoTransform[ii] = adfGeoTransform[ii];
* }
* }
*
* updateLandtypeImg8(ltimgFile);
*
* ltimgFile = null;
*
* }*/
Console.WriteLine("\nEnvironment map Updated.");
landunit SECLog_use = gl_landUnits.first();
using (StreamWriter fpLogFileSEC = new StreamWriter(fpLogFileSEC_name))
{
fpLogFileSEC.Write("Year: {0}\n", i);
int ii_count = 0;
for (; ii_count < gl_landUnits.Num_Landunits; ii_count++)
{
fpLogFileSEC.Write("Landtype{0}:\n", ii_count);
for (int jj_count = 1; jj_count <= specAtNum; jj_count++)
{
fpLogFileSEC.Write("spec{0}: {1:N6}, ", jj_count, SECLog_use.probRepro(jj_count));
}
SECLog_use = gl_landUnits.next();
fpLogFileSEC.Write("\n");
}
}
}
}
}
}//end if
Console.WriteLine("Processing succession at Year {0}", i);
singularLandisIteration(i, pPDP);
if (i % gl_sites.SuccessionTimeStep == 0 || i == numbOfIter * gl_sites.SuccessionTimeStep)
{
int[] frequency = new int[6] {
1, 1, 1, 1, 1, 1
};
if (i % (gl_sites.SuccessionTimeStep * freq[0]) == 0 && i_d_timestep <= numbOfIter)
{
In_Output.putOutput_Landis70Pro(0, i_d_timestep, freq);
}
if (i == gl_sites.SuccessionTimeStep * numbOfIter)
{
In_Output.putOutput_Landis70Pro(0, numbOfIter, frequency);
}
if (i % (gl_sites.SuccessionTimeStep * freq[4]) == 0 && i_d_timestep <= numbOfIter)
{
In_Output.putOutput(0, i_d_timestep, freq);
}
if (i == gl_sites.SuccessionTimeStep * numbOfIter)
{
In_Output.putOutput(0, numbOfIter, frequency);
}
In_Output.putOutput_AgeDistStat(i_d_timestep);
}
//}
//Simulation loops end/////////////////////////////////////////////////
}//end Run()
public override void Initialize()
{
DateTime now = DateTime.Now;
//Gdal.AllRegister();
Console.WriteLine("Start Landis Pro at {0}", now);
Timestep = gl_param.SuccessionTimestep;
gl_sites.SuccessionTimeStep = gl_param.SuccessionTimestep;
In_Output.Init_IO();
envOn = 0;
int reclYear = 0;
Console.Write("Beginning Landis 7.0 Pro Run\n Initializing...\n");
if (envOn > 0)
{
Console.Write("Environment will be updated every {0} iterarion\n", envOn);
}
Landispro.Timestep.timestep = (uint)gl_param.SuccessionTimestep;
numbOfIter = gl_param.Num_Iteration;
//gl_sites.Stocking_x_value = gl_param.Stocking_x_value;
//gl_sites.Stocking_y_value = gl_param.Stocking_y_value;
//gl_sites.Stocking_z_value = gl_param.Stocking_z_value;
gl_sites.CellSize = gl_param.CellSize;
#if !LANDISPRO_ONLY_SUCCESSION
gl_sites.TimeStep_BDA = gl_param.Timestep_BDA;
gl_sites.TimeStep_Fire = gl_param.Timestep_Fire;
gl_sites.TimeStep_Fuel = gl_param.Timestep_Fuel;
gl_sites.TimeStep_Harvest = gl_param.Timestep_Harvest; //Harvest Module
gl_sites.TimeStep_Wind = gl_param.Timestep_Wind;
#endif
for (int x = 0; x < 5; x++)
{
freq[x] = 1;
}
#if !LANDISPRO_ONLY_SUCCESSION
if ((Init_Output.GDLLMode & defines.G_HARVEST) != 0) //Harvest Module
{
freq[5] = 1;
}
if ((Init_Output.GDLLMode & defines.G_BDA) != 0)
{
Console.Write("BDA ");
}
if ((Init_Output.GDLLMode & defines.G_WIND) != 0)
{
Console.Write("Wind ");
}
if ((Init_Output.GDLLMode & defines.G_HARVEST) != 0)
{
Console.Write("Harvest ");
}
if ((Init_Output.GDLLMode & defines.G_FUEL) != 0)
{
Console.Write("Fuel ");
}
if ((Init_Output.GDLLMode & defines.G_FUELMANAGEMENT) != 0)
{
Console.Write("Fuel management ");
}
if ((Init_Output.GDLLMode & defines.G_FIRE) != 0)
{
Console.Write("Fire ");
}
if (Init_Output.GDLLMode != 0)
{
Console.Write("are(is) on\n");
}
#endif
//In_Output.getInput(freq, ageMaps, pPDP, BDANo, wAdfGeoTransform);
SiteVars.Initialize();
In_Output.getInput(freq, ageMaps, pPDP, wAdfGeoTransform);
#if !LANDISPRO_ONLY_SUCCESSION
if ((gDLLMode & DEFINES.G_HARVEST) != 0)
{
Console.WriteLine("Harvest Dll loaded in...");
GlobalFunctions.HarvestPass(sites, speciesAttrs);
sites.Harvest70outputdim();
}
#endif
Console.Write("Finish getting input\n");
In_Output.OutputScenario();
In_Output.initiateOutput_landis70Pro();
snr = gl_sites.numRows;
snc = gl_sites.numColumns;
specAtNum = gl_spe_Attrs.NumAttrs;
gl_sites.GetSeedDispersalProbability(null, gl_param.SeedRainFlag);
gl_sites.GetSpeciesGrowthRates(gl_param.GrowthFlagFile, gl_param.GrowthFlag);
gl_sites.GetSpeciesMortalityRates(gl_param.MortalityFile, gl_param.MortalityFlag);
gl_sites.GetVolumeRead(gl_param.VolumeFile, gl_param.VolumeFlag);
initiateRDofSite_Landis70();
if (reclYear != 0)
{
int local_num = reclYear / gl_sites.SuccessionTimeStep;
reclass3.reclassify(reclYear, ageMaps);
In_Output.putOutput(local_num, local_num, freq);
In_Output.putOutput_Landis70Pro(local_num, local_num, freq);
In_Output.putOutput_AgeDistStat(local_num);
Console.Write("Ending Landispro Succession.\n");
}
else
{
In_Output.putOutput(0, 0, freq);
In_Output.putOutput_Landis70Pro(0, 0, freq);
In_Output.putOutput_AgeDistStat(0);
}
if (gl_param.RandSeed == 0) //random
{
DateTime startTime = new DateTime(1970, 1, 1);
gl_param.RandSeed = (int)Convert.ToUInt32(Math.Abs((DateTime.Now - startTime).TotalSeconds));
}
system1.fseed(gl_param.RandSeed);
Console.WriteLine("gl_param.RandSeed = {0}", gl_param.RandSeed);
if (envOn > gl_param.Num_Iteration)
{
throw new Exception("Invalid year of interpretation for updating environment");
}
fpforTimeBU_name = gl_param.OutputDir + "/Running_Time_Stat.txt";
fpLogFileSEC_name = gl_param.OutputDir + "/SECLog.txt";
var now2 = DateTime.Now;
Console.WriteLine("\nFinish the initilization at {0}", now2);
var ltimeDiff = now2 - now;
Console.Write("it took {0} seconds\n", ltimeDiff);
gl_landUnits.initiateVariableVector(gl_param.Num_Iteration, gl_param.SuccessionTimestep, specAtNum, gl_param.FlagforSECFile);
using (StreamWriter fpforTimeBU = new StreamWriter(fpforTimeBU_name))
{
fpforTimeBU.Write("Initilization took: {0} seconds\n", ltimeDiff);
}
}