CLASS

/************************************************************/

//              DESCRIPTION

// 

// 

// Simple Scenario with one PWR UOX and one PWR MOX

// 

//                                   

//   _______     ____    _______    ___________

//  |       |   |    |  |       |  |	       |

//  |Reactor| =>|Pool|=>|Storage|=>|Fabrication| 

//  | UOx   |	|UOX |	| UOX 	|  |   Plant   |

//  |_______|   |____|  |_______|  |___________|

//  				      	||

//				      	\/

//               _______     ____     _______

//              |       |   |    |   |       |

//              |Storage|<= |Pool|<==|Reactor|

//              | MOX   |   |MOX |   | MOX   |

//              |_______|   |____|   |_______|

//

//

//@author FrA

/***********************************************************/

#include "CLASSHeaders.hxx"

#include <sstream>

#include <iomanip>

#include <math.h>

#include <string>

#include "XS/XSM_MLP.hxx"			//Load the include for Neural network cross section predictor

#include "Irradiation/IM_RK4.hxx"		//Load the include for Runge Kutta 4 resolution

#include "Equivalence/EQM_PWR_MLP_MOX.hxx"	//Load the include for Neural Network Equivalence Model (PWRMOX)

using namespace std;

int main(int argc, char** argv)

{

  //seconds in one year

  cSecond year = 3600*24.*365.25; 

  /******LOG MANAGEMENT**********************************/

  //Definition of the Log file : CLASS messages output 

  int Std_output_level 	= 0; // Only error are shown in terminal

  int File_output_level = 2; // Error + Warning + Info are shown in the file CLASS_OUTPUT.log

  CLASSLogger *Logger 	= new CLASSLogger("CLASS_OUTPUT.log",Std_output_level,File_output_level);

  /******SCENARIO****************************************/

  // The scenario start at year 2015

  Scenario *gCLASS=new Scenario(2015*year,Logger);

  gCLASS->SetStockManagement(true); //If false all the IsotopicVector in stocks are mixed together.

  gCLASS->SetTimeStep(year/4.); //the scenario calculation is updated every 3 months

  cSecond EndOfScenarioTime=2060*year; //Scenario ends in year 2060

  gCLASS->SetOutputFileName("Scenario6.root");//Set the name of the output file

  /******DATA BASES**************************************/

  //Geting CLASS to path

  string CLASS_PATH = getenv("CLASS_PATH");

  if (CLASS_PATH=="")

    {

      cout<<" Please setenv CLASS_PATH to your CLASS installation folder in your .bashs or .tcshrc"<<endl;

      exit(0);

    }

  string PATH_TO_DATA = CLASS_PATH + "/DATA_BASES/";

  /*===Decay data base===*/

  //The decay data base is taken from the file Decay.idx

  DecayDataBank* DecayDB = new DecayDataBank(gCLASS->GetLog(), PATH_TO_DATA + "DECAY/ALL/Decay.idx"); 

  //you may have to open this file and do the proper changes according your path

  gCLASS->SetDecayDataBase(DecayDB);

  //This decay data base will be used for all the decay calculations in this Scenario

  /*===Reactor data base===*/

  // Reprocessed fuel PWR MOX

  XSM_MLP* XSMOX = new XSM_MLP(gCLASS->GetLog(), PATH_TO_DATA + "PWR/MOX/XSModel/30Wg_FullMOX");

  //Defining the XS Predictor

  IM_RK4 *IMRK4 = new IM_RK4(gCLASS->GetLog());	//Bateman's equation solver method (RungeKutta4)

                                         //NO PATH FOR THE RK SOLVER???

  EQM_PWR_MLP_MOX* EQMMLPPWRMOX = new EQM_PWR_MLP_MOX(gCLASS->GetLog(),PATH_TO_DATA + "PWR/MOX/EQModel/MLP/EQM_MLP_PWR_MOX_3batch.xml");//Defining the EquivalenceModel

  PhysicsModels* PHYMOD = new PhysicsModels(XSMOX, EQMMLPPWRMOX, IMRK4);//The PhysicsModels containing the 3 object previously defined

  //Fixed fuel : PWR UOX

  EvolutionData *CYCLADE =new EvolutionData(gCLASS->GetLog(), PATH_TO_DATA + "PWR/UOX/FixedFuel/CYCLADES.dat");

  EvolutionData *GARANCE = new EvolutionData(gCLASS->GetLog(), PATH_TO_DATA + "PWR/UOX/FixedFuel/GARANCE.dat");

  EvolutionData *STD900 = new EvolutionData(gCLASS->GetLog(), PATH_TO_DATA + "PWR/UOX/FixedFuel/STD900.dat");

  /******FACILITIES*************************************/

  /*=== Stock===*/

  //Storage for infinite U and Pu

  Storage *Stock = new Storage(gCLASS->GetLog());    //Definition of the stock

  Stock->SetName("Stock");   //Its name

  Stock->AddToStock(ZAI(92,238,0) *  3e31);

  Stock->AddToStock(ZAI(92,235,0) * 0.04375* 3e31);

  /*  Stock->AddToStock(ZAI(94,238,0) *  10e28);

  Stock->AddToStock(ZAI(94,239,0) *  10e28);

  Stock->AddToStock(ZAI(94,240,0) *  10e28);

  Stock->AddToStock(ZAI(94,241,0) *  10e28);

  Stock->AddToStock(ZAI(94,242,0) *  10e28);*/

  gCLASS->Add(Stock);	//Adding the stock to the Scenario

  //Storage for MOX

  Storage *StockMOX = new Storage(gCLASS->GetLog());//Definition of the stock

  StockMOX->SetName("StockMOX");       	//Its name

  gCLASS->Add(StockMOX);	      	//Adding the stock to the Scenario

  //Storage for UOX

  Storage *StockUOX = new Storage(gCLASS->GetLog());//Definition of the stock

  StockUOX->SetName("StockUOX");	//Its name

  gCLASS->Add(StockUOX);	       	//Adding the stock to the Scenario

  /*===Pool===*/

  //Pool for UOX

  Pool *Cooling_UOX = new Pool(gCLASS->GetLog(),StockUOX, 5*year); //After 5 years of cooling, the pool sends its content to "StockUOX"

  Cooling_UOX->SetName("Pool_UOX");

  gCLASS->Add(Cooling_UOX);

  //Pool for MOX

  Pool *Cooling_MOX = new Pool(gCLASS->GetLog(),StockMOX, 5*year); //After 5 years of cooling, the pool sends its content to "StockMOX"

  Cooling_MOX->SetName("Pool_MOX");

  gCLASS->Add(Cooling_MOX);

  /*===A FabricationPlant===*/

  FabricationPlant *FP_MOX = new FabricationPlant(gCLASS->GetLog(), 3*year); //Declare a FabricationPlant. After the build of the fuel, it decays during 3years before to be loaded in Reactor

  FP_MOX->SetFiFo(false); //The latest isotopicVector to enter in "Stock" will be used to build the fuel (Opposite of First In First Out)

  FP_MOX->SetName("Fab_MOX");

  FP_MOX->AddFissileStorage(StockUOX);	//Tell the FP to look in Stock for fissionable material 

  FP_MOX->AddFertileStorage(Stock);//Tell the FP to look in Stock for fertile material 

  //If fertile stock is not defined (like here), CLASS get fertile from nature (OUTCOMING vector)

  //FP_MOX->SetReUsableStorage(wastestock);//By default the fabricationplant get the list of nuclei defined in the EquivalenceModel (here EQM_MLP_MOX) from stock and send the others nuclei in WASTE. If user want these nuclei to go in another stock  he can use the SetReUsableStorage function

  gCLASS->AddFabricationPlant(FP_MOX);

  /*=== Reactor===*/

   double	Power_CP0 = 2.66e10;	    //Thermal power (in W)

  //Combustibles type : CYCLADE

  double  BurnUp_Cyclade    = 47; 		// GWd/tHM

  double  HMMass_Cyclade    = 72.3;		//heavy metal mass (in tons)

  //REACTOR PWR UOX

  cSecond StartingTime =  2016*year;

  cSecond LifeTime     =  (EndOfScenarioTime - StartingTime);

  Reactor* PWR_UOX = new Reactor(gCLASS->GetLog(),// Log

   CYCLADE,	// The DataBase used

   Cooling_UOX,	// Connected Backend facility : The reactor discharge its fuel into the Pool "Cooling_UOX"

   StartingTime,// Starting time

   LifeTime,	// time of reactor life time

   Power_CP0,	// Power

   HMMass_Cyclade,// HM mass

   BurnUp_Cyclade,// BurnUp

   0.8);	// Load Factor

  PWR_UOX->SetName("PWR_UOX");// name of the reactor (as it will show up in the CLASSGui)

  gCLASS->AddReactor(PWR_UOX);//Add this reactor to the scenario

  //POWER reactor MOX

  double Power_CPY = 2.785e9;	    //Thermal power (in W)

  //Combustibles type : STANDARD 900

  double  BurnUpMOX = 35;             //GWd/tHM

  //double  BurnUp_STD900    = 33;      // GWd/tHM

  double  HMMass_STD900    = 72.3;    //heavy metal mass (in tons)

  //the PWR MOX

  StartingTime = 2022*year;

  LifeTime     =  EndOfScenarioTime - StartingTime;

  Reactor* PWR_MOX = new Reactor(gCLASS->GetLog(),// Log

   PHYMOD,			// The models used to build the fuel & to calculate its evolution

   FP_MOX,			// The FabricationPlant

   Cooling_MOX,			// Connected Backend

   StartingTime,		// Starting time

   LifeTime,			// time of reactor l

   Power_CPY,			// Power

   HMMass_STD900,		// HM mass

   BurnUpMOX,			// BurnUp

   0.8);			// Load Factor

  PWR_MOX->SetName("PWR_MOX");// name of the reactor (as it will appear in the CLASSGui)

  gCLASS->AddReactor(PWR_MOX);//Add this reactor to the scenario

  gCLASS->Evolution((double)EndOfScenarioTime);//Perform the calculation from year 2015(defined in Scenario declaration) to year 2060

	delete gCLASS;

}

//==========================================================================================

// Compilation

//==========================================================================================

/*

 \rm CLASS* ; g++ -o CLASS_Exec Scenario6.cxx -I $CLASS_include -L $CLASS_lib -lCLASSpkg `root-config --cflags` `root-config --libs` -fopenmp -lgomp -Wunused-result

 */