diff options
Diffstat (limited to 'Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m')
| -rwxr-xr-x | Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m | 516 |
1 files changed, 516 insertions, 0 deletions
diff --git a/Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m b/Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m new file mode 100755 index 0000000..8c131d9 --- /dev/null +++ b/Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m @@ -0,0 +1,516 @@ +*DECK BibCreaG2.c2m +*---------- +* Nom : BibCreaG2.c2m +* Type : procedure DRAGON +* Usage : Creer une bibliotheque pour cellule G-2 +* a partir des parametres de cellules pour +* du combustible neuf +* Auteur : G. Marleau +* Date : 2002/12/03 +* +* Description de la procedure : +* Biblio := BibCreaG2 :: +* <<RefLib>> +* <<TempComb>> <<EnriComb>> <<DensComb>> +* <<TempCalo>> <<PurtCalo>> <<DensCalo>> +* <<TempMod>> <<PurtMod>> <<DensMod>> <<BoreMod>> +* <<Impression>> +* ; +* Structures de donnees: +* Biblio : Bibliotheque de section efficace (L_LIBRARY). +* Parametres de controle: +* RefLib : Type de bibliotheque +* "Winfrith" -> Bibliotheque Winfrith WIMS-AECL +* "ENDFB5" -> Bibliotheque ENDF/B-V WIMS-AECL +* (defaut) +* "ENDFB6" -> Bibliotheque ENDF/B-VI WIMS-AECL +* "IGE227W" -> Bibliotheque Winfrith WIMS-AECL +* isotopes de IGE227 +* "IGE227E5" -> Bibliotheque ENDF/B-V WIMS-AECL +* isotopes de IGE227 +* TempComb : Temperature combustible (C) +* EnriComb : Enrichissement combustible (%) +* DensComb : Densite combustible (g/cc) +* TempCalo : Temperature caloporteur (C) +* PurtCalo : % D2O caloporteur +* DensCalo : Densite caloporteur (g/cc) +* A calculer si 0.0 +* TempMod : Temperature moderateur (C) +* PurtMod : % D2O moderateur +* DensMod : Densite moderateur (g/cc) +* A calculer si 0.0 +* BoreMod : Bore dans le moderateur (ppm) +* Impression : Niveau d'impression +* Modules DRAGON requis: +* INFO: : Gestion des densites et concentration isotopique. +* LIB: : Gestion des bibliotheque microscopiques. +* +*---------- +* Definition des parametres de la procedure +* +PARAMETER Biblio :: + EDIT 1 + ::: LINKED_LIST Biblio ; ; +STRING RefLib ; +REAL TempComb EnriComb DensComb + TempCalo PurtCalo DensCalo + TempMod PurtMod DensMod BoreMod ; +INTEGER Impression ; +*---- +* Reccuperer information transferee a la procedure +*---- +:: >>RefLib<< + >>TempComb<< >>EnriComb<< >>DensComb<< + >>TempCalo<< >>PurtCalo<< >>DensCalo<< + >>TempMod<< >>PurtMod<< >>DensMod<< >>BoreMod<< + >>Impression<< ; +*---- +* Imprimer donnees d'entree si requis +*---- +IF Impression 10 > THEN + ECHO "Bibliotheque =" RefLib ; + ECHO "Temperature combustible =" TempComb "C" ; + ECHO "Enrichissement combustible =" EnriComb "%" ; + ECHO "Densite combustible =" DensComb "g/cc" ; + ECHO "Temperature caloporteur =" TempCalo "C" ; + ECHO "Purete caloporteur =" PurtCalo "%" ; + IF DensCalo 0.0 > THEN + ECHO "Densite caloporteur =" DensCalo "g/cc" ; + ENDIF ; + ECHO "Temperature moderateur =" TempMod "C" ; + ECHO "Purete moderateur =" PurtMod "%" ; + IF DensMod 0.0 > THEN + ECHO "Densite moderateur =" DensMod "g/cc" ; + ENDIF ; + ECHO "Bore moderateur =" BoreMod "ppm" ; +ENDIF ; +*---- +* Modules +*-- +MODULE INFO: LIB: ; +*---- +* Local variables +*---- +STRING LibName ; +STRING IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoCr IsoNi IsoNb + IsoC12 IsoHe ; +*---- +* Nom de la bibliotheque et des isotopes +*---- +IF RefLib "Winfrith" = THEN +*---- +* Winfrith +*---- + EVALUATE LibName := "WIMSLIB" ; + EVALUATE IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoCr IsoNi IsoNb + IsoC12 IsoHe := + "U235" "U238" "U236" "PU239" "O16" + "H1" "D2" "B10" " " "B-NAT" + "ZR91" "FE56" "CR52" "NI58" "NB93" + "C12" "HE4" ; +ELSEIF RefLib "IGE227W" = THEN +*---- +* Winfrith +* IGE227 +*---- + EVALUATE LibName := "WIMSLIB" ; + EVALUATE IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoCr IsoNi IsoNb + IsoC12 IsoHe := + "U235" "U238" "U236" "PU239" "O16" + "H1" "D2" "B10" " " "B-NAT" + "ZR91" "FE56" "CR52" "NI58" "NB93" + "C12" "HE4" ; +ELSEIF RefLib "ENDFB6" = THEN +*---- +* ENDF/B-VI +*---- + EVALUATE LibName := "E6MLIB" ; + EVALUATE IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoCr IsoNi IsoNb + IsoC12 IsoHe := + "U235" "U238" "U236" "PU239" "O16" + "H1H2O" "DD2O" "B10" "B11" " " + " " "FE56" "CR52" "NI58" "NB93" + "C" " " ; +ELSEIF RefLib "IGE227E5" = THEN +*---- +* ENDF/B-V +* IGE227 +*---- + EVALUATE LibName := "E5WLIB" ; + EVALUATE IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoCr IsoNi IsoNb + IsoC12 IsoHe := + "U235" "U238" "U236" "PU239" "O16" + "H-H2O" "D-D2O" "B10" "B11" "B10" + "ZR-NAT" "FE" "CR" "NI" "NB93" + "C" "HE4" ; +ELSE +*---- +* ENDF/B-V +*---- + EVALUATE LibName := "E5WLIB" ; + EVALUATE IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoCr IsoNi IsoNb + IsoC12 IsoHe := + "U235" "U238" "U236" "PU239" "O16" + "H-H2O" "D-D2O" "B10" "B11" " " + "ZR-NAT" "FE" "CR" "NI" "NB93" + "C" "HE4" ; + ENDIF ; +*---- +* Parametres pour densite calo et moderateur +* ainsi que concentration de bore +*---- +REAL LDensCalo LDensMod ; +REAL BNZrNb BNZr BNMod := + 0.00021 0.00031 BoreMod 1.0E-4 * ; +REAL DZrNPt DZrNCt := + 97.5 99.71 ; +REAL B10ZrNb B10Zr + B11ZrNb B11Zr + B10Mod B11Mod := + 0.199 BNZrNb * 0.199 BNZr * + 0.801 BNZrNb * 0.801 BNZr * + 0.199 BNMod * 0.801 BNMod * ; +REAL DZr90Pt DZr91Pt + DZr92Pt DZr94Pt + DZr96Pt := + DZrNPt 51.45 * DZrNPt 11.22 * + DZrNPt 17.15 * DZrNPt 17.38 * + DZrNPt 2.80 * ; +REAL DZr90Ct DZr91Ct + DZr92Ct DZr94Ct + DZr96Ct := + DZrNCt 51.45 * DZrNCt 11.22 * + DZrNCt 17.15 * DZrNCt 17.38 * + DZrNCt 2.80 * ; +REAL + FwgtU8 FwgtU5 FwgtO16 + CwgtH1 CwgtD2 CwgtO16 Cdens + MwgtH1 MwgtD2 MwgtO16 Mdens ; +*---- +* Calcul densite et WGT% pour caloporteur et moderateur +* et WGT% combustible +*---- +IF DensCalo 0.0 > THEN + INFO: :: + TMP: <<TempCalo>> K + PUR: <<PurtCalo>> WGT% + LIB: WIMSAECL FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>CwgtH1<< >>CwgtD2<< >>CwgtO16<< ; + EVALUATE LDensCalo := DensCalo ; + ELSE + INFO: :: + TMP: <<TempCalo>> K + PUR: <<PurtCalo>> WGT% + CALC DENS WATER >>LDensCalo<< + LIB: WIMSAECL FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>CwgtH1<< >>CwgtD2<< >>CwgtO16<< ; + ENDIF ; +IF DensMod 0.0 > THEN + INFO: :: + TMP: <<TempMod>> K + PUR: <<PurtMod>> WGT% + LIB: WIMSAECL FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ; + EVALUATE LDensMod := DensMod ; + ELSE + INFO: :: + TMP: <<TempMod>> K + PUR: <<PurtMod>> WGT% + CALC DENS WATER >>LDensMod<< + LIB: WIMSAECL FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ; + ENDIF ; +INFO: :: + ENR: <<EnriComb>> WGT% + LIB: WIMSAECL FIL: <<LibName>> + ISO: 3 <<IsoU5>> <<IsoU8>> <<IsoO16>> + CALC WGT% UO2 >>FwgtU5<< >>FwgtU8<< >>FwgtO16<< ; +*---- +* donnees d'evolution de <<LibName>> +* section efficace microscopiques de <<LibName>> +*---- +*---- +* Nom de la bibliotheque et des isotopes +*---- +IF RefLib "Winfrith" = THEN +*---- +* Winfrith +*---- + Biblio := LIB: :: + EDIT <<Impression>> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <<LibName>> + MIXS LIB: WIMSAECL FIL: <<LibName>> + MIX 1 <<TempCalo>> <<LDensCalo>> + H1 = <<IsoH1>> <<CwgtH1>> + D2 = <<IsoD2>> <<CwgtD2>> + O16 = <<IsoO16>> <<CwgtO16>> + MIX 2 <<TempCalo>> 6.57 + ZR = <<IsoZr>> <<DZrNPt>> + NB = <<IsoNb>> 2.5 + BN = <<IsoBn>> <<BNZrNb>> + MIX 3 <<TempMod>> 0.0016987 + O16 = <<IsoO16>> 72.89 + C12 = <<IsoC12>> 27.11 + MIX 4 <<TempMod>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + BN = <<IsoBn>> <<BNZr>> + MIX 5 <<TempMod>> <<LDensMod>> + H1 = <<IsoH1>> <<MwgtH1>> + D2 = <<IsoD2>> <<MwgtD2>> + O16 = <<IsoO16>> <<MwgtO16>> + BN = <<IsoBn>> <<BNMod>> + MIX 6 <<TempComb>> <<DensComb>> + O16 = <<IsoO16>> <<FwgtO16>> + U235 = <<IsoU5>> <<FwgtU5>> 1 + U236 = <<IsoU6>> 0.0 1 + U238 = <<IsoU8>> <<FwgtU8>> 1 + PU239 = <<IsoPu9>> 0.0 1 + PU240 = PU240 0.0 1 + PU241 = PU241 0.0 1 + MIX 7 COMB 6 100.0 + MIX 8 COMB 6 100.0 + MIX 9 COMB 6 100.0 + MIX 10 <<TempCalo>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + BN = <<IsoBn>> <<BNZr>> ; +*---- +ELSEIF RefLib "IGE227W" = THEN +*---- +* Winfrith +* IGE227 +*---- + Biblio := LIB: :: + EDIT <<Impression>> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <<LibName>> + MIXS LIB: WIMSAECL FIL: <<LibName>> + MIX 1 <<TempCalo>> <<LDensCalo>> + H1 = <<IsoH1>> <<CwgtH1>> + D2 = <<IsoD2>> <<CwgtD2>> + O16 = <<IsoO16>> <<CwgtO16>> + MIX 2 <<TempCalo>> 6.57 + ZR = <<IsoZr>> <<DZrNPt>> + NB = <<IsoNb>> 2.5 + BN = <<IsoBn>> <<BNZrNb>> + MIX 3 <<TempMod>> 0.0014 + HE4 = <<IsoHe>> 100.0 + MIX 4 <<TempMod>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + BN = <<IsoBn>> <<BNZr>> + MIX 5 <<TempMod>> <<LDensMod>> + H1 = <<IsoH1>> <<MwgtH1>> + D2 = <<IsoD2>> <<MwgtD2>> + O16 = <<IsoO16>> <<MwgtO16>> + BN = <<IsoBn>> <<BNMod>> + MIX 6 <<TempComb>> <<DensComb>> + O16 = <<IsoO16>> <<FwgtO16>> + U235 = <<IsoU5>> <<FwgtU5>> 1 + U236 = <<IsoU6>> 0.0 1 + U238 = <<IsoU8>> <<FwgtU8>> 1 + PU239 = <<IsoPu9>> 0.0 1 + PU240 = PU240 0.0 1 + PU241 = PU241 0.0 1 + MIX 7 COMB 6 100.0 + MIX 8 COMB 6 100.0 + MIX 9 COMB 6 100.0 + MIX 10 <<TempCalo>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + BN = <<IsoBn>> <<BNZr>> ; +ELSEIF RefLib "ENDFB6" = THEN +*---- +* ENDF/B-VI +*---- + Biblio := LIB: :: + EDIT <<Impression>> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <<LibName>> + MIXS LIB: WIMSAECL FIL: <<LibName>> + MIX 1 <<TempCalo>> <<LDensCalo>> + H1 = <<IsoH1>> <<CwgtH1>> + D2 = <<IsoD2>> <<CwgtD2>> + O16 = <<IsoO16>> <<CwgtO16>> + MIX 2 <<TempCalo>> 6.57 + Zr90 = ZR90PT <<DZr90Pt>> + Zr91 = ZR91PT <<DZr91Pt>> + Zr92 = ZR92PT <<DZr92Pt>> + Zr94 = ZR94PT <<DZr94Pt>> + Zr96 = ZR96PT <<DZr96Pt>> + NB = <<IsoNb>> 2.5 + B10 = <<IsoB10>> <<B10ZrNb>> + B11 = <<IsoB11>> <<B11ZrNb>> + MIX 3 <<TempMod>> 0.0016987 + O16 = <<IsoO16>> 72.89 + C12 = <<IsoC12>> 27.11 + MIX 4 <<TempMod>> 6.44 + Zr90 = ZR90CT <<DZr90Ct>> + Zr91 = ZR91CT <<DZr91Ct>> + Zr92 = ZR92CT <<DZr92Ct>> + Zr94 = ZR94CT <<DZr94Ct>> + Zr96 = ZR96CT <<DZr96Ct>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + B10 = <<IsoB10>> <<B10Zr>> + B11 = <<IsoB11>> <<B11Zr>> + MIX 5 <<TempMod>> <<LDensMod>> + H1 = <<IsoH1>> <<MwgtH1>> + D2 = <<IsoD2>> <<MwgtD2>> + O16 = <<IsoO16>> <<MwgtO16>> + B10 = <<IsoB10>> <<B10Mod>> + B11 = <<IsoB11>> <<B11Mod>> + MIX 6 <<TempComb>> <<DensComb>> + O16 = <<IsoO16>> <<FwgtO16>> + U235 = <<IsoU5>> <<FwgtU5>> 1 + U236 = <<IsoU6>> 0.0 1 + U238 = <<IsoU8>> <<FwgtU8>> 1 + PU239 = <<IsoPu9>> 0.0 1 + PU240 = PU240 0.0 1 + PU241 = PU241 0.0 1 + MIX 7 COMB 6 100.0 + MIX 8 COMB 6 100.0 + MIX 9 COMB 6 100.0 + MIX 10 <<TempCalo>> 6.44 + Zr90 = ZR90CL <<DZr90Ct>> + Zr91 = ZR91CL <<DZr91Ct>> + Zr92 = ZR92CL <<DZr92Ct>> + Zr94 = ZR94CL <<DZr94Ct>> + Zr96 = ZR96CL <<DZr96Ct>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + B10 = <<IsoB10>> <<B10Zr>> + B11 = <<IsoB11>> <<B11Zr>> ; +ELSEIF RefLib "IGE227E5" = THEN +*---- +* ENDF/B-V +* IGE227 +*---- + Biblio := LIB: :: + EDIT <<Impression>> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <<LibName>> + MIXS LIB: WIMSAECL FIL: <<LibName>> + MIX 1 <<TempCalo>> <<LDensCalo>> + H1 = <<IsoH1>> <<CwgtH1>> + D2 = <<IsoD2>> <<CwgtD2>> + O16 = <<IsoO16>> <<CwgtO16>> + MIX 2 <<TempCalo>> 6.57 + ZR = <<IsoZr>> <<DZrNPt>> + NB = <<IsoNb>> 2.5 + BN = <<IsoBn>> <<BNZrNb>> + MIX 3 <<TempMod>> 0.0014 + HE4 = <<IsoHe>> 100.0 + MIX 4 <<TempMod>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + BN = <<IsoBn>> <<BNZr>> + MIX 5 <<TempMod>> <<LDensMod>> + H1 = <<IsoH1>> <<MwgtH1>> + D2 = <<IsoD2>> <<MwgtD2>> + O16 = <<IsoO16>> <<MwgtO16>> + BN = <<IsoBn>> <<BNMod>> + MIX 6 <<TempComb>> <<DensComb>> + O16 = <<IsoO16>> <<FwgtO16>> + U235 = <<IsoU5>> <<FwgtU5>> 1 + U236 = <<IsoU6>> 0.0 1 + U238 = <<IsoU8>> <<FwgtU8>> 1 + PU239 = <<IsoPu9>> 0.0 1 + PU240 = PU240 0.0 1 + PU241 = PU241 0.0 1 + MIX 7 COMB 6 100.0 + MIX 8 COMB 6 100.0 + MIX 9 COMB 6 100.0 + MIX 10 <<TempCalo>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + BN = <<IsoBn>> <<BNZr>> ; +ELSE +*---- +* ENDF/B-V +*---- + Biblio := LIB: :: + EDIT <<Impression>> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <<LibName>> + MIXS LIB: WIMSAECL FIL: <<LibName>> + MIX 1 <<TempCalo>> <<LDensCalo>> + H1 = <<IsoH1>> <<CwgtH1>> + D2 = <<IsoD2>> <<CwgtD2>> + O16 = <<IsoO16>> <<CwgtO16>> + MIX 2 <<TempCalo>> 6.57 + ZR = <<IsoZr>> <<DZrNPt>> + NB = <<IsoNb>> 2.5 + B10 = <<IsoB10>> <<B10ZrNb>> + B11 = <<IsoB11>> <<B11ZrNb>> + MIX 3 <<TempMod>> 0.0016987 + O16 = <<IsoO16>> 72.89 + C12 = <<IsoC12>> 27.11 + MIX 4 <<TempMod>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + B10 = <<IsoB10>> <<B10Zr>> + B11 = <<IsoB11>> <<B11Zr>> + MIX 5 <<TempMod>> <<LDensMod>> + H1 = <<IsoH1>> <<MwgtH1>> + D2 = <<IsoD2>> <<MwgtD2>> + O16 = <<IsoO16>> <<MwgtO16>> + B10 = <<IsoB10>> <<B10Mod>> + B11 = <<IsoB11>> <<B11Mod>> + MIX 6 <<TempComb>> <<DensComb>> + O16 = <<IsoO16>> <<FwgtO16>> + U235 = <<IsoU5>> <<FwgtU5>> 1 + U236 = <<IsoU6>> 0.0 1 + U238 = <<IsoU8>> <<FwgtU8>> 1 + PU239 = <<IsoPu9>> 0.0 1 + PU240 = PU240 0.0 1 + PU241 = PU241 0.0 1 + MIX 7 COMB 6 100.0 + MIX 8 COMB 6 100.0 + MIX 9 COMB 6 100.0 + MIX 10 <<TempCalo>> 6.44 + ZR = <<IsoZr>> <<DZrNCt>> + FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 + NI = <<IsoNi>> 0.06 + B10 = <<IsoB10>> <<B10Zr>> + B11 = <<IsoB11>> <<B11Zr>> ; + ENDIF ; +*---- +* terminer execution +*---- +QUIT "LIST" . |