From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001 From: stainer_t Date: Mon, 8 Sep 2025 13:48:49 +0200 Subject: Initial commit from Polytechnique Montreal --- Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m | 516 ++++++++++++++++++++++ 1 file changed, 516 insertions(+) create mode 100755 Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m (limited to 'Donjon/data/Test4x4coupl_nompi_proc/BibCreaG2.c2m') 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 :: +* <> +* <> <> <> +* <> <> <> +* <> <> <> <> +* <> +* ; +* 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: <> K + PUR: <> WGT% + LIB: WIMSAECL FIL: <> + ISO: 3 <> <> <> + CALC WGT% D2O >>CwgtH1<< >>CwgtD2<< >>CwgtO16<< ; + EVALUATE LDensCalo := DensCalo ; + ELSE + INFO: :: + TMP: <> K + PUR: <> WGT% + CALC DENS WATER >>LDensCalo<< + LIB: WIMSAECL FIL: <> + ISO: 3 <> <> <> + CALC WGT% D2O >>CwgtH1<< >>CwgtD2<< >>CwgtO16<< ; + ENDIF ; +IF DensMod 0.0 > THEN + INFO: :: + TMP: <> K + PUR: <> WGT% + LIB: WIMSAECL FIL: <> + ISO: 3 <> <> <> + CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ; + EVALUATE LDensMod := DensMod ; + ELSE + INFO: :: + TMP: <> K + PUR: <> WGT% + CALC DENS WATER >>LDensMod<< + LIB: WIMSAECL FIL: <> + ISO: 3 <> <> <> + CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ; + ENDIF ; +INFO: :: + ENR: <> WGT% + LIB: WIMSAECL FIL: <> + ISO: 3 <> <> <> + CALC WGT% UO2 >>FwgtU5<< >>FwgtU8<< >>FwgtO16<< ; +*---- +* donnees d'evolution de <> +* section efficace microscopiques de <> +*---- +*---- +* Nom de la bibliotheque et des isotopes +*---- +IF RefLib "Winfrith" = THEN +*---- +* Winfrith +*---- + Biblio := LIB: :: + EDIT <> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <> + MIXS LIB: WIMSAECL FIL: <> + MIX 1 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + MIX 2 <> 6.57 + ZR = <> <> + NB = <> 2.5 + BN = <> <> + MIX 3 <> 0.0016987 + O16 = <> 72.89 + C12 = <> 27.11 + MIX 4 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + BN = <> <> + MIX 5 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + BN = <> <> + MIX 6 <> <> + O16 = <> <> + U235 = <> <> 1 + U236 = <> 0.0 1 + U238 = <> <> 1 + PU239 = <> 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 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + BN = <> <> ; +*---- +ELSEIF RefLib "IGE227W" = THEN +*---- +* Winfrith +* IGE227 +*---- + Biblio := LIB: :: + EDIT <> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <> + MIXS LIB: WIMSAECL FIL: <> + MIX 1 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + MIX 2 <> 6.57 + ZR = <> <> + NB = <> 2.5 + BN = <> <> + MIX 3 <> 0.0014 + HE4 = <> 100.0 + MIX 4 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + BN = <> <> + MIX 5 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + BN = <> <> + MIX 6 <> <> + O16 = <> <> + U235 = <> <> 1 + U236 = <> 0.0 1 + U238 = <> <> 1 + PU239 = <> 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 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + BN = <> <> ; +ELSEIF RefLib "ENDFB6" = THEN +*---- +* ENDF/B-VI +*---- + Biblio := LIB: :: + EDIT <> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <> + MIXS LIB: WIMSAECL FIL: <> + MIX 1 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + MIX 2 <> 6.57 + Zr90 = ZR90PT <> + Zr91 = ZR91PT <> + Zr92 = ZR92PT <> + Zr94 = ZR94PT <> + Zr96 = ZR96PT <> + NB = <> 2.5 + B10 = <> <> + B11 = <> <> + MIX 3 <> 0.0016987 + O16 = <> 72.89 + C12 = <> 27.11 + MIX 4 <> 6.44 + Zr90 = ZR90CT <> + Zr91 = ZR91CT <> + Zr92 = ZR92CT <> + Zr94 = ZR94CT <> + Zr96 = ZR96CT <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + B10 = <> <> + B11 = <> <> + MIX 5 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + B10 = <> <> + B11 = <> <> + MIX 6 <> <> + O16 = <> <> + U235 = <> <> 1 + U236 = <> 0.0 1 + U238 = <> <> 1 + PU239 = <> 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 <> 6.44 + Zr90 = ZR90CL <> + Zr91 = ZR91CL <> + Zr92 = ZR92CL <> + Zr94 = ZR94CL <> + Zr96 = ZR96CL <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + B10 = <> <> + B11 = <> <> ; +ELSEIF RefLib "IGE227E5" = THEN +*---- +* ENDF/B-V +* IGE227 +*---- + Biblio := LIB: :: + EDIT <> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <> + MIXS LIB: WIMSAECL FIL: <> + MIX 1 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + MIX 2 <> 6.57 + ZR = <> <> + NB = <> 2.5 + BN = <> <> + MIX 3 <> 0.0014 + HE4 = <> 100.0 + MIX 4 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + BN = <> <> + MIX 5 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + BN = <> <> + MIX 6 <> <> + O16 = <> <> + U235 = <> <> 1 + U236 = <> 0.0 1 + U238 = <> <> 1 + PU239 = <> 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 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + BN = <> <> ; +ELSE +*---- +* ENDF/B-V +*---- + Biblio := LIB: :: + EDIT <> + NMIX 10 CTRA WIMS + DEPL LIB: WIMSAECL FIL: <> + MIXS LIB: WIMSAECL FIL: <> + MIX 1 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + MIX 2 <> 6.57 + ZR = <> <> + NB = <> 2.5 + B10 = <> <> + B11 = <> <> + MIX 3 <> 0.0016987 + O16 = <> 72.89 + C12 = <> 27.11 + MIX 4 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + B10 = <> <> + B11 = <> <> + MIX 5 <> <> + H1 = <> <> + D2 = <> <> + O16 = <> <> + B10 = <> <> + B11 = <> <> + MIX 6 <> <> + O16 = <> <> + U235 = <> <> 1 + U236 = <> 0.0 1 + U238 = <> <> 1 + PU239 = <> 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 <> 6.44 + ZR = <> <> + FE = <> 0.16 + CR = <> 0.11 + NI = <> 0.06 + B10 = <> <> + B11 = <> <> ; + ENDIF ; +*---- +* terminer execution +*---- +QUIT "LIST" . -- cgit v1.2.3