From 07efb4b866e67a8500c8a5f9949d9ff7e70fd84b Mon Sep 17 00:00:00 2001
From: michael <michael@localhost.localdomain>
Date: Fri, 10 Dec 2021 14:43:20 +0100
Subject: [PATCH] Mesh should be set befor giving input fields

---
 .../C/DiffusionEquation_1DHeatedRod.cxx       | 17 +++++++--------
 .../C/DiffusionEquation_1DHeatedRod_FE.cxx    | 21 +++++++++----------
 .../DiffusionEquation_1DHeatedRod_FE_MPI.cxx  | 17 +++++++--------
 .../DiffusionEquation_1DHeatedRod_FV_MPI.cxx  | 15 +++++++------
 4 files changed, 33 insertions(+), 37 deletions(-)

diff --git a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod.cxx b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod.cxx
index 0c275b2..9f181df 100755
--- a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod.cxx
+++ b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod.cxx
@@ -38,7 +38,13 @@ int main(int argc, char** argv)
  
     bool FEcalculation=false;
 	DiffusionEquation  myProblem(spaceDim,FEcalculation,rho_ur,cp_ur,lambda_ur);
-	Field VV("Solid temperature", CELLS, M, 1);
+
+	//Setting initial field
+	Vector VV_Constant(1);
+	VV_Constant(0) = 623;//Rod clad temperature
+
+	cout << "Building initial data" << endl;
+	myProblem.setInitialFieldConstant(M,VV_Constant);
 
 	//Set fluid temperature (temperature du fluide)
 	double fluidTemp=573;//fluid mean temperature
@@ -51,13 +57,6 @@ int main(int argc, char** argv)
 	myProblem.setHeatPowerField(Phi);
 	Phi.writeVTK("1DheatPowerField");
 
-	//Initial field creation
-	Vector VV_Constant(1);
-	VV_Constant(0) = 623;//Rod clad temperature
-
-	cout << "Building initial data" << endl;
-	myProblem.setInitialFieldConstant(M,VV_Constant);
-
 	//set the boundary conditions
 	myProblem.setNeumannBoundaryCondition("Neumann");
 
@@ -94,7 +93,7 @@ int main(int argc, char** argv)
 	else
 		cout << "Simulation "<<fileName<<"  failed ! " << endl;
 
-	cout << "------------ End of calculation -----------" << endl;
+	cout << "------------ End of simulation -----------" << endl;
 	myProblem.terminate();
 
 	return EXIT_SUCCESS;
diff --git a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE.cxx b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE.cxx
index ddb8a4c..7e258ba 100755
--- a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE.cxx
+++ b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE.cxx
@@ -38,7 +38,13 @@ int main(int argc, char** argv)
  
     bool FEcalculation=true;
 	DiffusionEquation  myProblem(spaceDim,FEcalculation,rho_ur,cp_ur,lambda_ur);
-	Field VV("Solid temperature", NODES, M, 1);
+
+	//Set initial field
+	Vector VV_Constant(1);
+	VV_Constant(0) = 623;//Rod clad temperature
+
+	cout << "Building initial data" << endl;
+	myProblem.setInitialFieldConstant(M,VV_Constant,NODES);
 
 	//Set fluid temperature (temperature du fluide)
 	double fluidTemp=573;//fluid mean temperature
@@ -51,13 +57,6 @@ int main(int argc, char** argv)
 	myProblem.setHeatPowerField(Phi);
 	Phi.writeVTK("1DheatPowerField");
 
-	//Initial field creation
-	Vector VV_Constant(1);
-	VV_Constant(0) = 623;//Rod clad temperature
-
-	cout << "Building initial data" << endl;
-	myProblem.setInitialFieldConstant(M,VV_Constant,NODES);
-
 	//set the boundary conditions
 	myProblem.setNeumannBoundaryCondition("Neumann");
 
@@ -90,11 +89,11 @@ int main(int argc, char** argv)
 	myProblem.initialize();
 	bool ok = myProblem.run();
 	if (ok)
-		cout << "Simulation "<<fileName<<" is successful !" << endl;
+		cout << "Simulation of "<<fileName<<" is successful !" << endl;
 	else
-		cout << "Simulation "<<fileName<<"  failed ! " << endl;
+		cout << "Simulation of "<<fileName<<"  failed ! " << endl;
 
-	cout << "------------ End of calculation -----------" << endl;
+	cout << "------------ End of simulation -----------" << endl;
 	myProblem.terminate();
 
 	return EXIT_SUCCESS;
diff --git a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE_MPI.cxx b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE_MPI.cxx
index ddb8a4c..61e88b8 100755
--- a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE_MPI.cxx
+++ b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FE_MPI.cxx
@@ -38,7 +38,13 @@ int main(int argc, char** argv)
  
     bool FEcalculation=true;
 	DiffusionEquation  myProblem(spaceDim,FEcalculation,rho_ur,cp_ur,lambda_ur);
-	Field VV("Solid temperature", NODES, M, 1);
+
+	//Set initial field 
+	Vector VV_Constant(1);
+	VV_Constant(0) = 623;//Rod clad temperature
+
+	cout << "Building initial data" << endl;
+	myProblem.setInitialFieldConstant(M,VV_Constant,NODES);
 
 	//Set fluid temperature (temperature du fluide)
 	double fluidTemp=573;//fluid mean temperature
@@ -51,13 +57,6 @@ int main(int argc, char** argv)
 	myProblem.setHeatPowerField(Phi);
 	Phi.writeVTK("1DheatPowerField");
 
-	//Initial field creation
-	Vector VV_Constant(1);
-	VV_Constant(0) = 623;//Rod clad temperature
-
-	cout << "Building initial data" << endl;
-	myProblem.setInitialFieldConstant(M,VV_Constant,NODES);
-
 	//set the boundary conditions
 	myProblem.setNeumannBoundaryCondition("Neumann");
 
@@ -94,7 +93,7 @@ int main(int argc, char** argv)
 	else
 		cout << "Simulation "<<fileName<<"  failed ! " << endl;
 
-	cout << "------------ End of calculation -----------" << endl;
+	cout << "------------ End of simulation -----------" << endl;
 	myProblem.terminate();
 
 	return EXIT_SUCCESS;
diff --git a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FV_MPI.cxx b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FV_MPI.cxx
index 0c275b2..4de4255 100755
--- a/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FV_MPI.cxx
+++ b/CoreFlows/examples/C/DiffusionEquation_1DHeatedRod_FV_MPI.cxx
@@ -38,7 +38,13 @@ int main(int argc, char** argv)
  
     bool FEcalculation=false;
 	DiffusionEquation  myProblem(spaceDim,FEcalculation,rho_ur,cp_ur,lambda_ur);
-	Field VV("Solid temperature", CELLS, M, 1);
+
+	// Set initial field
+	Vector VV_Constant(1);
+	VV_Constant(0) = 623;//Rod clad temperature
+
+	cout << "Building initial data" << endl;
+	myProblem.setInitialFieldConstant(M,VV_Constant);
 
 	//Set fluid temperature (temperature du fluide)
 	double fluidTemp=573;//fluid mean temperature
@@ -51,13 +57,6 @@ int main(int argc, char** argv)
 	myProblem.setHeatPowerField(Phi);
 	Phi.writeVTK("1DheatPowerField");
 
-	//Initial field creation
-	Vector VV_Constant(1);
-	VV_Constant(0) = 623;//Rod clad temperature
-
-	cout << "Building initial data" << endl;
-	myProblem.setInitialFieldConstant(M,VV_Constant);
-
 	//set the boundary conditions
 	myProblem.setNeumannBoundaryCondition("Neumann");
 
-- 
2.39.2