diff options
Diffstat (limited to 'Dragon/data/rep900_mco_proc')
20 files changed, 5012 insertions, 0 deletions
diff --git a/Dragon/data/rep900_mco_proc/Autop1.c2m b/Dragon/data/rep900_mco_proc/Autop1.c2m new file mode 100644 index 0000000..3571273 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Autop1.c2m @@ -0,0 +1,137 @@ +*DECK Autop1 +*---- +* Name : Autop1.c2m +* Type : DRAGON procedure +* Use : Realize the first self shielding calculation +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBRARY2 := Autop1 LIBRARY TRACKSS :: <<Fuel>> <<NbRegions>> ; +* +* Input data : +* LIBRARY : Library with SUBG key word +* TRACKSS : Tracking of the SS geometry +* Fuel : Type of fuel (UOX/MOX) +* NbRegions : Number of regions for the calculation (20/32) +* Output data : +* LIBRARY2 : Self shielded Library + + +PARAMETER LIBRARY2 LIBRARY TRACKSS :: + EDIT 1 + ::: LINKED_LIST LIBRARY2 LIBRARY TRACKSS ; + ; +* ; +MODULE USS: END: ; + +STRING Fuel ; + :: >>Fuel<< ; +INTEGER NbRegions ; + :: >>NbRegions<< ; + +IF Fuel "MOX" = THEN +LIBRARY2 := USS: LIBRARY TRACKSS :: + EDIT 0 + ARM + TRAN + PASS 2 + GRMIN 31 GRMAX 93 + CALC + REGI W1 U235 ALL + REGI W1 U236 ALL + REGI W1 Pu239 ALL + REGI W1 Pu240 ALL + REGI W1 Pu241 ALL + REGI W1 Pu242 ALL + REGI W1 Zr90g ALL + REGI W1 Zr91g ALL + REGI W1 Zr92g ALL + REGI W1 Zr94g ALL + REGI W1 Zr96g ALL + REGI W1 Zr90f ALL + REGI W1 Zr91f ALL + REGI W1 Zr92f ALL + REGI W1 Zr94f ALL + REGI W1 Zr96f ALL + + ENDC +; +ELSEIF Fuel "UOX" = THEN + IF NbRegions 20 = THEN +LIBRARY2 := USS: LIBRARY TRACKSS :: + EDIT 0 + ARM + TRAN + PASS 2 + GRMIN 31 GRMAX 93 + CALC + REGI W1 U235 ALL + REGI W1 U236 ALL + REGI W1 Pu239 ALL + REGI W1 Pu240 ALL + REGI W1 Pu241 ALL + REGI W1 Pu242 ALL + REGI W1 Zr90g ALL + REGI W1 Zr91g ALL + REGI W1 Zr92g ALL + REGI W1 Zr94g ALL + REGI W1 Zr96g ALL + REGI W1 Zr90f ALL + REGI W1 Zr91f ALL + REGI W1 Zr92f ALL + REGI W1 Zr94f ALL + REGI W1 Zr96f ALL + + REGI W1 U238 9 + REGI W2 U238 10 + REGI W3 U238 11 + REGI W4 U238 12 + REGI W5 U238 13 17 21 25 + REGI W6 U238 14 18 22 26 + REGI W7 U238 15 19 23 27 + REGI W8 U238 16 20 24 28 + ENDC +; + ELSEIF NbRegions 32 = THEN +LIBRARY2 := USS: LIBRARY TRACKSS :: + EDIT 0 + ARM + TRAN + PASS 2 + GRMIN 31 GRMAX 93 + CALC + REGI W1 U235 ALL + REGI W1 U236 ALL + REGI W1 Pu239 ALL + REGI W1 Pu240 ALL + REGI W1 Pu241 ALL + REGI W1 Pu242 ALL + REGI W1 Zr90g ALL + REGI W1 Zr91g ALL + REGI W1 Zr92g ALL + REGI W1 Zr94g ALL + REGI W1 Zr96g ALL + REGI W1 Zr90f ALL + REGI W1 Zr91f ALL + REGI W1 Zr92f ALL + REGI W1 Zr94f ALL + REGI W1 Zr96f ALL + + REGI W1 U238 9 21 + REGI W2 U238 10 22 + REGI W3 U238 11 23 + REGI W4 U238 12 24 + REGI W5 U238 13 17 25 29 33 37 + REGI W6 U238 14 18 26 30 34 38 + REGI W7 U238 15 19 27 31 35 39 + REGI W8 U238 16 20 28 32 36 40 + ENDC +; + ENDIF ; +ELSE + ECHO "Mauvaise donnée de Fuel" ; +ENDIF ; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Autop2.c2m b/Dragon/data/rep900_mco_proc/Autop2.c2m new file mode 100644 index 0000000..6b702df --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Autop2.c2m @@ -0,0 +1,138 @@ +*DECK Autop2 +*---- +* Name : Autop2.c2m +* Type : DRAGON procedure +* Use : Realize the rerun of self shielding calculation +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBRARY2 := Autop2 LIBRARY2 LIBRARY TRACKSS :: <<Fuel>> <<NbRegions>> ; +* +* Input data : +* LIBRARY : Library with SUBG key word +* LIBRARY2 : Library already self shielded +* TRACKSS : Tracking of the SS geometry +* Fuel : Type of fuel (UOX/MOX) +* NbRegions : Number of regions for the calculation (20/32) +* Output data : +* LIBRARY2 : Self shielded Library + + +PARAMETER LIBRARY2 LIBRARY TRACKSS :: + EDIT 1 + ::: LINKED_LIST LIBRARY2 LIBRARY TRACKSS ; + ; +* ; +MODULE USS: END: ; + +STRING Fuel ; + :: >>Fuel<< ; +INTEGER NbRegions ; + :: >>NbRegions<< ; + +IF Fuel "MOX" = THEN +LIBRARY2 := USS: LIBRARY LIBRARY2 TRACKSS :: + EDIT 0 + ARM + TRAN + PASS 1 + GRMIN 31 GRMAX 93 + CALC + REGI W1 U235 ALL + REGI W1 U236 ALL + REGI W1 Pu239 ALL + REGI W1 Pu240 ALL + REGI W1 Pu241 ALL + REGI W1 Pu242 ALL + REGI W1 Zr90g ALL + REGI W1 Zr91g ALL + REGI W1 Zr92g ALL + REGI W1 Zr94g ALL + REGI W1 Zr96g ALL + REGI W1 Zr90f ALL + REGI W1 Zr91f ALL + REGI W1 Zr92f ALL + REGI W1 Zr94f ALL + REGI W1 Zr96f ALL + + ENDC +; +ELSEIF Fuel "UOX" = THEN + IF NbRegions 20 = THEN +LIBRARY2 := USS: LIBRARY LIBRARY2 TRACKSS :: + EDIT 0 + ARM + TRAN + PASS 1 + GRMIN 31 GRMAX 93 + CALC + REGI W1 U235 ALL + REGI W1 U236 ALL + REGI W1 Pu239 ALL + REGI W1 Pu240 ALL + REGI W1 Pu241 ALL + REGI W1 Pu242 ALL + REGI W1 Zr90g ALL + REGI W1 Zr91g ALL + REGI W1 Zr92g ALL + REGI W1 Zr94g ALL + REGI W1 Zr96g ALL + REGI W1 Zr90f ALL + REGI W1 Zr91f ALL + REGI W1 Zr92f ALL + REGI W1 Zr94f ALL + REGI W1 Zr96f ALL + + REGI W1 U238 9 + REGI W2 U238 10 + REGI W3 U238 11 + REGI W4 U238 12 + REGI W5 U238 13 17 21 25 + REGI W6 U238 14 18 22 26 + REGI W7 U238 15 19 23 27 + REGI W8 U238 16 20 24 28 + ENDC +; + ELSEIF NbRegions 32 = THEN +LIBRARY2 := USS: LIBRARY LIBRARY2 TRACKSS :: + EDIT 0 + ARM + TRAN + PASS 1 + GRMIN 31 GRMAX 93 + CALC + REGI W1 U235 ALL + REGI W1 U236 ALL + REGI W1 Pu239 ALL + REGI W1 Pu240 ALL + REGI W1 Pu241 ALL + REGI W1 Pu242 ALL + REGI W1 Zr90g ALL + REGI W1 Zr91g ALL + REGI W1 Zr92g ALL + REGI W1 Zr94g ALL + REGI W1 Zr96g ALL + REGI W1 Zr90f ALL + REGI W1 Zr91f ALL + REGI W1 Zr92f ALL + REGI W1 Zr94f ALL + REGI W1 Zr96f ALL + + REGI W1 U238 9 21 + REGI W2 U238 10 22 + REGI W3 U238 11 23 + REGI W4 U238 12 24 + REGI W5 U238 13 17 25 29 33 37 + REGI W6 U238 14 18 26 30 34 38 + REGI W7 U238 15 19 27 31 35 39 + REGI W8 U238 16 20 28 32 36 40 + ENDC +; + ENDIF ; +ELSE + ECHO "Mauvaise donnée de Fuel" ; +ENDIF ; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/BranchFlux.c2m b/Dragon/data/rep900_mco_proc/BranchFlux.c2m new file mode 100644 index 0000000..1573312 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/BranchFlux.c2m @@ -0,0 +1,167 @@ +*DECK BranchFlux +*---- +* Name : BranchFlux.c2m +* Type : DRAGON procedure +* Use : Two-level branch flux calculation using the extended +* Le Tellier scheme. +* Author : A. Hebert, based on a procedure from R. Vallerent +* +* Procedure called as: +* +* LIBRARY2 FLUX1 FLUX2 LIBEQ := BranchFlux LIBRARY2 FLUX1 FLUX2 +* BURN164 LIBRARY TRACKSS TRACKN1 TRACKN2 TRKFILN2 :: (data) ; +* +* Input data : +* LIBRARY : Microlib with subgroups +* LIBRARY2 : Self-shielded microlib +* TRACKSS : Tracking of the self-shielding geometry +* TRACKN1 : Tracking of the level 1 geometry +* TRACKN2 : Tracking of the level 2 geometry +* TRKFILN2 : Sequential tracking file of the level 2 geometry +* FLUX1 : Multigroup flux of the first level +* FLUX2 : Multigroup flux of the second level +* BURN164 : Burnup object of level 2 +* Output data : +* LIBRARY2 : Self-shielded microlib +* FLUX1 : Multigroup flux of the first level +* FLUX2 : Multigroup flux of the second level +* LIBEQ : 26-group microlib +* Parameters : +* istep : Burnup step index corresponding to the calculation +* Fuel : Type of fuel ("MOX" or "UOX") +* NbRegions : Number of regions (32 for MOX, 20 or 32 for UOX) +* BUend : Burnup corresponding to the calculation +* BUnextSS : Burnup corresponding to the next self-shielding +* calculation +* Fuelpwr : Fuel power in MW/tonne + +PARAMETER LIBRARY2 FLUX1 FLUX2 LIBEQ BURN164 LIBRARY TRACKSS + TRACKN1 TRACKN2 TRKFILN2 :: + EDIT 1 + ::: LINKED_LIST LIBRARY2 FLUX1 FLUX2 BURN164 LIBEQ + LIBRARY TRACKSS TRACKN1 TRACKN2 ; + ::: SEQ_BINARY TRKFILN2 ; + ; +* ; +*---- +* Modules used in this procedure +*---- +MODULE LIB: ASM: FLU: EDI: EVO: GREP: ABORT: DELETE: END: ; +PROCEDURE Autop1 Autop2 ConcLIBEQ MultLIBEQ_20 MultLIBEQ_32 ; +LINKED_LIST EDITION SYS ; +REAL B2 Kinf1 Kinf2 Keff2 Tend ; + +STRING Fuel ; +INTEGER istep NbRegions ; +REAL BUend BUnextSS Fuelpwr ; + :: >>istep<< >>Fuel<< >>NbRegions<< >>BUend<< >>BUnextSS<< + >>Fuelpwr<< ; + + EVALUATE Tend := BUend Fuelpwr / ; + ECHO "-----------------------------------------------------------" ; + ECHO "Two-level flux calculation at step=" istep ; + +*---- +* Resonance self-shielding +*---- + ECHO "BURNUP_sortie = " BUend "Next SS = " BUnextSS ; + IF BUend BUnextSS - ABS 1E-5 < THEN + ECHO "Perform a self-shielding calculation" ; + ECHO "-----------------------------------------------------------" ; + IF istep 1 = THEN + LIBRARY2 := Autop1 LIBRARY TRACKSS :: <<Fuel>> <<NbRegions>> ; + ELSE + LIBRARY2 := Autop2 LIBRARY2 LIBRARY TRACKSS + :: <<Fuel>> <<NbRegions>> ; + ENDIF ; + ELSE + ECHO "No self-shielding calculation" ; + ECHO "-----------------------------------------------------------" ; + ENDIF ; + +************************************************************************ +* BEGIN FLUX CALCULATION * +************************************************************************ + SYS := ASM: LIBRARY2 TRACKN1 :: + EDIT 0 ARM ; + + IF istep 1 = THEN + FLUX1 := FLU: SYS LIBRARY2 TRACKN1 + :: EDIT 0 TYPE B B1 SIGS ; + ELSE + FLUX1 := FLU: FLUX1 SYS LIBRARY2 TRACKN1 + :: EDIT 0 TYPE B B1 SIGS ; + ENDIF ; + + GREP: FLUX1 :: GETVAL 'B2 B1HOM ' 1 1 1 >>B2<< ; + ECHO "B2-1 = " B2 " at BU = " BUend ; + + GREP: FLUX1 :: GETVAL 'K-INFINITY ' 1 1 1 >>Kinf1<< ; + ECHO "K-inf1_sortie is " Kinf1 " at BU = " BUend ; + +*---- +* 26 groups energy condensation +*---- + + EDITION := EDI: FLUX1 LIBRARY2 TRACKN1 :: + EDIT 0 + MICR ALL + MERG MIX + COND 10 14 18 26 33 40 49 56 65 74 86 196 227 + 230 233 238 241 244 247 254 259 263 267 272 277 + SAVE ON COND26 + ; + + LIBEQ := EDITION :: + STEP UP COND26 + ; + + SYS := DELETE: SYS ; + +*---- +* Second level flux calculation (MOC) +*---- + IF NbRegions 20 = THEN + LIBEQ := MultLIBEQ_20 LIBEQ ; + ELSEIF NbRegions 32 = THEN + LIBEQ := MultLIBEQ_32 LIBEQ ; + ENDIF ; + +*---- +* Recover number densities for the current burnup step +*---- + LIBEQ := ConcLIBEQ LIBEQ BURN164 :: <<istep>> ; + + LIBEQ := LIB: LIBEQ :: EDIT 0 MACR MIXS ; + + SYS := ASM: LIBEQ TRACKN2 TRKFILN2 :: + EDIT 0 ARM ; + + IF istep 1 = THEN + FLUX2 := FLU: SYS LIBEQ TRACKN2 TRKFILN2 :: + EDIT 0 + TYPE K B1 PNL BUCK <<B2>> + THER 5 5E-5 EXTE 30 1E-5 UNKT 1E-4 ; + ELSE + FLUX2 := FLU: FLUX2 SYS LIBEQ TRACKN2 TRKFILN2 :: + EDIT 0 + TYPE K B1 PNL BUCK <<B2>> + THER 5 5E-5 EXTE 30 1E-5 UNKT 1E-4 ; + ENDIF ; + + GREP: FLUX2 :: GETVAL 'K-INFINITY ' 1 1 1 >>Kinf2<< ; + ECHO "K-inf2_sortie is " Kinf2 " at BU = " BUend ; + + GREP: FLUX2 :: GETVAL 'K-EFFECTIVE ' 1 1 1 >>Keff2<< ; + ECHO "K-eff2_sortie is " Keff2 " at BU = " BUend ; + + GREP: FLUX2 :: GETVAL 'B2 B1HOM ' 1 1 1 >>B2<< ; + ECHO "B2-2 = " B2 " at BU = " BUend ; + + SYS EDITION := DELETE: SYS EDITION ; +************************************************************************ +* END FLUX CALCULATION * +************************************************************************ + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/CalcFlux.c2m b/Dragon/data/rep900_mco_proc/CalcFlux.c2m new file mode 100644 index 0000000..eed3d79 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/CalcFlux.c2m @@ -0,0 +1,179 @@ +*DECK CalcFlux +*---- +* Name : CalcFlux.c2m +* Type : DRAGON procedure +* Use : Two-level flux calculation using the extended +* Le Tellier scheme. +* Author : A. Hebert, based on a procedure from R. Vallerent +* +* Procedure called as: +* +* LIBRARY2 FLUX1 FLUX2 BURN164 LIBEQ := CalcFlux LIBRARY2 FLUX1 +* FLUX2 BURN164 LIBEQ LIBRARY TRACKSS TRACKN1 TRACKN2 TRKFILN2 +* :: (data) ; +* +* Input data : +* LIBRARY : Microlib with subgroups +* LIBRARY2 : Self-shielded microlib +* TRACKSS : Tracking of the self-shielding geometry +* TRACKN1 : Tracking of the level 1 geometry +* TRACKN2 : Tracking of the level 2 geometry +* TRKFILN2 : Sequential tracking file of the level 2 geometry +* FLUX1 : Multigroup flux of the first level +* FLUX2 : Multigroup flux of the second level +* BURN164 : Burnup object of level 2 +* Output data : +* LIBRARY2 : Self-shielded microlib +* FLUX1 : Multigroup flux of the first level +* FLUX2 : Multigroup flux of the second level +* BURN164 : Burnup object of level 2 +* LIBEQ : 26-group microlib +* Parameters : +* istep : Burnup step index corresponding to the calculation +* Fuel : Type of fuel ("MOX" or "UOX") +* NbRegions : Number of regions (32 for MOX, 20 or 32 for UOX) +* BUend : Burnup corresponding to the calculation +* BUnextSS : Burnup corresponding to the next self-shielding +* calculation +* Fuelpwr : Fuel power in MW/tonne + +PARAMETER LIBRARY2 FLUX1 FLUX2 BURN164 LIBEQ LIBRARY TRACKSS TRACKN1 + TRACKN2 TRKFILN2 :: + EDIT 1 + ::: LINKED_LIST LIBRARY2 FLUX1 FLUX2 BURN164 LIBEQ LIBRARY + TRACKSS TRACKN1 TRACKN2 ; + ::: SEQ_BINARY TRKFILN2 ; + ; +* ; +*---- +* Modules used in this procedure +*---- +MODULE LIB: ASM: FLU: EDI: EVO: GREP: ABORT: DELETE: END: ; +PROCEDURE Autop1 Autop2 ConcLIBEQ MultLIBEQ_20 MultLIBEQ_32 ; +LINKED_LIST EDITION SYS ; +REAL B2 Kinf1 Kinf2 Keff2 Tend ; + +STRING Fuel ; +INTEGER istep NbRegions ; +REAL BUend BUnextSS Fuelpwr ; + :: >>istep<< >>Fuel<< >>NbRegions<< >>BUend<< >>BUnextSS<< + >>Fuelpwr<< ; + + EVALUATE Tend := BUend Fuelpwr / ; + ECHO "-----------------------------------------------------------" ; + ECHO "Two-level flux calculation at step=" istep ; + ECHO "BURNUP_sortie = " BUend "Next SS = " BUnextSS ; + IF BUend BUnextSS - ABS 1E-5 < THEN + ECHO "Perform a self-shielding calculation" ; + ECHO "-----------------------------------------------------------" ; + IF istep 1 = THEN + LIBRARY2 := Autop1 LIBRARY TRACKSS :: <<Fuel>> <<NbRegions>> ; + ELSE + LIBRARY2 := Autop2 LIBRARY2 LIBRARY TRACKSS + :: <<Fuel>> <<NbRegions>> ; + ENDIF ; + ELSE + ECHO "No self-shielding calculation" ; + ECHO "-----------------------------------------------------------" ; + ENDIF ; +************************************************************************ +* BEGIN FLUX CALCULATION * +************************************************************************ + SYS := ASM: LIBRARY2 TRACKN1 :: + EDIT 0 ARM ; + + IF istep 1 = THEN + FLUX1 := FLU: SYS LIBRARY2 TRACKN1 + :: EDIT 0 TYPE B B1 SIGS ; + ELSE + FLUX1 := FLU: FLUX1 SYS LIBRARY2 TRACKN1 + :: EDIT 0 TYPE B B1 SIGS ; + ENDIF ; + + GREP: FLUX1 :: GETVAL 'B2 B1HOM ' 1 1 1 >>B2<< ; + ECHO "B2-1 = " B2 " at BU = " BUend ; + + GREP: FLUX1 :: GETVAL 'K-INFINITY ' 1 1 1 >>Kinf1<< ; + ECHO "K-inf1_sortie is " Kinf1 " at BU = " BUend ; + +*---- +* 26 groups energy condensation +*---- + + EDITION := EDI: FLUX1 LIBRARY2 TRACKN1 :: + EDIT 0 + MICR ALL + MERG MIX + COND 10 14 18 26 33 40 49 56 65 74 86 196 227 + 230 233 238 241 244 247 254 259 263 267 272 277 + SAVE ON COND26 + ; + + LIBEQ := EDITION :: + STEP UP COND26 + ; + + SYS := DELETE: SYS ; + +*---- +* Second level flux calculation (MOC) +*---- + +* added by remi( + + IF NbRegions 20 = THEN + LIBEQ := MultLIBEQ_20 LIBEQ ; + ELSEIF NbRegions 32 = THEN + LIBEQ := MultLIBEQ_32 LIBEQ ; + ENDIF ; + + IF istep 1 > THEN + LIBEQ := ConcLIBEQ LIBEQ BURN164 :: <<istep>> ; + ENDIF ; + + LIBEQ := LIB: LIBEQ :: EDIT 0 MACR MIXS ; +* )added by remi + + SYS := ASM: LIBEQ TRACKN2 TRKFILN2 :: + EDIT 0 ARM ; + + IF istep 1 = THEN + FLUX2 := FLU: SYS LIBEQ TRACKN2 TRKFILN2 :: + EDIT 0 + TYPE K B1 PNL BUCK <<B2>> + THER 5 5E-5 EXTE 30 1E-5 UNKT 1E-4 ; + ELSE + FLUX2 := FLU: FLUX2 SYS LIBEQ TRACKN2 TRKFILN2 :: + EDIT 0 + TYPE K B1 PNL BUCK <<B2>> + THER 5 5E-5 EXTE 30 1E-5 UNKT 1E-4 ; + ENDIF ; + + GREP: FLUX2 :: GETVAL 'K-INFINITY ' 1 1 1 >>Kinf2<< ; + ECHO "K-inf2_sortie is " Kinf2 " at BU = " BUend ; + + GREP: FLUX2 :: GETVAL 'K-EFFECTIVE ' 1 1 1 >>Keff2<< ; + ECHO "K-eff2_sortie is " Keff2 " at BU = " BUend ; + + GREP: FLUX2 :: GETVAL 'B2 B1HOM ' 1 1 1 >>B2<< ; + ECHO "B2-2 = " B2 " at BU = " BUend ; + + SYS EDITION := DELETE: SYS EDITION ; +************************************************************************ +* END FLUX CALCULATION * +************************************************************************ + +IF istep 1 = THEN + BURN164 LIBEQ := EVO: LIBEQ FLUX2 TRACKN2 :: + EDIT 0 + SAVE <<Tend>> DAY POWR <<Fuelpwr>> + ; +ELSE + BURN164 LIBEQ := EVO: BURN164 LIBEQ FLUX2 TRACKN2 :: + EDIT 0 + SAVE <<Tend>> DAY POWR <<Fuelpwr>> + ; +ENDIF ; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/ConcLIBEQ.c2m b/Dragon/data/rep900_mco_proc/ConcLIBEQ.c2m new file mode 100644 index 0000000..939307f --- /dev/null +++ b/Dragon/data/rep900_mco_proc/ConcLIBEQ.c2m @@ -0,0 +1,60 @@ +*DECK ConcLIBEQ +*---- +* Name : ConcLIBEQ.c2m +* Type : DRAGON procedure +* Use : Update of the LIBEQ concentrations +* with the ones of BURN164 +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBEQ := ConcLIBEQ LIBEQ BURN164 :: <<istep>> ; +* +* Input data : +* LIBEQ : Library which concentrations are to be updated +* BURN164 : Burnup object which contains the 156-fuel mixes +* concentrations +* istep : step of the calculation (where the concentrations +* will be read in BURN164) +* Output data : +* LIBEQ : Library which concentrations are updated + + +PARAMETER LIBEQ BURN164 :: + EDIT 1 + ::: LINKED_LIST LIBEQ BURN164 ; + ; +* ; +MODULE LIB: END: ; + +INTEGER istep ; + :: >>istep<< ; + +LIBEQ := LIB: LIBEQ BURN164 :: + EDIT 0 + BURN <<istep>> + + MIX 9 MIX 10 MIX 11 MIX 12 MIX 13 MIX 14 MIX 15 MIX 16 + MIX 17 MIX 18 MIX 19 MIX 20 MIX 21 MIX 22 MIX 23 MIX 24 + MIX 25 MIX 26 MIX 27 MIX 28 MIX 29 MIX 30 MIX 31 MIX 32 + MIX 33 MIX 34 MIX 35 MIX 36 MIX 37 MIX 38 MIX 39 MIX 40 + MIX 41 MIX 42 MIX 43 MIX 44 MIX 45 MIX 46 MIX 47 MIX 48 + MIX 49 MIX 50 MIX 51 MIX 52 MIX 53 MIX 54 MIX 55 MIX 56 + MIX 57 MIX 58 MIX 59 MIX 60 MIX 61 MIX 62 MIX 63 MIX 64 + MIX 65 MIX 66 MIX 67 MIX 68 MIX 69 MIX 70 MIX 71 MIX 72 + MIX 73 MIX 74 MIX 75 MIX 76 MIX 77 MIX 78 MIX 79 MIX 80 + MIX 81 MIX 82 MIX 83 MIX 84 MIX 85 MIX 86 MIX 87 MIX 88 + MIX 89 MIX 90 MIX 91 MIX 92 MIX 93 MIX 94 MIX 95 MIX 96 + MIX 97 MIX 98 MIX 99 MIX 100 MIX 101 MIX 102 MIX 103 MIX 104 + MIX 105 MIX 106 MIX 107 MIX 108 MIX 109 MIX 110 MIX 111 MIX 112 + MIX 113 MIX 114 MIX 115 MIX 116 MIX 117 MIX 118 MIX 119 MIX 120 + MIX 121 MIX 122 MIX 123 MIX 124 MIX 125 MIX 126 MIX 127 MIX 128 + MIX 129 MIX 130 MIX 131 MIX 132 MIX 133 MIX 134 MIX 135 MIX 136 + MIX 137 MIX 138 MIX 139 MIX 140 MIX 141 MIX 142 MIX 143 MIX 144 + MIX 145 MIX 146 MIX 147 MIX 148 MIX 149 MIX 150 MIX 151 MIX 152 + MIX 153 MIX 154 MIX 155 MIX 156 MIX 157 MIX 158 MIX 159 MIX 160 + MIX 161 MIX 162 MIX 163 MIX 164 +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/ConcLIBRARY2.c2m b/Dragon/data/rep900_mco_proc/ConcLIBRARY2.c2m new file mode 100644 index 0000000..45f10fc --- /dev/null +++ b/Dragon/data/rep900_mco_proc/ConcLIBRARY2.c2m @@ -0,0 +1,54 @@ +*DECK ConcLIBRARY2 +*---- +* Name : ConcLIBRARY2.c2m +* Type : DRAGON procedure +* Use : Update of the LIBRARY2 concentrations +* with the ones of LIBHOM +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBRARY2 := ConcLIBRARY2 LIBRARY2 LIBHOM :: <<NbRegions>> ; +* +* Input data : +* LIBRARY2 : Library which concentrations are to be updated +* LIBHOM : One group library after homogeneization and 2nd level +* evolution +* NbRegions : Number of fuel regions for the calculation (20/32) +* Output data : +* LIBRARY2 : Library which concentrations are updated + + +PARAMETER LIBRARY2 LIBHOM :: + EDIT 1 + ::: LINKED_LIST LIBRARY2 LIBHOM ; + ; +* ; +MODULE LIB: END: ; + +INTEGER NbRegions ; + :: >>NbRegions<< ; + +IF NbRegions 20 = THEN + LIBRARY2 := LIB: LIBRARY2 LIBHOM :: + EDIT 0 + MAXS + MIX 9 MIX 10 MIX 11 MIX 12 MIX 13 MIX 14 MIX 15 MIX 16 + MIX 17 MIX 18 MIX 19 MIX 20 MIX 21 MIX 22 MIX 23 MIX 24 + MIX 25 MIX 26 MIX 27 MIX 28 + ; +ELSEIF NbRegions 32 = THEN + LIBRARY2 := LIB: LIBRARY2 LIBHOM :: + EDIT 0 + MAXS + MIX 9 MIX 10 MIX 11 MIX 12 MIX 13 MIX 14 MIX 15 MIX 16 + MIX 17 MIX 18 MIX 19 MIX 20 MIX 21 MIX 22 MIX 23 MIX 24 + MIX 25 MIX 26 MIX 27 MIX 28 MIX 29 MIX 30 MIX 31 MIX 32 + MIX 33 MIX 34 MIX 35 MIX 36 MIX 37 MIX 38 MIX 39 MIX 40 + ; +ELSE + ECHO "Mauvaise donnee de NbRegions" ; +ENDIF ; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Condense_2gr.c2m b/Dragon/data/rep900_mco_proc/Condense_2gr.c2m new file mode 100644 index 0000000..719a647 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Condense_2gr.c2m @@ -0,0 +1,137 @@ +*DECK Condense_2gr +*---- +* Name : Condense_2gr.c2m +* Type : DRAGON procedure +* Use : Realize the deux groups condensation +* and pin by pin or assembly homogeneization +* Author : R. Vallerent +* +* Procedure called as: +* +* EDIOBJ := Condense_2gr FLUX2 LIBEQ TRACKN2 :: <<Library>> ; +* +* Input data : +* FLUX2 : Multigroup Flux +* LIBEQ : Library xs +* TRACKN2 : Tracking of the level 2 geometry +* Output data : +* EDIOBJ : Pin by pin homogeneization +* Parameters : +* Library : Name of input cross-section library +* TypeCompo : Type of edition: PinByPin or Assembly + +PARAMETER EDIOBJ FLUX2 LIBEQ TRACKN2 :: + EDIT 1 + ::: LINKED_LIST EDIOBJ FLUX2 LIBEQ TRACKN2 ; + ; +* ; +MODULE EDI: ABORT: END: ; + +STRING Library TypeCompo ; + :: >>Library<< >>TypeCompo<< ; + +STRING Am242M Pm148M ; +IF Library "CEA05V4_281" = THEN + EVALUATE Am242M Pm148M := 'Am242M' 'Pm148M' ; +ELSEIF Library "DLIBJ3_281" = THEN + EVALUATE Am242M Pm148M := 'Am242m' 'Pm148m' ; +ELSE + ECHO "Invalid parameter : Library" ; + ABORT: ; +ENDIF ; + +IF TypeCompo "PinByPin" = THEN + EDIOBJ := EDI: FLUX2 LIBEQ TRACKN2 :: + EDIT 0 + MICR 35 U234 U235 U236 U237 U238 Np237 Np238 Np239 Pu238 Pu239 + Pu240 Pu241 Pu242 Cm242 Cm243 Cm244 Cm245 Am241 <<Am242M>> Am243 + Pm147 Pm148 <<Pm148M>> Pm149 Sm147 Sm148 Sm149 Sm150 Nd146 Nd147 + Nd148 B10 B11 Xe135 I135 + MERG REGI + 1 1 1 1 1 1 1 1 1 1 1 1 2 2 + 2 2 2 2 2 2 2 2 2 2 2 2 2 2 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 4 4 4 4 4 4 4 4 4 4 4 4 + 4 4 4 4 4 4 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 6 6 6 6 6 6 + 6 6 6 6 6 6 6 6 6 6 7 7 7 7 + 7 7 7 7 7 7 7 7 7 7 7 7 7 7 + 8 8 8 8 8 8 8 8 8 8 8 8 8 8 + 8 8 9 9 9 9 9 9 9 9 9 9 9 9 + 9 9 9 9 9 9 10 10 10 10 10 10 10 10 + 10 10 10 10 10 10 10 10 11 11 11 11 11 11 + 11 11 11 11 11 11 11 11 11 11 11 11 11 11 + 11 11 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 13 13 13 13 + 13 13 13 13 13 13 13 13 13 13 13 13 13 13 + 13 13 13 13 14 14 14 14 14 14 14 14 14 14 + 14 14 14 14 14 14 14 14 14 14 14 14 15 15 + 15 15 15 15 15 15 15 15 15 15 15 15 15 15 + 15 15 15 15 15 15 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 17 17 17 17 17 17 17 17 17 17 17 17 17 17 + 17 17 17 17 17 17 17 17 17 17 17 18 18 18 + 18 18 18 18 18 18 18 18 18 18 18 18 18 19 + 19 19 19 19 19 19 19 19 19 19 19 19 19 19 + 19 19 19 19 19 19 19 20 20 20 20 20 20 20 + 20 20 20 20 20 20 20 20 20 20 20 20 20 20 + 20 21 21 21 21 21 21 21 21 21 21 21 21 21 + 21 21 21 21 21 21 21 21 21 22 22 22 22 22 + 22 22 22 22 22 22 22 22 22 22 22 22 22 22 + 22 22 22 23 23 23 23 23 23 23 23 23 23 23 + 23 23 23 23 23 23 23 23 23 23 23 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 25 25 25 25 25 25 + 25 25 25 25 25 25 25 25 25 25 25 25 26 26 + 26 26 26 26 26 26 26 26 26 26 26 26 26 26 + 26 26 26 26 26 26 27 27 27 27 27 27 27 27 + 27 27 27 27 27 27 27 27 27 27 27 27 27 27 + 28 28 28 28 28 28 28 28 28 28 28 28 28 28 + 28 28 28 28 28 28 28 28 28 28 28 28 28 29 + 29 29 29 29 29 29 29 29 29 29 29 29 29 29 + 29 29 29 29 29 29 29 30 30 30 30 30 30 30 + 30 30 30 30 30 30 30 30 30 30 30 30 30 30 + 30 30 30 30 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 32 32 32 32 32 32 32 32 + 32 32 32 32 32 32 32 32 32 32 32 32 32 32 + 33 33 33 33 33 33 33 33 33 33 33 33 33 33 + 33 33 33 33 33 33 33 33 34 34 34 34 34 34 + 34 34 34 34 34 34 34 34 34 34 34 34 34 34 + 34 34 35 35 35 35 35 35 35 35 35 35 35 35 + 35 35 35 35 35 35 35 35 35 35 35 35 35 36 + 36 36 36 36 36 36 36 36 36 36 36 36 36 36 + 36 36 36 37 37 37 37 37 37 37 37 37 37 37 + 37 37 37 37 37 37 37 37 37 37 37 38 38 38 + 38 38 38 38 38 38 38 38 38 38 38 38 38 38 + 38 38 38 38 38 39 39 39 39 39 39 39 39 39 + 39 39 39 39 39 39 39 39 39 39 39 39 39 39 + 39 39 40 40 40 40 40 40 40 40 40 40 40 40 + 40 40 40 40 41 41 41 41 41 41 41 41 41 41 + 41 41 41 41 41 41 41 41 41 41 41 41 42 42 + 42 42 42 42 42 42 42 42 42 42 42 42 42 42 + 42 42 42 42 42 42 42 42 42 43 43 43 43 43 + 43 43 43 43 43 43 43 43 43 43 43 44 44 44 + 44 44 44 44 44 44 44 44 44 44 44 44 44 44 + 44 44 44 44 44 44 44 44 45 45 45 45 45 45 + 45 45 45 45 45 45 45 45 45 45 45 45 45 45 + COND 19 + SAVE ON EDI2A + ; +ELSEIF TypeCompo "Assembly" = THEN + EDIOBJ := EDI: FLUX2 LIBEQ TRACKN2 :: + EDIT 0 + MICR 35 U234 U235 U236 U237 U238 Np237 Np238 Np239 Pu238 Pu239 + Pu240 Pu241 Pu242 Cm242 Cm243 Cm244 Cm245 Am241 <<Am242M>> Am243 + Pm147 Pm148 <<Pm148M>> Pm149 Sm147 Sm148 Sm149 Sm150 Nd146 Nd147 + Nd148 B10 B11 Xe135 I135 + MERG COMP + COND 19 + SAVE ON EDI2B + ; +ELSE + ECHO "TypeCompo incorrect" ; + ABORT: ; +ENDIF ; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Geo_N1_20.c2m b/Dragon/data/rep900_mco_proc/Geo_N1_20.c2m new file mode 100644 index 0000000..3ce8892 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Geo_N1_20.c2m @@ -0,0 +1,139 @@ +*DECK Geo_N1 +*---- +* Name : Geo_N1.c2m +* Type : DRAGON procedure +* Use : Geometry generation for 1st Level Flux Calculation +* with 20 fuel regions +* Author : R. Vallerent +* +* Procedure called as: +* +*GEON1 := Geo_N1 :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> +* <<R_int_TG>> <<R_ext_TG>> <<R_int_TI>> +* <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>> +* <<CoteL>> <<RmodeN1>> ; +* +* Input data : +* Rcomb1 : 50/100 of outer radius of fuel (cm) +* Rcomb2 : 80/100 of outer radius of fuel (cm) +* Rcomb3 : 95/100 of outer radius of fuel (cm) +* Rcomb4 : 100/100 of outer radius of fuel (cm) +* R_int_TG : Inner radius of cladding of guide tube (cm) +* R_ext_TG : Outer radius of cladding of guide tube (cm) +* R_int_TI : Inner radius of cladding of guide tube (cm) +* R_ext_TI : Outer radius of cladding of guide tube (cm) +* R_int_G : Inner radius of cladding of fuel tube (cm) +* R_ext_G : Outer radius of cladding of fuel tube (cm) +* Cote : Lattice pitch (cm) +* CoteL : Lattice pitch + Water space (cm) +* RmodeN1 : Parameter for flux calculation in Level 1 (cm) +* +* Output data : +* GEON1 : Geometry for 1st Level Flux Calculation + + +PARAMETER GEON1 :: + EDIT 0 + ::: LINKED_LIST GEON1 ; + ; +*---- +* Modules used in this procedure +*---- +MODULE GEO: END: ; + +*---- +* Input data recovery +*---- +* ; +REAL Rcomb1 Rcomb2 Rcomb3 Rcomb4 ; +:: >>Rcomb1<< >>Rcomb2<< >>Rcomb3<< >>Rcomb4<< ; +REAL R_int_TG R_ext_TG R_int_TI R_ext_TI ; +:: >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< ; +REAL R_int_G R_ext_G Cote CoteL ; +:: >>R_int_G<< >>R_ext_G<< >>Cote<< >>CoteL<< ; +REAL RmodeN1 ; +:: >>RmodeN1<< ; + +GEON1 := GEO: :: CAR2D 9 9 + EDIT 0 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + CELL TI C1 C1 T1 C1 C1 T2 C1 C4 + C2 C2 C1 C2 C2 C1 C2 C4 + C2 C1 C2 C2 C1 C2 C4 + T1 C1 C1 T2 C1 C4 + C2 C1 C1 C2 C4 + T2 C1 C3 C4 + C2 C3 C4 + C3 C4 + C5 + + MERGE 8 1 1 7 1 1 6 1 4 + 2 2 1 2 2 1 2 4 + 2 1 2 2 1 2 4 + 7 1 1 6 1 4 + 2 1 1 2 4 + 6 1 3 4 + 2 3 4 + 3 4 + 5 + +* T2 -> 6 , T1 -> 7 , TI -> 8 + + TURN A A E A A E A A A + A E F A E D A A + C B G C H G A + A G C A G A + A B D A A + A A A A + A G A + A A + A + + + ::: C1 := GEO: CARCEL 7 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> <<RmodeN1>> + MIX 9 10 11 12 3 2 4 5 + ; + ::: C2 := GEO: C1 + MIX 13 14 15 16 3 2 4 5 + ; + ::: C3 := GEO: C1 + MIX 17 18 19 20 3 2 4 5 + ; + ::: C4 := GEO: CARCEL 7 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> <<RmodeN1>> + MIX 21 22 23 24 3 2 4 5 + ; + ::: C5 := GEO: CARCEL 7 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<CoteL>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> <<RmodeN1>> + MIX 25 26 27 28 3 2 4 5 + ; + ::: T2 := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TG>> <<R_ext_TG>> + MIX 6 1 6 + ; + ::: T1 := GEO: T2 + MIX 7 1 7 + ; + ::: TI := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TI>> <<R_ext_TI>> + MIX 8 1 8 + ; +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Geo_N1_32.c2m b/Dragon/data/rep900_mco_proc/Geo_N1_32.c2m new file mode 100644 index 0000000..a648442 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Geo_N1_32.c2m @@ -0,0 +1,148 @@ +*DECK Geo_N1 +*---- +* Name : Geo_N1.c2m +* Type : DRAGON procedure +* Use : Geometry generation for 1st Level Flux Calculation +* with 32 fuel regions +* Author : R. Vallerent +* +* Procedure called as: +* +*GEON1 := Geo_N1 :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> +* <<R_int_TG>> <<R_ext_TG>> <<R_int_TI>> +* <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>> +* <<CoteL>> <<RmodeN1>> ; +* +* Input data : +* Rcomb1 : 50/100 of outer radius of fuel (cm) +* Rcomb2 : 80/100 of outer radius of fuel (cm) +* Rcomb3 : 95/100 of outer radius of fuel (cm) +* Rcomb4 : 100/100 of outer radius of fuel (cm) +* R_int_TG : Inner radius of cladding of guide tube (cm) +* R_ext_TG : Outer radius of cladding of guide tube (cm) +* R_int_TI : Inner radius of cladding of guide tube (cm) +* R_ext_TI : Outer radius of cladding of guide tube (cm) +* R_int_G : Inner radius of cladding of fuel tube (cm) +* R_ext_G : Outer radius of cladding of fuel tube (cm) +* Cote : Lattice pitch (cm) +* CoteL : Lattice pitch + Water space (cm) +* RmodeN1 : Parameter for flux calculation in Level 1 (cm) +* +* Output data : +* GEON1 : Geometry for 1st Level Flux Calculation + + +PARAMETER GEON1 :: + EDIT 0 + ::: LINKED_LIST GEON1 ; + ; +*---- +* Modules used in this procedure +*---- +MODULE GEO: END: ; + +*---- +* Input data recovery +*---- +* ; +REAL Rcomb1 Rcomb2 Rcomb3 Rcomb4 ; +:: >>Rcomb1<< >>Rcomb2<< >>Rcomb3<< >>Rcomb4<< ; +REAL R_int_TG R_ext_TG R_int_TI R_ext_TI ; +:: >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< ; +REAL R_int_G R_ext_G Cote CoteL ; +:: >>R_int_G<< >>R_ext_G<< >>Cote<< >>CoteL<< ; +REAL RmodeN1 ; +:: >>RmodeN1<< ; + +GEON1 := GEO: :: CAR2D 9 9 + EDIT 0 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + CELL TI C1 C1 T1 C1 C1 T2 C4 C6 + C2 C2 C1 C2 C2 C1 C2 C6 + C2 C1 C2 C2 C1 C2 C6 + T1 C1 C1 T2 C4 C6 + C2 C1 C1 C2 C6 + T2 C1 C3 C6 + C2 C3 C6 + C5 C7 + C8 + + MERGE 11 1 1 10 1 1 9 4 6 + 2 2 1 2 2 1 2 6 + 2 1 2 2 1 2 6 + 10 1 1 9 4 6 + 2 1 1 2 6 + 9 1 3 6 + 2 3 6 + 5 7 + 8 + +* T2 -> 9 , T1 -> 10 , TI -> 11 + + TURN A A E A A E A A A + A E F A E D A A + C B G C H G A + A G C A G A + A B D A A + A A A A + A G A + A A + A + + + ::: C1 := GEO: CARCEL 7 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> <<RmodeN1>> + MIX 9 10 11 12 3 2 4 5 + ; + ::: C2 := GEO: C1 + MIX 13 14 15 16 3 2 4 5 + ; + ::: C3 := GEO: C1 + MIX 17 18 19 20 3 2 4 5 + ; + ::: C4 := GEO: C1 + MIX 21 22 23 24 3 2 4 5 + ; + ::: C5 := GEO: C1 + MIX 25 26 27 28 3 2 4 5 + ; + ::: C6 := GEO: CARCEL 7 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> <<RmodeN1>> + MIX 29 30 31 32 3 2 4 5 + ; + ::: C7 := GEO: C6 + MIX 33 34 35 36 3 2 4 5 + ; + ::: C8 := GEO: CARCEL 7 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<CoteL>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> <<RmodeN1>> + MIX 37 38 39 40 3 2 4 5 + ; + ::: T2 := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TG>> <<R_ext_TG>> + MIX 6 1 6 + ; + ::: T1 := GEO: T2 + MIX 7 1 7 + ; + ::: TI := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TI>> <<R_ext_TI>> + MIX 8 1 8 + ; +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Geo_N2_20.c2m b/Dragon/data/rep900_mco_proc/Geo_N2_20.c2m new file mode 100644 index 0000000..8ed81a1 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Geo_N2_20.c2m @@ -0,0 +1,797 @@ +*DECK Geo_N2 +*---- +* Name : Geo_N2.c2m +* Type : DRAGON procedure +* Use : Geometry generation for 1st Level Flux Calculation +* Author : R. Vallerent +* +* Procedure called as: +* +*GEON2 := Geo_N2 :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> +* <<R_int_TG>> <<R_ext_TG>> <<R_int_TI>> +* <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>> +* <<Lame>> <<RmodeN2>> +* <<R_ext_G2>> <<R_ext_TGm>> <<R_ext_TIm>> ; +* +* Input data : +* Rcomb1 : 50/100 of outer radius of fuel (cm) +* Rcomb2 : 80/100 of outer radius of fuel (cm) +* Rcomb3 : 95/100 of outer radius of fuel (cm) +* Rcomb4 : 100/100 of outer radius of fuel (cm) +* R_int_TG : Inner radius of cladding of guide tube (cm) +* R_ext_TG : Outer radius of cladding of guide tube (cm) +* R_int_TI : Inner radius of cladding of guide tube (cm) +* R_ext_TI : Outer radius of cladding of guide tube (cm) +* R_int_G : Inner radius of cladding of fuel tube (cm) +* R_ext_G : Outer radius of cladding of fuel tube (cm) +* Cote : Lattice pitch (cm) +* Lame : Water space on the side of the assembly (cm) +* RmodeN2 : Parameter for flux calculation in Level 2 (cm) +* R_ext_G2 : R_ext_G * 2 (cm) +* R_ext_TGm : -R_ext_TGm (cm) +* R_ext_TIm : -R_ext_TIm (cm) +* +* Output data : +* GEON2 : Geometry for Self-Shielding Calculation + + +PARAMETER GEON2 :: + EDIT 0 + ::: LINKED_LIST GEON2 ; + ; +*---- +* Modules used in this procedure +*---- +MODULE GEO: END: ; + +*---- +* Input data recovery +*---- +* ; +REAL Rcomb1 Rcomb2 Rcomb3 Rcomb4 ; +:: >>Rcomb1<< >>Rcomb2<< >>Rcomb3<< >>Rcomb4<< ; +REAL R_int_TG R_ext_TG R_int_TI R_ext_TI ; +:: >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< ; +REAL R_int_G R_ext_G Cote Lame ; +:: >>R_int_G<< >>R_ext_G<< >>Cote<< >>Lame<< ; +REAL RmodeN2 ; +:: >>RmodeN2<< ; + +REAL Angle := 0.785398163 ; +REAL meshc1 meshc2 := 0.37 0.89 ; + +REAL R_ext_G2 R_ext_TGm R_ext_TIm ; +:: >>R_ext_G2<< >>R_ext_TGm<< >>R_ext_TIm<< ; + +GEON2 := GEO: :: CAR2D 10 10 + EDIT 0 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + CELL TI C01 C03 T1 C08 C13 T2 C23 C31 L1 + C02 C04 C06 C09 C14 C18 C24 C32 L1 + C05 C07 C10 C15 C19 C25 C33 L1 + T1 C11 C16 T2 C26 C34 L1 + C12 C17 C20 C27 C35 L1 + T2 C21 C28 C36 L1 + C22 C29 C37 L1 + C30 C38 L1 + C39 L1 + L2 + + ::: C01 := GEO: CARCEL 2 3 3 + MESHX 0.0 <<meshc1>> <<meshc2>> <<Cote>> + MESHY 0.0 <<meshc1>> <<meshc2>> <<Cote>> + SPLITX 1 1 1 + SPLITY 1 1 1 + RADIUS 0.0 <<R_ext_G>> <<RmodeN2>> + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + CLUSTER PIN + ::: PIN := GEO: TUBE 6 + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MESHX 0.0 <<R_ext_G2>> SPLITX 1 + MESHY 0.0 <<R_ext_G2>> SPLITY 1 + MIX 9 10 11 12 3 2 + NPIN 1 RPIN 0.0 APIN <<Angle>> + ; + ; + + + ::: C02 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 13 14 15 16 3 2 + ; + ; + + + ::: C03 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 41 42 43 44 3 2 + ; + ; + + + ::: C04 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 89 90 91 92 3 2 + ; + ; + + + ::: C05 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 93 94 95 96 3 2 + ; + ; + + + ::: C06 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 45 46 47 48 3 2 + ; + ; + + + ::: C07 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 49 50 51 52 3 2 + ; + ; + + + ::: C08 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 53 54 55 56 3 2 + ; + ; + + + ::: C09 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 97 98 99 100 3 2 + ; + ; + + + ::: C10 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 101 102 103 104 3 2 + ; + ; + + + ::: C11 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 57 58 59 60 3 2 + ; + ; + + + ::: C12 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 105 106 107 108 3 2 + ; + ; + + + ::: C13 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 61 62 63 64 3 2 + ; + ; + + + ::: C14 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 109 110 111 112 3 2 + ; + ; + + + ::: C15 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 113 114 115 116 3 2 + ; + ; + + + ::: C16 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 65 66 67 68 3 2 + ; + ; + + + ::: C17 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 69 70 71 72 3 2 + ; + ; + + + ::: C18 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 73 74 75 76 3 2 + ; + ; + + + ::: C19 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 77 78 79 80 3 2 + ; + ; + + + ::: C20 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 81 82 83 84 3 2 + ; + ; + + + ::: C21 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 85 86 87 88 3 2 + ; + ; + + + ::: C22 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 117 118 119 120 3 2 + ; + ; + + + ::: C23 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 33 34 35 36 3 2 + ; + ; + + + ::: C24 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 121 122 123 124 3 2 + ; + ; + + + ::: C25 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 125 126 127 128 3 2 + ; + ; + + + ::: C26 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 137 138 139 140 3 2 + ; + ; + + + ::: C27 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 129 130 131 132 3 2 + ; + ; + + + ::: C28 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 17 18 19 20 3 2 + ; + ; + + + ::: C29 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 133 134 135 136 3 2 + ; + ; + + + ::: C30 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 29 30 31 32 3 2 + ; + ; + + + ::: C31 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 21 22 23 24 3 2 + ; + ; + + + ::: C32 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 141 142 143 144 3 2 + ; + ; + + + ::: C33 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 145 146 147 148 3 2 + ; + ; + + + ::: C34 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 149 150 151 152 3 2 + ; + ; + + + ::: C35 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 153 154 155 156 3 2 + ; + ; + + + ::: C36 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 157 158 159 160 3 2 + ; + ; + + + ::: C37 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 161 162 163 164 3 2 + ; + ; + + + ::: C38 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 37 38 39 40 3 2 + ; + ; + + + ::: C39 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 25 26 27 28 3 2 + ; + ; + + + ::: T2 := GEO: CARCEL 1 1 1 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_ext_TG>> + MIX 1 6 + CLUSTER PIN + ::: PIN := GEO: TUBE 2 3 3 + RADIUS 0.0 <<R_int_TG>> <<R_ext_TG>> + SPLITR 3 1 + MESHX <<R_ext_TGm>> -0.26 0.26 <<R_ext_TG>> + MESHY <<R_ext_TGm>> -0.26 0.26 <<R_ext_TG>> + MIX 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + NPIN 1 RPIN 0.0 APIN <<Angle>> + ; + ; + + + ::: T1 := GEO: T2 + MIX 1 7 + ::: PIN := GEO: PIN + MIX 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + ; + ; + + + ::: TI := GEO: CARCEL 1 1 1 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_ext_TI>> + MIX 1 8 + CLUSTER PIN + ::: PIN := GEO: TUBE 2 3 3 + RADIUS 0.0 <<R_int_TI>> <<R_ext_TI>> + SPLITR 3 1 + MESHX <<R_ext_TIm>> -0.26 0.26 <<R_ext_TI>> + MESHY <<R_ext_TIm>> -0.26 0.26 <<R_ext_TI>> + MIX 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + NPIN 1 RPIN 0.0 APIN <<Angle>> + ; + ; + + + ::: L1 := GEO: CARCEL 0 1 3 + MESHX 0.0 <<Lame>> + MESHY 0.0 <<meshc1>> <<meshc2>> <<Cote>> + MIX 5 5 5 + ; + + + ::: L2 := GEO: CARCEL 0 1 1 + MESHX 0.0 <<Lame>> + MESHY 0.0 <<Lame>> + MIX 5 + ; +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Geo_N2_32.c2m b/Dragon/data/rep900_mco_proc/Geo_N2_32.c2m new file mode 100644 index 0000000..69b9e2f --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Geo_N2_32.c2m @@ -0,0 +1,797 @@ +*DECK Geo_N2 +*---- +* Name : Geo_N2.c2m +* Type : DRAGON procedure +* Use : Geometry generation for 1st Level Flux Calculation +* Author : R. Vallerent +* +* Procedure called as: +* +*GEON2 := Geo_N2 :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> +* <<R_int_TG>> <<R_ext_TG>> <<R_int_TI>> +* <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>> +* <<Lame>> <<RmodeN2>> +* <<R_ext_G2>> <<R_ext_TGm>> <<R_ext_TIm>> ; +* +* Input data : +* Rcomb1 : 50/100 of outer radius of fuel (cm) +* Rcomb2 : 80/100 of outer radius of fuel (cm) +* Rcomb3 : 95/100 of outer radius of fuel (cm) +* Rcomb4 : 100/100 of outer radius of fuel (cm) +* R_int_TG : Inner radius of cladding of guide tube (cm) +* R_ext_TG : Outer radius of cladding of guide tube (cm) +* R_int_TI : Inner radius of cladding of guide tube (cm) +* R_ext_TI : Outer radius of cladding of guide tube (cm) +* R_int_G : Inner radius of cladding of fuel tube (cm) +* R_ext_G : Outer radius of cladding of fuel tube (cm) +* Cote : Lattice pitch (cm) +* Lame : Water space on the side of the assembly (cm) +* RmodeN2 : Parameter for flux calculation in Level 2 (cm) +* R_ext_G2 : R_ext_G * 2 (cm) +* R_ext_TGm : -R_ext_TGm (cm) +* R_ext_TIm : -R_ext_TIm (cm) +* +* Output data : +* GEON2 : Geometry for Self-Shielding Calculation + + +PARAMETER GEON2 :: + EDIT 0 + ::: LINKED_LIST GEON2 ; + ; +*---- +* Modules used in this procedure +*---- +MODULE GEO: END: ; + +*---- +* Input data recovery +*---- +* ; +REAL Rcomb1 Rcomb2 Rcomb3 Rcomb4 ; +:: >>Rcomb1<< >>Rcomb2<< >>Rcomb3<< >>Rcomb4<< ; +REAL R_int_TG R_ext_TG R_int_TI R_ext_TI ; +:: >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< ; +REAL R_int_G R_ext_G Cote Lame ; +:: >>R_int_G<< >>R_ext_G<< >>Cote<< >>Lame<< ; +REAL RmodeN2 ; +:: >>RmodeN2<< ; + +REAL Angle := 0.785398163 ; +REAL meshc1 meshc2 := 0.37 0.89 ; + +REAL R_ext_G2 R_ext_TGm R_ext_TIm ; +:: >>R_ext_G2<< >>R_ext_TGm<< >>R_ext_TIm<< ; + +GEON2 := GEO: :: CAR2D 10 10 + EDIT 0 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + CELL TI C01 C03 T1 C08 C13 T2 C23 C31 L1 + C02 C04 C06 C09 C14 C18 C24 C32 L1 + C05 C07 C10 C15 C19 C25 C33 L1 + T1 C11 C16 T2 C26 C34 L1 + C12 C17 C20 C27 C35 L1 + T2 C21 C28 C36 L1 + C22 C29 C37 L1 + C30 C38 L1 + C39 L1 + L2 + + ::: C01 := GEO: CARCEL 2 3 3 + MESHX 0.0 <<meshc1>> <<meshc2>> <<Cote>> + MESHY 0.0 <<meshc1>> <<meshc2>> <<Cote>> + SPLITX 1 1 1 + SPLITY 1 1 1 + RADIUS 0.0 <<R_ext_G>> <<RmodeN2>> + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + CLUSTER PIN + ::: PIN := GEO: TUBE 6 + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MESHX 0.0 <<R_ext_G2>> SPLITX 1 + MESHY 0.0 <<R_ext_G2>> SPLITY 1 + MIX 9 10 11 12 3 2 + NPIN 1 RPIN 0.0 APIN <<Angle>> + ; + ; + + + ::: C02 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 13 14 15 16 3 2 + ; + ; + + + ::: C03 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 41 42 43 44 3 2 + ; + ; + + + ::: C04 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 89 90 91 92 3 2 + ; + ; + + + ::: C05 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 93 94 95 96 3 2 + ; + ; + + + ::: C06 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 45 46 47 48 3 2 + ; + ; + + + ::: C07 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 49 50 51 52 3 2 + ; + ; + + + ::: C08 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 53 54 55 56 3 2 + ; + ; + + + ::: C09 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 97 98 99 100 3 2 + ; + ; + + + ::: C10 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 101 102 103 104 3 2 + ; + ; + + + ::: C11 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 57 58 59 60 3 2 + ; + ; + + + ::: C12 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 105 106 107 108 3 2 + ; + ; + + + ::: C13 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 61 62 63 64 3 2 + ; + ; + + + ::: C14 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 109 110 111 112 3 2 + ; + ; + + + ::: C15 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 113 114 115 116 3 2 + ; + ; + + + ::: C16 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 65 66 67 68 3 2 + ; + ; + + + ::: C17 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 69 70 71 72 3 2 + ; + ; + + + ::: C18 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 73 74 75 76 3 2 + ; + ; + + + ::: C19 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 77 78 79 80 3 2 + ; + ; + + + ::: C20 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 81 82 83 84 3 2 + ; + ; + + + ::: C21 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 85 86 87 88 3 2 + ; + ; + + + ::: C22 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 117 118 119 120 3 2 + ; + ; + + + ::: C23 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 21 22 23 24 3 2 + ; + ; + + + ::: C24 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 121 122 123 124 3 2 + ; + ; + + + ::: C25 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 125 126 127 128 3 2 + ; + ; + + + ::: C26 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 137 138 139 140 3 2 + ; + ; + + + ::: C27 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 129 130 131 132 3 2 + ; + ; + + + ::: C28 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 17 18 19 20 3 2 + ; + ; + + + ::: C29 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 133 134 135 136 3 2 + ; + ; + + + ::: C30 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 25 26 27 28 3 2 + ; + ; + + + ::: C31 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 29 30 31 32 3 2 + ; + ; + + + ::: C32 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 141 142 143 144 3 2 + ; + ; + + + ::: C33 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 145 146 147 148 3 2 + ; + ; + + + ::: C34 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 149 150 151 152 3 2 + ; + ; + + + ::: C35 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 153 154 155 156 3 2 + ; + ; + + + ::: C36 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 157 158 159 160 3 2 + ; + ; + + + ::: C37 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 161 162 163 164 3 2 + ; + ; + + + ::: C38 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 33 34 35 36 3 2 + ; + ; + + + ::: C39 := GEO: C01 + MIX 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + 2 4 5 + ::: PIN := GEO: PIN + MIX 37 38 39 40 3 2 + ; + ; + + + ::: T2 := GEO: CARCEL 1 1 1 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_ext_TG>> + MIX 1 6 + CLUSTER PIN + ::: PIN := GEO: TUBE 2 3 3 + RADIUS 0.0 <<R_int_TG>> <<R_ext_TG>> + SPLITR 3 1 + MESHX <<R_ext_TGm>> -0.26 0.26 <<R_ext_TG>> + MESHY <<R_ext_TGm>> -0.26 0.26 <<R_ext_TG>> + MIX 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + 6 1 + NPIN 1 RPIN 0.0 APIN <<Angle>> + ; + ; + + + ::: T1 := GEO: T2 + MIX 1 7 + ::: PIN := GEO: PIN + MIX 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + 7 1 + ; + ; + + + ::: TI := GEO: CARCEL 1 1 1 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_ext_TI>> + MIX 1 8 + CLUSTER PIN + ::: PIN := GEO: TUBE 2 3 3 + RADIUS 0.0 <<R_int_TI>> <<R_ext_TI>> + SPLITR 3 1 + MESHX <<R_ext_TIm>> -0.26 0.26 <<R_ext_TI>> + MESHY <<R_ext_TIm>> -0.26 0.26 <<R_ext_TI>> + MIX 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + 8 1 + NPIN 1 RPIN 0.0 APIN <<Angle>> + ; + ; + + + ::: L1 := GEO: CARCEL 0 1 3 + MESHX 0.0 <<Lame>> + MESHY 0.0 <<meshc1>> <<meshc2>> <<Cote>> + MIX 5 5 5 + ; + + + ::: L2 := GEO: CARCEL 0 1 1 + MESHX 0.0 <<Lame>> + MESHY 0.0 <<Lame>> + MIX 5 + ; +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Geo_SS_20.c2m b/Dragon/data/rep900_mco_proc/Geo_SS_20.c2m new file mode 100644 index 0000000..ca5c9ad --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Geo_SS_20.c2m @@ -0,0 +1,138 @@ +*DECK Geo_SS +*---- +* Name : Geo_SS.c2m +* Type : DRAGON procedure +* Use : Geometry generation for Self-Shielding Calculation +* with 20 fuel regions +* Author : R. Vallerent +* +* Procedure called as: +* +*GEOSS := Geo_SS :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> +* <<R_int_TG>> <<R_ext_TG>> <<R_int_TI>> +* <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>> +* <<CoteL>> ; +* +* Input data : +* Rcomb1 : 50/100 of outer radius of fuel (cm) +* Rcomb2 : 80/100 of outer radius of fuel (cm) +* Rcomb3 : 95/100 of outer radius of fuel (cm) +* Rcomb4 : 100/100 of outer radius of fuel (cm) +* R_int_TG : Inner radius of cladding of guide tube (cm) +* R_ext_TG : Outer radius of cladding of guide tube (cm) +* R_int_TI : Inner radius of cladding of guide tube (cm) +* R_ext_TI : Outer radius of cladding of guide tube (cm) +* R_int_G : Inner radius of cladding of fuel tube (cm) +* R_ext_G : Outer radius of cladding of fuel tube (cm) +* Cote : Lattice pitch (cm) +* CoteL : Lattice pitch + Water space (cm) +* +* Output data : +* GEOSS : Geometry for Self-Shielding Calculation + + +PARAMETER GEOSS :: + EDIT 0 + ::: LINKED_LIST GEOSS ; + ; +*---- +* Modules used in this procedure +*---- +MODULE GEO: END: ; + +*---- +* Input data recovery +*---- +* ; +REAL Rcomb1 Rcomb2 Rcomb3 Rcomb4 ; +:: >>Rcomb1<< >>Rcomb2<< >>Rcomb3<< >>Rcomb4<< ; +REAL R_int_TG R_ext_TG R_int_TI R_ext_TI ; +:: >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< ; +REAL R_int_G R_ext_G Cote CoteL ; +:: >>R_int_G<< >>R_ext_G<< >>Cote<< >>CoteL<< ; + + + +GEOSS := GEO: :: CAR2D 9 9 + EDIT 0 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + CELL TI C1 C1 T1 C1 C1 T2 C1 C4 + C2 C2 C1 C2 C2 C1 C2 C4 + C2 C1 C2 C2 C1 C2 C4 + T1 C1 C1 T2 C1 C4 + C2 C1 C1 C2 C4 + T2 C1 C3 C4 + C2 C3 C4 + C3 C4 + C5 + + MERGE 8 1 1 7 1 1 6 1 4 + 2 2 1 2 2 1 2 4 + 2 1 2 2 1 2 4 + 7 1 1 6 1 4 + 2 1 1 2 4 + 6 1 3 4 + 2 3 4 + 3 4 + 5 + +* T2 -> 6 , T1 -> 7 , TI -> 8 + + TURN A A E A A E A A A + A E F A E D A A + C B G C H G A + A G C A G A + A B D A A + A A A A + A G A + A A + A + + + ::: C1 := GEO: CARCEL 6 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MIX 9 10 11 12 3 2 4 + ; + ::: C2 := GEO: C1 + MIX 13 14 15 16 3 2 4 + ; + ::: C3 := GEO: C1 + MIX 17 18 19 20 3 2 4 + ; + ::: C4 := GEO: CARCEL 6 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MIX 21 22 23 24 3 2 4 + ; + ::: C5 := GEO: CARCEL 6 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<CoteL>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MIX 25 26 27 28 3 2 4 + ; + ::: T2 := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TG>> <<R_ext_TG>> + MIX 6 1 6 + ; + ::: T1 := GEO: T2 + MIX 7 1 7 + ; + ::: TI := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TI>> <<R_ext_TI>> + MIX 8 1 8 + ; +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Geo_SS_32.c2m b/Dragon/data/rep900_mco_proc/Geo_SS_32.c2m new file mode 100644 index 0000000..cc8b9cd --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Geo_SS_32.c2m @@ -0,0 +1,147 @@ +*DECK Geo_SS +*---- +* Name : Geo_SS.c2m +* Type : DRAGON procedure +* Use : Geometry generation for Self-Shielding Calculation +* with 32 fuel regions +* Author : R. Vallerent +* +* Procedure called as: +* +*GEOSS := Geo_SS :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> +* <<R_int_TG>> <<R_ext_TG>> <<R_int_TI>> +* <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>> +* <<CoteL>> ; +* +* Input data : +* Rcomb1 : 50/100 of outer radius of fuel (cm) +* Rcomb2 : 80/100 of outer radius of fuel (cm) +* Rcomb3 : 95/100 of outer radius of fuel (cm) +* Rcomb4 : 100/100 of outer radius of fuel (cm) +* R_int_TG : Inner radius of cladding of guide tube (cm) +* R_ext_TG : Outer radius of cladding of guide tube (cm) +* R_int_TI : Inner radius of cladding of guide tube (cm) +* R_ext_TI : Outer radius of cladding of guide tube (cm) +* R_int_G : Inner radius of cladding of fuel tube (cm) +* R_ext_G : Outer radius of cladding of fuel tube (cm) +* Cote : Lattice pitch (cm) +* CoteL : Lattice pitch + Water space (cm) +* +* Output data : +* GEOSS : Geometry for Self-Shielding Calculation + + +PARAMETER GEOSS :: + EDIT 0 + ::: LINKED_LIST GEOSS ; + ; +*---- +* Modules used in this procedure +*---- +MODULE GEO: END: ; + +*---- +* Input data recovery +*---- +* ; +REAL Rcomb1 Rcomb2 Rcomb3 Rcomb4 ; +:: >>Rcomb1<< >>Rcomb2<< >>Rcomb3<< >>Rcomb4<< ; +REAL R_int_TG R_ext_TG R_int_TI R_ext_TI ; +:: >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< ; +REAL R_int_G R_ext_G Cote CoteL ; +:: >>R_int_G<< >>R_ext_G<< >>Cote<< >>CoteL<< ; + + + +GEOSS := GEO: :: CAR2D 9 9 + EDIT 0 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + CELL TI C1 C1 T1 C1 C1 T2 C4 C6 + C2 C2 C1 C2 C2 C1 C2 C6 + C2 C1 C2 C2 C1 C2 C6 + T1 C1 C1 T2 C4 C6 + C2 C1 C1 C2 C6 + T2 C1 C3 C6 + C2 C3 C6 + C5 C7 + C8 + + MERGE 11 1 1 10 1 1 9 4 6 + 2 2 1 2 2 1 2 6 + 2 1 2 2 1 2 6 + 10 1 1 9 4 6 + 2 1 1 2 6 + 9 1 3 6 + 2 3 6 + 5 7 + 8 + +* T2 -> 9 , T1 -> 10 , TI -> 11 + + TURN A A E A A E A A A + A E F A E D A A + C B G C H G A + A G C A G A + A B D A A + A A A A + A G A + A A + A + + + ::: C1 := GEO: CARCEL 6 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MIX 9 10 11 12 3 2 4 + ; + ::: C2 := GEO: C1 + MIX 13 14 15 16 3 2 4 + ; + ::: C3 := GEO: C1 + MIX 17 18 19 20 3 2 4 + ; + ::: C4 := GEO: C1 + MIX 21 22 23 24 3 2 4 + ; + ::: C5 := GEO: C1 + MIX 25 26 27 28 3 2 4 + ; + ::: C6 := GEO: CARCEL 6 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MIX 29 30 31 32 3 2 4 + ; + ::: C7 := GEO: C6 + MIX 33 34 35 36 3 2 4 + ; + ::: C8 := GEO: CARCEL 6 + MESHX 0.0 <<CoteL>> + MESHY 0.0 <<CoteL>> + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> + <<R_int_G>> <<R_ext_G>> + MIX 37 38 39 40 3 2 4 + ; + ::: T2 := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TG>> <<R_ext_TG>> + MIX 6 1 6 + ; + ::: T1 := GEO: T2 + MIX 7 1 7 + ; + ::: TI := GEO: CARCEL 2 + MESHX 0.0 <<Cote>> + MESHY 0.0 <<Cote>> + RADIUS 0.0 <<R_int_TI>> <<R_ext_TI>> + MIX 8 1 8 + ; +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/GetDim.c2m b/Dragon/data/rep900_mco_proc/GetDim.c2m new file mode 100644 index 0000000..9f31f4f --- /dev/null +++ b/Dragon/data/rep900_mco_proc/GetDim.c2m @@ -0,0 +1,120 @@ +*DECK GetDim.c2m +*---- +* Name : GetDim.c2m +* Type : DRAGON procedure +* Use : Public dimensions for assembly calculation +* Author : R. Vallerent +* +* Procedure called as: +* +* GetDim :: +* >>Rcomb1<< >>Rcomb2<< >>Rcomb3<< >>Rcomb4<< +* >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< +* >>R_int_G<< >>R_ext_G<< >>Cote<< >>Lame<< >>CoteL<< +* >>RmodeN1<< >>RmodeN2<< >>R_ext_G2<< >>R_ext_TGm<< >>R_ext_TIm<< +* >>N_U5_UOX<< >>N_U5_MOX<< >>N_Pu8<< >>N_Pu9<< >>N_Pu0<< +* >>N_Pu1<< >>N_Pu2<< >>N_Am1<< >>densU_UOX<< >>densU_MOX<< +* >>densPu_MOX<< >>teneur1_Pu<< >>teneur2_Pu<< >>teneur3_Pu<< +* >>temp_comb<< >>temp_mode<< >>dens_mod<< ; +* +* Input data : none +* +* Output data: +* Rcomb1 : 50/100 of outer radius of fuel (cm) +* Rcomb2 : 80/100 of outer radius of fuel (cm) +* Rcomb3 : 95/100 of outer radius of fuel (cm) +* Rcomb4 : 100/100 of outer radius of fuel (cm) +* R_int_TG : Inner radius of cladding of guide tube (cm) +* R_ext_TG : Outer radius of cladding of guide tube (cm) +* R_int_TI : Inner radius of cladding of guide tube (cm) +* R_ext_TI : Outer radius of cladding of guide tube (cm) +* R_int_G : Inner radius of cladding of fuel tube (cm) +* R_ext_G : Outer radius of cladding of fuel tube (cm) +* Cote : Lattice pitch (cm) +* Lame : Water space on the side of the assembly (cm) +* CoteL : Lattice pitch + Water space (cm) +* RmodeN1 : Parameter for flux calculation in Level 1 (cm) +* RmodeN2 : Parameter for flux calculation in Level 2 (cm) +* R_ext_G2 : R_ext_G * 2 (cm) +* R_ext_TGm : -R_ext_TGm (cm) +* R_ext_TIm : -R_ext_TIm (cm) +* N_U5_UOX : U235 enrichment in UOX-Fuel +* N_U5_MOX : U235 enrichment in MOX-Fuel +* N_Pu8 : Pu238 percentage in Pu-vector +* N_Pu9 : Pu239 percentage in Pu-vector +* N_Pu0 : Pu240 percentage in Pu-vector +* N_Pu1 : Pu241 percentage in Pu-vector +* N_Pu2 : Pu242 percentage in Pu-vector +* N_Am1 : Am241 percentage in Pu-vector +* densU_UOX : Theoretical density of Uranium in UOX-Fuel +* densU_MOX : Theoretical density of Uranium in MOX-Fuel +* densPu_MOX : Theoretical density of Plutonium in MOX-Fuel +* teneur1_Pu : Plutonium content in MOX-high enrichment region +* teneur2_Pu : Plutonium content in MOX-medium enrichment region +* teneur3_Pu : Plutonium content in MOX-low enrichment region +* temp_comb : Theoretical Fuel temperature (K) +* temp_mode : Theoretical moderator temperature (K) +* dens_mod : Theoretical moderator density (g/cc) + +MODULE + END: ; + +REAL Rcomb1 Rcomb2 Rcomb3 Rcomb4 + R_int_TG R_ext_TG R_int_TI R_ext_TI R_int_G R_ext_G + Cote Lame CoteL + RmodeN1 RmodeN2 + R_ext_G2 R_ext_TGm R_ext_TIm + N_U5_UOX N_U5_MOX + N_Pu8 N_Pu9 N_Pu0 N_Pu1 N_Pu2 N_Am1 + densU_UOX densU_MOX densPu_MOX + teneur1_Pu teneur2_Pu teneur3_Pu + temp_comb temp_mode dens_mod ; + + +EVALUATE Rcomb4 := 0.41 ; +EVALUATE Rcomb1 := 0.5 SQRT Rcomb4 * ; +EVALUATE Rcomb2 := 0.8 SQRT Rcomb4 * ; +EVALUATE Rcomb3 := 0.95 SQRT Rcomb4 * ; + +EVALUATE R_int_TG := 0.56 ; +EVALUATE R_ext_TG := 0.62 ; +EVALUATE R_int_TI := 0.56 ; +EVALUATE R_ext_TI := 0.62 ; +EVALUATE R_int_G := 0.418 ; +EVALUATE R_ext_G := 0.48 ; + +EVALUATE Cote := 1.26 ; +EVALUATE Lame := 0.04 ; +EVALUATE CoteL := Cote Lame + ; + +EVALUATE RmodeN1 := 0.60 ; +EVALUATE RmodeN2 := 0.627 ; + +EVALUATE R_ext_G2 := R_ext_G 2.0 * ; +EVALUATE R_ext_TGm := R_ext_TG -1.0 * ; +EVALUATE R_ext_TIm := R_ext_TI -1.0 * ; + +EVALUATE N_U5_UOX := 3.70E-02 ; +EVALUATE N_U5_MOX := 0.25E-02 ; + +EVALUATE N_Pu8 N_Pu9 N_Pu0 N_Pu1 N_Pu2 N_Am1 := + 1.20E-02 59.50E-02 24.10E-02 9.00E-02 4.70E-02 1.50E-02 ; + +EVALUATE densU_UOX densU_MOX densPu_MOX := 10.0 10.0 10.5 ; + +EVALUATE teneur1_Pu teneur2_Pu teneur3_Pu := + 7.50E-02 4.90E-02 2.80E-02 ; + +EVALUATE temp_comb temp_mode dens_mod := 800.0 600.0 0.659 ; + + +:: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> <<R_int_TG>> <<R_ext_TG>> + <<R_int_TI>> <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>> <<Lame>> + <<CoteL>> <<RmodeN1>> <<RmodeN2>> <<R_ext_G2>> <<R_ext_TGm>> + <<R_ext_TIm>> <<N_U5_UOX>> <<N_U5_MOX>> <<N_Pu8>> <<N_Pu9>> <<N_Pu0>> + <<N_Pu1>> <<N_Pu2>> <<N_Am1>> <<densU_UOX>> <<densU_MOX>> + <<densPu_MOX>> <<teneur1_Pu>> <<teneur2_Pu>> <<teneur3_Pu>> + <<temp_comb>> <<temp_mode>> <<dens_mod>> ; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Homogenize.c2m b/Dragon/data/rep900_mco_proc/Homogenize.c2m new file mode 100644 index 0000000..52140c4 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Homogenize.c2m @@ -0,0 +1,82 @@ +*DECK Homogenize +*---- +* Name : Homogenize.c2m +* Type : DRAGON procedure +* Use : Update of the LIBRARY2 concentrations +* with the ones of LIBHOM +* Author : R. Vallerent and A. Hebert +* +* Procedure called as: +* +* EDIHOM := Homogenize FLUX2 LIBEQ TRACKN2 :: <<NbRegions>> ; +* +* Input data : +* FLUX2 : Second level Flux +* LIBEQ : 26 groups library for second level calculation +* TRACKN2 : Tracking of the second level geometry +* NbRegions : Number of fuel regions for the calculation (20/32) +* Output data : +* LIBHOM : Microlib object with homogenized concentrations + + +PARAMETER LIBHOM FLUX2 LIBEQ TRACKN2 :: + EDIT 1 + ::: LINKED_LIST LIBHOM FLUX2 LIBEQ TRACKN2 ; + ; +* ; +MODULE EDI: DELETE: END: ; +LINKED_LIST EDIHOM ; + +INTEGER NbRegions ; + :: >>NbRegions<< ; + +IF NbRegions 20 = THEN + EDIHOM := EDI: FLUX2 LIBEQ TRACKN2 :: + EDIT 0 + MICR ALL + COND + MERG MIX + 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 + 16 17 18 19 20 21 22 23 24 25 26 27 28 17 18 + 19 20 9 10 11 12 21 22 23 24 9 10 11 12 9 + 10 11 12 9 10 11 12 9 10 11 12 9 10 11 12 + 9 10 11 12 9 10 11 12 9 10 11 12 9 10 11 + 12 9 10 11 12 9 10 11 12 9 10 11 12 13 14 + 15 16 13 14 15 16 13 14 15 16 13 14 15 16 13 + 14 15 16 13 14 15 16 13 14 15 16 13 14 15 16 + 13 14 15 16 13 14 15 16 13 14 15 16 17 18 19 + 20 9 10 11 12 21 22 23 24 21 22 23 24 21 22 + 23 24 21 22 23 24 21 22 23 24 21 22 23 24 + SAVE ON HOMOGENE + ; +ELSEIF NbRegions 32 = THEN + EDIHOM := EDI: FLUX2 LIBEQ TRACKN2 :: + EDIT 0 + MICR ALL + COND + MERG MIX + 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 + 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 + 31 32 33 34 35 36 37 38 39 40 9 10 11 12 9 + 10 11 12 9 10 11 12 9 10 11 12 9 10 11 12 + 9 10 11 12 9 10 11 12 9 10 11 12 9 10 11 + 12 9 10 11 12 9 10 11 12 9 10 11 12 13 14 + 15 16 13 14 15 16 13 14 15 16 13 14 15 16 13 + 14 15 16 13 14 15 16 13 14 15 16 13 14 15 16 + 13 14 15 16 13 14 15 16 13 14 15 16 17 18 19 + 20 21 22 23 24 29 30 31 32 29 30 31 32 29 30 + 31 32 29 30 31 32 29 30 31 32 29 30 31 32 + SAVE ON HOMOGENE + ; +ELSE + ECHO "Wrong value for NbRegions" ; +ENDIF ; + +LIBHOM := EDIHOM :: + STEP UP HOMOGENE +; + +EDIHOM := DELETE: EDIHOM ; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Mix_MOX_32.c2m b/Dragon/data/rep900_mco_proc/Mix_MOX_32.c2m new file mode 100644 index 0000000..df76323 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Mix_MOX_32.c2m @@ -0,0 +1,453 @@ +*DECK Mix_Assb_DRA +*---- +* Name : Mix_Assb_DRA.c2m +* Type : DRAGON procedure +* Use : Microlib generation with Draglibs for +* MOX calculation with 32 fuel regions +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBRARY := Mix_Assb_DRA :: +* <<dens_mod>> <<pbore>> <<N_U5_MOX>> <<N_Pu8>> <<N_Pu9>> <<N_Pu0>> +* <<N_Pu1>> <<N_Pu2>> <<N_Am1>> <<densU_MOX>> <<densPu_MOX>> +* <<teneur1_Pu>> <<teneur2_Pu>> <<teneur3_Pu>> <<temp_comb>> +* <<temp_mode>> ; +* ; +* +* Input data : +* dens_mod : Theoretical moderator density +* pbore : Natural boron concentration (ppm) +* N_U5_MOX : U235 enrichment in MOX-Fuel +* N_Pu8 : Pu238 percentage in Pu-vector +* N_Pu9 : Pu239 percentage in Pu-vector +* N_Pu0 : Pu240 percentage in Pu-vector +* N_Pu1 : Pu241 percentage in Pu-vector +* N_Pu2 : Pu242 percentage in Pu-vector +* N_Am1 : Am241 percentage in Pu-vector +* densU_MOX : Theoretical density of Uranium in MOX-Fuel +* densPu_MOX : Theoretical density of Plutonium in MOX-Fuel +* teneur1_Pu : Plutonium content in MOX-high enrichment region +* teneur2_Pu : Plutonium content in MOX-medium enrichment region +* teneur3_Pu : Plutonium content in MOX-low enrichment region +* temp_comb : Theoretical Fuel temperature (K) +* temp_mode : Theoretical moderator temperature (K) +* +* Output data : +* LIBRARY : Microlib with table of probabilities + + +PARAMETER LIBRARY :: + EDIT 0 + ::: LINKED_LIST LIBRARY ; + ; + +*---- +* Input data recovery +*---- +STRING Library ; +:: >>Library<< ; + +REAL dens_mod pbore N_U5_MOX N_Pu8 N_Pu9 N_Pu0 + N_Pu1 N_Pu2 N_Am1 densU_MOX densPu_MOX + teneur1_Pu teneur2_Pu teneur3_Pu temp_comb + temp_mode ; +:: >>dens_mod<< >>pbore<< >>N_U5_MOX<< >>N_Pu8<< >>N_Pu9<< >>N_Pu0<< + >>N_Pu1<< >>N_Pu2<< >>N_Am1<< >>densU_MOX<< >>densPu_MOX<< + >>teneur1_Pu<< >>teneur2_Pu<< >>teneur3_Pu<< >>temp_comb<< + >>temp_mode<< ; + +*---- +* Modules used in this procedure +*---- +MODULE LIB: UTL: DELETE: END: ; + +REAL TfuelK := temp_comb ; +REAL TmodeK := temp_mode ; + +REAL TgainK := TmodeK ; + +* Calculation of moderator data +* NOMBRE D AVOGADRO +REAL A := 6.022094E-1 ; ! Normalisation des concentrations par 10E-24 + +REAL M_H2O := 15.9994 2.0 1.00794 * + ; +REAL N_H2O := dens_mod 0.6022 * M_H2O / ; +REAL b10 b11 := 0.199 0.801 ; +REAL coeff := 0.199 A 1.0E-06 * * 10.8110164 / ; +REAL N_B10 := coeff 0.659 pbore * * ; ! densite nominale +REAL N_2H2O N_B11 := 2.0 N_H2O * N_B10 b11 * b10 / ; + +* Data for all enrichment MOX +REAL den_uo2t den_puo2t := densU_MOX densPu_MOX ; + +REAL N_U5 := N_U5_MOX ; +REAL N_U4 := N_U5 54.0E-04 * 0.72 / ; +REAL N_U6 := 1.00E-14 ; +REAL N_U8 := 1. N_U5 - N_U4 - N_U6 - ; + +REAL N_Np7 := 1.00E-14 ; + +REAL M_U4 M_U5 M_U8 M_U6 := + 234.040952088 235.043929918 238.050788247 236.045568006 ; +REAL M_Pu8 M_Pu9 M_Pu0 M_Pu1 := + 238.049559894 239.052163381 240.053813545 241.056851456 ; +REAL M_Pu2 M_Am1 M_Np7 M_O16 := + 242.058742611 241.056829144 237.048173444 15.99491461956 ; +*taken from : http://amdc.in2p3.fr/masstables/Ame2003/mass.mas03 +* ; +*---- +* MOX-high enrichment region +*---- +REAL teneur1_U := 1. teneur1_Pu - ; +REAL xu1 xp1 := teneur1_U den_uo2t / teneur1_Pu den_puo2t / ; +REAL den_mox1 := 1. xu1 xp1 + / ; + +REAL X1_U4 X1_U5 := N_U4 teneur1_U * N_U5 teneur1_U * ; +REAL X1_U6 X1_U8 := N_U6 teneur1_U * N_U8 teneur1_U * ; +REAL X1_Pu8 X1_Pu9 := N_Pu8 teneur1_Pu * N_Pu9 teneur1_Pu * ; +REAL X1_Pu0 X1_Pu1 := N_Pu0 teneur1_Pu * N_Pu1 teneur1_Pu * ; +REAL X1_Pu2 X1_Am1 := N_Pu2 teneur1_Pu * N_Am1 teneur1_Pu * ; +REAL X1_Np7 := N_Np7 teneur1_Pu * ; + +REAL X1_O16 := X1_U4 M_U4 / X1_U5 M_U5 / + ; +EVALUATE X1_O16 := X1_O16 X1_U8 M_U8 / + X1_U6 M_U6 / + ; +EVALUATE X1_O16 := X1_O16 X1_Pu8 M_Pu8 / + X1_Pu9 M_Pu9 / + ; +EVALUATE X1_O16 := X1_O16 X1_Pu0 M_Pu0 / + X1_Pu1 M_Pu1 / + ; +EVALUATE X1_O16 := X1_O16 X1_Pu2 M_Pu2 / + X1_Am1 M_Am1 / + ; +EVALUATE X1_O16 := X1_O16 X1_Np7 M_Np7 / + 2.0 * M_O16 * ; + +REAL X1_TOT := X1_U4 X1_U5 + X1_U6 + X1_U8 + ; +EVALUATE X1_TOT := X1_TOT X1_Pu8 + X1_Pu9 + X1_Pu0 + ; +EVALUATE X1_TOT := X1_TOT X1_Pu1 + X1_Pu2 + X1_Am1 + X1_Np7 + X1_O16 + ; + +EVALUATE X1_U4 X1_U5 := X1_U4 X1_TOT / X1_U5 X1_TOT / ; +EVALUATE X1_U6 X1_U8 := X1_U6 X1_TOT / X1_U8 X1_TOT / ; +EVALUATE X1_Pu8 X1_Pu9 := X1_Pu8 X1_TOT / X1_Pu9 X1_TOT / ; +EVALUATE X1_Pu0 X1_Pu1 := X1_Pu0 X1_TOT / X1_Pu1 X1_TOT / ; +EVALUATE X1_Pu2 X1_Am1 := X1_Pu2 X1_TOT / X1_Am1 X1_TOT / ; +EVALUATE X1_Np7 X1_O16 := X1_Np7 X1_TOT / X1_O16 X1_TOT / ; + +*---- +* MOX-medium enrichment region +*---- +REAL teneur2_U := 1. teneur2_Pu - ; +REAL xu2 xp2 := teneur2_U den_uo2t / teneur2_Pu den_puo2t / ; +REAL den_mox2 := 1. xu2 xp2 + / ; + +REAL X2_U4 X2_U5 := N_U4 teneur2_U * N_U5 teneur2_U * ; +REAL X2_U6 X2_U8 := N_U6 teneur2_U * N_U8 teneur2_U * ; +REAL X2_Pu8 X2_Pu9 := N_Pu8 teneur2_Pu * N_Pu9 teneur2_Pu * ; +REAL X2_Pu0 X2_Pu1 := N_Pu0 teneur2_Pu * N_Pu1 teneur2_Pu * ; +REAL X2_Pu2 X2_Am1 := N_Pu2 teneur2_Pu * N_Am1 teneur2_Pu * ; +REAL X2_Np7 := N_Np7 teneur2_Pu * ; + +REAL X2_O16 := X2_U4 M_U4 / X2_U5 M_U5 / + ; +EVALUATE X2_O16 := X2_O16 X2_U8 M_U8 / + X2_U6 M_U6 / + ; +EVALUATE X2_O16 := X2_O16 X2_Pu8 M_Pu8 / + X2_Pu9 M_Pu9 / + ; +EVALUATE X2_O16 := X2_O16 X2_Pu0 M_Pu0 / + X2_Pu1 M_Pu1 / + ; +EVALUATE X2_O16 := X2_O16 X2_Pu2 M_Pu2 / + X2_Am1 M_Am1 / + ; +EVALUATE X2_O16 := X2_O16 X2_Np7 M_Np7 / + 2.0 * M_O16 * ; + +REAL X2_TOT := X2_U4 X2_U5 + X2_U6 + X2_U8 + ; +EVALUATE X2_TOT := X2_TOT X2_Pu8 + X2_Pu9 + X2_Pu0 + ; +EVALUATE X2_TOT := X2_TOT X2_Pu1 + X2_Pu2 + X2_Am1 + X2_Np7 + X2_O16 + ; + +EVALUATE X2_U4 X2_U5 := X2_U4 X2_TOT / X2_U5 X2_TOT / ; +EVALUATE X2_U6 X2_U8 := X2_U6 X2_TOT / X2_U8 X2_TOT / ; +EVALUATE X2_Pu8 X2_Pu9 := X2_Pu8 X2_TOT / X2_Pu9 X2_TOT / ; +EVALUATE X2_Pu0 X2_Pu1 := X2_Pu0 X2_TOT / X2_Pu1 X2_TOT / ; +EVALUATE X2_Pu2 X2_Am1 := X2_Pu2 X2_TOT / X2_Am1 X2_TOT / ; +EVALUATE X2_Np7 X2_O16 := X2_Np7 X2_TOT / X2_O16 X2_TOT / ; + + +*---- +* MOX-low enrichment region +*---- +REAL teneur3_U := 1. teneur3_Pu - ; +REAL xu3 xp3 := teneur3_U den_uo2t / teneur3_Pu den_puo2t / ; +REAL den_mox3 := 1. xu3 xp3 + / ; + +REAL X3_U4 X3_U5 := N_U4 teneur3_U * N_U5 teneur3_U * ; +REAL X3_U6 X3_U8 := N_U6 teneur3_U * N_U8 teneur3_U * ; +REAL X3_Pu8 X3_Pu9 := N_Pu8 teneur3_Pu * N_Pu9 teneur3_Pu * ; +REAL X3_Pu0 X3_Pu1 := N_Pu0 teneur3_Pu * N_Pu1 teneur3_Pu * ; +REAL X3_Pu2 X3_Am1 := N_Pu2 teneur3_Pu * N_Am1 teneur3_Pu * ; +REAL X3_Np7 := N_Np7 teneur3_Pu * ; + +REAL X3_O16 := X3_U4 M_U4 / X3_U5 M_U5 / + ; +EVALUATE X3_O16 := X3_O16 X3_U8 M_U8 / + X3_U6 M_U6 / + ; +EVALUATE X3_O16 := X3_O16 X3_Pu8 M_Pu8 / + X3_Pu9 M_Pu9 / + ; +EVALUATE X3_O16 := X3_O16 X3_Pu0 M_Pu0 / + X3_Pu1 M_Pu1 / + ; +EVALUATE X3_O16 := X3_O16 X3_Pu2 M_Pu2 / + X3_Am1 M_Am1 / + ; +EVALUATE X3_O16 := X3_O16 X3_Np7 M_Np7 / + 2.0 * M_O16 * ; + +REAL X3_TOT := X3_U4 X3_U5 + X3_U6 + X3_U8 + ; +EVALUATE X3_TOT := X3_TOT X3_Pu8 + X3_Pu9 + X3_Pu0 + ; +EVALUATE X3_TOT := X3_TOT X3_Pu1 + X3_Pu2 + X3_Am1 + X3_Np7 + X3_O16 + ; + +EVALUATE X3_U4 X3_U5 := X3_U4 X3_TOT / X3_U5 X3_TOT / ; +EVALUATE X3_U6 X3_U8 := X3_U6 X3_TOT / X3_U8 X3_TOT / ; +EVALUATE X3_Pu8 X3_Pu9 := X3_Pu8 X3_TOT / X3_Pu9 X3_TOT / ; +EVALUATE X3_Pu0 X3_Pu1 := X3_Pu0 X3_TOT / X3_Pu1 X3_TOT / ; +EVALUATE X3_Pu2 X3_Am1 := X3_Pu2 X3_TOT / X3_Am1 X3_TOT / ; +EVALUATE X3_Np7 X3_O16 := X3_Np7 X3_TOT / X3_O16 X3_TOT / ; + +* -------------------------------- +* Cladding calculation +* -------------------------------- +REAL P_Zr P_Fe P_Cr P_O16 P_Sn := + 98.115E-02 0.210E-02 0.100E-02 0.125E-02 1.450E-02 ; + +REAL P_Zr90 P_Zr91 P_Zr92 P_Zr94 P_Zr96 := + 51.45E-02 11.22E-02 17.15E-02 17.38E-02 2.80E-02 ; +REAL P_Fe54 P_Fe56 P_Fe57 P_Fe58 := + 5.602E-02 91.865E-02 2.243E-02 0.290E-02 ; +*isot 5.8 % 91.72% 2.20% 0.28% +REAL P_Cr50 P_Cr52 P_Cr53 P_Cr54 := + 4.178E-02 83.701E-02 9.673E-02 2.448E-02 ; +*isot 4.35% 83.79% 9.50% 2.36% +REAL P_Sn112 P_Sn114 P_Sn115 P_Sn116 P_Sn117 + P_Sn118 P_Sn119 P_Sn120 P_Sn122 P_Sn124 := + 0.92E-02 0.62E-02 0.33E-02 14.21E-02 7.56E-02 + 24.06E-02 8.59E-02 32.92E-02 4.75E-02 6.04E-02 ; + +REAL M_Zr := P_Zr90 89.904711 * P_Zr91 90.905643 * P_Zr92 91.905039 * + P_Zr94 93.906320 * P_Zr96 95.90829 * + + + + ; +REAL M_Fe := P_Fe54 53.9396148 * P_Fe56 55.9349421 * + P_Fe57 56.9353987 * P_Fe58 57.9332805 * + + + ; +REAL M_Cr := P_Cr50 49.9460496 * P_Cr52 51.9405119 * + P_Cr53 52.9406538 * P_Cr54 53.9388849 * + + + ; +REAL M_Sn := P_Sn112 111.90483 * P_Sn114 113.90278 * P_Sn115 114.90335 * + P_Sn116 115.901748 * P_Sn117 116.902961 * + P_Sn118 117.901613 * P_Sn119 118.903316 * + P_Sn120 119.902207 * P_Sn122 121.903451 * + P_Sn124 123.905283 * + + + + + + + + + ; + +REAL rho_clad := 6.55 ; +REAL N_Zr_clad := P_Zr rho_clad A * M_Zr / * ; +REAL N_Zr90 N_Zr91 N_Zr92 N_Zr94 N_Zr96 := + N_Zr_clad P_Zr90 * N_Zr_clad P_Zr91 * N_Zr_clad P_Zr92 * + N_Zr_clad P_Zr94 * N_Zr_clad P_Zr96 * ; + +REAL N_Fe_clad := P_Fe rho_clad A * M_Fe / * ; +REAL N_Fe54 N_Fe56 N_Fe57 N_Fe58 := + N_Fe_clad P_Fe54 * N_Fe_clad P_Fe56 * N_Fe_clad P_Fe57 * + N_Fe_clad P_Fe58 * ; + +REAL N_Cr_clad := P_Cr rho_clad A * M_Cr / * ; +REAL N_Cr50 N_Cr52 N_Cr53 N_Cr54 := + N_Cr_clad P_Cr50 * N_Cr_clad P_Cr52 * N_Cr_clad P_Cr53 * + N_Cr_clad P_Cr54 * ; + +REAL N_O16G := P_O16 rho_clad A * 15.9949150 / * ; + +REAL N_Sn_clad := P_Sn rho_clad A * M_Sn / * ; +REAL N_Sn112 N_Sn114 N_Sn115 N_Sn116 + N_Sn117 N_Sn118 N_Sn119 N_Sn120 + N_Sn122 N_Sn124 := + N_Sn_clad P_Sn112 * N_Sn_clad P_Sn114 * N_Sn_clad P_Sn115 * + N_Sn_clad P_Sn116 * N_Sn_clad P_Sn117 * N_Sn_clad P_Sn118 * + N_Sn_clad P_Sn119 * N_Sn_clad P_Sn120 * N_Sn_clad P_Sn122 * + N_Sn_clad P_Sn124 * ; + +* ; +LIBRARY := LIB: :: + EDIT 0 + DEPL LIB: DRAGON FIL: <<Library>> + NMIX 40 CTRA APOL + SUBG + MIXS LIB: DRAGON FIL: <<Library>> + +*---- +* Cladding TG et TI +*---- + MIX 1 <<TgainK>> NOEV + Zr90g = Zr90 <<N_Zr90>> 2 + Zr91g = Zr91 <<N_Zr91>> 2 + Zr92g = Zr92 <<N_Zr92>> 2 + Zr94g = Zr94 <<N_Zr94>> 2 + Zr96g = Zr96 <<N_Zr96>> 2 + Fe54g = Fe54 <<N_Fe54>> + Fe56g = Fe56 <<N_Fe56>> + Fe57g = Fe57 <<N_Fe57>> + Fe58g = Fe58 <<N_Fe58>> + Cr50g = Cr50 <<N_Cr50>> + Cr52g = Cr52 <<N_Cr52>> + Cr53g = Cr53 <<N_Cr53>> + Cr54g = Cr54 <<N_Cr54>> + O16Gg = O16 <<N_O16G>> + Sn112g = Sn112 <<N_Sn112>> + Sn114g = Sn114 <<N_Sn114>> + Sn115g = Sn115 <<N_Sn115>> + Sn116g = Sn116 <<N_Sn116>> + Sn117g = Sn117 <<N_Sn117>> + Sn118g = Sn118 <<N_Sn118>> + Sn119g = Sn119 <<N_Sn119>> + Sn120g = Sn120 <<N_Sn120>> + Sn122g = Sn122 <<N_Sn122>> + Sn124g = Sn124 <<N_Sn124>> + +*---- +* Cladding fuel +*---- + MIX 2 <<TgainK>> NOEV + Zr90f = Zr90 <<N_Zr90>> 3 + Zr91f = Zr91 <<N_Zr91>> 3 + Zr92f = Zr92 <<N_Zr92>> 3 + Zr94f = Zr94 <<N_Zr94>> 3 + Zr96f = Zr96 <<N_Zr96>> 3 + Fe54f = Fe54 <<N_Fe54>> + Fe56f = Fe56 <<N_Fe56>> + Fe57f = Fe57 <<N_Fe57>> + Fe58f = Fe58 <<N_Fe58>> + Cr50f = Cr50 <<N_Cr50>> + Cr52f = Cr52 <<N_Cr52>> + Cr53f = Cr53 <<N_Cr53>> + Cr54f = Cr54 <<N_Cr54>> + O16Gf = O16 <<N_O16G>> + Sn112f = Sn112 <<N_Sn112>> + Sn114f = Sn114 <<N_Sn114>> + Sn115f = Sn115 <<N_Sn115>> + Sn116f = Sn116 <<N_Sn116>> + Sn117f = Sn117 <<N_Sn117>> + Sn118f = Sn118 <<N_Sn118>> + Sn119f = Sn119 <<N_Sn119>> + Sn120f = Sn120 <<N_Sn120>> + Sn122f = Sn122 <<N_Sn122>> + Sn124f = Sn124 <<N_Sn124>> + +*---- +* Gap +*---- + MIX 3 <<TgainK>> NOEV + Al27 = Al27 1.00E-08 + +*---- +* Moderator +*---- + MIX 4 <<TmodeK>> NOEV ! pastille d'U + H1_H2O = H1_H2O <<N_2H2O>> + O16_H2O = O16 <<N_H2O>> + B10 = B10 <<N_B10>> + B11 = B11 <<N_B11>> + MIX 5 + COMB 4 1.0 ! pastille d'U (2) + MIX 6 + COMB 4 1.0 ! gaines T2 + MIX 7 + COMB 4 1.0 ! gaines T1 + MIX 8 + COMB 4 1.0 ! gaine TI + +*---- +* Fuel -> 8 * 4 evolving mixes +*---- + MIX 9 <<TfuelK>> <<den_mox1>> + U234 = U234 <<X1_U4>> + U235 = U235 <<X1_U5>> 1 + U236 = U236 <<X1_U6>> 1 + U238 = U238 <<X1_U8>> 1 + Pu238 = Pu238 <<X1_Pu8>> + Pu239 = Pu239 <<X1_Pu9>> 1 + Pu240 = Pu240 <<X1_Pu0>> 1 + Pu241 = Pu241 <<X1_Pu1>> 1 + Pu242 = Pu242 <<X1_Pu2>> 1 + Am241 = Am241 <<X1_Am1>> + Np237 = Np237 <<X1_Np7>> + O16 = O16 <<X1_O16>> + MIX 10 + COMB 9 1.0 + MIX 11 + COMB 9 1.0 + MIX 12 + COMB 9 1.0 + + MIX 13 + COMB 9 1.0 + MIX 14 + COMB 9 1.0 + MIX 15 + COMB 9 1.0 + MIX 16 + COMB 9 1.0 + + MIX 17 + COMB 9 1.0 + MIX 18 + COMB 9 1.0 + MIX 19 + COMB 9 1.0 + MIX 20 + COMB 9 1.0 + + MIX 21 <<TfuelK>> <<den_mox2>> + U234 = U234 <<X2_U4>> + U235 = U235 <<X2_U5>> 1 + U236 = U236 <<X2_U6>> 1 + U238 = U238 <<X2_U8>> 1 + Pu238 = Pu238 <<X2_Pu8>> + Pu239 = Pu239 <<X2_Pu9>> 1 + Pu240 = Pu240 <<X2_Pu0>> 1 + Pu241 = Pu241 <<X2_Pu1>> 1 + Pu242 = Pu242 <<X2_Pu2>> 1 + Am241 = Am241 <<X2_Am1>> + Np237 = Np237 <<X2_Np7>> + O16 = O16 <<X2_O16>> + MIX 22 + COMB 21 1.0 + MIX 23 + COMB 21 1.0 + MIX 24 + COMB 21 1.0 + + MIX 25 + COMB 21 1.0 + MIX 26 + COMB 21 1.0 + MIX 27 + COMB 21 1.0 + MIX 28 + COMB 21 1.0 + + MIX 29 + COMB 21 1.0 + MIX 30 + COMB 21 1.0 + MIX 31 + COMB 21 1.0 + MIX 32 + COMB 21 1.0 + + MIX 33 <<TfuelK>> <<den_mox3>> + U234 = U234 <<X3_U4>> + U235 = U235 <<X3_U5>> 1 + U236 = U236 <<X3_U6>> 1 + U238 = U238 <<X3_U8>> 1 + Pu238 = Pu238 <<X3_Pu8>> + Pu239 = Pu239 <<X3_Pu9>> 1 + Pu240 = Pu240 <<X3_Pu0>> 1 + Pu241 = Pu241 <<X3_Pu1>> 1 + Pu242 = Pu242 <<X3_Pu2>> 1 + Am241 = Am241 <<X3_Am1>> + Np237 = Np237 <<X3_Np7>> + O16 = O16 <<X3_O16>> + MIX 34 + COMB 33 1.0 + MIX 35 + COMB 33 1.0 + MIX 36 + COMB 33 1.0 + + MIX 37 + COMB 33 1.0 + MIX 38 + COMB 33 1.0 + MIX 39 + COMB 33 1.0 + MIX 40 + COMB 33 1.0 +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Mix_UOX_20.c2m b/Dragon/data/rep900_mco_proc/Mix_UOX_20.c2m new file mode 100644 index 0000000..bd209a7 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Mix_UOX_20.c2m @@ -0,0 +1,311 @@ +*DECK Mix_Assb_DRA +*---- +* Name : Mix_Assb_DRA.c2m +* Type : DRAGON procedure +* Use : Microlib generation with Draglibs for +* UOX calculation with 20 fuel regions +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBRARY := Mix_Assb_DRA :: +* <<dens_mod>> <<pbore>> <<N_U5_UOX>> <<densU_UOX>> +* <<temp_comb>> <<temp_mode>> ; +* +* Input data : +* dens_mod : Theoretical moderator density +* pbore : Natural boron concentration (ppm) +* N_U5_UOX : U235 enrichment in UOX-Fuel +* densU_UOX : Theoretical density of Uranium in UOX-Fuel +* temp_comb : Theoretical Fuel temperature (K) +* temp_mode : Theoretical moderator temperature (K) +* +* Output data : +* LIBRARY : Microlib with table of probabilities + + +PARAMETER LIBRARY :: + EDIT 0 + ::: LINKED_LIST LIBRARY ; + ; + +*---- +* Input data recovery +*---- +STRING Library ; +:: >>Library<< ; + +REAL dens_mod pbore N_U5_UOX densU_UOX ; +:: >>dens_mod<< >>pbore<< >>N_U5_UOX<< >>densU_UOX<< ; +REAL temp_comb temp_mode ; +:: >>temp_comb<< >>temp_mode<< ; + +*---- +* Modules used in this procedure +*---- +MODULE LIB: UTL: DELETE: END: ; + +REAL TfuelK := temp_comb ; +REAL TmodeK := temp_mode ; + +REAL TgainK := TmodeK ; + +* Calculation of moderator data +* NOMBRE D AVOGADRO +REAL A := 6.022094E-1 ; ! Normalisation des concentrations par 10E-24 + +REAL M_H2O := 15.9994 2.0 1.00794 * + ; +REAL N_H2O := dens_mod 0.6022 * M_H2O / ; +REAL b10 b11 := 0.199 0.801 ; ! Bore naturel +REAL coeff := 0.199 A 1.0E-06 * * 10.8110164 / ; +REAL N_B10 := coeff 0.659 pbore * * ; ! densite nominale +REAL N_2H2O N_B11 := 2.0 N_H2O * N_B10 b11 * b10 / ; + +* Calculation of fuel data +REAL den_uo2 := densU_UOX ; +REAL M_U4 M_U5 M_U8 M_U6 := + 234.040952088 235.043929918 238.050788247 236.045568006 ; +REAL M_Pu8 M_Pu9 M_Pu0 M_Pu1 := + 238.049559894 239.052163381 240.053813545 241.056851456 ; +REAL M_Pu2 M_Am1 M_O16 := + 242.058742611 241.056829144 15.99491461956 ; + +REAL N_U5 := N_U5_UOX ; +REAL N_U4 := N_U5 54.0E-04 * 0.72 / ; +REAL N_U6 N_Pu8 N_Pu9 N_Pu0 N_Pu1 N_Pu2 N_Am1 := + 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 ; +REAL N_U8 := 1.0 N_U5 - N_U4 - N_Pu8 7.0 * - ; +REAL N_O16F := N_U4 M_U4 / N_U5 M_U5 / + N_U8 M_U8 / + N_Pu8 M_Pu8 / + ; +EVALUATE N_O16F := N_O16F N_U6 M_U6 / + ; +EVALUATE N_O16F := N_O16F N_Pu9 M_Pu9 / + N_Pu0 M_Pu0 / + ; +EVALUATE N_O16F := N_O16F N_Pu1 M_Pu1 / + N_Pu2 M_Pu2 / + ; +EVALUATE N_O16F := N_O16F N_Am1 M_Am1 / + 2.0 * M_O16 * ; + +REAL N_TOTF := N_U4 N_U5 + N_U6 + N_U8 + N_Pu8 + N_Pu9 + N_Pu0 + ; +EVALUATE N_TOTF := N_TOTF N_Pu1 + N_Pu2 + N_Am1 + N_O16F + ; +EVALUATE N_U4 := N_U4 N_TOTF / ; +EVALUATE N_U5 := N_U5 N_TOTF / ; +EVALUATE N_U6 := N_U6 N_TOTF / ; +EVALUATE N_U8 := N_U8 N_TOTF / ; +EVALUATE N_Pu8 := N_Pu8 N_TOTF / ; +EVALUATE N_Pu9 := N_Pu9 N_TOTF / ; +EVALUATE N_Pu0 := N_Pu0 N_TOTF / ; +EVALUATE N_Pu1 := N_Pu1 N_TOTF / ; +EVALUATE N_Pu2 := N_Pu2 N_TOTF / ; +EVALUATE N_Am1 := N_Am1 N_TOTF / ; +EVALUATE N_O16F := N_O16F N_TOTF / ; + +* -------------------------------- +* Cladding calculation +* -------------------------------- +REAL P_Zr P_Fe P_Cr P_O16 P_Sn := + 98.115E-02 0.210E-02 0.100E-02 0.125E-02 1.450E-02 ; + +REAL P_Zr90 P_Zr91 P_Zr92 P_Zr94 P_Zr96 := + 51.45E-02 11.22E-02 17.15E-02 17.38E-02 2.80E-02 ; +REAL P_Fe54 P_Fe56 P_Fe57 P_Fe58 := + 5.602E-02 91.865E-02 2.243E-02 0.290E-02 ; +*isot 5.8 % 91.72% 2.20% 0.28% +REAL P_Cr50 P_Cr52 P_Cr53 P_Cr54 := + 4.178E-02 83.701E-02 9.673E-02 2.448E-02 ; +*isot 4.35% 83.79% 9.50% 2.36% +REAL P_Sn112 P_Sn114 P_Sn115 P_Sn116 P_Sn117 + P_Sn118 P_Sn119 P_Sn120 P_Sn122 P_Sn124 := + 0.92E-02 0.62E-02 0.33E-02 14.21E-02 7.56E-02 + 24.06E-02 8.59E-02 32.92E-02 4.75E-02 6.04E-02 ; + +REAL M_Zr := P_Zr90 89.904711 * P_Zr91 90.905643 * P_Zr92 91.905039 * + P_Zr94 93.906320 * P_Zr96 95.90829 * + + + + ; +REAL M_Fe := P_Fe54 53.9396148 * P_Fe56 55.9349421 * + P_Fe57 56.9353987 * P_Fe58 57.9332805 * + + + ; +REAL M_Cr := P_Cr50 49.9460496 * P_Cr52 51.9405119 * + P_Cr53 52.9406538 * P_Cr54 53.9388849 * + + + ; +REAL M_Sn := P_Sn112 111.90483 * P_Sn114 113.90278 * P_Sn115 114.90335 * + P_Sn116 115.901748 * P_Sn117 116.902961 * + P_Sn118 117.901613 * P_Sn119 118.903316 * + P_Sn120 119.902207 * P_Sn122 121.903451 * + P_Sn124 123.905283 * + + + + + + + + + ; + +REAL rho_clad := 6.55 ; +REAL N_Zr_clad := P_Zr rho_clad A * M_Zr / * ; +REAL N_Zr90 N_Zr91 N_Zr92 N_Zr94 N_Zr96 := + N_Zr_clad P_Zr90 * N_Zr_clad P_Zr91 * N_Zr_clad P_Zr92 * + N_Zr_clad P_Zr94 * N_Zr_clad P_Zr96 * ; + +REAL N_Fe_clad := P_Fe rho_clad A * M_Fe / * ; +REAL N_Fe54 N_Fe56 N_Fe57 N_Fe58 := + N_Fe_clad P_Fe54 * N_Fe_clad P_Fe56 * N_Fe_clad P_Fe57 * + N_Fe_clad P_Fe58 * ; + +REAL N_Cr_clad := P_Cr rho_clad A * M_Cr / * ; +REAL N_Cr50 N_Cr52 N_Cr53 N_Cr54 := + N_Cr_clad P_Cr50 * N_Cr_clad P_Cr52 * N_Cr_clad P_Cr53 * + N_Cr_clad P_Cr54 * ; + +REAL N_O16G := P_O16 rho_clad A * 15.9949150 / * ; + +REAL N_Sn_clad := P_Sn rho_clad A * M_Sn / * ; +REAL N_Sn112 N_Sn114 N_Sn115 N_Sn116 + N_Sn117 N_Sn118 N_Sn119 N_Sn120 + N_Sn122 N_Sn124 := + N_Sn_clad P_Sn112 * N_Sn_clad P_Sn114 * N_Sn_clad P_Sn115 * + N_Sn_clad P_Sn116 * N_Sn_clad P_Sn117 * N_Sn_clad P_Sn118 * + N_Sn_clad P_Sn119 * N_Sn_clad P_Sn120 * N_Sn_clad P_Sn122 * + N_Sn_clad P_Sn124 * ; + +* ; +LIBRARY := LIB: :: + EDIT 0 + DEPL LIB: DRAGON FIL: <<Library>> + + NMIX 28 CTRA APOL + SUBG + MIXS LIB: DRAGON FIL: <<Library>> + +*---- +* Cladding TG et TI +*---- + MIX 1 <<TgainK>> NOEV + Zr90g = Zr90 <<N_Zr90>> 2 + Zr91g = Zr91 <<N_Zr91>> 2 + Zr92g = Zr92 <<N_Zr92>> 2 + Zr94g = Zr94 <<N_Zr94>> 2 + Zr96g = Zr96 <<N_Zr96>> 2 + Fe54g = Fe54 <<N_Fe54>> + Fe56g = Fe56 <<N_Fe56>> + Fe57g = Fe57 <<N_Fe57>> + Fe58g = Fe58 <<N_Fe58>> + Cr50g = Cr50 <<N_Cr50>> + Cr52g = Cr52 <<N_Cr52>> + Cr53g = Cr53 <<N_Cr53>> + Cr54g = Cr54 <<N_Cr54>> + O16Gg = O16 <<N_O16G>> + Sn112g = Sn112 <<N_Sn112>> + Sn114g = Sn114 <<N_Sn114>> + Sn115g = Sn115 <<N_Sn115>> + Sn116g = Sn116 <<N_Sn116>> + Sn117g = Sn117 <<N_Sn117>> + Sn118g = Sn118 <<N_Sn118>> + Sn119g = Sn119 <<N_Sn119>> + Sn120g = Sn120 <<N_Sn120>> + Sn122g = Sn122 <<N_Sn122>> + Sn124g = Sn124 <<N_Sn124>> + +*---- +* Cladding fuel +*---- + MIX 2 <<TgainK>> NOEV + Zr90f = Zr90 <<N_Zr90>> 3 + Zr91f = Zr91 <<N_Zr91>> 3 + Zr92f = Zr92 <<N_Zr92>> 3 + Zr94f = Zr94 <<N_Zr94>> 3 + Zr96f = Zr96 <<N_Zr96>> 3 + Fe54f = Fe54 <<N_Fe54>> + Fe56f = Fe56 <<N_Fe56>> + Fe57f = Fe57 <<N_Fe57>> + Fe58f = Fe58 <<N_Fe58>> + Cr50f = Cr50 <<N_Cr50>> + Cr52f = Cr52 <<N_Cr52>> + Cr53f = Cr53 <<N_Cr53>> + Cr54f = Cr54 <<N_Cr54>> + O16Gf = O16 <<N_O16G>> + Sn112f = Sn112 <<N_Sn112>> + Sn114f = Sn114 <<N_Sn114>> + Sn115f = Sn115 <<N_Sn115>> + Sn116f = Sn116 <<N_Sn116>> + Sn117f = Sn117 <<N_Sn117>> + Sn118f = Sn118 <<N_Sn118>> + Sn119f = Sn119 <<N_Sn119>> + Sn120f = Sn120 <<N_Sn120>> + Sn122f = Sn122 <<N_Sn122>> + Sn124f = Sn124 <<N_Sn124>> + +*---- +* Gap +*---- + MIX 3 <<TgainK>> NOEV + Al27 = Al27 1.00E-08 + +*---- +* Moderator +*---- + MIX 4 <<TmodeK>> NOEV ! pastille d'U + H1_H2O = H1_H2O <<N_2H2O>> + O16_H2O = O16 <<N_H2O>> + B10 = B10 <<N_B10>> + B11 = B11 <<N_B11>> + MIX 5 + COMB 4 1.0 ! pastille d'U (2) + MIX 6 + COMB 4 1.0 ! gaines T2 + MIX 7 + COMB 4 1.0 ! gaines T1 + MIX 8 + COMB 4 1.0 ! gaine TI + +*---- +* Fuel -> 5 * 4 evolving mixes +*---- +* C1 Cell + MIX 9 <<TfuelK>> <<den_uo2>> + U234 = U234 <<N_U4>> + U235 = U235 <<N_U5>> 1 + U236 = U236 <<N_U6>> 1 + U238 = U238 <<N_U8>> 1 + Pu238 = Pu238 <<N_Pu8>> + Pu239 = Pu239 <<N_Pu9>> 1 + Pu240 = Pu240 <<N_Pu0>> 1 + Pu241 = Pu241 <<N_Pu1>> 1 + Pu242 = Pu242 <<N_Pu2>> 1 + Am241 = Am241 <<N_Am1>> + O16F = O16 <<N_O16F>> + MIX 10 + COMB 9 1.0 + MIX 11 + COMB 9 1.0 + MIX 12 + COMB 9 1.0 + +* C2 Cell + MIX 13 + COMB 9 1.0 + MIX 14 + COMB 9 1.0 + MIX 15 + COMB 9 1.0 + MIX 16 + COMB 9 1.0 + +* C3 Cell + MIX 17 + COMB 9 1.0 + MIX 18 + COMB 9 1.0 + MIX 19 + COMB 9 1.0 + MIX 20 + COMB 9 1.0 + +* C4 Cell + MIX 21 + COMB 9 1.0 + MIX 22 + COMB 9 1.0 + MIX 23 + COMB 9 1.0 + MIX 24 + COMB 9 1.0 + +* C5 Cell + MIX 25 + COMB 9 1.0 + MIX 26 + COMB 9 1.0 + MIX 27 + COMB 9 1.0 + MIX 28 + COMB 9 1.0 +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/Mix_UOX_32.c2m b/Dragon/data/rep900_mco_proc/Mix_UOX_32.c2m new file mode 100644 index 0000000..95ef3e1 --- /dev/null +++ b/Dragon/data/rep900_mco_proc/Mix_UOX_32.c2m @@ -0,0 +1,341 @@ +*DECK Mix_Assb_DRA +*---- +* Name : Mix_Assb_DRA.c2m +* Type : DRAGON procedure +* Use : Microlib generation with Draglibs for +* UOX calculation with 32 fuel regions +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBRARY := Mix_Assb_DRA :: +* <<dens_mod>> <<pbore>> <<N_U5_UOX>> <<densU_UOX>> +* <<temp_comb>> <<temp_mode>> ; +* +* Input data : +* dens_mod : Theoretical moderator density +* pbore : Natural boron concentration (ppm) +* N_U5_UOX : U235 enrichment in UOX-Fuel +* densU_UOX : Theoretical density of Uranium in UOX-Fuel +* temp_comb : Theoretical Fuel temperature (K) +* temp_mode : Theoretical moderator temperature (K) +* +* Output data : +* LIBRARY : Microlib with table of probabilities + + +PARAMETER LIBRARY :: + EDIT 0 + ::: LINKED_LIST LIBRARY ; + ; + +*---- +* Input data recovery +*---- +STRING Library ; +:: >>Library<< ; + +REAL dens_mod pbore N_U5_UOX densU_UOX ; +:: >>dens_mod<< >>pbore<< >>N_U5_UOX<< >>densU_UOX<< ; +REAL temp_comb temp_mode ; +:: >>temp_comb<< >>temp_mode<< ; + +*---- +* Modules used in this procedure +*---- +MODULE LIB: UTL: DELETE: END: ; + +REAL TfuelK := temp_comb ; +REAL TmodeK := temp_mode ; + +REAL TgainK := TmodeK ; + +* Calculation of moderator data +* NOMBRE D AVOGADRO +REAL A := 6.022094E-1 ; ! Normalisation des concentrations par 10E-24 + +REAL M_H2O := 15.9994 2.0 1.00794 * + ; +REAL N_H2O := dens_mod A * M_H2O / ; +REAL b10 b11 := 0.199 0.801 ; ! Bore naturel +REAL coeff := 0.199 A 1.0E-06 * * 10.8110164 / ; +REAL N_B10 := coeff 0.659 pbore * * ; ! densite nominale +REAL N_2H2O N_B11 := 2.0 N_H2O * N_B10 b11 * b10 / ; + +* Calculation of fuel data +REAL den_uo2 := densU_UOX ; +REAL M_U4 M_U5 M_U8 M_U6 := + 234.040952088 235.043929918 238.050788247 236.045568006 ; +REAL M_Pu8 M_Pu9 M_Pu0 M_Pu1 := + 238.049559894 239.052163381 240.053813545 241.056851456 ; +REAL M_Pu2 M_Am1 M_O16 := + 242.058742611 241.056829144 15.99491461956 ; + +REAL N_U5 := N_U5_UOX ; +REAL N_U4 := N_U5 54.0E-04 * 0.72 / ; +REAL N_U6 N_Pu8 N_Pu9 N_Pu0 N_Pu1 N_Pu2 N_Am1 := + 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 ; +REAL N_U8 := 1.0 N_U5 - N_U4 - N_Pu8 7.0 * - ; +REAL N_O16F := N_U4 M_U4 / N_U5 M_U5 / + N_U8 M_U8 / + N_Pu8 M_Pu8 / + ; +EVALUATE N_O16F := N_O16F N_U6 M_U6 / + ; +EVALUATE N_O16F := N_O16F N_Pu9 M_Pu9 / + N_Pu0 M_Pu0 / + ; +EVALUATE N_O16F := N_O16F N_Pu1 M_Pu1 / + N_Pu2 M_Pu2 / + ; +EVALUATE N_O16F := N_O16F N_Am1 M_Am1 / + 2.0 * M_O16 * ; + +REAL N_TOTF := N_U4 N_U5 + N_U6 + N_U8 + N_Pu8 + N_Pu9 + N_Pu0 + ; +EVALUATE N_TOTF := N_TOTF N_Pu1 + N_Pu2 + N_Am1 + N_O16F + ; +EVALUATE N_U4 := N_U4 N_TOTF / ; +EVALUATE N_U5 := N_U5 N_TOTF / ; +EVALUATE N_U6 := N_U6 N_TOTF / ; +EVALUATE N_U8 := N_U8 N_TOTF / ; +EVALUATE N_Pu8 := N_Pu8 N_TOTF / ; +EVALUATE N_Pu9 := N_Pu9 N_TOTF / ; +EVALUATE N_Pu0 := N_Pu0 N_TOTF / ; +EVALUATE N_Pu1 := N_Pu1 N_TOTF / ; +EVALUATE N_Pu2 := N_Pu2 N_TOTF / ; +EVALUATE N_Am1 := N_Am1 N_TOTF / ; +EVALUATE N_O16F := N_O16F N_TOTF / ; + +* -------------------------------- +* Cladding calculation +* -------------------------------- +REAL P_Zr P_Fe P_Cr P_O16 P_Sn := + 98.115E-02 0.210E-02 0.100E-02 0.125E-02 1.450E-02 ; + +REAL P_Zr90 P_Zr91 P_Zr92 P_Zr94 P_Zr96 := + 51.45E-02 11.22E-02 17.15E-02 17.38E-02 2.80E-02 ; +REAL P_Fe54 P_Fe56 P_Fe57 P_Fe58 := + 5.602E-02 91.865E-02 2.243E-02 0.290E-02 ; +*isot 5.8 % 91.72% 2.20% 0.28% +REAL P_Cr50 P_Cr52 P_Cr53 P_Cr54 := + 4.178E-02 83.701E-02 9.673E-02 2.448E-02 ; +*isot 4.35% 83.79% 9.50% 2.36% +REAL P_Sn112 P_Sn114 P_Sn115 P_Sn116 P_Sn117 + P_Sn118 P_Sn119 P_Sn120 P_Sn122 P_Sn124 := + 0.92E-02 0.62E-02 0.33E-02 14.21E-02 7.56E-02 + 24.06E-02 8.59E-02 32.92E-02 4.75E-02 6.04E-02 ; + +REAL M_Zr := P_Zr90 89.904711 * P_Zr91 90.905643 * P_Zr92 91.905039 * + P_Zr94 93.906320 * P_Zr96 95.90829 * + + + + ; +REAL M_Fe := P_Fe54 53.9396148 * P_Fe56 55.9349421 * + P_Fe57 56.9353987 * P_Fe58 57.9332805 * + + + ; +REAL M_Cr := P_Cr50 49.9460496 * P_Cr52 51.9405119 * + P_Cr53 52.9406538 * P_Cr54 53.9388849 * + + + ; +REAL M_Sn := P_Sn112 111.90483 * P_Sn114 113.90278 * P_Sn115 114.90335 * + P_Sn116 115.901748 * P_Sn117 116.902961 * + P_Sn118 117.901613 * P_Sn119 118.903316 * + P_Sn120 119.902207 * P_Sn122 121.903451 * + P_Sn124 123.905283 * + + + + + + + + + ; + +REAL rho_clad := 6.55 ; +REAL N_Zr_clad := P_Zr rho_clad A * M_Zr / * ; +REAL N_Zr90 N_Zr91 N_Zr92 N_Zr94 N_Zr96 := + N_Zr_clad P_Zr90 * N_Zr_clad P_Zr91 * N_Zr_clad P_Zr92 * + N_Zr_clad P_Zr94 * N_Zr_clad P_Zr96 * ; + +REAL N_Fe_clad := P_Fe rho_clad A * M_Fe / * ; +REAL N_Fe54 N_Fe56 N_Fe57 N_Fe58 := + N_Fe_clad P_Fe54 * N_Fe_clad P_Fe56 * N_Fe_clad P_Fe57 * + N_Fe_clad P_Fe58 * ; + +REAL N_Cr_clad := P_Cr rho_clad A * M_Cr / * ; +REAL N_Cr50 N_Cr52 N_Cr53 N_Cr54 := + N_Cr_clad P_Cr50 * N_Cr_clad P_Cr52 * N_Cr_clad P_Cr53 * + N_Cr_clad P_Cr54 * ; + +REAL N_O16G := P_O16 rho_clad A * 15.9949150 / * ; + +REAL N_Sn_clad := P_Sn rho_clad A * M_Sn / * ; +REAL N_Sn112 N_Sn114 N_Sn115 N_Sn116 + N_Sn117 N_Sn118 N_Sn119 N_Sn120 + N_Sn122 N_Sn124 := + N_Sn_clad P_Sn112 * N_Sn_clad P_Sn114 * N_Sn_clad P_Sn115 * + N_Sn_clad P_Sn116 * N_Sn_clad P_Sn117 * N_Sn_clad P_Sn118 * + N_Sn_clad P_Sn119 * N_Sn_clad P_Sn120 * N_Sn_clad P_Sn122 * + N_Sn_clad P_Sn124 * ; + +* ; +LIBRARY := LIB: :: + EDIT 0 + DEPL LIB: DRAGON FIL: <<Library>> + + NMIX 40 CTRA APOL + SUBG + MIXS LIB: DRAGON FIL: <<Library>> + +*---- +* Cladding TG et TI +*---- + MIX 1 <<TgainK>> NOEV + Zr90g = Zr90 <<N_Zr90>> 2 + Zr91g = Zr91 <<N_Zr91>> 2 + Zr92g = Zr92 <<N_Zr92>> 2 + Zr94g = Zr94 <<N_Zr94>> 2 + Zr96g = Zr96 <<N_Zr96>> 2 + Fe54g = Fe54 <<N_Fe54>> + Fe56g = Fe56 <<N_Fe56>> + Fe57g = Fe57 <<N_Fe57>> + Fe58g = Fe58 <<N_Fe58>> + Cr50g = Cr50 <<N_Cr50>> + Cr52g = Cr52 <<N_Cr52>> + Cr53g = Cr53 <<N_Cr53>> + Cr54g = Cr54 <<N_Cr54>> + O16Gg = O16 <<N_O16G>> + Sn112g = Sn112 <<N_Sn112>> + Sn114g = Sn114 <<N_Sn114>> + Sn115g = Sn115 <<N_Sn115>> + Sn116g = Sn116 <<N_Sn116>> + Sn117g = Sn117 <<N_Sn117>> + Sn118g = Sn118 <<N_Sn118>> + Sn119g = Sn119 <<N_Sn119>> + Sn120g = Sn120 <<N_Sn120>> + Sn122g = Sn122 <<N_Sn122>> + Sn124g = Sn124 <<N_Sn124>> + +*---- +* Cladding fuel +*---- + MIX 2 <<TgainK>> NOEV + Zr90f = Zr90 <<N_Zr90>> 3 + Zr91f = Zr91 <<N_Zr91>> 3 + Zr92f = Zr92 <<N_Zr92>> 3 + Zr94f = Zr94 <<N_Zr94>> 3 + Zr96f = Zr96 <<N_Zr96>> 3 + Fe54f = Fe54 <<N_Fe54>> + Fe56f = Fe56 <<N_Fe56>> + Fe57f = Fe57 <<N_Fe57>> + Fe58f = Fe58 <<N_Fe58>> + Cr50f = Cr50 <<N_Cr50>> + Cr52f = Cr52 <<N_Cr52>> + Cr53f = Cr53 <<N_Cr53>> + Cr54f = Cr54 <<N_Cr54>> + O16Gf = O16 <<N_O16G>> + Sn112f = Sn112 <<N_Sn112>> + Sn114f = Sn114 <<N_Sn114>> + Sn115f = Sn115 <<N_Sn115>> + Sn116f = Sn116 <<N_Sn116>> + Sn117f = Sn117 <<N_Sn117>> + Sn118f = Sn118 <<N_Sn118>> + Sn119f = Sn119 <<N_Sn119>> + Sn120f = Sn120 <<N_Sn120>> + Sn122f = Sn122 <<N_Sn122>> + Sn124f = Sn124 <<N_Sn124>> + +*---- +* Gap +*---- + MIX 3 <<TgainK>> NOEV + Al27 = Al27 1.00E-08 + +*---- +* Moderator +*---- + MIX 4 <<TmodeK>> NOEV ! pastille d'U + H1_H2O = H1_H2O <<N_2H2O>> + O16_H2O = O16 <<N_H2O>> + B10 = B10 <<N_B10>> + B11 = B11 <<N_B11>> + MIX 5 + COMB 4 1.0 ! pastille d'U (2) + MIX 6 + COMB 4 1.0 ! gaines T2 + MIX 7 + COMB 4 1.0 ! gaines T1 + MIX 8 + COMB 4 1.0 ! gaine TI + +*---- +* Fuel -> 8 * 4 evolving mixes +*---- +* C1 Cell + MIX 9 <<TfuelK>> <<den_uo2>> + U234 = U234 <<N_U4>> + U235 = U235 <<N_U5>> 1 + U236 = U236 <<N_U6>> 1 + U238 = U238 <<N_U8>> 1 + Pu238 = Pu238 <<N_Pu8>> + Pu239 = Pu239 <<N_Pu9>> 1 + Pu240 = Pu240 <<N_Pu0>> 1 + Pu241 = Pu241 <<N_Pu1>> 1 + Pu242 = Pu242 <<N_Pu2>> 1 + Am241 = Am241 <<N_Am1>> + O16F = O16 <<N_O16F>> + MIX 10 + COMB 9 1.0 + MIX 11 + COMB 9 1.0 + MIX 12 + COMB 9 1.0 + +* C2 Cell + MIX 13 + COMB 9 1.0 + MIX 14 + COMB 9 1.0 + MIX 15 + COMB 9 1.0 + MIX 16 + COMB 9 1.0 + +* C3 Cell + MIX 17 + COMB 9 1.0 + MIX 18 + COMB 9 1.0 + MIX 19 + COMB 9 1.0 + MIX 20 + COMB 9 1.0 + +* C4 Cell + MIX 21 + COMB 9 1.0 + MIX 22 + COMB 9 1.0 + MIX 23 + COMB 9 1.0 + MIX 24 + COMB 9 1.0 + +* C5 Cell + MIX 25 + COMB 9 1.0 + MIX 26 + COMB 9 1.0 + MIX 27 + COMB 9 1.0 + MIX 28 + COMB 9 1.0 + +* C6 Cell + MIX 29 + COMB 9 1.0 + MIX 30 + COMB 9 1.0 + MIX 31 + COMB 9 1.0 + MIX 32 + COMB 9 1.0 + +* C7 Cell + MIX 33 + COMB 9 1.0 + MIX 34 + COMB 9 1.0 + MIX 35 + COMB 9 1.0 + MIX 36 + COMB 9 1.0 + +* C8 Cell + MIX 37 + COMB 9 1.0 + MIX 38 + COMB 9 1.0 + MIX 39 + COMB 9 1.0 + MIX 40 + COMB 9 1.0 +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/MultLIBEQ_20.c2m b/Dragon/data/rep900_mco_proc/MultLIBEQ_20.c2m new file mode 100644 index 0000000..6b1733e --- /dev/null +++ b/Dragon/data/rep900_mco_proc/MultLIBEQ_20.c2m @@ -0,0 +1,349 @@ +*DECK MultLIBEQ_20 +*---- +* Name : MultLIBEQ_20.c2m +* Type : DRAGON procedure +* Use : Increase the number of mixes in the microlib +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBEQ := MultLIBEQ_20 LIBEQ ; +* +* Input data : +* LIBEQ : Microlib with the number of mixs of the 1st level +* +* Output data : +* LIBEQ : Microlib with the number of mixs of the 2nd level + + +PARAMETER LIBEQ :: + EDIT 0 + ::: LINKED_LIST LIBEQ ; + ; +* ; +MODULE LIB: END: ; + +LIBEQ := LIB: LIBEQ :: + EDIT 0 + NMIX 164 + + MIXS LIB: DRAGON FIL: LIBEQ + +* Multiplication pour C1 +MIX 41 + COMB 9 1.0 +MIX 42 + COMB 10 1.0 +MIX 43 + COMB 11 1.0 +MIX 44 + COMB 12 1.0 + +MIX 45 + COMB 9 1.0 +MIX 46 + COMB 10 1.0 +MIX 47 + COMB 11 1.0 +MIX 48 + COMB 12 1.0 + +MIX 49 + COMB 9 1.0 +MIX 50 + COMB 10 1.0 +MIX 51 + COMB 11 1.0 +MIX 52 + COMB 12 1.0 + +MIX 53 + COMB 9 1.0 +MIX 54 + COMB 10 1.0 +MIX 55 + COMB 11 1.0 +MIX 56 + COMB 12 1.0 + +MIX 57 + COMB 9 1.0 +MIX 58 + COMB 10 1.0 +MIX 59 + COMB 11 1.0 +MIX 60 + COMB 12 1.0 + +MIX 61 + COMB 9 1.0 +MIX 62 + COMB 10 1.0 +MIX 63 + COMB 11 1.0 +MIX 64 + COMB 12 1.0 + +MIX 65 + COMB 9 1.0 +MIX 66 + COMB 10 1.0 +MIX 67 + COMB 11 1.0 +MIX 68 + COMB 12 1.0 + +MIX 69 + COMB 9 1.0 +MIX 70 + COMB 10 1.0 +MIX 71 + COMB 11 1.0 +MIX 72 + COMB 12 1.0 + +MIX 73 + COMB 9 1.0 +MIX 74 + COMB 10 1.0 +MIX 75 + COMB 11 1.0 +MIX 76 + COMB 12 1.0 + +MIX 77 + COMB 9 1.0 +MIX 78 + COMB 10 1.0 +MIX 79 + COMB 11 1.0 +MIX 80 + COMB 12 1.0 + +MIX 81 + COMB 9 1.0 +MIX 82 + COMB 10 1.0 +MIX 83 + COMB 11 1.0 +MIX 84 + COMB 12 1.0 + +MIX 85 + COMB 9 1.0 +MIX 86 + COMB 10 1.0 +MIX 87 + COMB 11 1.0 +MIX 88 + COMB 12 1.0 + +* Multiplication pour C2 +MIX 89 + COMB 13 1.0 +MIX 90 + COMB 14 1.0 +MIX 91 + COMB 15 1.0 +MIX 92 + COMB 16 1.0 + +MIX 93 + COMB 13 1.0 +MIX 94 + COMB 14 1.0 +MIX 95 + COMB 15 1.0 +MIX 96 + COMB 16 1.0 + +MIX 97 + COMB 13 1.0 +MIX 98 + COMB 14 1.0 +MIX 99 + COMB 15 1.0 +MIX 100 + COMB 16 1.0 + +MIX 101 + COMB 13 1.0 +MIX 102 + COMB 14 1.0 +MIX 103 + COMB 15 1.0 +MIX 104 + COMB 16 1.0 + +MIX 105 + COMB 13 1.0 +MIX 106 + COMB 14 1.0 +MIX 107 + COMB 15 1.0 +MIX 108 + COMB 16 1.0 + +MIX 109 + COMB 13 1.0 +MIX 110 + COMB 14 1.0 +MIX 111 + COMB 15 1.0 +MIX 112 + COMB 16 1.0 + +MIX 113 + COMB 13 1.0 +MIX 114 + COMB 14 1.0 +MIX 115 + COMB 15 1.0 +MIX 116 + COMB 16 1.0 + +MIX 117 + COMB 13 1.0 +MIX 118 + COMB 14 1.0 +MIX 119 + COMB 15 1.0 +MIX 120 + COMB 16 1.0 + +MIX 121 + COMB 13 1.0 +MIX 122 + COMB 14 1.0 +MIX 123 + COMB 15 1.0 +MIX 124 + COMB 16 1.0 + +MIX 125 + COMB 13 1.0 +MIX 126 + COMB 14 1.0 +MIX 127 + COMB 15 1.0 +MIX 128 + COMB 16 1.0 + +MIX 129 + COMB 13 1.0 +MIX 130 + COMB 14 1.0 +MIX 131 + COMB 15 1.0 +MIX 132 + COMB 16 1.0 + +* Multiplication pour C3 +MIX 133 + COMB 17 1.0 +MIX 134 + COMB 18 1.0 +MIX 135 + COMB 19 1.0 +MIX 136 + COMB 20 1.0 + +* Multiplication pour C4 -> devient type C1 -> C26 +MIX 137 + COMB 9 1.0 +MIX 138 + COMB 10 1.0 +MIX 139 + COMB 11 1.0 +MIX 140 + COMB 12 1.0 + +* Multiplication pour C6 -> devient du C4 +MIX 141 + COMB 21 1.0 +MIX 142 + COMB 22 1.0 +MIX 143 + COMB 23 1.0 +MIX 144 + COMB 24 1.0 + +MIX 145 + COMB 21 1.0 +MIX 146 + COMB 22 1.0 +MIX 147 + COMB 23 1.0 +MIX 148 + COMB 24 1.0 + +MIX 149 + COMB 21 1.0 +MIX 150 + COMB 22 1.0 +MIX 151 + COMB 23 1.0 +MIX 152 + COMB 24 1.0 + +MIX 153 + COMB 21 1.0 +MIX 154 + COMB 22 1.0 +MIX 155 + COMB 23 1.0 +MIX 156 + COMB 24 1.0 + +MIX 157 + COMB 21 1.0 +MIX 158 + COMB 22 1.0 +MIX 159 + COMB 23 1.0 +MIX 160 + COMB 24 1.0 + +MIX 161 + COMB 21 1.0 +MIX 162 + COMB 22 1.0 +MIX 163 + COMB 23 1.0 +MIX 164 + COMB 24 1.0 + + +* Ajout des milieux C5 (=C3) -> cellule C30 +MIX 29 + COMB 17 1.0 +MIX 30 + COMB 18 1.0 +MIX 31 + COMB 19 1.0 +MIX 32 + COMB 20 1.0 + +* Ajout des milieux C4 (=C1) -> cellule C23 +MIX 33 + COMB 9 1.0 +MIX 34 + COMB 10 1.0 +MIX 35 + COMB 11 1.0 +MIX 36 + COMB 12 1.0 + +* Ajout des milieux C7 (=C4) -> cellule C38 +MIX 37 + COMB 21 1.0 +MIX 38 + COMB 22 1.0 +MIX 39 + COMB 23 1.0 +MIX 40 + COMB 24 1.0 +; + +END: ; +QUIT . diff --git a/Dragon/data/rep900_mco_proc/MultLIBEQ_32.c2m b/Dragon/data/rep900_mco_proc/MultLIBEQ_32.c2m new file mode 100644 index 0000000..b63dd9c --- /dev/null +++ b/Dragon/data/rep900_mco_proc/MultLIBEQ_32.c2m @@ -0,0 +1,318 @@ +*DECK MultLIBEQ_32 +*---- +* Name : MultLIBEQ_32.c2m +* Type : DRAGON procedure +* Use : Increase the number of mixes in the microlib +* Author : R. Vallerent +* +* Procedure called as: +* +* LIBEQ := MultLIBEQ_32 LIBEQ ; +* +* Input data : +* LIBEQ : Microlib with the number of mixs of the 1st level +* +* Output data : +* LIBEQ : Microlib with the number of mixs of the 2nd level + + +PARAMETER LIBEQ :: + EDIT 0 + ::: LINKED_LIST LIBEQ ; + ; +* ; +MODULE LIB: END: ; + +LIBEQ := LIB: LIBEQ :: + EDIT 0 + NMIX 164 + + MIXS LIB: DRAGON FIL: LIBEQ + +* Multiplication pour C1 +MIX 41 + COMB 9 1.0 +MIX 42 + COMB 10 1.0 +MIX 43 + COMB 11 1.0 +MIX 44 + COMB 12 1.0 + +MIX 45 + COMB 9 1.0 +MIX 46 + COMB 10 1.0 +MIX 47 + COMB 11 1.0 +MIX 48 + COMB 12 1.0 + +MIX 49 + COMB 9 1.0 +MIX 50 + COMB 10 1.0 +MIX 51 + COMB 11 1.0 +MIX 52 + COMB 12 1.0 + +MIX 53 + COMB 9 1.0 +MIX 54 + COMB 10 1.0 +MIX 55 + COMB 11 1.0 +MIX 56 + COMB 12 1.0 + +MIX 57 + COMB 9 1.0 +MIX 58 + COMB 10 1.0 +MIX 59 + COMB 11 1.0 +MIX 60 + COMB 12 1.0 + +MIX 61 + COMB 9 1.0 +MIX 62 + COMB 10 1.0 +MIX 63 + COMB 11 1.0 +MIX 64 + COMB 12 1.0 + +MIX 65 + COMB 9 1.0 +MIX 66 + COMB 10 1.0 +MIX 67 + COMB 11 1.0 +MIX 68 + COMB 12 1.0 + +MIX 69 + COMB 9 1.0 +MIX 70 + COMB 10 1.0 +MIX 71 + COMB 11 1.0 +MIX 72 + COMB 12 1.0 + +MIX 73 + COMB 9 1.0 +MIX 74 + COMB 10 1.0 +MIX 75 + COMB 11 1.0 +MIX 76 + COMB 12 1.0 + +MIX 77 + COMB 9 1.0 +MIX 78 + COMB 10 1.0 +MIX 79 + COMB 11 1.0 +MIX 80 + COMB 12 1.0 + +MIX 81 + COMB 9 1.0 +MIX 82 + COMB 10 1.0 +MIX 83 + COMB 11 1.0 +MIX 84 + COMB 12 1.0 + +MIX 85 + COMB 9 1.0 +MIX 86 + COMB 10 1.0 +MIX 87 + COMB 11 1.0 +MIX 88 + COMB 12 1.0 + +* Multiplication pour C2 +MIX 89 + COMB 13 1.0 +MIX 90 + COMB 14 1.0 +MIX 91 + COMB 15 1.0 +MIX 92 + COMB 16 1.0 + +MIX 93 + COMB 13 1.0 +MIX 94 + COMB 14 1.0 +MIX 95 + COMB 15 1.0 +MIX 96 + COMB 16 1.0 + +MIX 97 + COMB 13 1.0 +MIX 98 + COMB 14 1.0 +MIX 99 + COMB 15 1.0 +MIX 100 + COMB 16 1.0 + +MIX 101 + COMB 13 1.0 +MIX 102 + COMB 14 1.0 +MIX 103 + COMB 15 1.0 +MIX 104 + COMB 16 1.0 + +MIX 105 + COMB 13 1.0 +MIX 106 + COMB 14 1.0 +MIX 107 + COMB 15 1.0 +MIX 108 + COMB 16 1.0 + +MIX 109 + COMB 13 1.0 +MIX 110 + COMB 14 1.0 +MIX 111 + COMB 15 1.0 +MIX 112 + COMB 16 1.0 + +MIX 113 + COMB 13 1.0 +MIX 114 + COMB 14 1.0 +MIX 115 + COMB 15 1.0 +MIX 116 + COMB 16 1.0 + +MIX 117 + COMB 13 1.0 +MIX 118 + COMB 14 1.0 +MIX 119 + COMB 15 1.0 +MIX 120 + COMB 16 1.0 + +MIX 121 + COMB 13 1.0 +MIX 122 + COMB 14 1.0 +MIX 123 + COMB 15 1.0 +MIX 124 + COMB 16 1.0 + +MIX 125 + COMB 13 1.0 +MIX 126 + COMB 14 1.0 +MIX 127 + COMB 15 1.0 +MIX 128 + COMB 16 1.0 + +MIX 129 + COMB 13 1.0 +MIX 130 + COMB 14 1.0 +MIX 131 + COMB 15 1.0 +MIX 132 + COMB 16 1.0 + +* Multiplication pour C3 +MIX 133 + COMB 17 1.0 +MIX 134 + COMB 18 1.0 +MIX 135 + COMB 19 1.0 +MIX 136 + COMB 20 1.0 + +* Multiplication pour C4 +MIX 137 + COMB 21 1.0 +MIX 138 + COMB 22 1.0 +MIX 139 + COMB 23 1.0 +MIX 140 + COMB 24 1.0 + +* Multiplication pour C6 +MIX 141 + COMB 29 1.0 +MIX 142 + COMB 30 1.0 +MIX 143 + COMB 31 1.0 +MIX 144 + COMB 32 1.0 + +MIX 145 + COMB 29 1.0 +MIX 146 + COMB 30 1.0 +MIX 147 + COMB 31 1.0 +MIX 148 + COMB 32 1.0 + +MIX 149 + COMB 29 1.0 +MIX 150 + COMB 30 1.0 +MIX 151 + COMB 31 1.0 +MIX 152 + COMB 32 1.0 + +MIX 153 + COMB 29 1.0 +MIX 154 + COMB 30 1.0 +MIX 155 + COMB 31 1.0 +MIX 156 + COMB 32 1.0 + +MIX 157 + COMB 29 1.0 +MIX 158 + COMB 30 1.0 +MIX 159 + COMB 31 1.0 +MIX 160 + COMB 32 1.0 + +MIX 161 + COMB 29 1.0 +MIX 162 + COMB 30 1.0 +MIX 163 + COMB 31 1.0 +MIX 164 + COMB 32 1.0 +; + +END: ; +QUIT . |