From ff47bcc40eaf79909d83645b6aea033ece861781 Mon Sep 17 00:00:00 2001
From: Baptiste Mouginot <mouginot.baptiste@gmail.com>
Date: Tue, 27 Jan 2015 16:28:22 +0000
Subject: [PATCH] Adding Variable Equivalence model for UOX (Use for REP 900 GW
 with 1/3 fuel loading pattern)

git-svn-id: svn+ssh://svn.in2p3.fr/class@484 0e7d625b-0364-4367-a6be-d5be4a48d228
---
 .../Model/Equivalence/EQM_POL_PWR_UO2.cxx     | 34 +++++++++++++++++++
 .../Model/Equivalence/EQM_POL_PWR_UO2.hxx     | 26 ++++++++++++++
 source/trunk/src/Makefile                     |  2 +-
 3 files changed, 61 insertions(+), 1 deletion(-)
 create mode 100644 source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.cxx
 create mode 100644 source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.hxx

diff --git a/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.cxx b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.cxx
new file mode 100644
index 000000000..ebf2dd1f7
--- /dev/null
+++ b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.cxx
@@ -0,0 +1,34 @@
+#include "EquivalenceModel.hxx" 
+#include "EQM_POL_PWR_UO2.hxx" 
+#include "CLASSLogger.hxx"
+/*Whatever include you need */
+// ________________________________________________________________________ 
+// EQM_POL_PWR_UO2
+//
+// Brief description
+// ________________________________________________________________________ 
+//Constructor(s)
+EQM_POL_PWR_UO2::EQM_POL_PWR_UO2 ( /*parameters*/ ) 
+{
+//.....Do whatever you want with your parameters 
+/*
+Fill the two isotopic vectors fFissileList and fFertileList
+see explanation in the manual */
+// Fertile
+ZAI U8(92 ,238 ,0) ;
+fFertileList = U8*1;
+// Fissile
+ZAI U5(92 ,235 ,0) ;
+// ...
+fFissileList = U5*1;
+}
+// _______________________________________________________________________
+double EQM_POL_PWR_UO2::GetFissileMolarFraction ( IsotopicVector Fissil , IsotopicVector Fertil , double BurnUp )
+{
+/*Code your Equivalence Model : This function has to return the molar fraction of fissile in 
+the fuel needed to reach the BurnUp(GWd/tHM) according to
+the composition of the Fissil and Fertil vectors
+*/
+return 0.0125575 + 0.00053602*BurnUp + 1.80023e-06*BurnUp*BurnUp ;
+
+}
\ No newline at end of file
diff --git a/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.hxx b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.hxx
new file mode 100644
index 000000000..90add1cbd
--- /dev/null
+++ b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.hxx
@@ -0,0 +1,26 @@
+#ifndef _EQM_POL_PWR_UO2_HXX
+#define _EQM_POL_PWR_UO2_HXX
+#include "EquivalenceModel.hxx" 
+using namespace std;
+//−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−//
+/*!
+Define a EQM_POL_PWR_UO2
+Explain briefly what is it . @author YourName
+@version 3.0
+*/
+//________________________________________________________________________
+class EQM_POL_PWR_UO2 : public EquivalenceModel
+{
+  public:
+/*Constructor*/
+EQM_POL_PWR_UO2(/*parameters*/ ); //!< Explain what is the parameters (if any)
+
+/**This function IS the equivalence model**/
+double GetFissileMolarFraction(IsotopicVector Fissil, IsotopicVector Fertil,double BurnUp) ; // !<Return the molar fraction of fissile element
+
+  private:
+/*Your private variables */ 
+
+};
+
+#endif
\ No newline at end of file
diff --git a/source/trunk/src/Makefile b/source/trunk/src/Makefile
index 8c2be94fe..68387ed5a 100755
--- a/source/trunk/src/Makefile
+++ b/source/trunk/src/Makefile
@@ -40,7 +40,7 @@ OBJS		= 	CLASSLogger.o \
 			CLASSFuelPlan.o\
 			Scenario.o
 
-OBJMODEL	=	$(EQM)/EQM_MLP_PWR_MOX.o $(EQM)/EQM_QUAD_PWR_MOX.o $(EQM)/EQM_LIN_PWR_MOX.o \
+OBJMODEL	=	$(EQM)/EQM_MLP_PWR_MOX.o $(EQM)/EQM_QUAD_PWR_MOX.o $(EQM)/EQM_LIN_PWR_MOX.o $(EQM)/EQM_POL_PWR_UO2.o \
 			$(XSM)/XSM_MLP.o $(XSM)/XSM_CLOSEST.o \
 			$(IM)/IM_RK4.o $(IM)/IM_Matrix.o
 
-- 
GitLab