Initial commit from Polytechnique Montreal

1 files changed, 1548 insertions, 0 deletions

diff --git a/Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m b/Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m
new file mode 100644
index 0000000..607e3a6
--- /dev/null
+++ b/Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m
@@ -0,0 +1,1548 @@
+*DECK MetalMix.c2m
+*======================================================
+* Name:   MetalMix.c2m
+* Type:   DRAGON procedure
+* Use:    Create isotopic mixtures of metals.
+* Author: A. Bruneau 10.2013
+*         [CS-SI | IRSN/PSN-EXP/SNC/LNR]
+* Mod:    A. Bruneau 12.2013
+*         [CS-SI | IRSN/PSN-EXP/SNC/LNR]
+*         Added AIC with AgNat, CdNat & In for AIC rod description.
+* Mod:    A. Bruneau 02.2014
+*         [CS-SI | IRSN/PSN-EXP/SNC/LNR]
+*         Added GdNat with Gd for Gd2O3
+* Mod:    G. Tixier 09.2016
+*         [CS-SI | IRSN/PSN-EXP/SNC/LNR]
+*         Ajouts S32, S33 et S34 pour materiau M5
+*------------------------------------------------------
+* Procedure called as:
+*------------------------------------------------------
+* Input Required:
+*  > Molar Masses
+* Output Provided:
+*  < Molar Masses
+*  < Metals Concentration
+*======================================================
+PARAMETER   FeNat CrNat SnNat NiNat BNat
+            ZrNat TiNat SiNat KNat
+            AgNat CdNat SNat InNat Zirc4 Incon SS304
+            M5 AIC GdNat Gd2O3 BP MolarMasses ::
+        ::: LINKED_LIST
+            FeNat CrNat SnNat NiNat BNat
+            ZrNat TiNat SiNat KNat
+            AgNat CdNat SNat InNat Zirc4 Incon SS304
+            M5 AIC GdNat Gd2O3 BP MolarMasses ; ;
+MODULE      GREP: UTL: END: ABORT: ;
+STRING      LibType ; ! ENDFVII1 JEFF31 JEFF311 JEFF312...
+STRING      Gestion ;
+REAL        avogadro := 6.022094E+23 ;
+REAL        avo_conc := avogadro 1.E-24 * ;
+REAL        ReNorm ; ! Renormalization Factor
+REAL        SumFr  ; ! Fraction Sum (needed when &IMPU)
+:: >>LibType<< >>Gestion<< ;
+
+*------------------------------------------------------
+REAL am241 ;
+REAL al27 ;
+REAL b10 ;
+REAL b10b ;
+REAL b11 ;
+REAL b11b ;
+REAL co59 ;
+REAL cr50 ;
+REAL cr52 ;
+REAL cr53 ;
+REAL cr54 ;
+REAL fe54 ;
+REAL fe56 ;
+REAL fe57 ;
+REAL fe58 ;
+REAL h2o h1 ;
+REAL mn55 ;
+REAL mo95 ;
+REAL nb93 ;
+REAL ni58 ;
+REAL ni60 ;
+REAL ni61 ;
+REAL ni62 ;
+REAL ni64 ;
+REAL o16 ;
+REAL pu238 ;
+REAL pu239 ;
+REAL pu240 ;
+REAL pu241 ;
+REAL pu242 ;
+REAL siNat si28 si29 si30 ;
+REAL sn112 ;
+REAL sn114 ;
+REAL sn115 ;
+REAL sn116 ;
+REAL sn117 ;
+REAL sn118 ;
+REAL sn119 ;
+REAL sn120 ;
+REAL sn122 ;
+REAL sn124 ;
+REAL tiNat ti46 ti47 ti48 ti49 ti50 ;
+REAL u234 ;
+REAL u235 ;
+REAL u236 ;
+REAL u238 ;
+REAL zrNat zr90 zr91 zr92 zr93 zr94 zr95 zr96 ;
+REAL cnat ;
+*---- added for AIC
+REAL ag107 ;
+REAL ag109 ;
+REAL cd106 ;
+REAL cd108 ;
+REAL cd110 ;
+REAL cd111 ;
+REAL cd112 ;
+REAL cd113 ;
+REAL cd114 ;
+REAL cd116 ;
+REAL in115 in113 ;
+*---- added for Gd2O3
+REAL gd154 ;
+REAL gd155 ;
+REAL gd156 ;
+REAL gd157 ;
+REAL gd158 ;
+REAL gd160 ;
+*---- added for Pyrex
+REAL kNat k39 k40 k41 ;
+REAL na23 ;
+*---- added for M5
+REAL s32  ;
+REAL s33  ;
+REAL s34  ;
+
+*------------------------------------------------------
+*------------------------------------------------------
+* Recuperation des fractions massiques moyennes
+*------------------------------------------------------
+
+* Pour les isotopes non conflictuels
+
+GREP: MolarMasses :: GETVAL 'B10  ' 1 >>b10<<   ;
+GREP: MolarMasses :: GETVAL 'B10B ' 1 >>b10b<<  ;
+GREP: MolarMasses :: GETVAL 'B11  ' 1 >>b11<<   ;
+GREP: MolarMasses :: GETVAL 'B11B ' 1 >>b11b<<  ;
+GREP: MolarMasses :: GETVAL 'CNat ' 1 >>cnat<<  ;
+GREP: MolarMasses :: GETVAL 'O16  ' 1 >>o16<<   ;
+GREP: MolarMasses :: GETVAL 'Al27 ' 1 >>al27<<  ;
+GREP: MolarMasses :: GETVAL 'Cr50 ' 1 >>cr50<<  ;
+GREP: MolarMasses :: GETVAL 'Cr52 ' 1 >>cr52<<  ;
+GREP: MolarMasses :: GETVAL 'Cr53 ' 1 >>cr53<<  ;
+GREP: MolarMasses :: GETVAL 'Cr54 ' 1 >>cr54<<  ;
+GREP: MolarMasses :: GETVAL 'Mn55 ' 1 >>mn55<<  ;
+GREP: MolarMasses :: GETVAL 'Fe54 ' 1 >>fe54<<  ;
+GREP: MolarMasses :: GETVAL 'Fe56 ' 1 >>fe56<<  ;
+GREP: MolarMasses :: GETVAL 'Fe57 ' 1 >>fe57<<  ;
+GREP: MolarMasses :: GETVAL 'Fe58 ' 1 >>fe58<<  ;
+GREP: MolarMasses :: GETVAL 'Co59 ' 1 >>co59<<  ;
+GREP: MolarMasses :: GETVAL 'Ni58 ' 1 >>ni58<<  ;
+GREP: MolarMasses :: GETVAL 'Ni60 ' 1 >>ni60<<  ;
+GREP: MolarMasses :: GETVAL 'Ni61 ' 1 >>ni61<<  ;
+GREP: MolarMasses :: GETVAL 'Ni62 ' 1 >>ni62<<  ;
+GREP: MolarMasses :: GETVAL 'Ni64 ' 1 >>ni64<<  ;
+GREP: MolarMasses :: GETVAL 'Nb93 ' 1 >>nb93<<  ;
+GREP: MolarMasses :: GETVAL 'Mo95 ' 1 >>mo95<<  ;
+GREP: MolarMasses :: GETVAL 'Ag107' 1 >>ag107<< ;
+GREP: MolarMasses :: GETVAL 'Ag109' 1 >>ag109<< ;
+GREP: MolarMasses :: GETVAL 'Cd110' 1 >>cd110<< ;
+GREP: MolarMasses :: GETVAL 'Cd111' 1 >>cd111<< ;
+GREP: MolarMasses :: GETVAL 'Cd112' 1 >>cd112<< ;
+GREP: MolarMasses :: GETVAL 'Cd113' 1 >>cd113<< ;
+GREP: MolarMasses :: GETVAL 'Cd114' 1 >>cd114<< ;
+GREP: MolarMasses :: GETVAL 'Cd116' 1 >>cd116<< ;
+GREP: MolarMasses :: GETVAL 'Sn115' 1 >>sn115<< ;
+GREP: MolarMasses :: GETVAL 'Sn116' 1 >>sn116<< ;
+GREP: MolarMasses :: GETVAL 'Sn117' 1 >>sn117<< ;
+GREP: MolarMasses :: GETVAL 'Sn118' 1 >>sn118<< ;
+GREP: MolarMasses :: GETVAL 'Sn119' 1 >>sn119<< ;
+GREP: MolarMasses :: GETVAL 'Sn120' 1 >>sn120<< ;
+GREP: MolarMasses :: GETVAL 'Sn122' 1 >>sn122<< ;
+GREP: MolarMasses :: GETVAL 'Sn124' 1 >>sn124<< ;
+GREP: MolarMasses :: GETVAL 'Gd154' 1 >>gd154<< ;
+GREP: MolarMasses :: GETVAL 'Gd155' 1 >>gd155<< ;
+GREP: MolarMasses :: GETVAL 'Gd156' 1 >>gd156<< ;
+GREP: MolarMasses :: GETVAL 'Gd157' 1 >>gd157<< ;
+GREP: MolarMasses :: GETVAL 'Gd158' 1 >>gd158<< ;
+GREP: MolarMasses :: GETVAL 'Gd160' 1 >>gd160<< ;
+GREP: MolarMasses :: GETVAL 'U234 ' 1 >>u234<<  ;
+GREP: MolarMasses :: GETVAL 'U235 ' 1 >>u235<<  ;
+GREP: MolarMasses :: GETVAL 'U236 ' 1 >>u236<<  ;
+GREP: MolarMasses :: GETVAL 'U238 ' 1 >>u238<<  ;
+GREP: MolarMasses :: GETVAL 'Pu238' 1 >>pu238<< ;
+GREP: MolarMasses :: GETVAL 'Pu239' 1 >>pu239<< ;
+GREP: MolarMasses :: GETVAL 'Pu240' 1 >>pu240<< ;
+GREP: MolarMasses :: GETVAL 'Pu241' 1 >>pu241<< ;
+GREP: MolarMasses :: GETVAL 'Pu242' 1 >>pu242<< ;
+GREP: MolarMasses :: GETVAL 'Am241' 1 >>am241<< ;
+GREP: MolarMasses :: GETVAL 'Na23 ' 1 >>na23<< ;
+
+* Pour les isotopes conflictuels (absent de certaines bibli.)
+
+IF LibType 'ENDFVI8'    = THEN
+  GREP: MolarMasses :: GETVAL 'H1   ' 1 >>h1<<    ;
+  GREP: MolarMasses :: GETVAL 'SiNat' 1 >>siNat<< ;
+  GREP: MolarMasses :: GETVAL 'TiNat' 1 >>tiNat<< ;
+  GREP: MolarMasses :: GETVAL 'ZrNat' 1 >>zrNat<< ;
+  GREP: MolarMasses :: GETVAL 'InNat' 1 >>in115<< ;
+      ! In115 non defini dans ENDFVI8
+  GREP: MolarMasses :: GETVAL 'Knat ' 1 >>kNat<< ;
+  GREP: MolarMasses :: GETVAL 'S32  ' 1 >>s32<< ;
+
+ELSEIF LibType 'JEF22'       = THEN
+  GREP: MolarMasses :: GETVAL 'H1   ' 1 >>h1<<    ;
+  GREP: MolarMasses :: GETVAL 'SiNat' 1 >>siNat<< ;
+  GREP: MolarMasses :: GETVAL 'TiNat' 1 >>tiNat<< ;
+  GREP: MolarMasses :: GETVAL 'ZrNat' 1 >>zrNat<< ;
+  GREP: MolarMasses :: GETVAL 'In115' 1 >>in115<< ;
+  GREP: MolarMasses :: GETVAL 'Knat ' 1 >>kNat<< ;
+  GREP: MolarMasses :: GETVAL 'S32  ' 1 >>s32<< ;
+  GREP: MolarMasses :: GETVAL 'S33  ' 1 >>s33<< ;
+  GREP: MolarMasses :: GETVAL 'S34  ' 1 >>s34<< ;
+
+ELSEIF LibType 'JEFF31'   =   LibType 'JEFF311'  = +
+       LibType 'JEFF312'  = +
+       LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+  GREP: MolarMasses :: GETVAL 'H1   ' 1 >>h1<<    ;
+  GREP: MolarMasses :: GETVAL 'Si28 ' 1 >>si28<<  ;
+  GREP: MolarMasses :: GETVAL 'Si29 ' 1 >>si29<<  ;
+  GREP: MolarMasses :: GETVAL 'Si30 ' 1 >>si30<<  ;
+  GREP: MolarMasses :: GETVAL 'Ti46 ' 1 >>ti46<<  ;
+  GREP: MolarMasses :: GETVAL 'Ti47 ' 1 >>ti47<<  ;
+  GREP: MolarMasses :: GETVAL 'Ti48 ' 1 >>ti48<<  ;
+  GREP: MolarMasses :: GETVAL 'Ti49 ' 1 >>ti49<<  ;
+  GREP: MolarMasses :: GETVAL 'Ti50 ' 1 >>ti50<<  ;
+  GREP: MolarMasses :: GETVAL 'Zr90 ' 1 >>zr90<<  ;
+  GREP: MolarMasses :: GETVAL 'Zr91 ' 1 >>zr91<<  ;
+  GREP: MolarMasses :: GETVAL 'Zr92 ' 1 >>zr92<<  ;
+  GREP: MolarMasses :: GETVAL 'Zr93 ' 1 >>zr93<<  ;
+  GREP: MolarMasses :: GETVAL 'Zr94 ' 1 >>zr94<<  ;
+  GREP: MolarMasses :: GETVAL 'Zr95 ' 1 >>zr95<<  ;
+  GREP: MolarMasses :: GETVAL 'Zr96 ' 1 >>zr96<<  ;
+  GREP: MolarMasses :: GETVAL 'Cd106' 1 >>cd106<< ;
+  GREP: MolarMasses :: GETVAL 'Cd108' 1 >>cd108<< ;
+  GREP: MolarMasses :: GETVAL 'In115' 1 >>in115<< ;
+  GREP: MolarMasses :: GETVAL 'Sn112' 1 >>sn112<< ;
+  GREP: MolarMasses :: GETVAL 'Sn114' 1 >>sn114<< ;
+  GREP: MolarMasses :: GETVAL 'K39  ' 1 >>k39<< ;
+  GREP: MolarMasses :: GETVAL 'K40  ' 1 >>k40<< ;
+  GREP: MolarMasses :: GETVAL 'K41  ' 1 >>k41<< ;
+  GREP: MolarMasses :: GETVAL 'S32  ' 1 >>s32<< ;
+  GREP: MolarMasses :: GETVAL 'S33  ' 1 >>s33<< ;
+  GREP: MolarMasses :: GETVAL 'S34  ' 1 >>s34<< ;
+
+ELSE ! CEA99v9p5, CEA514T2, CLA99CEA93
+  GREP: MolarMasses :: GETVAL 'H1'    1 >>h1<<    ;
+  GREP: MolarMasses :: GETVAL 'SiNat' 1 >>siNat<< ;
+  GREP: MolarMasses :: GETVAL 'TiNat' 1 >>tiNat<< ;
+  GREP: MolarMasses :: GETVAL 'ZrNat' 1 >>zrNat<< ;
+  GREP: MolarMasses :: GETVAL 'In115' 1 >>in115<< ;
+  GREP: MolarMasses :: GETVAL 'Knat ' 1 >>kNat<< ;
+  GREP: MolarMasses :: GETVAL 'S32  ' 1 >>s32<< ;
+  GREP: MolarMasses :: GETVAL 'S33  ' 1 >>s33<< ;
+  GREP: MolarMasses :: GETVAL 'S34  ' 1 >>s34<< ;
+
+ENDIF ;
+
+*------------------------------------------------------
+*  ZrNat : Natural Zirconium
+*  Isotopic Composition : Zr90 .5145
+*                         Zr91 .1123
+*                         Zr92 .1716
+*                         Zr93 .
+*                         Zr94 .1738
+*                         Zr95 .
+*                         Zr96 .0281
+*------------------------------------------------------
+REAL Zr90f :=  .0 ;
+REAL Zr91f :=  .0 ;
+REAL Zr92f :=  .0 ;
+REAL Zr93f :=  .0 ;
+REAL Zr94f :=  .0 ;
+REAL Zr95f :=  .0 ;
+REAL Zr96f :=  .0 ;
+REAL MZrNat ; ! Average Weight Ratio
+
+IF LibType 'JEFF31'   =   LibType 'JEFF311' = +
+   LibType 'JEFF312'  = +
+   LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+
+  EVALUATE Zr90f :=  .5145 ;
+  EVALUATE Zr91f :=  .1122 ;
+  EVALUATE Zr92f :=  .1715 ;
+  EVALUATE Zr93f := 0.     ;
+  EVALUATE Zr94f :=  .1738 ;
+  EVALUATE Zr95f := 0.     ;
+  EVALUATE Zr96f :=  .0280 ;
+
+  *---- ZrNat Average Weight Ratio ----
+  EVALUATE MZrNat :=        zr90  Zr90f  *   zr91  Zr91f  * +
+                            zr92  Zr92f  * + zr93  Zr93f  * +
+                            zr94  Zr94f  * + zr95  Zr95f  * +
+                            zr96  Zr96f  * + ;
+
+  ECHO "Masse molaire du zirconium " MZrNat ;
+
+ELSE ! JEF22
+
+  EVALUATE MZrNat := zrNat ;
+
+ENDIF ;
+
+*---- ZrNat Table Creation ----
+ZrNat := UTL:       :: CREA 'AvgMolM  ' 1 = <<MZrNat>> ;
+ZrNat := UTL: ZrNat :: CREA 'Zr90frIso' 1 = <<Zr90f>>  ;
+ZrNat := UTL: ZrNat :: CREA 'Zr91frIso' 1 = <<Zr91f>>  ;
+ZrNat := UTL: ZrNat :: CREA 'Zr92frIso' 1 = <<Zr92f>>  ;
+ZrNat := UTL: ZrNat :: CREA 'Zr93frIso' 1 = <<Zr93f>>  ;
+ZrNat := UTL: ZrNat :: CREA 'Zr94frIso' 1 = <<Zr94f>>  ;
+ZrNat := UTL: ZrNat :: CREA 'Zr95frIso' 1 = <<Zr95f>>  ;
+ZrNat := UTL: ZrNat :: CREA 'Zr96frIso' 1 = <<Zr96f>>  ;
+
+*------------------------------------------------------
+*  FeNat : Natural Iron
+*  Isotopic Composition : Fe54 .0580
+*                         Fe56 .9172
+*                         Fe57 .0220
+*                         Fe58 .0028
+*------------------------------------------------------
+REAL Fe54f := .0580 ;
+REAL Fe56f := .9172 ;
+REAL Fe57f := .0220 ;
+REAL Fe58f := .0028 ;
+
+*---- FeNat Average Weight Ratio ----
+REAL     MFeNat ; ! Average Weight Ratio
+EVALUATE MFeNat :=        fe54  Fe54f  *   fe56  Fe56f  * + ;
+EVALUATE MFeNat := MFeNat fe57  Fe57f  * + fe58  Fe58f  * + ;
+
+*---- FeNat Table Creation ----
+FeNat := UTL:       :: CREA 'AvgMolM  ' 1 = <<MFeNat>> ;
+FeNat := UTL: FeNat :: CREA 'Fe54frIso' 1 = <<Fe54f>>  ;
+FeNat := UTL: FeNat :: CREA 'Fe56frIso' 1 = <<Fe56f>>  ;
+FeNat := UTL: FeNat :: CREA 'Fe57frIso' 1 = <<Fe57f>>  ;
+FeNat := UTL: FeNat :: CREA 'Fe58frIso' 1 = <<Fe58f>>  ;
+
+*------------------------------------------------------
+*  TiNat : Natural Titanium
+*  Isotopic Composition : Ti46 .0825
+*                         Ti47 .0744
+*                         Ti48 .7372
+*                         Ti49 .0541
+*                         Ti50 .0518
+*------------------------------------------------------
+REAL Ti46f := .0 ;
+REAL Ti47f := .0 ;
+REAL Ti48f := .0 ;
+REAL Ti49f := .0 ;
+REAL Ti50f := .0 ;
+REAL MTiNat ; ! Average Weight Ratio
+
+IF LibType 'JEFF31'   =   LibType 'JEFF311' = +
+   LibType 'JEFF312'  = +
+   LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+
+  EVALUATE Ti46f := .0825 ;
+  EVALUATE Ti47f := .0744 ;
+  EVALUATE Ti48f := .7372 ;
+  EVALUATE Ti49f := .0541 ;
+  EVALUATE Ti50f := .0518 ;
+
+  *---- TiNat Average Weight Ratio ----
+  EVALUATE MTiNat := ti46  Ti46f  *   ti47  Ti47f  * +
+                     ti48  Ti48f  * +
+                     ti49  Ti49f  * + ti50  Ti50f  * + ;
+
+  ECHO "Masse molaire du titane " MTiNat ;
+
+ELSE ! JEF22 ENDFVI8
+
+  EVALUATE MTiNat := tiNat ;
+
+ENDIF ;
+
+*---- TiNat Table Creation ----
+TiNat := UTL:       :: CREA 'AvgMolM  ' 1 = <<MTiNat>> ;
+TiNat := UTL: TiNat :: CREA 'Ti46frIso' 1 = <<Ti46f>>  ;
+TiNat := UTL: TiNat :: CREA 'Ti47frIso' 1 = <<Ti47f>>  ;
+TiNat := UTL: TiNat :: CREA 'Ti48frIso' 1 = <<Ti48f>>  ;
+TiNat := UTL: TiNat :: CREA 'Ti49frIso' 1 = <<Ti49f>>  ;
+TiNat := UTL: TiNat :: CREA 'Ti50frIso' 1 = <<Ti50f>>  ;
+
+*------------------------------------------------------
+*  SiNat : Natural Silicium
+*  Isotopic Composition : Si28 .922298
+*                         Si29 .046833
+*                         Si30 .030873
+*------------------------------------------------------
+REAL Si28f := .0 ;
+REAL Si29f := .0 ;
+REAL Si30f := .0 ;
+REAL MSiNat ; ! Average Weight Ratio
+
+
+IF LibType 'JEFF31'   =   LibType 'JEFF311' = +
+   LibType 'JEFF312'  = +
+   LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+
+  EVALUATE Si28f := .9222 ;
+  EVALUATE Si29f := .0468 ;
+  EVALUATE Si30f := .0309 ;
+
+  *---- SiNat Average Weight Ratio ----
+  EVALUATE MSiNat := si28  Si28f  *   si29  Si29f  * +
+                     si30  Si30f  * + ;
+
+  ECHO "Masse molaire du silicium " MSiNat ;
+
+ELSE ! JEF22 ENDFVI8
+
+  EVALUATE MSiNat := siNat ;
+
+ENDIF ;
+
+*---- SiNat Table Creation ----
+SiNat := UTL:       :: CREA 'AvgMolM  ' 1 = <<MSiNat>> ;
+SiNat := UTL: SiNat :: CREA 'Si28frIso' 1 = <<Si28f>>  ;
+SiNat := UTL: SiNat :: CREA 'Si29frIso' 1 = <<Si29f>>  ;
+SiNat := UTL: SiNat :: CREA 'Si30frIso' 1 = <<Si30f>>  ;
+
+*------------------------------------------------------
+*  CrNat : Natural Chromium
+*  Isotopic Composition : Cr50 .0435
+*                         Cr52 .8379
+*                         Cr53 .0950
+*                         Cr54 .0236
+*------------------------------------------------------
+REAL Cr50f := .0435 ;
+REAL Cr52f := .8379 ;
+REAL Cr53f := .0950 ;
+REAL Cr54f := .0236 ;
+
+*---- CrNat Average Weight Ratio ----
+REAL     MCrNat ; ! Average Weight Ratio
+EVALUATE MCrNat :=        cr50  Cr50f  *   cr52  Cr52f  * + ;
+EVALUATE MCrNat := MCrNat cr53  Cr53f  * + cr54  Cr54f  * + ;
+
+*---- CrNat Table Creation ----
+CrNat := UTL:       :: CREA 'AvgMolM  ' 1 = <<MCrNat>> ;
+CrNat := UTL: CrNat :: CREA 'Cr50frIso' 1 = <<Cr50f>>  ;
+CrNat := UTL: CrNat :: CREA 'Cr52frIso' 1 = <<Cr52f>>  ;
+CrNat := UTL: CrNat :: CREA 'Cr53frIso' 1 = <<Cr53f>>  ;
+CrNat := UTL: CrNat :: CREA 'Cr54frIso' 1 = <<Cr54f>>  ;
+
+*------------------------------------------------------
+*  SnNat : Natural Tin
+*  Isotopic Composition : Sn112 .010  Sauf JEF22 et ENDFVI8
+*                         Sn114 .007  Sauf JEF22 et ENDFVI8
+*                         Sn115 .004
+*                         Sn116 .147 / Sn116 .164 (JEF22, ENDFVI8)
+*                         Sn117 .077
+*                         Sn118 .243
+*                         Sn119 .086
+*                         Sn120 .324
+*                         Sn122 .046
+*                         Sn124 .056
+* Pour Sn112 et Sn114 avec JEF22 ou ENDFVI8, on reporte les
+* fractions atomiques sur Sn116 (premier isotope pair defini
+* dans ces bibliotheques)
+*------------------------------------------------------
+REAL Sn112f := .0  ;
+REAL Sn114f := .0  ;
+REAL Sn115f := .0  ;
+REAL Sn116f := .0  ;
+REAL Sn117f := .0  ;
+REAL Sn118f := .0  ;
+REAL Sn119f := .0  ;
+REAL Sn120f := .0  ;
+REAL Sn122f := .0  ;
+REAL Sn124f := .0  ;
+REAL MSnNat ; ! Average Weight Ratio
+
+IF LibType 'JEFF31'   =   LibType 'JEFF311'  = +
+   LibType 'JEFF312'  = +
+   LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+
+  EVALUATE Sn112f := .010  ;
+  EVALUATE Sn114f := .007  ;
+  EVALUATE Sn115f := .004  ;
+  EVALUATE Sn116f := .147  ;
+  EVALUATE Sn117f := .077  ;
+  EVALUATE Sn118f := .243  ;
+  EVALUATE Sn119f := .086  ;
+  EVALUATE Sn120f := .324  ;
+  EVALUATE Sn122f := .046  ;
+  EVALUATE Sn124f := .056  ;
+
+  *---- SnNat Average Weight Ratio ----
+  EVALUATE MSnNat  :=        sn112 Sn112f *   sn114 Sn114f * + ;
+  EVALUATE MSnNat  := MSnNat sn115 Sn115f * + sn116 Sn116f * + ;
+  EVALUATE MSnNat  := MSnNat sn117 Sn117f * + sn118 Sn118f * + ;
+  EVALUATE MSnNat  := MSnNat sn119 Sn119f * + sn120 Sn120f * + ;
+  EVALUATE MSnNat  := MSnNat sn122 Sn122f * + sn124 Sn124f * + ;
+
+ELSE ! JEF22, ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93
+
+  EVALUATE Sn115f := .004  ;
+  EVALUATE Sn116f := .164  ;
+  EVALUATE Sn117f := .077  ;
+  EVALUATE Sn118f := .243  ;
+  EVALUATE Sn119f := .086  ;
+  EVALUATE Sn120f := .324  ;
+  EVALUATE Sn122f := .046  ;
+  EVALUATE Sn124f := .056  ;
+
+  *---- SnNat Average Weight Ratio ----
+  EVALUATE MSnNat  :=        sn115 Sn115f *   sn116 Sn116f * + ;
+  EVALUATE MSnNat  := MSnNat sn117 Sn117f * + sn118 Sn118f * + ;
+  EVALUATE MSnNat  := MSnNat sn119 Sn119f * + sn120 Sn120f * + ;
+  EVALUATE MSnNat  := MSnNat sn122 Sn122f * + sn124 Sn124f * + ;
+
+ENDIF ;
+
+*---- SnNat Table Creation ----
+SnNat := UTL:       :: CREA 'AvgMolM   ' 1 = <<MSnNat>> ;
+SnNat := UTL: SnNat :: CREA 'Sn112frIso' 1 = <<Sn112f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn114frIso' 1 = <<Sn114f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn115frIso' 1 = <<Sn115f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn116frIso' 1 = <<Sn116f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn117frIso' 1 = <<Sn117f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn118frIso' 1 = <<Sn118f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn119frIso' 1 = <<Sn119f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn120frIso' 1 = <<Sn120f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn122frIso' 1 = <<Sn122f>> ;
+SnNat := UTL: SnNat :: CREA 'Sn124frIso' 1 = <<Sn124f>> ;
+
+*------------------------------------------------------
+*  NiNat : Natural Nickel
+*  Isotopic Composition : Ni58 .6827
+*                         Ni60 .2610
+*                         Ni61 .0113
+*                         Ni62 .0359
+*                         Ni64 .0091
+*------------------------------------------------------
+REAL Ni58f := .6827 ;
+REAL Ni60f := .2610 ;
+REAL Ni61f := .0113 ;
+REAL Ni62f := .0359 ;
+REAL Ni64f := .0091 ;
+
+*---- NiNat Average Weight Ratio ----
+REAL     MNiNat ; ! Average Weight Ratio
+EVALUATE MNiNat :=        ni58  Ni58f  *   ni60  Ni60f  * +
+                          ni61  Ni61f  * + ni62  Ni62f  * +
+                          ni64  Ni64f  * + ;
+
+*---- NiNat Table Creation ----
+NiNat := UTL:       :: CREA 'AvgMolM  ' 1 = <<MNiNat>> ;
+NiNat := UTL: NiNat :: CREA 'Ni58frIso' 1 = <<Ni58f>>  ;
+NiNat := UTL: NiNat :: CREA 'Ni60frIso' 1 = <<Ni60f>>  ;
+NiNat := UTL: NiNat :: CREA 'Ni61frIso' 1 = <<Ni61f>>  ;
+NiNat := UTL: NiNat :: CREA 'Ni62frIso' 1 = <<Ni62f>>  ;
+NiNat := UTL: NiNat :: CREA 'Ni64frIso' 1 = <<Ni64f>>  ;
+
+*------------------------------------------------------
+*  BNat : Natural Boron
+*  Isotopic Composition : B10 .2
+*                         B11 .8
+*------------------------------------------------------
+REAL B10f := .2 ;
+REAL B11f := .8 ;
+IF
+Gestion "BEAV" = THEN
+EVALUATE B10f := 0.199 ;
+EVALUATE B11f := 0.801 ;
+ENDIF ;
+
+*---- BNat Average Weight Ratio ----
+REAL     MBNat ; ! Average Weight Ratio
+EVALUATE MBNat := b10 B10f * b11 B11f * + ;
+IF Gestion "BEAV" = THEN
+EVALUATE MBNat := 10.8135 ;
+ENDIF ;
+
+*---- BNat Table Creation ----
+BNat := UTL:      :: CREA 'AvgMolM ' 1 = <<MBNat>> ;
+BNat := UTL: BNat :: CREA 'B10frIso' 1 = <<B10f>>  ;
+BNat := UTL: BNat :: CREA 'B11frIso' 1 = <<B11f>>  ;
+
+*------------------------------------------------------
+*  AgNat : Natural Silver
+*  Isotopic Composition :  Ag107 .5183
+*                          Ag109 .4817
+*------------------------------------------------------
+REAL Ag107f := .5183 ;
+REAL Ag109f := .4817 ;
+
+*---- AgNat Average Weight Ratio ----
+REAL     MAgNat ;
+EVALUATE MAgNat := ag107 Ag107f * ag109 Ag109f * + ;
+
+*---- AgNat Table Creation ----
+AgNat := UTL:       :: CREA 'AvgMolM    ' 1 = <<MAgNat>> ;
+AgNat := UTL: AgNat :: CREA 'Ag107frIso ' 1 = <<Ag107f>> ;
+AgNat := UTL: AgNat :: CREA 'Ag109frIso ' 1 = <<Ag109f>> ;
+
+*------------------------------------------------------
+*  CdNat : Natural Cadmium
+*  Isotopic Composition : Cd106 .0125 *
+*                         Cd108 .0089 **
+*                         Cd110 .1251
+*                         Cd111 .1281
+*                         Cd112 .2413 **
+*                         Cd113 .1222
+*                         Cd114 .2872
+*                         Cd116 .0747
+*
+*  *  : Cd106 is not available on JEF22 and ENDFVI8. It is replaced
+*       by Cd112.
+*       For JEF22 and ENDFVI8 :
+*                     ==> Cd106 .0000
+*                     ==> Cd108 .0000
+*                     ==> Cd112 .2627
+*  ** : Cd108 is not available on JEF22 and ENDFVI8.
+*       It is replaced by Cd112 which has relatively close enough
+*       neutronic properties
+*
+*                     ==> Cd108 .0000
+*                     ==> Cd112 .2502
+*------------------------------------------------------
+REAL Cd106f := .0 ;
+REAL Cd108f := .0 ;
+REAL Cd110f := .0 ;
+REAL Cd111f := .0 ;
+REAL Cd112f := .0 ;
+REAL Cd113f := .0 ;
+REAL Cd114f := .0 ;
+REAL Cd116f := .0 ;
+REAL MCdNat ;
+
+
+IF LibType 'JEFF31'   =   LibType 'JEFF311'  = +
+   LibType 'JEFF312'  = +
+   LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+
+  EVALUATE Cd106f := .0125 ;
+  EVALUATE Cd108f := .0089 ;
+  EVALUATE Cd110f := .1251 ;
+  EVALUATE Cd111f := .1281 ;
+  EVALUATE Cd112f := .2413 ;
+  EVALUATE Cd113f := .1222 ;
+  EVALUATE Cd114f := .2872 ;
+  EVALUATE Cd116f := .0747 ;
+
+*---- CdNat Average Weight Ratio ----
+  EVALUATE MCdNat :=        cd106 Cd106f *   cd108 Cd108f * + ;
+  EVALUATE MCdNat := MCdNat cd110 Cd110f * + cd111 Cd111f * + ;
+  EVALUATE MCdNat := MCdNat cd112 Cd112f * + cd113 Cd113f * + ;
+  EVALUATE MCdNat := MCdNat cd114 Cd114f * + cd116 Cd116f * + ;
+
+ELSE ! JEF22, ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93
+
+  EVALUATE Cd110f := .1251 ;
+  EVALUATE Cd111f := .1281 ;
+  EVALUATE Cd112f := .2627 ;
+  EVALUATE Cd113f := .1222 ;
+  EVALUATE Cd114f := .2872 ;
+  EVALUATE Cd116f := .0747 ;
+
+*---- CdNat Average Weight Ratio ----
+  EVALUATE MCdNat :=        cd110 Cd110f * ;
+  EVALUATE MCdNat := MCdNat cd111 Cd111f * + cd112 Cd112f * + ;
+  EVALUATE MCdNat := MCdNat cd113 Cd113f * + cd114 Cd114f * + ;
+  EVALUATE MCdNat := MCdNat cd116 Cd116f * + ;
+
+ENDIF ;
+
+*---- CdNat Table Creation ----
+CdNat := UTL:       :: CREA 'AvgMolM    ' 1 = <<MCdNat>> ;
+CdNat := UTL: CdNat :: CREA 'Cd106frIso ' 1 = <<Cd106f>> ;
+CdNat := UTL: CdNat :: CREA 'Cd108frIso ' 1 = <<Cd108f>> ;
+CdNat := UTL: CdNat :: CREA 'Cd110frIso ' 1 = <<Cd110f>> ;
+CdNat := UTL: CdNat :: CREA 'Cd111frIso ' 1 = <<Cd111f>> ;
+CdNat := UTL: CdNat :: CREA 'Cd112frIso ' 1 = <<Cd112f>> ;
+CdNat := UTL: CdNat :: CREA 'Cd113frIso ' 1 = <<Cd113f>> ;
+CdNat := UTL: CdNat :: CREA 'Cd114frIso ' 1 = <<Cd114f>> ;
+CdNat := UTL: CdNat :: CREA 'Cd116frIso ' 1 = <<Cd116f>> ;
+
+*------------------------------------------------------
+*  KNat : Natural Potassium
+*  Reference : www.wikipedia.org/wiki/Potassium
+*  Isotopic Composition : K39 .9326
+*                         K40 .0001167
+*                         K41 .0673
+*------------------------------------------------------
+REAL K39f := .0 ;
+REAL K40f := .0 ;
+REAL K41f := .0 ;
+REAL MKNat ; ! Average Weight Ratio
+
+
+IF LibType 'JEFF31'   =   LibType 'JEFF311' = +
+   LibType 'JEFF312'  = +
+   LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+
+  EVALUATE K39f := .9326 ;
+  EVALUATE K40f := .0001167 ;
+  EVALUATE K41f := .0673 ;
+
+*---- KNat Average Weight Ratio ----
+  EVALUATE MKNat := k39  K39f  *   k40  K40f  * + k41  K41f  * + ;
+
+  ECHO "Masse molaire du potassium " MKNat ;
+
+ELSE ! JEF22, ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93
+
+  EVALUATE MKNat := kNat ;
+
+ENDIF ;
+
+*---- KNat Table Creation ----
+KNat := UTL:      :: CREA 'AvgMolM  ' 1 = <<MKNat>> ;
+KNat := UTL: KNat :: CREA 'K39frIso' 1 = <<K39f>>  ;
+KNat := UTL: KNat :: CREA 'K40frIso' 1 = <<K40f>>  ;
+KNat := UTL: KNat :: CREA 'K41frIso' 1 = <<K41f>>  ;
+
+
+*------------------------------------------------------
+*  SNat : Natural Sulfur
+*  Reference : www.wikipedia.org/wiki/Soufre
+*  Isotopic Composition :  S32 0.9502
+*                          S33 0.0076
+*                          S34 0.0422
+*  NB: pas de S33 et S34 dans ENDFBVI8
+* Pour S34 dans ENDFVI8, on reporte les
+* fractions atomiques sur S32 (premier isotope pair defini
+* dans ces bibliotheques)
+*------------------------------------------------------
+
+REAL S32f := .0  ;
+REAL S33f := .0  ;
+REAL S34f := .0  ;
+
+REAL MSNat ; ! Average Weight Ratio
+ECHO "libtype:" LibType ;
+IF LibType 'JEFF31'   =   LibType 'JEFF311'  = +
+   LibType 'JEFF312'  = + LibType 'JEFF22'  = +
+   LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN
+
+  EVALUATE S32f := .9502  ;
+  EVALUATE S33f := .0076  ;
+  EVALUATE S34f := .0422  ;
+
+*---- SNat Average Weight Ratio ----
+  EVALUATE MSNat := S32f s32 * S33f s33 * + ;
+  EVALUATE MSNat := MSNat S34f s34 * + ;
+
+ELSE ! ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93
+
+  EVALUATE S32f := .9924  ;
+
+  *---- SnNat Average Weight Ratio ----
+  EVALUATE MSNat  := s32 S32f * ;
+
+ENDIF ;
+
+*---- SNat Table Creation ----
+SNat := UTL:      :: CREA 'AvgMolM ' 1 = <<MSNat>> ;
+SNat := UTL: SNat :: CREA 'S32frIso' 1 = <<S32f>>  ;
+SNat := UTL: SNat :: CREA 'S33frIso' 1 = <<S33f>>  ;
+SNat := UTL: SNat :: CREA 'S34frIso' 1 = <<S34f>>  ;
+
+*------------------------------------------------------
+*  InNat : Natural Indium
+*  Reference : www.webelements.com/indium/isotopes.html
+*  Isotopic Composition : In115 0.9571
+*                         In113 0.0429
+* NB : In113 non disponible. On considere que
+* l'InNat correspond a 100% de In115
+*------------------------------------------------------
+REAL In113f := 0. ;
+REAL In115f := 1. ;
+REAL MInNat ; ! Average Weight Ratio
+
+*---- SNat Average Weight Ratio ----
+EVALUATE MInNat := In115f in115 * ;
+
+*---- InNat Table Creation ----
+InNat := UTL:       :: CREA 'AvgMolM   ' 1 = <<MInNat>> ;
+InNat := UTL: InNat :: CREA 'In115frIso' 1 = <<In115f>>  ;
+InNat := UTL: InNat :: CREA 'In113frIso' 1 = <<In113f>>  ;
+
+*------------------------------------------------------
+*  Zirc4 : Zircalloy4
+*  Massic Composition : FeNat .002
+*                       CrNat .001
+*                       016   .001
+*                       ZrNat .981
+*                       SnNat .015
+*------------------------------------------------------
+REAL FeZ4f  := .002 ;
+REAL CrZ4f  := .001 ;
+REAL O16Z4f := .001 ;
+REAL ZrZ4f  := .981 ;
+REAL SnZ4f  := .015 ;
+
+REAL Zr90Z4f  Zr91Z4f  Zr92Z4f  Zr93Z4f
+     Zr94Z4f  Zr95Z4f  Zr96Z4f ;
+REAL Fe54Z4f  Fe56Z4f  Fe57Z4f  Fe58Z4f ;
+REAL Cr50Z4f  Cr52Z4f  Cr53Z4f  Cr54Z4f ;
+REAL Sn112Z4f Sn114Z4f Sn115Z4f Sn116Z4f Sn117Z4f
+     Sn118Z4f Sn119Z4f Sn120Z4f Sn122Z4f Sn124Z4f ;
+
+*---- Mass Ratios -> Isotopic Ratios ----
+EVALUATE FeZ4f  := FeZ4f  avo_conc * MFeNat / ;
+EVALUATE CrZ4f  := CrZ4f  avo_conc * MCrNat / ;
+EVALUATE O16Z4f := O16Z4f avo_conc * o16    / ;
+EVALUATE ZrZ4f  := ZrZ4f  avo_conc * MZrNat / ;
+EVALUATE SnZ4f  := SnZ4f  avo_conc * MSnNat / ;
+
+EVALUATE ReNorm := FeZ4f CrZ4f + O16Z4f + ZrZ4f + SnZ4f + ;
+
+EVALUATE FeZ4f  := FeZ4f  ReNorm / ;
+EVALUATE CrZ4f  := CrZ4f  ReNorm / ;
+EVALUATE O16Z4f := O16Z4f ReNorm / ;
+EVALUATE ZrZ4f  := ZrZ4f  ReNorm / ;
+EVALUATE SnZ4f  := SnZ4f  ReNorm / ;
+
+*---- Zirc4 Average Weight Ratio ---
+REAL     MZirc4 ; ! Average Weight Ratio
+EVALUATE MZirc4 := FeZ4f  MFeNat *   CrZ4f MCrNat * +
+                   O16Z4f o16    * + ZrZ4f MZrNat  * +
+                   SnZ4f  MSnNat * + ;
+
+ECHO "Masse molaire du zircaloy " MZirc4 ;
+
+*---- Isotopic Ratios for all Zr, Fe, Cr & Sn ----
+* Rq : si les isotopes n'existent pas ce n'est pas grave,
+*      leur valeur sera egale a 0. dans tous les calculs
+*      et leur appel sera commente dans l'appel a LIB de
+*      la procedure CreaMix.c2m. Leur presence n'a aucune
+*      influence sur les calculs du fait de leur valeur = 0.
+*-------------------------------------------------
+EVALUATE Zr90Z4f  := ZrZ4f Zr90f  * ;
+EVALUATE Zr91Z4f  := ZrZ4f Zr91f  * ;
+EVALUATE Zr92Z4f  := ZrZ4f Zr92f  * ;
+EVALUATE Zr93Z4f  := ZrZ4f Zr93f  * ;
+EVALUATE Zr94Z4f  := ZrZ4f Zr94f  * ;
+EVALUATE Zr95Z4f  := ZrZ4f Zr95f  * ;
+EVALUATE Zr96Z4f  := ZrZ4f Zr96f  * ;
+
+EVALUATE Fe54Z4f  := FeZ4f Fe54f  * ;
+EVALUATE Fe56Z4f  := FeZ4f Fe56f  * ;
+EVALUATE Fe57Z4f  := FeZ4f Fe57f  * ;
+EVALUATE Fe58Z4f  := FeZ4f Fe58f  * ;
+
+EVALUATE Cr50Z4f  := CrZ4f Cr50f  * ;
+EVALUATE Cr52Z4f  := CrZ4f Cr52f  * ;
+EVALUATE Cr53Z4f  := CrZ4f Cr53f  * ;
+EVALUATE Cr54Z4f  := CrZ4f Cr54f  * ;
+
+EVALUATE Sn112Z4f := SnZ4f Sn112f * ;
+EVALUATE Sn114Z4f := SnZ4f Sn114f * ;
+EVALUATE Sn115Z4f := SnZ4f Sn115f * ;
+EVALUATE Sn116Z4f := SnZ4f Sn116f * ;
+EVALUATE Sn117Z4f := SnZ4f Sn117f * ;
+EVALUATE Sn118Z4f := SnZ4f Sn118f * ;
+EVALUATE Sn119Z4f := SnZ4f Sn119f * ;
+EVALUATE Sn120Z4f := SnZ4f Sn120f * ;
+EVALUATE Sn122Z4f := SnZ4f Sn122f * ;
+EVALUATE Sn124Z4f := SnZ4f Sn124f * ;
+
+*---- Zirc4 Table Creation ----
+
+Zirc4 := UTL:       :: CREA 'AvgMolM     ' 1 = <<MZirc4>>   ;
+Zirc4 := UTL: Zirc4 :: CREA 'O16Z4frIso  ' 1 = <<O16Z4f>>   ;
+Zirc4 := UTL: Zirc4 :: CREA 'Fe54Z4frIso ' 1 = <<Fe54Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Fe56Z4frIso ' 1 = <<Fe56Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Fe57Z4frIso ' 1 = <<Fe57Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Fe58Z4frIso ' 1 = <<Fe58Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Cr50Z4frIso ' 1 = <<Cr50Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Cr52Z4frIso ' 1 = <<Cr52Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Cr53Z4frIso ' 1 = <<Cr53Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Cr54Z4frIso ' 1 = <<Cr54Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'ZrNatfrIso  ' 1 = <<ZrZ4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Zr90Z4frIso ' 1 = <<Zr90Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Zr91Z4frIso ' 1 = <<Zr91Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Zr92Z4frIso ' 1 = <<Zr92Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Zr93Z4frIso ' 1 = <<Zr93Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Zr94Z4frIso ' 1 = <<Zr94Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Zr95Z4frIso ' 1 = <<Zr95Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Zr96Z4frIso ' 1 = <<Zr96Z4f>>  ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn112Z4frIso' 1 = <<Sn112Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn114Z4frIso' 1 = <<Sn114Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn115Z4frIso' 1 = <<Sn115Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn116Z4frIso' 1 = <<Sn116Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn117Z4frIso' 1 = <<Sn117Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn118Z4frIso' 1 = <<Sn118Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn119Z4frIso' 1 = <<Sn119Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn120Z4frIso' 1 = <<Sn120Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn122Z4frIso' 1 = <<Sn122Z4f>> ;
+Zirc4 := UTL: Zirc4 :: CREA 'Sn124Z4frIso' 1 = <<Sn124Z4f>> ;
+
+*------------------------------------------------------
+*  Incon : Inconel - AFA 2G
+*  Massic Composition : NiNat .5250
+*                       CrNat .1900
+*                       FeNat .1823
+*                       BNatB .00003
+*                       Al27  .0050
+*                       Mo95  .0305
+*                       TiNat .0088
+*                       SiNat .0018
+*                       Co59  .0005
+*                       Mn55  .0018
+*                       Nb93  .05175
+*                      (&IMPU .00252)
+*------------------------------------------------------
+REAL NiIncf := .5250  ;
+REAL CrIncf := .1900  ;
+REAL FeIncf := .1823  ;
+REAL BIncf  := .00003 ;
+REAL AlIncf := .0050  ;
+REAL MoIncf := .0305  ;
+REAL TiIncf := .0088  ;
+REAL SiIncf := .0018  ;
+REAL CoIncf := .0005  ;
+REAL MnIncf := .0018  ;
+REAL NbIncf := .05175 ;
+
+IF Gestion "BEAV" = THEN
+EVALUATE NiIncf := .5119  ;
+EVALUATE CrIncf := .1896  ;
+EVALUATE FeIncf := .2863  ;
+EVALUATE BIncf  := 1E-15 ;
+EVALUATE AlIncf := 1E-15  ;
+EVALUATE MoIncf := 1E-15  ;
+EVALUATE TiIncf := 1E-15  ;
+EVALUATE SiIncf := .0035  ;
+EVALUATE CoIncf := .0005  ;
+EVALUATE MnIncf := .0087  ;
+EVALUATE NbIncf := 1E-15 ;
+ENDIF ;
+
+REAL Ni58Inf Ni60Inf Ni61Inf Ni62Inf Ni64Inf ;
+REAL Cr50Inf Cr52Inf Cr53Inf Cr54Inf ;
+REAL Fe54Inf Fe56Inf Fe57Inf Fe58Inf ;
+REAL Ti46Inf Ti47Inf Ti48Inf Ti49Inf Ti50Inf ;
+REAL Si28Inf Si29Inf Si30Inf ;
+REAL B10Inf  B11Inf  ;
+
+*---- Fraction Sum ----
+EVALUATE SumFr  := NiIncf CrIncf + FeIncf + BIncf  + AlIncf + MoIncf + ;
+EVALUATE SumFr  := SumFr  TiIncf + SiIncf + CoIncf + MnIncf + NbIncf + ;
+
+*---- Mass Ratios -> Isotopic Ratios ----
+EVALUATE NiIncf := NiIncf avo_conc * MNiNat / ;
+EVALUATE CrIncf := CrIncf avo_conc * MCrNat / ;
+EVALUATE FeIncf := FeIncf avo_conc * MFeNat / ;
+EVALUATE BIncf  := BIncf  avo_conc * MBNat  / ;
+EVALUATE AlIncf := AlIncf avo_conc * al27   / ;
+EVALUATE MoIncf := MoIncf avo_conc * mo95   / ;
+EVALUATE TiIncf := TiIncf avo_conc * MTiNat  / ;
+EVALUATE SiIncf := SiIncf avo_conc * MSiNat  / ;
+EVALUATE CoIncf := CoIncf avo_conc * co59   / ;
+EVALUATE MnIncf := MnIncf avo_conc * mn55   / ;
+EVALUATE NbIncf := NbIncf avo_conc * nb93   / ;
+
+EVALUATE ReNorm := NiIncf CrIncf + FeIncf + BIncf  + AlIncf + MoIncf + ;
+EVALUATE ReNorm := ReNorm TiIncf + SiIncf + CoIncf + MnIncf + NbIncf + ;
+EVALUATE ReNorm := ReNorm SumFr  / ;
+
+EVALUATE NiIncf := NiIncf ReNorm / ;
+EVALUATE CrIncf := CrIncf ReNorm / ;
+EVALUATE FeIncf := FeIncf ReNorm / ;
+EVALUATE BIncf  := BIncf  ReNorm / ;
+EVALUATE AlIncf := AlIncf ReNorm / ;
+EVALUATE MoIncf := MoIncf ReNorm / ;
+EVALUATE TiIncf := TiIncf ReNorm / ;
+EVALUATE SiIncf := SiIncf ReNorm / ;
+EVALUATE CoIncf := CoIncf ReNorm / ;
+EVALUATE MnIncf := MnIncf ReNorm / ;
+EVALUATE NbIncf := NbIncf ReNorm / ;
+
+*---- Isotopic Ratios for all Ni, Cr, Fe & B ----
+EVALUATE Ni58Inf := Ni58f NiIncf * ;
+EVALUATE Ni60Inf := Ni60f NiIncf * ;
+EVALUATE Ni61Inf := Ni61f NiIncf * ;
+EVALUATE Ni62Inf := Ni62f NiIncf * ;
+EVALUATE Ni64Inf := Ni64f NiIncf * ;
+
+EVALUATE Ti46Inf := Ti46f TiIncf * ;
+EVALUATE Ti47Inf := Ti47f TiIncf * ;
+EVALUATE Ti48Inf := Ti48f TiIncf * ;
+EVALUATE Ti49Inf := Ti49f TiIncf * ;
+EVALUATE Ti50Inf := Ti50f TiIncf * ;
+
+EVALUATE Si28Inf := Si28f SiIncf * ;
+EVALUATE Si29Inf := Si29f SiIncf * ;
+EVALUATE Si30Inf := Si30f SiIncf * ;
+
+EVALUATE Cr50Inf := Cr50f CrIncf * ;
+EVALUATE Cr52Inf := Cr52f CrIncf * ;
+EVALUATE Cr53Inf := Cr53f CrIncf * ;
+EVALUATE Cr54Inf := Cr54f CrIncf * ;
+
+EVALUATE Fe54Inf := Fe54f FeIncf * ;
+EVALUATE Fe56Inf := Fe56f FeIncf * ;
+EVALUATE Fe57Inf := Fe57f FeIncf * ;
+EVALUATE Fe58Inf := Fe58f FeIncf * ;
+
+EVALUATE B10Inf  := B10f  BIncf  * ;
+EVALUATE B11Inf  := B11f  BIncf  * ;
+
+*---- Inconel Average Weight Ratio ---
+REAL     MIncon ; ! Average Weight Ratio
+EVALUATE MIncon :=  NiIncf MNiNat *   CrIncf MCrNat * +
+                    FeIncf MFeNat * + BIncf  MBNat  * +
+                    AlIncf al27   * + MoIncf mo95   * +
+                    TiIncf MTiNat * + SiIncf MSiNat * +
+                    CoIncf co59   * + MnIncf mn55   * +
+                    NbIncf nb93   * + ;
+
+ECHO "Masse molaire de l'inconel " MIncon ;
+
+*---- Due to &IMPU, fractions are summing to "1-&IMPU".
+*     These fractions are stored, as they represent the correct proportion.
+*     MIncon however has to be renormalized so that fractions are summing to 1. ----
+EVALUATE MIncon := MIncon SumFr  / ;
+
+*---- Incon Table Creation ----
+
+Incon := UTL:       :: CREA 'AvgMolM     ' 1 = <<MIncon>>  ;
+Incon := UTL: Incon :: CREA 'B10InfrIso  ' 1 = <<B10Inf>>  ;
+Incon := UTL: Incon :: CREA 'B11InfrIso  ' 1 = <<B11Inf>>  ;
+Incon := UTL: Incon :: CREA 'Al27InfrIso ' 1 = <<AlIncf>>  ;
+Incon := UTL: Incon :: CREA 'SiNatInfrIso' 1 = <<SiIncf>>  ;
+Incon := UTL: Incon :: CREA 'Si28InfrIso ' 1 = <<Si28Inf>> ;
+Incon := UTL: Incon :: CREA 'Si29InfrIso ' 1 = <<Si29Inf>> ;
+Incon := UTL: Incon :: CREA 'Si30InfrIso ' 1 = <<Si30Inf>> ;
+Incon := UTL: Incon :: CREA 'TiNatInfrIso' 1 = <<TiIncf>>  ;
+Incon := UTL: Incon :: CREA 'Ti46InfrIso ' 1 = <<Ti46Inf>> ;
+Incon := UTL: Incon :: CREA 'Ti47InfrIso ' 1 = <<Ti47Inf>> ;
+Incon := UTL: Incon :: CREA 'Ti48InfrIso ' 1 = <<Ti48Inf>> ;
+Incon := UTL: Incon :: CREA 'Ti49InfrIso ' 1 = <<Ti49Inf>> ;
+Incon := UTL: Incon :: CREA 'Ti50InfrIso ' 1 = <<Ti50Inf>> ;
+Incon := UTL: Incon :: CREA 'Cr50InfrIso ' 1 = <<Cr50Inf>> ;
+Incon := UTL: Incon :: CREA 'Cr52InfrIso ' 1 = <<Cr52Inf>> ;
+Incon := UTL: Incon :: CREA 'Cr53InfrIso ' 1 = <<Cr53Inf>> ;
+Incon := UTL: Incon :: CREA 'Cr54InfrIso ' 1 = <<Cr54Inf>> ;
+Incon := UTL: Incon :: CREA 'Mn55InfrIso ' 1 = <<MnIncf>>  ;
+Incon := UTL: Incon :: CREA 'Fe54InfrIso ' 1 = <<Fe54Inf>> ;
+Incon := UTL: Incon :: CREA 'Fe56InfrIso ' 1 = <<Fe56Inf>> ;
+Incon := UTL: Incon :: CREA 'Fe57InfrIso ' 1 = <<Fe57Inf>> ;
+Incon := UTL: Incon :: CREA 'Fe58InfrIso ' 1 = <<Fe58Inf>> ;
+Incon := UTL: Incon :: CREA 'Co59InfrIso ' 1 = <<CoIncf>>  ;
+Incon := UTL: Incon :: CREA 'Ni58InfrIso ' 1 = <<Ni58Inf>> ;
+Incon := UTL: Incon :: CREA 'Ni60InfrIso ' 1 = <<Ni60Inf>> ;
+Incon := UTL: Incon :: CREA 'Ni61InfrIso ' 1 = <<Ni61Inf>> ;
+Incon := UTL: Incon :: CREA 'Ni62InfrIso ' 1 = <<Ni62Inf>> ;
+Incon := UTL: Incon :: CREA 'Ni64InfrIso ' 1 = <<Ni64Inf>> ;
+Incon := UTL: Incon :: CREA 'Nb93InfrIso ' 1 = <<NbIncf>>  ;
+Incon := UTL: Incon :: CREA 'Mo95InfrIso ' 1 = <<MoIncf>>  ;
+
+*------------------------------------------------------
+*  SS304 : control rod clad -- Stainless Steel
+*  Massic Composition : FeNat .706
+*                       NiNat .100
+*                       CrNat .180
+*                       Mn55  .010
+*                       SiNat .004
+*------------------------------------------------------
+REAL FeSSf  := .706 ;
+REAL NiSSf  := .100 ;
+REAL CrSSf  := .180 ;
+REAL MnSSf  := .010 ;
+REAL SiSSf  := .004 ;
+
+IF Gestion "BEAV" = THEN
+ EVALUATE FeSSf  := .684 ;
+ EVALUATE NiSSf  := .100 ;
+ EVALUATE CrSSf  := .190 ;
+ EVALUATE MnSSf  := .020 ;
+ EVALUATE SiSSf  := .006 ;
+ENDIF ;
+
+REAL Fe54SSf  Fe56SSf  Fe57SSf  Fe58SSf  ;
+REAL Ni58SSf  Ni60SSf  Ni61SSf  Ni62SSf  Ni64SSf ;
+REAL Cr50SSf  Cr52SSf  Cr53SSf  Cr54SSf  ;
+REAL Si28SSf  Si29SSf  Si30SSf ;
+
+*---- Mass Ratios -> Isotopic Ratios ----
+EVALUATE FeSSf  := FeSSf  avo_conc * MFeNat / ;
+EVALUATE NiSSf  := NiSSf  avo_conc * MNiNat / ;
+EVALUATE CrSSf  := CrSSf  avo_conc * MCrNat / ;
+EVALUATE MnSSf  := MnSSf  avo_conc * mn55   / ;
+EVALUATE SiSSf  := SiSSf  avo_conc * MSiNat  / ;
+
+EVALUATE ReNorm := FeSSf NiSSf + CrSSf + MnSSf + SiSSf + ;
+
+EVALUATE FeSSf  := FeSSf  ReNorm / ;
+EVALUATE NiSSf  := NiSSf  ReNorm / ;
+EVALUATE CrSSf  := CrSSf  ReNorm / ;
+EVALUATE MnSSf  := MnSSf  ReNorm / ;
+EVALUATE SiSSf  := SiSSf  ReNorm / ;
+
+*---- SS304 Average Weight Ratio ---
+REAL     MSS304 ; ! Average Weight Ratio
+EVALUATE MSS304 :=        FeSSf  MFeNat *   NiSSf MNiNat * +
+                          CrSSf  MCrNat * + MnSSf mn55   * +
+                          SiSSf  MSiNat * + ;
+
+ECHO "Masse molaire du SS304 " MSS304 ;
+
+*---- Isotopic Ratios for all Fe, Ni & Cr ----
+EVALUATE Fe54SSf := FeSSf Fe54f * ;
+EVALUATE Fe56SSf := FeSSf Fe56f * ;
+EVALUATE Fe57SSf := FeSSf Fe57f * ;
+EVALUATE Fe58SSf := FeSSf Fe58f * ;
+
+EVALUATE Ni58SSf := NiSSf Ni58f * ;
+EVALUATE Ni60SSf := NiSSf Ni60f * ;
+EVALUATE Ni61SSf := NiSSf Ni61f * ;
+EVALUATE Ni62SSf := NiSSf Ni62f * ;
+EVALUATE Ni64SSf := NiSSf Ni64f * ;
+
+EVALUATE Cr50SSf := CrSSf Cr50f * ;
+EVALUATE Cr52SSf := CrSSf Cr52f * ;
+EVALUATE Cr53SSf := CrSSf Cr53f * ;
+EVALUATE Cr54SSf := CrSSf Cr54f * ;
+
+EVALUATE Si28SSf  := SiSSf Si28f * ;
+EVALUATE Si29SSf  := SiSSf Si29f * ;
+EVALUATE Si30SSf  := SiSSf Si30f * ;
+
+*---- SS304 Table Creation ----
+
+* Isotopes non conflictuels
+SS304 := UTL:       :: CREA 'AvgMolM     ' 1 = <<MSS304>>  ;
+SS304 := UTL: SS304 :: CREA 'SiNatSSfrIso' 1 = <<SiSSf>> ;
+SS304 := UTL: SS304 :: CREA 'Si28SSfrIso ' 1 = <<Si28SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Si29SSfrIso ' 1 = <<Si29SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Si30SSfrIso ' 1 = <<Si30SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Fe54SSfrIso ' 1 = <<Fe54SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Fe56SSfrIso ' 1 = <<Fe56SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Fe57SSfrIso ' 1 = <<Fe57SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Fe58SSfrIso ' 1 = <<Fe58SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Ni58SSfrIso ' 1 = <<Ni58SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Ni60SSfrIso ' 1 = <<Ni60SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Ni61SSfrIso ' 1 = <<Ni61SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Ni62SSfrIso ' 1 = <<Ni62SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Ni64SSfrIso ' 1 = <<Ni64SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Cr50SSfrIso ' 1 = <<Cr50SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Cr52SSfrIso ' 1 = <<Cr52SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Cr53SSfrIso ' 1 = <<Cr53SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Cr54SSfrIso ' 1 = <<Cr54SSf>> ;
+SS304 := UTL: SS304 :: CREA 'Mn55SSfrIso ' 1 = <<MnSSf>>   ;
+
+*------------------------------------------------------
+*  M5
+*  Ref: C. Lemaignan "Zirconium Alloys: Properties
+*                     and Characteristics"
+*  Massic Composition : FeNat .00050
+*                       O16   .00135
+*                       ZrNat .98813
+*                       Nb93  .01000
+*                       SNat  .00002
+*------------------------------------------------------
+REAL FeM5f   := 0.00050 ;
+REAL O16M5f  := 0.00135 ;
+REAL ZrM5f   := 0.98813 ;
+REAL Nb93M5f := 0.01000 ;
+REAL SM5f    := 0.00002 ;
+
+REAL Zr90M5f Zr91M5f Zr92M5f Zr93M5f Zr94M5f Zr95M5f Zr96M5f ;
+REAL Fe54M5f Fe56M5f Fe57M5f Fe58M5f ;
+REAL S32M5f  S33M5f  S34M5f          ;
+
+*---- Mass Ratios -> Isotopic Ratios ----
+EVALUATE FeM5f   := FeM5f   avo_conc * MFeNat / ;
+EVALUATE O16M5f  := O16M5f  avo_conc * o16    / ;
+EVALUATE ZrM5f   := ZrM5f   avo_conc * MZrNat  / ;
+EVALUATE Nb93M5f := Nb93M5f avo_conc * nb93   / ;
+EVALUATE SM5f    := SM5f    avo_conc * MSNat  / ;
+
+EVALUATE ReNorm := FeM5f O16M5f + ZrM5f + Nb93M5f + SM5f + ;
+
+EVALUATE FeM5f   := FeM5f   ReNorm / ;
+EVALUATE O16M5f  := O16M5f  ReNorm / ;
+EVALUATE ZrM5f   := ZrM5f   ReNorm / ;
+EVALUATE Nb93M5f := Nb93M5f ReNorm / ;
+EVALUATE SM5f    := SM5f    ReNorm / ;
+
+*---- M5 Average Weight Ratio ----
+REAL MM5 ; ! Average Weight Ratio
+EVALUATE MM5 :=     FeM5f MFeNat *   O16M5f  o16  * + ;
+EVALUATE MM5 := MM5 ZrM5f MZrNat  * + Nb93M5f nb93 * + ;
+EVALUATE MM5 := MM5 SM5f  MSNat  * + ;
+
+*---- Isotopic Ratios for all Zr & Fe & S ----
+EVALUATE Zr90M5f := ZrM5f Zr90f * ;
+EVALUATE Zr91M5f := ZrM5f Zr91f * ;
+EVALUATE Zr92M5f := ZrM5f Zr92f * ;
+EVALUATE Zr93M5f := ZrM5f Zr93f * ;
+EVALUATE Zr94M5f := ZrM5f Zr94f * ;
+EVALUATE Zr95M5f := ZrM5f Zr95f * ;
+EVALUATE Zr96M5f := ZrM5f Zr96f * ;
+
+EVALUATE Fe54M5f := FeM5f Fe54f * ;
+EVALUATE Fe56M5f := FeM5f Fe56f * ;
+EVALUATE Fe57M5f := FeM5f Fe57f * ;
+EVALUATE Fe58M5f := FeM5f Fe58f * ;
+
+EVALUATE S32M5f  := SM5f  S32f  * ;
+EVALUATE S33M5f  := SM5f  S33f  * ;
+EVALUATE S34M5f  := SM5f  S34f  * ;
+
+*---- M5 Table Creation ----
+M5 := UTL:    :: CREA 'AvgMolM    ' 1 = <<MM5>>     ;
+M5 := UTL: M5 :: CREA 'Fe54M5frIso' 1 = <<Fe54M5f>> ;
+M5 := UTL: M5 :: CREA 'Fe56M5frIso' 1 = <<Fe56M5f>> ;
+M5 := UTL: M5 :: CREA 'Fe57M5frIso' 1 = <<Fe57M5f>> ;
+M5 := UTL: M5 :: CREA 'Fe58M5frIso' 1 = <<Fe58M5f>> ;
+M5 := UTL: M5 :: CREA 'O16M5frIso ' 1 = <<O16M5f>>  ;
+M5 := UTL: M5 :: CREA 'ZrM5frIso  ' 1 = <<ZrM5f>>   ;
+M5 := UTL: M5 :: CREA 'Zr90M5frIso' 1 = <<Zr90M5f>> ;
+M5 := UTL: M5 :: CREA 'Zr91M5frIso' 1 = <<Zr91M5f>> ;
+M5 := UTL: M5 :: CREA 'Zr92M5frIso' 1 = <<Zr92M5f>> ;
+M5 := UTL: M5 :: CREA 'Zr93M5frIso' 1 = <<Zr93M5f>> ;
+M5 := UTL: M5 :: CREA 'Zr94M5frIso' 1 = <<Zr94M5f>> ;
+M5 := UTL: M5 :: CREA 'Zr95M5frIso' 1 = <<Zr95M5f>> ;
+M5 := UTL: M5 :: CREA 'Zr96M5frIso' 1 = <<Zr96M5f>> ;
+M5 := UTL: M5 :: CREA 'Nb93M5frIso' 1 = <<Nb93M5f>> ;
+M5 := UTL: M5 :: CREA 'S32M5frIso ' 1 = <<S32M5f>>  ;
+M5 := UTL: M5 :: CREA 'S33M5frIso ' 1 = <<S33M5f>>  ;
+M5 := UTL: M5 :: CREA 'S34M5frIso ' 1 = <<S34M5f>>  ;
+
+*------------------------------------------------------
+*  AIC : control rod neutron absorber -- AgInCd
+*  Massic Composition : AgNat .80611
+*                       CdNat .04909
+*                       In115 (1-(Ag+Cd))
+*  Dans le cas de ENDFVI8 : In115 est en realite InNat
+*------------------------------------------------------
+REAL AgAICf := .80611 ;
+REAL CdAICf := .04909 ;
+REAL InAICf := 1. AgAICf - CdAICf - ;
+
+REAL Ag107AICf Ag109AICf ;
+REAL Cd106AICf Cd108AICf Cd110AICf Cd111AICf Cd112AICf Cd113AICf ;
+REAL Cd114AICf Cd116AICf ;
+REAL In115AICf In113AICf ;
+
+*---- Mass Ratios -> Isotopic Ratios ----
+EVALUATE AgAICf := AgAICf avo_conc * MAgNat / ;
+EVALUATE CdAICf := CdAICf avo_conc * MCdNat / ;
+*EVALUATE InAICf := InAICf avo_conc * in115  / ;
+EVALUATE InAICf := InAICf avo_conc * MInNat  / ;
+
+EVALUATE ReNorm := AgAICf CdAICf + InAICf + ;
+
+EVALUATE AgAICf := AgAICf ReNorm / ;
+EVALUATE CdAICf := CdAICf ReNorm / ;
+EVALUATE InAICf := InAICf ReNorm / ;
+
+*---- AIC Average Weight Ratio ---
+REAL     MAIC   ; ! Average Weight Ratio
+EVALUATE MAIC   :=        AgAICf MAgNat *   CdAICf MCdNat * + ;
+*EVALUATE MAIC   := MAIC   InAICf in115  * + ;
+EVALUATE MAIC   := MAIC   InAICf MInNat  * + ;
+
+*---- Isotopic Ratios for all Ag & Cd ----
+EVALUATE Ag107AICf := AgAICf Ag107f * ;
+EVALUATE Ag109AICf := AgAICf Ag109f * ;
+
+EVALUATE Cd106AICf := CdAICf Cd106f * ;
+EVALUATE Cd108AICf := CdAICf Cd108f * ;
+EVALUATE Cd110AICf := CdAICf Cd110f * ;
+EVALUATE Cd111AICf := CdAICf Cd111f * ;
+EVALUATE Cd112AICf := CdAICf Cd112f * ;
+EVALUATE Cd113AICf := CdAICf Cd113f * ;
+EVALUATE Cd114AICf := CdAICf Cd114f * ;
+EVALUATE Cd116AICf := CdAICf Cd116f * ;
+
+EVALUATE In115AICf := InAICf In115f * ;
+EVALUATE In113AICf := InAICf In113f * ;
+
+*---- AIC Table Creation ----
+AIC := UTL:     :: CREA 'AvgMolM     ' 1 = <<MAIC>>      ;
+AIC := UTL: AIC :: CREA 'Ag107AICfIso' 1 = <<Ag107AICf>> ;
+AIC := UTL: AIC :: CREA 'Ag109AICfIso' 1 = <<Ag109AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd106AICfIso' 1 = <<Cd106AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd108AICfIso' 1 = <<Cd108AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd110AICfIso' 1 = <<Cd110AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd111AICfIso' 1 = <<Cd111AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd112AICfIso' 1 = <<Cd112AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd113AICfIso' 1 = <<Cd113AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd114AICfIso' 1 = <<Cd114AICf>> ;
+AIC := UTL: AIC :: CREA 'Cd116AICfIso' 1 = <<Cd116AICf>> ;
+AIC := UTL: AIC :: CREA 'In115AICfIso' 1 = <<In115AICf>> ;
+AIC := UTL: AIC :: CREA 'In113AICfIso' 1 = <<In113AICf>> ;
+
+*------------------------------------------------------
+*  GdNat : Natural Gadolinium
+*  Isotopic Composition : Gd152 .002 ***
+*                         Gd154 .021 ***
+*                         Gd155 .148
+*                         Gd156 .206
+*                         Gd157 .157
+*                         Gd158 .248
+*                         Gd160 .218
+*
+*  *** : Gd152 is not available. It is replaced by Gd154
+*        which has relatively close enough neutronic properties
+*
+*                      => Gd152 .000
+*                      => Gd154 .023
+*------------------------------------------------------
+*REAL Gd152f := .000 ;
+REAL Gd154f := .023 ;
+REAL Gd155f := .148 ;
+REAL Gd156f := .206 ;
+REAL Gd157f := .157 ;
+REAL Gd158f := .248 ;
+REAL Gd160f := .218 ;
+
+*---- GdNat Average Weight Ratio ----
+REAL     MGdNat ; ! Average Weight Ratio
+EVALUATE MGdNat :=        gd154 Gd154f *   gd155 Gd155f * + ;
+EVALUATE MGdNat := MGdNat gd156 Gd156f * + gd157 Gd157f * + ;
+EVALUATE MGdNat := MGdNat gd158 Gd158f * + gd160 Gd160f * + ;
+
+*---- GdNat Table Creation ----
+GdNat := UTL:       :: CREA 'AvgMolM   ' 1 = <<MGdNat>> ;
+GdNat := UTL: GdNat :: CREA 'Gd154frIso' 1 = <<Gd154f>> ;
+GdNat := UTL: GdNat :: CREA 'Gd155frIso' 1 = <<Gd155f>> ;
+GdNat := UTL: GdNat :: CREA 'Gd156frIso' 1 = <<Gd156f>> ;
+GdNat := UTL: GdNat :: CREA 'Gd157frIso' 1 = <<Gd157f>> ;
+GdNat := UTL: GdNat :: CREA 'Gd158frIso' 1 = <<Gd158f>> ;
+GdNat := UTL: GdNat :: CREA 'Gd160frIso' 1 = <<Gd160f>> ;
+
+*------------------------------------------------------
+*  Gd2O3 : Gadolinium Oxyde
+*  Isotopic Composition : GdNat 2
+*                         O16   3
+*------------------------------------------------------
+REAL GadOxGdf := 2. ;
+REAL GadOxOf  := 3. ;
+
+REAL GadOxG4f GadOxG5f GadOxG6f GadOxG7f GadOxG8f GadOxG0f ;
+
+EVALUATE GadOxG4f := GadOxGdf Gd154f * ;
+EVALUATE GadOxG5f := GadOxGdf Gd155f * ;
+EVALUATE GadOxG6f := GadOxGdf Gd156f * ;
+EVALUATE GadOxG7f := GadOxGdf Gd157f * ;
+EVALUATE GadOxG8f := GadOxGdf Gd158f * ;
+EVALUATE GadOxG0f := GadOxGdf Gd160f * ;
+
+REAL     MGd2O3 ; ! Average Weight Ratio
+EVALUATE MGd2O3 := MGdNat GadOxGdf * o16 GadOxOf * + ;
+
+*---- Gd2O3 Table Creation ----
+Gd2O3 := UTL:       :: CREA 'AvgMolM     ' 1 = <<MGd2O3>>   ;
+Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG4frIso' 1 = <<GadOxG4f>> ;
+Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG5frIso' 1 = <<GadOxG5f>> ;
+Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG6frIso' 1 = <<GadOxG6f>> ;
+Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG7frIso' 1 = <<GadOxG7f>> ;
+Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG8frIso' 1 = <<GadOxG8f>> ;
+Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG0frIso' 1 = <<GadOxG0f>> ;
+Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxOfrIso ' 1 = <<GadOxOf>>  ;
+
+*------------------------------------------------------
+*  B2O3 : Solid White Glass Diboron Trioxide # B
+*   Quantity of BNat in B2O3 : BNatBq = 2.
+*   Quantity of O16 in B2O3  : O16Bq  = 3.
+*------------------------------------------------------
+REAL BNatBq := 2. ;
+REAL O16Bq :=  3. ;
+*---- B2O3 Average Weight Ratio ----
+REAL     MB2O3 ; ! Average Weight Ratio
+EVALUATE MB2O3 := MBNat BNatBq * o16 O16Bq * + ;
+*------------------------------------------------------
+*  SiO2 : Silica Crystalline  # S
+*   Quantity of SiNat in SiO2 : SiNatSq = 1.
+*   Quantity of O16 in SiO2   : O16Sq   = 2.
+*------------------------------------------------------
+REAL SiNatSq := 1. ;
+REAL O16Sq := 2. ;
+*---- SiO2  Average Weight Ratio ----
+REAL     MSiO2  ; ! Average Weight Ratio
+EVALUATE MSiO2 := MSiNat SiNatSq * o16 O16Sq * + ;
+*------------------------------------------------------
+*  Al2O3 : Aluminium Oxide  # A
+*   Quantity of Aluminium in Al2O3 : Al27Aq = 2.
+*   Quantity of O16 in Al2O3       : O16Aq  = 3.
+*------------------------------------------------------
+REAL Al27Aq := 2. ;
+REAL O16Aq  := 3. ;
+*---- Al2O3  Average Weight Ratio ----
+REAL     MAl2O3  ; ! Average Weight Ratio
+EVALUATE MAl2O3 := al27 Al27Aq * o16 O16Aq * + ;
+*------------------------------------------------------
+*  K2O : Potassium Oxide  # K
+*   Quantity of KNat in K2O  : KNatKq = 2.
+*   Quantity of O16 in K2O : O16Kq  = 1.
+*------------------------------------------------------
+REAL KNatKq := 2. ;
+REAL O16Kq  := 1. ;
+*---- K2O Average Weight Ratio ----
+REAL     MK2O ; ! Average Weight Ratio
+EVALUATE MK2O := MKNat KNatKq * o16 O16Kq * + ;
+*------------------------------------------------------
+*  Na2O : Sodium oxide  # N
+*   Quantity of Na23 in Na2O  : Na23Nq = 2.
+*   Quantity of O16 in Na2O : O16Nq  = 1.
+*------------------------------------------------------
+REAL Na23Nq := 2. ;
+REAL O16Nq  := 1. ;
+*---- Na2O Average Weight Ratio ----
+REAL     MNa2O ; ! Average Weight Ratio
+EVALUATE MNa2O := na23 Na23Nq * o16 O16Nq * + ;
+*------------------------------------------------------
+
+*------------------------------------------------------
+*  BP : Burnable Poison Borosilicate Glass
+*  Tihange Specifications CEA-N-2092
+*  Massic Composition : B2O3  .181
+*                       SiO2  .648
+*                       Al2O3 .076
+*                       K2O   .032
+*                       Na2O  .031
+* The sum of massic composition %m is equal to 0.968 !!
+*------------------------------------------------------
+*               Massic Fraction : %m
+*------------------------------------------------------
+REAL B2O3BPf  ;
+REAL SiO2BPf  ;
+REAL Al2O3BPf ;
+REAL K2OBPf   ;
+REAL Na2OBPf  ;
+
+IF Gestion "BEAV" = THEN
+*  Weight Fraction    : Onat  0.5481
+*                       AlNat 0.0344
+*                       B10   0.0071
+*                       B11   0.0317
+*                       Sinat 0.3787
+  EVALUATE B2O3BPf   := .125 ;
+  EVALUATE SiO2BPf   := .810 ; ! .806
+  EVALUATE Al2O3BPf  := .065 ; ! .024
+  EVALUATE K2OBPf    := .0 ; ! .022
+  EVALUATE Na2OBPf   := .0 ; ! .021
+ELSE
+  EVALUATE B2O3BPf   := .181 ;
+  EVALUATE SiO2BPf   := .648 ;
+  EVALUATE Al2O3BPf  := .076 ;
+  EVALUATE K2OBPf    := .032 ;
+  EVALUATE Na2OBPf   := .031 ;
+ENDIF ;
+EVALUATE SumFr     := B2O3BPf SiO2BPf Al2O3BPf K2OBPf Na2OBPf + + + + ;
+
+*
+*-----------------------------------------------------
+*     Material Massic To Material Fraction
+*
+* %m_i : Massic   fraction of material i
+* %n_i : Material fraction of material i
+*
+*    %n_i = [%m_i / M_i] / [ Sum ( %m_k / M_k ) ]
+*-----------------------------------------------------
+*
+EVALUATE B2O3BPf := B2O3BPf MB2O3 / ;
+EVALUATE SiO2BPf := SiO2BPf MSiO2 / ;
+EVALUATE Al2O3BPf := Al2O3BPf MAl2O3  / ;
+EVALUATE K2OBPf := K2OBPf MK2O / ;
+EVALUATE Na2OBPf := Na2OBPf MNa2O / ;
+*
+EVALUATE ReNorm := B2O3BPf SiO2BPf Al2O3BPf K2OBPf Na2OBPf + + + + ;
+*
+EVALUATE B2O3BPf := B2O3BPf ReNorm / ;
+EVALUATE SiO2BPf := SiO2BPf ReNorm / ;
+EVALUATE Al2O3BPf := Al2O3BPf ReNorm / ;
+EVALUATE K2OBPf := K2OBPf ReNorm / ;
+EVALUATE Na2OBPf := Na2OBPf ReNorm / ;
+
+*
+*-----------------------------------------------------
+*         Molar Mass of BP Calculation : MBP
+*-----------------------------------------------------
+REAL     MBP := 0.0 ; ! BP Average Weight Ratio
+EVALUATE MBP := B2O3BPf  MB2O3 *   SiO2BPf MSiO2 * + ;
+EVALUATE MBP := MBP Al2O3BPf MAl2O3  * + K2OBPf  MK2O  * + ;
+EVALUATE MBP := MBP Na2OBPf  MNa2O  * +  ;
+*
+*----------------------------------------------------
+*  AMF : Ajusted Material Fraction
+*   because  Sum of %m is 0.968 and not 1
+*    --> AMF =  Material Fraction (%n_i) * SumFr
+*----------------------------------------------------
+EVALUATE B2O3BPf := B2O3BPf SumFr * ;
+EVALUATE SiO2BPf := SiO2BPf SumFr * ;
+EVALUATE Al2O3BPf := Al2O3BPf SumFr * ;
+EVALUATE K2OBPf := K2OBPf SumFr * ;
+EVALUATE Na2OBPf := Na2OBPf SumFr * ;
+*-----------------------------------------------------
+*             Isotopic Ratio Calculation
+*-----------------------------------------------------
+* %f = AMF x Qem x Fie
+*  Qem = Quantity of Element in Material
+*  Fie = Fraction of Isotope In Element
+*
+*  example :  B10 isotopic fraction
+*   We find 0,20% of B10 in element BNat -> Fie = 0.2
+*   There are 2 BNat in B203             -> Qem = 2.
+*   The molar fraction of B2O3 is  %n    -> B2O3BPf
+*
+*   -->  B10BPf = 0.2 * 2. * B2O3BPf
+*-----------------------------------------------------
+*------ Isotopic Ratios for all isotopes in BP -------
+REAL B10BPf ;
+REAL B11BPf ;
+REAL O16BPf ;
+REAL SiNatBPf ;
+REAL Al27BPf ;
+REAL KNatBPf ;
+REAL Na23BPf ;
+
+*
+EVALUATE B10BPf := B10f BNatBq B2O3BPf   * *   ;
+EVALUATE B11BPf := B11f BNatBq B2O3BPf   * *   ;
+EVALUATE O16BPf :=        O16Bq B2O3BPf * O16Sq SiO2BPf * + ;
+EVALUATE O16BPf := O16BPf O16Aq Al2O3BPf * + ;
+EVALUATE O16BPf := O16BPf O16Kq K2OBPf * + ;
+EVALUATE O16BPf := O16BPf O16Nq Na2OBPf * + ;
+EVALUATE Na23BPf := Na23Nq Na2OBPf   *   ;
+EVALUATE Al27BPf := Al27Aq Al2O3BPf   *   ;
+EVALUATE KNatBPf := KNatKq K2OBPf   *   ;
+REAL K39BPf := KNatKq K39f * K2OBPf * ;
+REAL K40BPf := KNatKq K40f * K2OBPf * ;
+REAL K41BPf := KNatKq K41f * K2OBPf * ;
+EVALUATE SiNatBPf := SiNatSq SiO2BPf   *   ;
+REAL Si28BPf := SiNatSq Si28f * SiO2BPf  * ;
+REAL Si29BPf := SiNatSq Si29f * SiO2BPf  * ;
+REAL Si30BPf := SiNatSq Si30f * SiO2BPf  * ;
+*
+*---- BP Table Creation ----
+BP := UTL:     :: CREA 'AvgMolM     ' 1 = <<MBP>>      ;
+BP := UTL: BP  :: CREA 'B10BPfrIso  ' 1 = <<B10BPf>>   ;
+BP := UTL: BP  :: CREA 'B11BPfrIso  ' 1 = <<B11BPf>>   ;
+BP := UTL: BP  :: CREA 'O16BPfrIso  ' 1 = <<O16BPf>>   ;
+BP := UTL: BP  :: CREA 'SiNatBPfrIso' 1 = <<SiNatBPf>>   ;
+BP := UTL: BP  :: CREA 'Si28BPfrIso' 1 = <<Si28BPf>>   ;
+BP := UTL: BP  :: CREA 'Si29BPfrIso' 1 = <<Si29BPf>>   ;
+BP := UTL: BP  :: CREA 'Si30BPfrIso' 1 = <<Si30BPf>>   ;
+BP := UTL: BP  :: CREA 'Al27BPfrIso ' 1 = <<Al27BPf>>   ;
+BP := UTL: BP  :: CREA 'KNatBPfrIso ' 1 = <<KNatBPf>>   ;
+BP := UTL: BP  :: CREA 'K39BPfrIso ' 1 = <<K39BPf>>   ;
+BP := UTL: BP  :: CREA 'K40BPfrIso ' 1 = <<K40BPf>>   ;
+BP := UTL: BP  :: CREA 'K41BPfrIso ' 1 = <<K41BPf>>   ;
+BP := UTL: BP  :: CREA 'Na23BPfrIso ' 1 = <<Na23BPf>>   ;
+
+END: ;
+QUIT "LIST" .
+