#5: Corrections and improvements to the Baff-Refl reflector model in Donjon

2 files changed, 9 insertions, 7 deletions

diff --git a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m
index 48db04e..9b23d60 100755
--- a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m
+++ b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m
@@ -36,6 +36,7 @@ REAL avo_conc := avogadro 1.E-24 * ;
 REAL dSS   :=  7.792815 ; ! Densite de l'acier inoxydable (SS304)
 REAL dZr4  :=  6.49012  ; ! Densite du zircaloy 4
 REAL dInc  :=  8.303    ; ! Densite inconel (a 20 degres C...)
+INTEGER nbZones         ; ! Nombre de zones dans le reflecteur
 INTEGER iMil            ; ! Numero de la zone reflecteur
 REAL TMil               ; ! Temperature du moderateur en degres C (suppose = T de tout le milieu)
 REAL dZone              ; ! Densite de la zone consideree
@@ -49,7 +50,7 @@ INTEGER LibId ;
 *--------------------------------------------------------------------
 * Fractions volumiques des differents milieux presents dans la zone
 *--------------------------------------------------------------------
-:: >>dMod<< >>iMil<< >>TMil<< >>CB<< >>fvMod<< >>fvZr4<< >>fvInc<< >>fvSS304<<
+:: >>dMod<< >>nbZones<< >>iMil<< >>TMil<< >>CB<< >>fvMod<< >>fvZr4<< >>fvInc<< >>fvSS304<<
    >>NuclData<< >>LibType<< ;
 
 *------------------------------------------------------
@@ -423,9 +424,10 @@ ECHO "Nb93   " Nb93 ;
 * Definition du melange correspondant dans la librairie
 *--------------------------------------------------------------------
 ECHO "CreaDilut: call LIB" ;
+INTEGER nmix := nbZones 1 + ;
 IF iMil 1 = THEN
   LibDilut := LIB: ::
-    MXIS 300 NMIX 8
+    MXIS 300 NMIX <<nmix>>
     MIXS LIB: <<LibGenre>> FIL: <<NuclData>>
 *   Milieu generique issu de la dilution de tous les constituants
     MIX 1 <<TMil>>
@@ -465,7 +467,7 @@ IF iMil 1 = THEN
   ;
 ELSE
   LibDilut := LIB: LibDilut ::
-    MXIS 300 NMIX 8
+    MXIS 300 NMIX <<nmix>>
     MIXS LIB: <<LibGenre>> FIL: <<NuclData>>
 *   Milieu generique issu de la dilution de tous les constituants
     MIX <<iMil>> <<TMil>>
diff --git a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m
index 6a2f2cc..8d55fb4 100755
--- a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m
+++ b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m
@@ -7,7 +7,7 @@
 * CALL      :                                                                  *
 *  ApexRefl := REFL_APEX Geom FVol Param :: <<nameFeed>> <<nameRefl>>          *
 *                            <<Palier>> <<htype>> <<LibType>> <<NuclData>>     *
-*                            <<OptSph>> <<OptNoal>> <<OptSS>> ;                *
+*                            <<OptSph>> <<OptNoal>> <<OptSS>> <<nbZones>> ;    *
 *                                                                              *
 *  Input objects:                                                              *
 *    Geom    : macro-geometry LCM object                                       *
@@ -25,6 +25,7 @@
 *    OptSph  : 'SPH'/'NOSP' transform DF into ADF/keep DF                      *
 *    OptNoal : 'ALBE'/'NOAL' compute albedo/force right albedo to 0            *
 *    OptSS   : 'TONE'/'NOSS' do/don't activate self-shielding of Fe56 and Cr52 *
+*    nbZones : number of zones in the reflector                                *
 *                                                                              *
 *  Output objects:                                                             *
 *    ApexRefl : reflector APEX file                                            *
@@ -38,6 +39,7 @@ STRING nameFeed nameRefl Palier htype ;
 INTEGER nlf ; :: >>nlf<< ;
 STRING LibType NuclData OptSph OptNoal OptSS ;
   ::  >>LibType<< >>NuclData<< >>OptSph<< >>OptNoal<< >>OptSS<< ;
+INTEGER nbZones ; :: >>nbZones<< ;
 PROCEDURE GetMolarMass MetalMix CreaDilut assertS ;
 LINKED_LIST
   DONNEES Geom_refl LibDilut MACRO LibFuel MolarMasses FeNat CrNat SnNat
@@ -166,8 +168,6 @@ TRACK_refl := SNT: Geom_refl ::
   MAXR 100000
   SN 16 SCAT <<iscat>> ;
 
-INTEGER nbZones := 7 ;
-
 REAL dModBas := 0.753 ; ! Densite du moderateur en entree du coeur
 REAL TMilBas := 286.4 ; ! Temp du reflecteur bas (soit Tmod,entree)
 EVALUATE TMilBas := TMilBas 273.15 + ; ! Conversion en Kelvin
@@ -251,7 +251,7 @@ WHILE IndexCB iBoreMax < DO
 
     ! Define MIX 1 in reflector macrolib
     LibDilut := CreaDilut LibDilut MolarMasses BNat Zirc4 Incon SS304 ::
-    <<dMod>> <<premierMIX>> <<TMil>> <<CB>> <<fvMod>> <<fvZr4>> <<fvInc>>
+    <<dMod>> <<nbZones>> <<premierMIX>> <<TMil>> <<CB>> <<fvMod>> <<fvZr4>> <<fvInc>>
     <<fvSS304>> <<NuclData>> <<LibType>> ;
 
     EVALUATE iZone := iZone 1 + ;