public void test_RefinementGenerator_write_the_bws_file_in_the_proper_Directory_aritm()
{
WSProject project = new WSProject(project_file.FullName);
string bws_filename = "test_bws_position";
RefinementGenerator.AritmeticalGrading test_grading = new RefinementGenerator.AritmeticalGrading(1.2,12,400.0);
test_grading.Rectangle(project, bws_filename, 50, 50, 50, 50);
string expected_file_path = project.file.Directory.FullName + "\\dtm\\" + "test_bws_position" + ".bws"; ;
FileInfo expected_file = new FileInfo(expected_file_path);
Assert.IsTrue(expected_file.Exists);
//
}
public void test_RefinementGenerator_AritmeticalGrading_heightFirstLayer()
{
RefinementGenerator.AritmeticalGrading test_grading = new RefinementGenerator.AritmeticalGrading(1.2, 10, 10);
Assert.AreEqual(Math.Round(0.5263157894737, 4), Math.Round(test_grading.heightFirstLayer, 4));
}
public void test_RefinementGenerator_AritmeticalGrading_additiveCostant()
{
RefinementGenerator.AritmeticalGrading test_grading = new RefinementGenerator.AritmeticalGrading(1.2, 10, 10);
Assert.AreEqual(Math.Round(0.1052633157894737, 4), Math.Round(test_grading.additiveCostant, 4));
}
public void test_RefinementGenerator_AritmeticalGrading_heightDistributionFactor()
{
RefinementGenerator.AritmeticalGrading test_grading = new RefinementGenerator.AritmeticalGrading(1.2,10,10);
Assert.AreEqual(Math.Round(0.3571428571429, 4), Math.Round(test_grading.heightDistributionFactor, 4));
}
public ExcelBatch(WSProject referenceProject, FileInfo excelBatchFile)
{
// controllo l-esistenza del file excel
if (verifyExcelFile(excelBatchFile) && verifyProjectFiles(excelBatchFile, referenceProject))
{
// setto la directory
// riempo l'experimentBatch
Workbook book = Workbook.Load(excelBatchFile.FullName);
Worksheet sheet = book.Worksheets[0];
string batteryName = sheet.Cells[0, 1].StringValue;
double numberOfEsperiments = Convert.ToInt32(sheet.Cells[1, 1].Value);
Worksheet results_geom = book.Worksheets[1];
results_geom.Cells[0, 0] = new Cell("UCRT Profiles for Geometrical Grading");
Worksheet results_ke_geom = book.Worksheets[2];
results_ke_geom.Cells[0, 0] = new Cell("KE Profiles for Geometrical Grading");
Worksheet results_zcen_geom = book.Worksheets[3];
results_zcen_geom.Cells[0, 0] = new Cell("Z Profiles for Geometrical Grading");
Worksheet results_arit = book.Worksheets[4];
results_arit.Cells[0, 0] = new Cell("UCRT Profiles for Aritmetical Grading");
Worksheet results_ke_arit = book.Worksheets[5];
results_ke_arit.Cells[0, 0] = new Cell("KE Profiles for Aritmetical Grading");
Worksheet results_zcen_arit = book.Worksheets[6];
results_zcen_arit.Cells[0, 0] = new Cell("Z Profiles for Aritmetical Grading");
Worksheet results_tem1_arit = book.Worksheets[7];
results_tem1_arit.Cells[0, 0] = new Cell("TEM1 Profiles for Aritmetical Grading");
Worksheet results_tem1_geom = book.Worksheets[8];
results_tem1_geom.Cells[0, 0] = new Cell("TEM1 Profiles for Geometrical Grading");
Console.WriteLine("Total simulations in the excel :"+ numberOfEsperiments);
for (int i = 0; i < numberOfEsperiments; i++)
{
int todo = Convert.ToInt32(sheet.Cells[6, i + 1].Value);
if (todo == 0)
{
Console.WriteLine("Skipping case : " + i);
}
if( todo == 1 || todo == 2)
{
Console.WriteLine("Running case : " + i + " todo " + todo);
WSProject testcase = null;
WSProject testcase_arit = null;
bool restart = false;
string testCasetargetDirectory = excelBatchFile.Directory.FullName + "\\" + batteryName + "_" + sheet.Cells[7, i + 1].StringValue;
double z0_convergence_treshold = Convert.ToDouble(sheet.Cells[14, i + 1].Value);
int z0_convergence_sweeps = Convert.ToInt32(sheet.Cells[15, i + 1].Value);
int z0_cycles_to_be_checked = Convert.ToInt32(sheet.Cells[16, i + 1].Value);
int z0_monitoring_x = Convert.ToInt32(sheet.Cells[17, i + 1].Value);
int z0_monioring_y = Convert.ToInt32(sheet.Cells[18, i + 1].Value);
int nx = Convert.ToInt32(sheet.Cells[8, i + 1].Value);
int nj = Convert.ToInt32(sheet.Cells[9, i + 1].Value);
int verticalResolution = Convert.ToInt32(sheet.Cells[10, i + 1].Value);
float totalHeight = (float)Convert.ToDouble(sheet.Cells[11, i + 1].Value);
float roughness = (float)Convert.ToDouble(sheet.Cells[12, i + 1].Value);
double heithFirstLayerGeometrical = Convert.ToDouble(sheet.Cells[21, i + 1].Value);
RefinementGenerator.GeometricalGrading geom_refinement = new RefinementGenerator.GeometricalGrading(heithFirstLayerGeometrical, totalHeight, verticalResolution);
if (todo == 1) {
testcase = new WSProject(referenceProject.file.FullName, testCasetargetDirectory+"_geom");
testcase.load_parameters_from_excel(sheet);
#region read batch case parameters
FileInfo gwsfile = new FileInfo(excelBatchFile.Directory.FullName + "\\gws_files\\" + sheet.Cells[12, i + 1].ToString());
testcase.parameters.CFD.Height = totalHeight;
testcase.parameters.DTM.Roughness = roughness;
testcase.save();
testcase.load_gws(gwsfile.FullName);
// run the case of geometrical grid..
geom_refinement.Rectangle(testcase,"autogenerated_geometrical_bws_"+nx+"_"+nj, testcase.parameters.DTM.XMax,testcase.parameters.DTM.YMax,nx,nj);
#endregion
testcase.run(1);
}
else if (todo == 2)
{
restart = true;
z0_convergence_treshold = Convert.ToDouble(sheet.Cells[37, i + 1].Value);
z0_convergence_sweeps = Convert.ToInt32(sheet.Cells[38, i + 1].Value);
z0_cycles_to_be_checked = Convert.ToInt32(sheet.Cells[39, i + 1].Value);
FileInfo geometrico = new FileInfo(testCasetargetDirectory + "_geom\\" + referenceProject.file.Name);
testcase = new WSProject(geometrico.FullName);
}
Console.WriteLine("...G");
double[] z0results_geom = testcase.run_windfield_z0_conv(z0_convergence_treshold,z0_convergence_sweeps,z0_cycles_to_be_checked,z0_monitoring_x,z0_monioring_y, restart);
int resultRowShift = 0;
if (todo == 2) { resultRowShift = 21;}
#region write results geometrical
// write results
sheet.Cells[20 + resultRowShift, i + 1] = new Cell(geom_refinement.expansionFactor);
//height first cell
sheet.Cells[21 + resultRowShift, i + 1] = new Cell(geom_refinement.heightFirstLayer);
//z0
//// { z0 , (u*/k) , sigma, r2 , sigmaz0 }
sheet.Cells[22 + resultRowShift, i + 1] = new Cell(z0results_geom[0]);
//sigma z0
sheet.Cells[23 + resultRowShift, i + 1] = new Cell(z0results_geom[4]);
//r2
sheet.Cells[24 + resultRowShift, i + 1] = new Cell(z0results_geom[3]);
// add ustar
double ustar_geom = MyMath.ustar_neutral(testcase.parameters.WindField.VelocityBoundaryLayer, testcase.parameters.WindField.HeightBoundaryLayer, testcase.parameters.DTM.Roughness, 0.4);
if (testcase.parameters.WindField.Temperature == 1)
{
ustar_geom = MyMath.ustar_neutral(testcase.parameters.WindField.WindspeedInReferenceHeight, testcase.parameters.WindField.ReferenceHeight, testcase.parameters.DTM.Roughness, 0.4);
}
sheet.Cells[25 + resultRowShift, i + 1] = new Cell(ustar_geom);
//add the project to the dictionary
book.Save(excelBatchFile.FullName);
#endregion
#region write profiles geometrical
// write the profile in the profiles file
double[] profile_geom = testcase.WField[270].ucrt(z0_monitoring_x, z0_monioring_y);
double[] zcen_geom = testcase.WField[270].zcen(z0_monitoring_x, z0_monioring_y);
double[] profile_geom_ke = testcase.WField[270].ke(z0_monitoring_x, z0_monioring_y);
double[] theoretical_mo_neural_geom = MyMath.mo_neutral_u(zcen_geom, 0.4, ustar_geom, roughness);
double rmse_geom = MyMath.rmse(profile_geom, theoretical_mo_neural_geom);
sheet.Cells[26 + resultRowShift, i + 1] = new Cell(rmse_geom);
results_geom.Cells[1, i] = new Cell(i + 1);
results_zcen_geom.Cells[1, i] = new Cell(i + 1);
results_ke_geom.Cells[1, i] = new Cell(i + 1);
for (int count = 0; count < profile_geom.Length; count++)
{
results_geom.Cells[count+3, i] = new Cell(profile_geom[count]);
results_zcen_geom.Cells[count+3, i] = new Cell(zcen_geom[count]);
results_ke_geom.Cells[count + 3, i] = new Cell(profile_geom_ke[count]);
}
if (testcase.parameters.WindField.Temperature == 1)
{
double[] profile_geom_tem1 = testcase.WField[270].tem1(z0_monitoring_x, z0_monioring_y);
for (int count = 0; count < profile_geom.Length; count++)
{
results_tem1_geom.Cells[count + 3, i] = new Cell(profile_geom_tem1[count]);
}
}
book.Save(excelBatchFile.FullName);
#endregion
Console.WriteLine("...A");
RefinementGenerator.AritmeticalGrading arit_refinement = new RefinementGenerator.AritmeticalGrading(geom_refinement.expansionFactor, geom_refinement.numbersOfVericalLayer, geom_refinement.totalHeight);
if (todo == 1)
{
// duplicate the project
testcase_arit = new WSProject(testcase.file.FullName, testCasetargetDirectory + "_arit");
// create the bws
// assign the bws
arit_refinement.Rectangle(testcase_arit, "autogenerated_aritmeical_bws_" + nx + "_" + nj, testcase_arit.parameters.DTM.XMax, testcase_arit.parameters.DTM.YMax, nx, nj);
// run terrain
testcase_arit.run(1);
// run the project
}
else if (todo == 2)
{
z0_convergence_treshold = Convert.ToDouble(sheet.Cells[37, i + 1].Value);
z0_convergence_sweeps = Convert.ToInt32(sheet.Cells[38, i + 1].Value);
z0_cycles_to_be_checked = Convert.ToInt32(sheet.Cells[39, i + 1].Value);
FileInfo aritmetico = new FileInfo(testCasetargetDirectory + "_arit\\" + referenceProject.file.Name);
testcase_arit = new WSProject(aritmetico.FullName);
}
double[] z0results_arit = testcase_arit.run_windfield_z0_conv(z0_convergence_treshold, z0_convergence_sweeps, z0_cycles_to_be_checked, z0_monitoring_x, z0_monioring_y,restart);
#region write resuls aritmetical
// write he results z0
sheet.Cells[28 + resultRowShift, i + 1] = new Cell(arit_refinement.heightDistributionFactor);
//height first cell
sheet.Cells[29 + resultRowShift, i + 1] = new Cell(arit_refinement.heightFirstLayer);
//z0
//// { z0 , (u*/k) , sigma, r2 , sigmaz0 }
sheet.Cells[30 + resultRowShift, i + 1] = new Cell(z0results_arit[0]);
//sigma z0
sheet.Cells[31 + resultRowShift, i + 1] = new Cell(z0results_arit[4]);
//r2
sheet.Cells[32 + resultRowShift, i + 1] = new Cell(z0results_arit[3]);
// add ustar
double ustar_arit = MyMath.ustar_neutral(testcase_arit.parameters.WindField.VelocityBoundaryLayer, testcase_arit.parameters.WindField.HeightBoundaryLayer, testcase_arit.parameters.DTM.Roughness, 0.4);
if (testcase.parameters.WindField.Temperature == 1)
{
ustar_arit = MyMath.ustar_neutral(testcase_arit.parameters.WindField.WindspeedInReferenceHeight, testcase_arit.parameters.WindField.ReferenceHeight, testcase_arit.parameters.DTM.Roughness, 0.4);
}
sheet.Cells[25 + resultRowShift, i + 1] = new Cell(ustar_geom);
sheet.Cells[33 + resultRowShift, i + 1] = new Cell(ustar_arit);
//add the project to the dictionary
book.Save(excelBatchFile.FullName);
#endregion
#region write profiles aritmetical
// write the profile in the profiles file
double[] profile_arit = testcase_arit.WField[270].ucrt(z0_monitoring_x, z0_monioring_y);
double[] zcen_arit = testcase_arit.WField[270].zcen(z0_monitoring_x, z0_monioring_y);
double[] profile_arit_ke = testcase_arit.WField[270].ke(z0_monitoring_x, z0_monioring_y);
// add rmse
double[] theoretical_mo_neural_arit = MyMath.mo_neutral_u(zcen_arit, 0.4, ustar_arit, roughness);
double rmse_arit = MyMath.rmse(profile_arit, theoretical_mo_neural_arit);
sheet.Cells[34 + resultRowShift, i + 1] = new Cell(rmse_arit);
results_arit.Cells[1, i] = new Cell(i + 1);
results_zcen_arit.Cells[1, i] = new Cell(i + 1);
results_ke_arit.Cells[1, i] = new Cell(i + 1);
for (int count = 0; count < profile_arit.Length; count++)
{
results_arit.Cells[count+3, i] = new Cell(profile_arit[count]);
results_zcen_arit.Cells[count+3, i] = new Cell(zcen_arit[count]);
results_ke_arit.Cells[count+3, i] = new Cell(profile_arit_ke[count]);
}
if (testcase.parameters.WindField.Temperature == 1)
{
double[] profile_arit_tem1 = testcase_arit.WField[270].tem1(z0_monitoring_x, z0_monioring_y);
for (int count = 0; count < profile_arit.Length; count++)
{
results_tem1_arit.Cells[count + 3, i] = new Cell(profile_arit_tem1[count]);
}
}
book.Save(excelBatchFile.FullName);
//add the project to the dictionary
#endregion
}
}
}
else
{
throw new System.ArgumentException("Some error occurred regarding the excel file");
}
}