CLASS

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

//              DESCRIPTION

// Simple scenario :

// This park is constituted by a simple PWR UOX

// Reactor and a Storage.

//   _______     _______    _______

//  |       |   |       |  |       |

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

//  |_______|   |_______|  |_______|

//

// The spent fuel goes to a Storage

//

//@author BaL

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

#include "CLASSHeaders.hxx"

#include <sstream>

#include <iomanip>

#include <math.h>

#include <string>

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*15/512);  //1024 timesteps (max number of columns in odc file)

	gCLASS->SetOutputFileName("SimpleReactorMT.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");

	gCLASS->SetDecayDataBase(DecayDB);//This decay data base will be used for all the decay calculations in this Scenario

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

	//The file STD900.dat correspond to a fuel evolution of a UOX PWR (see manual for details)

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

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

	/*===A Stock===*/

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

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

	//Stock->AddToStock(ZAI(92,238,0) * 10);// Just to illustrate, here we add ten nuclei of 238U in this stock

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

	/*===Pool===*/

	//Pool for UOX

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

	Cooling_UOX->SetName("Pool_UOX");

	gCLASS->Add(Cooling_UOX);

	/*===A Reactor : PWR_UOX===*/

	double  HMMass = 72.5;//heavy metal mass (in tons)

	double	Power_CP0 = 2660e6;//Thermal power (in W)

	double  BurnUp = 33; //33 GWd/tHM

	cSecond StartingTime =  2015*year;

	cSecond LifeTime     = 30*year;

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

							   STD900,				// Data base

							   Cooling_UOX,			// Connected Backend facility (here the Storage "Stock" previously declared)

							   StartingTime,		// Starting time

							   LifeTime,			// time of reactor life time

							   Power_CP0,			// Power

							   HMMass,// HM mass

							   BurnUp,				// BurnUp

							   0.8);			// Load Factor

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

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

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

	delete gCLASS;

}

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

// Compilation

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

/*

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

 */