Project

General

Profile

Too many steps in routine RungeKutta » Parc_201540ans.cpp

Jeu de données scénario - broida jeremy, 07/02/2014 09:42 AM

 
1
using namespace std;
2

    
3

    
4
#include "/produits/sec/CODES/CLASS/CLASS_V3.0_02_2014/source/include/CLASSHeaders.hxx"
5
#include <sstream>
6
#include <iomanip>
7
string dtoa(double num)
8
{
9
	ostringstream os(ostringstream::out);
10
	os<<setprecision(3)<<num;
11
	return os.str();
12
}
13

    
14
int main(int argc, char** argv)
15
{
16

    
17
	CLASS gCLASS(-2*365.25*3600*24);		//gCLASS = Parc
18
	gCLASS.SetStockManagement(true);
19
	cSecond year = 3600*8766; // 3600*24*365.25;//Year duration in seconds
20
	gCLASS.SetTimeStep(year/12);
21
	
22
	
23
	cout << "DB Definition... \t" << endl;
24
	cout << "\t -decay \t" ;
25
	DataBank<ZAI>* DecayDB = new DataBank<ZAI>(gCLASS.GetLog(), "/produits/sec/CODES/CLASS/CLASS_V3.0_02_2014/DataBase/Decay.idx");		//loi de décroissance radioactive
26
	gCLASS.SetDecayDataBase(DecayDB);
27
	cout << "\t...OK!" << endl;
28
	
29
	cout << "\t -MOX \t" ;
30
	DataBank<IsotopicVector>* MOX = new DataBank<IsotopicVector>(gCLASS.GetLog(), "/produits/sec/CODES/CLASS/CLASS_V3.0_02_2014/DataBase/Data_REP_MOX_IRSN/REP_MOX_BU_40.idx");		//base de données MOX Burn Up 40
31
	cout << "Calculation of Distance Parameters : " << endl;
32
	//    MOX->SetDistanceType(1);//Use coefficient calculated by CLASS to estimate the distance between an Isotopic Vector and an EvolutionData from the MOX DataBank
33
	cout << "\t...OK!" << endl;
34

    
35
	cout << "\t -UOx \t" ;
36
	EvolutionData DB_REP_UOX_3_7 = EvolutionData(gCLASS.GetLog(), "/produits/sec/CODES/CLASS/CLASS_V3.0_02_2014/DataBase/Data_REP_UOX_IRSN/REP_UOX_3_7_U5.dat");		//Base de données UOX //
37
	EvolutionData DB_REP_UOX_4_0 = EvolutionData(gCLASS.GetLog(), "/produits/sec/CODES/CLASS/CLASS_V3.0_02_2014/DataBase/Data_REP_UOX_IRSN/REP_UOX_4_0_U5.dat");		//Base de données avec différents enrichissement
38
	EvolutionData DB_REP_UOX_4_2 = EvolutionData(gCLASS.GetLog(), "/produits/sec/CODES/CLASS/CLASS_V3.0_02_2014/DataBase/Data_REP_UOX_IRSN/REP_UOX_4_2_U5.dat");		//
39
	//EvolutionData DB_REP_UOX = EvolutionData(gCLASS.GetLog(), "/produits/sec/CODES/CLASS/CLASS_V3.0_02_2014/DataBase/Data_REP_UOX_IRSN/REP_UOX_3_1_U5.dat");		
40
	cout << "\t...OK!" << endl;
41
	cout << "DB Definition...\t \t \t...Done!" << endl;
42

    
43
	cout << "Storage Definition... \t"<<endl  ;			//Définition des stocks
44
	Storage *Stock = new Storage(gCLASS.GetLog());
45
	Stock->SetName("Stock Pu REP MOX");
46
	Storage *ReUsable = new Storage(gCLASS.GetLog());
47
	ReUsable->SetName("Dechets");
48
	cout << "\t...OK!" << endl;
49
	
50
	cout <<"Creation d'un stock initial"<< endl;
51

    
52
	double PuQuantiteIsotopique[5] ; 
53
	PuQuantiteIsotopique[0]= 9.49e27 ;                 // quantité de 238 Pu
54
	PuQuantiteIsotopique[1]= 1.65e29 ;                // quantité de 239 Pu
55
	PuQuantiteIsotopique[2]= 1.20e29 ;                // quantité de 240 			
56
	PuQuantiteIsotopique[3]= 4.68e28 ;                // quantité de 241 Pu
57
	PuQuantiteIsotopique[4]= 3.50e28 ;                // quantité de 242 Pu
58

    
59

    
60
	IsotopicVector IVStock;
61
	int IsotopePu = 238 ;
62
	for(int i=0;i<6;i++)         
63
	{
64
	IVStock += ZAI(94,IsotopePu,0) * PuQuantiteIsotopique[i];
65
	IsotopePu = IsotopePu+1;
66
	}
67

    
68

    
69

    
70

    
71
	cout << "Pool Definition... \t" <<endl ;		//Définition des piscines
72
	Pool *Cooling_UOX = new Pool(gCLASS.GetLog(),Stock, gCLASS.GetAbsoluteTime(), (double)year*5 );
73
	Cooling_UOX->SetName("Cooling Pool for UOX");
74
	Pool *Cooling_MOX = new Pool(gCLASS.GetLog(),ReUsable, gCLASS.GetAbsoluteTime(), (double)year*5 );
75
	Cooling_MOX->SetName("Cooling Pool for MOX");
76
	cout << "\t...OK!" << endl;
77

    
78
	cout << "FabricationPlant Definition... \t" <<endl ;		//Usine de fabrication MOX
79
	FabricationPlant *FP_MOX = new FabricationPlant(gCLASS.GetLog(),Stock, ReUsable);
80
	FP_MOX->SetName("Fabrication Plant for MOX");
81
	cout << "\t...OK!" << endl;
82

    
83

    
84
	double UOX_Cycle = year*4;
85
	double MOX_Cycle = year*4;
86
	
87
	cout << "Reactor Definition... \t"<<endl ;				// Définition des réacteurs
88
	
89
	Reactor *Belleville_1;
90
	Reactor *Belleville_2;
91

    
92
	Reactor *Blayais_1_MOX;
93
	Reactor *Blayais_1_UOX;
94
	Reactor *Blayais_2_MOX;
95
	Reactor *Blayais_2_UOX;
96
	Reactor *Blayais_3;
97
	Reactor *Blayais_4;
98

    
99
	Reactor *Bugey_2;
100
	Reactor *Bugey_3;
101
	Reactor *Bugey_4;
102
	Reactor *Bugey_5;
103

    
104
	Reactor *Cattenom_1;
105
	Reactor *Cattenom_2;
106
	Reactor *Cattenom_3;
107
	Reactor *Cattenom_4;
108

    
109
	Reactor *Chinon_1_MOX;
110
	Reactor *Chinon_1_UOX;
111
	Reactor *Chinon_2_MOX;
112
	Reactor *Chinon_2_UOX;
113
	Reactor *Chinon_3_MOX;
114
	Reactor *Chinon_3_UOX;
115
	Reactor *Chinon_4_MOX;
116
	Reactor *Chinon_4_UOX;
117

    
118
	Reactor *Chooz_1;
119
	Reactor *Chooz_2;
120

    
121
	Reactor *Civaux_1;
122
	Reactor *Civaux_2;
123

    
124
	Reactor *Cruas_1;
125
	Reactor *Cruas_2;
126
	Reactor *Cruas_3;
127
	Reactor *Cruas_4;
128

    
129
	Reactor *Dampierre_1_MOX;
130
	Reactor *Dampierre_1_UOX;
131
	Reactor *Dampierre_2_MOX;
132
	Reactor *Dampierre_2_UOX;
133
	Reactor *Dampierre_3_MOX;
134
	Reactor *Dampierre_3_UOX;
135
	Reactor *Dampierre_4_MOX;
136
	Reactor *Dampierre_4_UOX;
137

    
138
	Reactor *Fessenheim_1;
139
	Reactor *Fessenheim_2;
140

    
141
	Reactor *Flamanville_1;
142
	Reactor *Flamanville_2;
143

    
144
	Reactor *Golfech_1;
145
	Reactor *Golfech_2;
146

    
147
	Reactor *Gravelines_1;
148
	Reactor *Gravelines_2_MOX;
149
	Reactor *Gravelines_2_UOX;
150
	Reactor *Gravelines_3_MOX;
151
	Reactor *Gravelines_3_UOX;
152
	Reactor *Gravelines_4_MOX;
153
	Reactor *Gravelines_4_UOX;
154
	Reactor *Gravelines_5_MOX;
155
	Reactor *Gravelines_5_UOX;
156
	Reactor *Gravelines_6;
157

    
158
	Reactor *Nogent_1;
159
	Reactor *Nogent_2;
160

    
161
	Reactor *Paluel_1;
162
	Reactor *Paluel_2;
163
	Reactor *Paluel_3;
164
	Reactor *Paluel_4;
165

    
166
	Reactor *Penly_1;
167
	Reactor *Penly_2;
168

    
169
	Reactor *StAlban_1;
170
	Reactor *StAlban_2;
171

    
172
	Reactor *StLaurent_1_MOX;
173
	Reactor *StLaurent_1_UOX;
174
	Reactor *StLaurent_2_MOX;
175
	Reactor *StLaurent_2_UOX;
176

    
177
	Reactor *Tricastin_1_MOX;
178
	Reactor *Tricastin_1_UOX;
179
	Reactor *Tricastin_2_MOX;
180
	Reactor *Tricastin_2_UOX;
181
	Reactor *Tricastin_3_MOX;
182
	Reactor *Tricastin_3_UOX;
183
	Reactor *Tricastin_4_MOX;
184
	Reactor *Tricastin_4_UOX;
185

    
186

    
187

    
188

    
189

    
190
	
191

    
192
	Belleville_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*13, 3.93e9,108.09,43,0.75) ;			
193
	Belleville_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*14, 3.93e9,108.09,43,0.75) ;
194

    
195
	Blayais_1_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*6, 0.819e9,21.612,35,0.75) ;
196
	Blayais_1_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.911e9,50.428,43,0.75) ;
197
	Blayais_2_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*8, 0.819e9,21.612,35,0.75) ;
198
	Blayais_2_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*8, 1.911e9,50.428,43,0.75) ;
199
	Blayais_3 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*8, 2.73e9,71.84,42,0.75) ;
200
	Blayais_4 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*8, 2.73e9,71.84,42,0.75) ;
201

    
202
	Bugey_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_2, Cooling_UOX,0, double(year)*4, 2.73e9,69.87,42,0.75) ;
203
	Bugey_3 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_2, Cooling_UOX,0, double(year)*4, 2.73e9,69.87,42,0.75) ;
204
	Bugey_4 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_2, Cooling_UOX,0, double(year)*4, 2.64e9,67.65,42,0.75) ;
205
	Bugey_5 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_2, Cooling_UOX,0, double(year)*5, 2.64e9,67.65,42,0.75) ;
206

    
207
	Cattenom_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*12, 3.9e9,107.34,43,0.75) ;
208
	Cattenom_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*13, 3.9e9,107.34,43,0.75) ;
209
	Cattenom_3 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*16, 3.9e9,107.34,43,0.75) ;
210
	Cattenom_4 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*17, 3.9e9,107.34,43,0.75) ;
211

    
212
	Chinon_1_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*9, 0.8145e9,21.504,35,0.75) ;
213
	Chinon_1_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*9, 1.9005e9,50.176,43,0.75) ;
214
	Chinon_2_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*9, 0.8145e9,21.504,35,0.75) ;
215
	Chinon_2_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*9, 1.9005e9,50.176,43,0.75) ;
216
	Chinon_3_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*12, 0.8145e9,21.504,35,0.75) ;
217
	Chinon_3_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*12, 1.9005e9,50.176,43,0.75) ;
218
	Chinon_4_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*13, 0.8145e9,21.504,35,0.75) ;
219
	Chinon_4_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*13, 1.9005e9,50.176,43,0.75) ;
220

    
221
	Chooz_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*25, 4.5e9,107.56,43,0.75) ;
222
	Chooz_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*25, 4.5e9,107.56,43,0.75) ;
223

    
224
	Civaux_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*27, 4.485e9,107.24,43,0.75) ;
225
	Civaux_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*27, 4.485e9,107.24,43,0.75) ;
226

    
227
	Cruas_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*9, 2.745e9,72.32,42,0.75) ;
228
	Cruas_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*10, 2.745e9,72.32,42,0.75) ;
229
	Cruas_3 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*9, 2.745e9,72.32,42,0.75) ;
230
	Cruas_4 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*10, 2.745e9,72.32,42,0.75) ;
231

    
232
	Dampierre_1_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*5, 0.801e9,21.144,35,0.75) ;
233
	Dampierre_1_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*5, 1.869e9,49.336,43,0.75) ;
234
	Dampierre_2_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*5, 0.801e9,21.144,35,0.75) ;
235
	Dampierre_2_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.869e9,49.336,43,0.75) ;
236
	Dampierre_3_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*5, 0.801e9,21.144,35,0.75) ;
237
	Dampierre_3_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.869e9,49.336,43,0.75) ;
238
	Dampierre_4_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*5, 0.801e9,21.144,35,0.75) ;
239
	Dampierre_4_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.869e9,49.336,43,0.75) ;
240

    
241
	Fessenheim_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_2, Cooling_UOX,0, double(year)*3, 2.64e9,67.65,42,0.75) ;
242
	Fessenheim_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_2, Cooling_UOX,0, double(year)*3, 2.64e9,67.65,42,0.75) ;
243

    
244
	Flamanville_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*11, 3.99e9,109.71,43,0.75) ;
245
	Flamanville_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*12, 3.99e9,109.71,43,0.75) ;
246

    
247
	Golfech_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*16, 3.93e9,108.09,43,0.75) ;
248
	Golfech_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*19, 3.93e9,108.09,43,0.75) ;
249

    
250
	Gravelines_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*5, 2.73e9,71.84,42,0.75) ;
251
	Gravelines_2_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*5, 0.819e9,21.612,35,0.75) ;
252
	Gravelines_2_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*5, 1.911e9,50.428,42,0.75) ;
253
	Gravelines_3_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*6, 0.819e9,21.612,35,0.75) ;
254
	Gravelines_3_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.911e9,50.428,42,0.75) ;
255
	Gravelines_4_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*6, 0.819e9,21.612,35,0.75) ;
256
	Gravelines_4_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.911e9,50.428,42,0.75) ;
257
	Gravelines_5_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*10, 0.819e9,21.612,35,0.75) ;
258
	Gravelines_5_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*10, 1.911e9,50.428,42,0.75) ;
259
	Gravelines_6 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*10, 2.73e9,71.84,42,0.75) ;
260

    
261
	Nogent_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*13, 3.93e9,108.09,43,0.75) ;
262
	Nogent_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*14, 3.93e9,108.09,43,0.75) ;
263

    
264
	Paluel_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*10, 3.99e9,109.71,43,0.75) ;
265
	Paluel_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*10, 3.99e9,109.71,43,0.75) ;
266
	Paluel_3 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*11, 3.99e9,109.71,43,0.75) ;
267
	Paluel_4 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*11, 3.99e9,109.71,43,0.75) ;
268

    
269
	Penly_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*15, 3.99e9,109.71,43,0.75) ;
270
	Penly_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*17, 3.99e9,109.71,43,0.75) ;
271

    
272
	StAlban_1 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*11, 4.005e9,110.22,43,0.75) ;
273
	StAlban_2 = new Reactor(gCLASS.GetLog(),DB_REP_UOX_4_0, Cooling_UOX,0, double(year)*12, 4.005e9,110.22,43,0.75) ;
274

    
275
	StLaurent_1_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*8, 0.8235e9,21.756,35,0.75) ;
276
	StLaurent_1_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*8, 1.9215e9,50.764,43,0.75) ;
277
	StLaurent_2_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*8, 0.8235e9,21.756,35,0.75);
278
	StLaurent_2_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*8, 1.9215e9,50.764,43,0.75) ;
279

    
280
	Tricastin_1_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*5, 0.8235e9,21.756,35,0.75) ;
281
	Tricastin_1_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*5, 1.9215e9,50.764,43,0.75) ;
282
	Tricastin_2_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*5, 0.8235e9,21.756,35,0.75) ;
283
	Tricastin_2_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*5, 1.9215e9,50.764,43,0.75) ;
284
	Tricastin_3_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*6, 0.8235e9,21.756,35,0.75) ;
285
	Tricastin_3_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.9215e9,50.764,43,0.75) ;
286
	Tricastin_4_MOX = new Reactor(gCLASS.GetLog(),MOX, FP_MOX, Cooling_MOX,0, double(year)*6, 0.8235e9,21.756,35,0.75) ;
287
	Tricastin_4_UOX = new Reactor(gCLASS.GetLog(),DB_REP_UOX_3_7, Cooling_UOX,0, double(year)*6, 1.9215e9,50.764,43,0.75) ;
288

    
289

    
290

    
291
	cout << "\t...OK!" << endl;
292
	cout << endl << endl;
293

    
294
	cout << "Addind Phase : " << endl;
295
	cout << "\t -FabricationPlant" ;
296
	gCLASS.AddFabricationPlant(FP_MOX);		//Ajout de tous les éléments à notre parc
297
	cout << "\t...OK!" << endl;
298

    
299
	cout << "\t -Pool \t\t" ;
300
	gCLASS.AddPool(Cooling_UOX);
301
	gCLASS.AddPool(Cooling_MOX);
302
	cout << "\t...OK!" << endl;
303

    
304
	cout << "\t -Reactor \t\t" ;
305
		
306

    
307
	gCLASS.AddReactor(Belleville_1);
308
	gCLASS.AddReactor(Belleville_2);
309
	gCLASS.AddReactor(Blayais_1_MOX);
310
	gCLASS.AddReactor(Blayais_1_UOX);
311
	gCLASS.AddReactor(Blayais_2_MOX);
312
	gCLASS.AddReactor(Blayais_2_UOX);
313
	gCLASS.AddReactor(Blayais_3);
314
	gCLASS.AddReactor(Blayais_4);
315
	gCLASS.AddReactor(Bugey_2);
316
	gCLASS.AddReactor(Bugey_3);
317
	gCLASS.AddReactor(Bugey_4);
318
	gCLASS.AddReactor(Bugey_5);
319
	gCLASS.AddReactor(Cattenom_1);
320
	gCLASS.AddReactor(Cattenom_2);
321
	gCLASS.AddReactor(Cattenom_3);
322
	gCLASS.AddReactor(Cattenom_4);
323
	gCLASS.AddReactor(Chinon_1_MOX);
324
	gCLASS.AddReactor(Chinon_1_UOX);
325
	gCLASS.AddReactor(Chinon_2_MOX);
326
	gCLASS.AddReactor(Chinon_2_UOX);
327
	gCLASS.AddReactor(Chinon_3_MOX);
328
	gCLASS.AddReactor(Chinon_3_UOX);
329
	gCLASS.AddReactor(Chinon_4_MOX);
330
	gCLASS.AddReactor(Chinon_4_UOX);
331
	gCLASS.AddReactor(Chooz_1);
332
	gCLASS.AddReactor(Chooz_2);
333
	gCLASS.AddReactor(Civaux_1);
334
	gCLASS.AddReactor(Civaux_2);
335
	gCLASS.AddReactor(Cruas_1);
336
	gCLASS.AddReactor(Cruas_2);
337
	gCLASS.AddReactor(Cruas_3);
338
	gCLASS.AddReactor(Cruas_4);
339
	gCLASS.AddReactor(Dampierre_1_MOX);
340
	gCLASS.AddReactor(Dampierre_1_UOX);
341
	gCLASS.AddReactor(Dampierre_2_MOX);
342
	gCLASS.AddReactor(Dampierre_2_UOX);
343
	gCLASS.AddReactor(Dampierre_3_MOX);
344
	gCLASS.AddReactor(Dampierre_3_UOX);
345
	gCLASS.AddReactor(Dampierre_4_MOX);
346
	gCLASS.AddReactor(Dampierre_4_UOX);
347
	gCLASS.AddReactor(Fessenheim_1);
348
	gCLASS.AddReactor(Fessenheim_2);
349
	gCLASS.AddReactor(Flamanville_1);
350
	gCLASS.AddReactor(Flamanville_2);
351
	gCLASS.AddReactor(Golfech_1);
352
	gCLASS.AddReactor(Golfech_2);
353
	gCLASS.AddReactor(Gravelines_1);
354
	gCLASS.AddReactor(Gravelines_2_MOX);
355
	gCLASS.AddReactor(Gravelines_2_UOX);
356
	gCLASS.AddReactor(Gravelines_3_MOX);
357
	gCLASS.AddReactor(Gravelines_3_UOX);
358
	gCLASS.AddReactor(Gravelines_4_MOX);
359
	gCLASS.AddReactor(Gravelines_4_UOX);
360
	gCLASS.AddReactor(Gravelines_5_MOX);
361
	gCLASS.AddReactor(Gravelines_5_UOX);
362
	gCLASS.AddReactor(Gravelines_6);
363
	gCLASS.AddReactor(Nogent_1);
364
	gCLASS.AddReactor(Nogent_2);
365
	gCLASS.AddReactor(Paluel_1);
366
	gCLASS.AddReactor(Paluel_2);
367
	gCLASS.AddReactor(Paluel_3);
368
	gCLASS.AddReactor(Paluel_4);
369
	gCLASS.AddReactor(Penly_1);
370
	gCLASS.AddReactor(Penly_2);
371
	gCLASS.AddReactor(StAlban_1);
372
	gCLASS.AddReactor(StAlban_2);
373
	gCLASS.AddReactor(StLaurent_1_MOX);
374
	gCLASS.AddReactor(StLaurent_1_UOX);
375
	gCLASS.AddReactor(StLaurent_2_MOX);
376
	gCLASS.AddReactor(StLaurent_2_UOX);
377
	gCLASS.AddReactor(Tricastin_1_MOX);
378
	gCLASS.AddReactor(Tricastin_1_UOX);
379
	gCLASS.AddReactor(Tricastin_2_MOX);
380
	gCLASS.AddReactor(Tricastin_2_UOX);
381
	gCLASS.AddReactor(Tricastin_3_MOX);
382
	gCLASS.AddReactor(Tricastin_3_UOX);
383
	gCLASS.AddReactor(Tricastin_4_MOX);
384
	gCLASS.AddReactor(Tricastin_4_UOX);
385

    
386
	
387
	cout << "...OK!" << endl;
388

    
389
	cout << "\t -Storage \t\t" ;
390
	gCLASS.AddStorage(Stock);
391
	Stock->AddToStock(IVStock);
392
	gCLASS.AddStorage(ReUsable);
393

    
394
	cout << "...OK!" << endl;
395

    
396
	cout << "Addind Phase...\t \t \t...Done!" << endl;
397
	cout << endl << endl;
398
	
399

    
400

    
401
	cout << "Beginning the Evolution" << endl;
402
	gCLASS.Evolution((double)year*50);
403

    
404
}
(1-1/2)