Initial commit from Polytechnique Montreal

1 files changed, 723 insertions, 0 deletions

diff --git a/doc/IGE351/SectDexcell.tex b/doc/IGE351/SectDexcell.tex
new file mode 100644
index 0000000..c3e28a4
--- /dev/null
+++ b/doc/IGE351/SectDexcell.tex
@@ -0,0 +1,723 @@
+\subsection{The \moc{excelt} dependent records on a
+\dir{tracking} directory}\label{sect:excelltrackingdir}
+
+When the \moc{EXCELT:} modules is used ($\mathsf{CDOOR}$={\tt 'EXCELL'}), the following elements in the vector
+$\mathcal{S}^{t}_{i}$ will also be defined.
+
+\begin{itemize}
+\item $\mathcal{S}^{t}_{6}$ is the number of Legendre orders used for the flux expansions, where
+\begin{displaymath}
+\mathcal{S}^{t}_{6} = \left\{
+\begin{array}{rl}
+ 1 & \textrm{isotropic} \\
+ 2 & \textrm{linearly anisotropic.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{7}$ is the specific EXCELL tracking procedure considered where
+\begin{displaymath}
+\mathcal{S}^{t}_{7} = \left\{
+\begin{array}{rl}
+ 1 & \textrm{Cartesian 2D/3D assembly using {\tt EXCELT:}} \\
+ 2 & \textrm{hexagonal 2D/3D assembly using {\tt EXCELT:}} \\
+ 3 & \textrm{2D cluster geometry using {\tt EXCELT:}} \\
+ 4 & \textrm{2D or 3D geometry using {\tt NXT:} or {\tt SALT:}} \\
+ 5 & \textrm{2D multicell surfacic approximation.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{8}$ is the track normalization flag where
+\begin{displaymath}
+\mathcal{S}^{t}_{8} = \left\{
+\begin{array}{rl}
+ -1 & \textrm{direction dependent track normalization to exact volumes;} \\
+ 0 & \textrm{global track normalization to exact volumes;} \\
+ 1 & \textrm{no normalization.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{9}$ is the tracking type where
+\begin{displaymath}
+\mathcal{S}^{t}_{9} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{means that a standard tracking procedure was considered (\moc{TISO});} \\
+ 1 & \textrm{means that a cyclic tracking procedure was considered  (\moc{TSPC}).} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{10}$ is the type of boundary conditions where
+\begin{displaymath}
+\mathcal{S}^{t}_{10} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{means that isotropic (white) boundary conditions will be considered (\moc{PISO});} \\
+ 1 & \textrm{means that mirror-like (specular) boundary conditions will be considered (\moc{PSPC}).} \\
+\end{array} \right.
+\end{displaymath}
+
+Note that mirror-like boundary conditions ($\mathcal{S}^{t}_{10} = 1$) can be used only if
+a cyclic tracking procedure was considered ($\mathcal{S}^{t}_{9} = 1$). 
+
+\item $\mathcal{S}^{t}_{11}=N_\Omega$ is the order of the azimuthal (2-D) or solid (3-D) angular quadrature. For 2-D geometry, the order 
+of the azimuthal quadrature represents:
+\begin{itemize}
+\item the number of equal sectors (trapezoidal quadrature) in the $[0,\pi]$ range when the {\tt EXCELT:} module is
+used for Cartesian assemblies; 
+\item the number of equal sectors (trapezoidal quadrature) in the $[0,2\pi/3]$ range when the {\tt EXCELT:} module 
+is used for hexagonal geometries;
+\item the number of equal sectors (trapezoidal quadrature) in the $[0,\max (\mathcal{S}^{t}_{12},2) \, \pi]$ range when the {\tt EXCELT:}
+module is used for cluster geometries; 
+\item the number of trapezoidal sectors in the $[0,\pi/2]$ range when the {\tt NXT:} module is used. 
+\end{itemize}
+
+For 3-D geometry, the order of the solid angle quadrature is: 
+\begin{itemize}
+\item the order $n$ of the $EQ_n$ quadrature in a quadrant ($0\le \varphi\le \pi/2$ and $0\le \theta \le \pi/2$) when the
+{\tt EXCELT:} module is used for Cartesian assemblies for $N_{\rm dir}=n(n+2)/8$ direction in each 
+quadrant; 
+\item the number of equal sectors (trapezoidal quadrature) in the $[0,2\pi/3]$ range when the {\tt EXCELT:} module 
+is used for hexagonal geometries; 
+\item not used for the {\tt EXCELT:} module in cluster geometries; 
+\item the order $n$ of the $EQ_n$ , $LC_n$ or $LT_n$ quadrature in a quadrant with $0\le \varphi\le \pi/2$ and $0\le \theta \le \pi/2$.
+When the {\tt NXT:} modules is used, the number $N_{\rm dir}$ of directions for the azimuthal (2-D) or solid (3-D) angle quadrature
+results in $N_{\rm dir}=n(n+2)/8$ directions for the $EQ_n$ quadrature, 
+$N_{\rm dir}=3n(n+2)/8$ for the $LC_n$ quadrature and $N_{\rm dir}=3n^2/2$ for the $LT_n$ quadrature.
+\end{itemize}
+
+\item $\mathcal{S}^{t}_{12}$ is the angular symmetry factor;
+
+\item $\mathcal{S}^{t}_{13}$ is the polar quadrature type:
+\begin{displaymath}
+\mathcal{S}^{t}_{13} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{Gauss-Legendre} \\
+ 1 & \textrm{CACTUS type 1} \\
+ 2 & \textrm{CACTUS type 2} \\
+ 3 & \textrm{McDaniel} \\
+ 4 & \textrm{McDaniel with $P_1$ constraint} \\
+ 5 & \textrm{Gauss optimized.}
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{14}$ is the order of the polar quadrature.
+
+\item $\mathcal{S}^{t}_{15}$ is the azimuthal (2-D) or solid (3-D) angle quadrature type where
+\begin{displaymath}
+\mathcal{S}^{t}_{15} = \left\{
+\begin{array}{rl}
+ 1 & \textrm{for a $EQ_n$ (3-D) or trapezoidal (2-D) quadrature;} \\
+ 2 & \textrm{for a Gauss quadrature (2-D hexagonal geometries);} \\
+ 3 & \textrm{for a median angle quadrature;} \\
+ 4 & \textrm{for a $LC_n$ 3-D quadrature;} \\
+ 5 & \textrm{for a $LT_n$ 3-D quadrature;} \\
+ 6 & \textrm{for a $\mu_1$--optimized level-symmetric 3-D quadrature;} \\
+ 7 & \textrm{for a quadrupole range (QR) 3-D quadrature.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{16}$ is the number of geometric dimensions (1, 2 or 3).
+
+\item $\mathcal{S}^{t}_{17}$ is the number of tracking points on a line.
+
+\item $\mathcal{S}^{t}_{18}$ is the maximum length of a track.
+
+\item $\mathcal{S}^{t}_{19}$ is the total number of tracks generated.
+
+\item $\mathcal{S}^{t}_{20}$ is the total number of track directions processed.
+
+\item $\mathcal{S}^{t}_{21}$ is the line format option for \moc{TLM:} module, where
+\begin{displaymath}
+\mathcal{S}^{t}_{21} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{short format} \\
+ 1 & \textrm{complete format.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{22}$ is the vectorization option for computing collision probability matrices where
+\begin{displaymath}
+\mathcal{S}^{t}_{22} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{scalar algorithm. The tracking is readed for each energy group.} \\
+ 1 & \textrm{vectorial algorithm. The tracking is readed once and the collision
+ probability} \\
+   & \textrm{matrices are computed in many energy groups.} \\
+ 2 & \textrm{vectorial algorithm of type \moc{EXCELL:}. The tracking is computed in \moc{DOORVP}} \\
+   & \textrm{and the collision probability matrices are computed in many energy groups.} \\
+ 3 & \textrm{vectorial algorithm of type \moc{NXT:}. The tracking is computed in \moc{DOORVP}} \\
+   & \textrm{and the collision probability matrices are computed in many energy groups.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{23}$ is the tracking flag, where
+\begin{displaymath}
+\mathcal{S}^{t}_{23} = \left\{
+\begin{array}{rl}
+ -1 & \textrm{the LCM object information is used for Monte-Carlo calculations. No} \\
+   & \textrm{tracking file produced.} \\
+ 0 & \textrm{this option de-activates the tracking file production if \moc{TISO} or \moc{TSPC} are} \\
+   & \textrm{specified.} \\
+ 1 & \textrm{this option produces the tracking file.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{24}$ is the multicell surfacic flag, where
+\begin{displaymath}
+\mathcal{S}^{t}_{24} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{do not use the multicell surfacic approximation.} \\
+ {\tt NMACRO} & \textrm{number of macro geometries for the multicell surfacic approximation.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{t}_{25}$ is the number of normal planes considered in NXT tracking.
+
+\item $\mathcal{S}^{t}_{26}$ is the {\tt MERGMIX} activation flag (0/1: off/on).
+
+\item $\mathcal{S}^{t}_{27}$ is the number of tracks assigned to each OpenMP core.
+
+\item $\mathcal{S}^{t}_{28}$ is the number of interface currents for the multicell surfacic approximation. These interface currents are additional unknowns
+components if a flux-current iteration solution is used.
+
+\item $\mathcal{S}^{t}_{29}$ is the sum $\sum_{i=1}^{\mathcal{S}^{t}_{24}} N_i$ where $N_i$ is the number of perimeter elements after unfolding in macro
+geometry $i$.
+
+\item $\mathcal{S}^{t}_{39}$ is the prismatic tracking activation flag along $z$ axis for 3D geometry (0/3: off/on-$z$ axis).
+
+\end{itemize}
+
+The following records will also be present on the main level of a \dir{tracking} directory.
+
+\begin{DescriptionEnregistrement}{The \moc{EXCELT:} records in \dir{tracking}}{8.0cm}
+\label{tabl:main-structure}
+\RealEnr
+  {EXCELTRACKOP}{$40$}{}
+  {array ${\cal R}_i$ containing additional EXCELL or NXT tracking parameters.} 
+\IntEnr
+  {ICODE\blank{7}}{$6$}
+  {array ${\cal I}_{\beta,k}$ containing the surface albedo index (geometric surface albedo
+  $\beta_{{\rm g},k}$ are used if ${\cal I}_{\beta,k}< 0$ while physical 
+  surface albedo $\beta_{{\rm p},k}$ are used if ${\cal I}_{\beta,k}> 0$).} 
+\RealEnr
+  {ALBEDO\blank{6}}{$6$}{}
+  {array $\beta_{{\rm g},k}$ containing the geometric surface albedo (used only if ${\cal I}_{\beta,k}\ge 0$).}
+\OptDirEnr
+  {EXCELL\blank{6}}{$\mathcal{S}^{t}_{7}< 4$}
+  {directory containing additional {\tt EXCELT:}  records for the cases where $\mathcal{S}^{t}_{7}=1$ or $\mathcal{S}^{t}_{7}=3$.}
+\OptDirEnr
+  {NXTRecords\blank{2}}{$\mathcal{S}^{t}_{7}=4$}
+  {directory containing additional {\tt NXT:} records.}
+\OptDirEnr
+  {PROJECTION\blank{2}}{$\mathcal{S}^{t}_{39}>0$}
+  {directory containing the analysis of the projection of a 3D prismatic geometry.}
+\end{DescriptionEnregistrement}
+
+The record ${\cal R}_i$ contains the following information:
+\begin{itemize}
+\item ${\cal R}_1$ is the maximum error allowed on the exponential function. 
+\item ${\cal R}_2$ is the user requested tracking density in cm$^{-1}$ and in cm$^{-2}$ respectively for 2D and 3D calculations. 
+\item ${\cal R}_3$ is the maximum distance in cm between an integration line and a surface. 
+\item ${\cal R}_4$ is the computed tracking density in cm$^{-1}$ and in cm$^{-2}$ respectively for 2D and 3D calculations (used 
+only if $\mathcal{S}^{t}_{7}=4$).
+\item ${\cal R}_5$ is the computed line spacing in cm (used only if $\mathcal{S}^{t}_{7}=4$). 
+\item ${\cal R}_6$ is the weight of the spatial quadrature (used only if $\mathcal{S}^{t}_{7}=4$). 
+\item ${\cal R}_7$ is the minimal radius of the circle (2-D) or sphere (3-D) containing the geometry (used only if $\mathcal{S}^{t}_{7}=4$). 
+\item ${\cal R}_8$ is the $x$ position of the center of the minimal circle (2-D) or sphere (3-D) containing the geometry (used 
+only if $\mathcal{S}^{t}_{7}=4$). 
+\item ${\cal R}_9$ is the $y$ position of the center of the minimal circle (2-D) or sphere (3-D) containing the geometry (used 
+only if $\mathcal{S}^{t}_{7}=4$). 
+\item ${\cal R}_{10}$ is the $z$ position of the center of the minimal circle (2-D) or sphere (3-D) containing the geometry (used 
+only if $\mathcal{S}^{t}_{7}=4$). 
+\item ${\cal R}_{11}$ is the spatial cutoff factor for tracking.
+\item ${\cal R}_{12}$ is the stopping criterion for flux-current iterations of the multicell surfacic approximation.
+\item ${\cal R}_{40}$ user requested tracking density in cm$^{-1}$ for inline contruction of 3D tracks when a prismatic tracking is considered (only used if $\mathcal{S}^{t}_{39}>0$).
+\end{itemize}
+
+Additional records are defined in the structure of Table~\ref{tabl:main-structure} in case where the multicell surfacic
+approximation is used (if $\mathcal{S}^{t}_{7}=5$):
+
+\begin{DescriptionEnregistrement}{The additional multicell surfacic records in \dir{tracking}}{8.0cm}
+\DirlEnr
+  {MACRO-TRACK\blank{1}}{$\mathcal{S}^{t}_{24}$}
+  {List of directories. Each component of this list contains the the tracking information for a macro geometry.
+  Each directory of the list follows the specification presented in Table~\ref{tabl:main-structure} for a case with
+  a unique macro geometry ($\mathcal{S}^{t}_{7}=4$ and $\mathcal{S}^{t}_{24}=1$).}
+\IntEnr
+  {NMC\_NODE\blank{4}}{$\mathcal{S}^{t}_{24}+1$}
+  {Offset of the first region in each macro geometry}
+\IntEnr
+  {NMC\_SURF\blank{4}}{$\mathcal{S}^{t}_{24}+1$}
+  {Offset of the first boundary surface in each macro geometry}
+\IntEnr
+  {IFR\blank{9}}{$\mathcal{S}^{t}_{29}$}
+  {Index numbers of incoming currents}
+\RealEnr
+  {ALB\blank{9}}{$\mathcal{S}^{t}_{29}$}{}
+  {Albedo or transmission factors corresponding to incoming currents}
+\IntEnr
+  {MIX\blank{9}}{$\mathcal{S}^{t}_{29}$}
+  {Index numbers of outgoing currents}
+\RealEnr
+  {DVX\blank{9}}{$\mathcal{S}^{t}_{29}$}{}
+  {Weights corresponding to outgoing currents}
+\RealEnr
+  {SUR\blank{9}}{$\mathcal{S}^{t}_{29}$}{cm}
+  {Interface surfaces corresponding to incoming currents}
+\IntEnr
+  {MERGE\_MACRO\blank{1}}{$\mathcal{S}^{t}_{1}$}
+  {Merge index assigned to each node in the global surfacic geometry}
+\end{DescriptionEnregistrement}
+
+The \namedir{NXTRecords} directory contains the information required to track the geometry using the \moc{NXT:} module module once it has been analyzed. The contents of this directory is presented in Table~\ref{tabl:NXTRecords}.
+
+\begin{DescriptionEnregistrement}{Global geometry records in \namedir{NXTRecords}}{8.0cm}\label{tabl:NXTRecords}
+\IntEnr
+  {G00000001DIM}{$40$}
+  {array $N^{\text{GG}}_{i}$ containing the dimensioning information required to rebuilt the assembly} 
+\IntEnr
+  {G00000001CUF}{$2,N^{\text{GG}}_{5}$}
+  {array $D^{\text{GG}}_{i,j}$ containing the assembly description of the geometry in terms of cells and rotations. The first element ($i=1$) identifies the cell number while the second element identifies the cell rotation} 
+\IntEnr
+  {G00000001CIS}{$4,N^{\text{GG}}_{4}$}
+  {array $S^{\text{GG}}_{i,j}$ containing the cell intrinsic symmetry properties. A value of $1$ indicates that a center cell reflexion symmetry is present while a value of $0$ indicates that the symmetry is not considered (see below for a more complete description of this array)} 
+\IntEnr
+  {G00000001CFE}{$0:8,N^{\text{GG}}_{4}$}
+  {array $F^{\text{GG}}_{i,j}$ containing the assembly external surface identification index (see below for a more complete description of this array)}
+\DbleEnr
+  {G00000001SMX}{$0:N^{\text{GG}}_{13}$}{cm}
+  {array $x^{\text{GG}}$ containing the $x$-directed mesh for the cell assembly in a Cartesian or Cylindrical geometry and the $x$ position of the cell center for an hexagonal assembly (see below for more explanations)}
+\DbleEnr
+  {G00000001SMY}{$0:N^{\text{GG}}_{14}$}{cm}
+  {array $y^{\text{GG}}$ containing the $y$-directed mesh for the cell assembly in a Cartesian or Cylindrical geometry and the $y$ position of the cell center for an hexagonal assembly (see below for more explanations)}
+\OptDbleEnr
+  {G00000001SMZ}{$0:N^{\text{GG}}_{15}$}{$N^{\text{GG}}_{1}=3$}{cm}
+  {array $z^{\text{GG}}$ containing the $z$-directed mesh for the cell assembly (see below for more explanations)}
+\OptDbleEnr
+  {G00000001SMR}{$0:1$}{$N^{\text{GG}}_{2}=1$}{cm}
+  {the radius $r^{\text{GG}}$ of the outer assembly boundary (see below for more explanations)}
+\IntEnr
+  {KEYMRG\blank{6}}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}
+  {array $\text{MRG}_{i}$ containing the merged surface and region number associated with each individual surfaces and regions in this geometry}
+\IntEnr
+  {MATALB\blank{6}}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}
+  {array containing the albedo number associated with each surface and the physical mixture number associated with each region in this geometry}
+\IntEnr
+  {HOMMATALB\blank{3}}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}
+  {array containing the albedo number associated with each surface and the virtual (homogenization) mixture number associated with each region in this geometry}
+\DbleEnr
+  {SAreaRvolume}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}{}
+  {array containing the area ($S_{\alpha}$ in cm for 2-D and cm$^{2}$ for 3-D problems) and volume ($V_{i}$ cm$^{2}$ for 2-D and
+cm$^{3}$ for 3-D problems) associated with each surface and region in this geometry}
+\end{DescriptionEnregistrement}
+
+The dimensioning vector for the global geometry contains the following information:
+\begin{itemize}
+\item $N^{\text{GG}}_{1}$ number of dimensions for the problem.
+\item $N^{\text{GG}}_{2}$ type of boundary, defined as follows:
+  \begin{itemize}
+  \item $i=0$: Cartesian geometry;
+  \item $i=1$: cylindrical geometry;
+  \item $i=2$: isocel geometry with specular reflection;
+  \item $i=3$: hexagonal geometry with translation;
+  \item $i=4$: isocel geometry with RA60 rotation and translation;
+  \item $i=5$: lozenge geometry with R120 rotation and translation.
+  \end{itemize}
+\item $N^{\text{GG}}_{3}$ first direction to process in the analysis. For cylinder, this is the direction of the first axis of the plane normal to the cylinder axis. For Cartesian and hexagonal geometries a value of 1 ($x$-axis) is selected by default. 
+\item $N^{\text{GG}}_{4}$ number of cells in the original geometry (before unfolding).
+\item $N^{\text{GG}}_{5}$ number of cells in the geometry after the original geometry is unfolded according to the symmetries.
+\item $N^{\text{GG}}_{6}$ diagonal symmetry flag. A value of $0$ indicates that this symmetry is not used. A value of $-1$ indicates that the symmetry is used for the $x_{-}=y_{+}$ plane and a value of $1$ that the symmetry is used for the $x_{+}=y_{}$ plane.
+\item $N^{\text{GG}}_{7}$ flag to identify symmetries with respect to the $x$-axis ($x_{-}$ or $x_{+}$). A value of $0$ indicates that no symmetry is present, $N^{\text{GG}}_{7}=\pm 1$ is for a \moc{SYME} symmetry at the $x_{\pm}$ plane, $N^{\text{GG}}_{7}=\pm 2$ represents a \moc{SSYM} symmetry at the $x_{\pm}$ plane and $N^{\text{GG}}_{7}= 3$ implies a translation symmetry is the $x$ direction ($x_{-}=x_{+}$).
+\item $N^{\text{GG}}_{8}$ flag to identify symmetries with respect to the $y$-axis ($y_{-}$ or $y_{+}$). A value of $0$ indicates that no symmetry is present, $N^{\text{GG}}_{7}=\pm 1$ is for a \moc{SYME} symmetry at the $y_{\pm}$ plane, $N^{\text{GG}}_{7}=\pm 2$ represents a \moc{SSYM} symmetry at the $y_{\pm}$ plane and $N^{\text{GG}}_{7}= 3$ implies a translation symmetry is the $y$ direction ($y_{-}=y_{+}$).
+\item $N^{\text{GG}}_{9}$ flag to identify symmetries with respect to the $z$-axis ($z_{-}$ or $z_{+}$). A value of $0$ indicates that no symmetry is present, $N^{\text{GG}}_{7}=\pm 1$ is for a \moc{SYME} symmetry at the $z_{\pm}$ plane, $N^{\text{GG}}_{7}=\pm 2$ represents a \moc{SSYM} symmetry at the $z_{\pm}$ plane and $N^{\text{GG}}_{7}= 3$ implies a translation symmetry is the $z$ direction ($z_{-}=z_{+}$).
+\item $N^{\text{GG}}_{10}$ number of $x$ mesh subdivisions or hexagons in the original geometry.
+\item $N^{\text{GG}}_{11}$ number of $y$ mesh subdivisions or hexagons in the original geometry.
+\item $N^{\text{GG}}_{12}$ number of $z$ mesh subdivisions in the original geometry.
+\item $N^{\text{GG}}_{13}$ number of $x$ mesh subdivisions or hexagons in the unfolded geometry.
+\item $N^{\text{GG}}_{14}$ number of $y$ mesh subdivisions or hexagons in the unfolded geometry.
+\item $N^{\text{GG}}_{15}$ number of $z$ mesh subdivisions in the unfolded geometry.
+\item $N^{\text{GG}}_{16}$ maximum number cells required to represent this geometry.
+\item $N^{\text{GG}}_{17}$ maximum number of region for this geometry.
+\item $N^{\text{GG}}_{18}$ total number of clusters in this geometry.
+\item $N^{\text{GG}}_{19}$ maximum number of pins in this geometry.
+\item $N^{\text{GG}}_{20}$ maximum dimensions of any mesh array for a cell in this geometry.
+\item $N^{\text{GG}}_{21}$ maximum dimensions of any mesh array for a pin in this geometry. 
+\item $N^{\text{GG}}_{22}$ number of external surfaces for this geometry.
+\item $N^{\text{GG}}_{23}$ number of regions for this geometry.
+\item $N^{\text{GG}}_{24}$ maximum number of external surfaces in a sub-geometry included in this geometry.
+\item $N^{\text{GG}}_{25}$ maximum number of regions in a sub-geometry included in this geometry.
+\item $N^{\text{GG}}_{26}$ MERGE flag defined as follows:
+  \begin{itemize}
+  \item $i=0$: no merge;
+  \item $i=1$: {\tt MERGE MIX} applied to regions.
+  \end{itemize}
+\end{itemize}
+
+The indexing of array $S^{\text{GG}}_{i,j}$ for the axis of symmetry is as follows
+\begin{enumerate}
+ \item Cartesian  assemblies: 
+  \begin{itemize}
+  \item $i=1$ refers to a reflexion of the geometry on a plane normal the $x$-axis; 
+  \item $i=2$ refers to a reflexion of the geometry on a plane normal the $y$-axis; 
+  \item $i=3$ refers to a reflexion of the geometry on the plane $x=y$; 
+  \item $i=4$ refers to a reflexion of the geometry on a plane normal the $z$-axis. 
+  \end{itemize}
+\item Hexagonal assemblies (symmetries not yet programmed).
+  \begin{itemize}
+  \item $i=1$ refers to a reflexion of the geometry on a plane normal the $u$-axis; 
+  \item $i=2$ refers to a reflexion of the geometry on a plane normal the $v$-axis; 
+  \item $i=3$ refers to a reflexion of the geometry on the plane $w$; 
+  \item $i=4$ refers to a reflexion of the geometry on a plane normal the $z$ axis. 
+  \end{itemize}
+\end{enumerate}
+
+\begin{figure}[htbp]  
+\begin{center} 
+\parbox{6cm}{\epsfxsize=6cm \epsffile{HexAssmbR.eps}}\hspace{1.0cm} \parbox{6cm}{\epsfxsize=6cm \epsffile{HexFaces.eps}}
+\caption{Example of an assembly of hexagons (left) and external faces identification for an hexagon}\label{fig:HexAssmbR}   
+\end{center}  
+\end{figure}
+
+The indexing of array $F^{\text{GG}}_{i,j}$ for external surface identification is as follows. First $F^{\text{GG}}_{0,j}$ represents the number of times the cell appears in the geometry after it has been unfolded. For $i>0$, $F^{\text{GG}}_{i,j}$ can take the following values
+$$
+F^{\text{GG}}_{i,j} =\cases{ 1 & \text{surface associated with direction $i$ of cell $j$ is an external boundary of the assembly}\cr
+                                   0 & \text{surface associated with direction $i$ of cell $j$ is not an external boundary of the assembly}\cr
+                      }
+$$
+with the following planes associated with different values of $i$:
+\begin{enumerate}
+ \item Cartesian assemblies: 
+  \begin{itemize}
+  \item $i=1$ surfaces on the $x_{-}$ plane for cell $j$; 
+  \item $i=2$ surfaces on the $x_{+}$ plane for cell $j$;
+  \item $i=3$ surfaces on the $y_{-}$ plane for cell $j$;
+  \item $i=4$ surfaces on the $y_{+}$ plane for cell $j$;
+  \item $i=5$ surfaces on the $z_{-}$ plane for cell $j$; 
+  \item $i=6$ surfaces on the $z_{+}$ plane for cell $j$.
+  \end{itemize}
+%\item Cylindrical assemblies (not yet programmed).
+%  \begin{itemize}
+%  \item $i=5$ identify the surfaces on the $z_{-}$ plane for cell $j$; 
+%  \item $i=6$ identify the surfaces on the $z_{+}$ plane for cell $j$.
+%  \item $i=8$ identify the surfaces on the $r_{-}$ plane for cell $j$; 
+%  \end{itemize}
+\item Hexagonal assemblies (see \Fig{HexAssmbR}):
+  \begin{itemize}
+  \item $i=1$ surfaces on the $u_{-}$ plane for cell $j$; 
+  \item $i=2$ surfaces on the $u_{+}$ plane for cell $j$;
+  \item $i=3$ surfaces on the $v_{-}$ plane for cell $j$;
+  \item $i=4$ surfaces on the $v_{+}$ plane for cell $j$;
+  \item $i=5$ surfaces on the $z_{-}$ plane for cell $j$; 
+  \item $i=6$ surfaces on the $z_{+}$ plane for cell $j$;
+  \item $i=7$ surfaces on the $w_{-}$ plane for cell $j$;
+  \item $i=8$ surfaces on the $w_{+}$ plane for cell $j$.
+  \end{itemize}
+\end{enumerate}
+
+The arrays $x^{\text{GG}}$, $y^{\text{GG}}$, $z^{\text{GG}}$ and $r^{\text{GG}}$ contain the following information:
+\begin{enumerate}
+\item Cartesian assemblies: 
+  \begin{itemize}
+  \item $x^{\text{GG}}_{i-1}$ and $x^{\text{GG}}_{i}$ are the lower and upper $x$ limits of mesh element $i$ ($i=1,n^{x}$);
+  \item $y^{\text{GG}}_{j-1}$ and $y^{\text{GG}}_{j}$ are the lower and upper $y$ limits of mesh element $j$ ($j=1,n^{y}$);
+  \item $z^{\text{GG}}_{k-1}$ and $z^{\text{GG}}_{k}$ are the lower and upper $z$ limits of mesh element $k$ ($k=1,n^{z}$).
+  \end{itemize}
+%\item Cylindrical assemblies: not yet programmed.
+\item Hexagonal assemblies (see \Fig{HexAssmbR}): 
+  \begin{itemize}
+  \item $x^{\text{GG}}_{0}=h$ is the width of one face of the hexagon and $x^{\text{GG}}_{i}$ is the position in $x$ of the center of cell $i$ in the assembly;
+  \item $y^{\text{GG}}_{0}=h$ is the width of one face of the hexagon and $y^{\text{GG}}_{j}$ is the position in $y$ of the center of cell $j$ in the assembly;
+  \item $z^{\text{GG}}_{k-1}$ and $z^{\text{GG}}_{k}$ are the lower and upper $z$ limits of mesh element $k$ ($k=1,n^{z}$).
+  \end{itemize}
+\end{enumerate}
+
+As we noted above, the global geometry is always an assembly containing cells. For each cell $i$ in this assembly, several records will be generated in the \namedir{NXTRecords} directory. These records are identified using a FORTRAN \verb|CHARACTER*12| variable as follows 
+\begin{quote}
+\begin{verbatim}
+INTEGER      I
+CHARACTER*12 NAMREC
+CHARACTER*3  NREC
+WRITE(NAMREC,'(A1,I8.8,A3)') 'C',I,NREC
+\end{verbatim}
+\end{quote}
+where the variable \verb|NREC| can take the following values:
+\begin{itemize}
+\item DIM for dimensioning information;
+\item SMR for the radial mesh description;
+\item SMX for the $x$-directed mesh description;
+\item SMY for the $y$-directed mesh description;
+\item SMZ for the $z$-directed mesh description;
+\item MIX for physical mixture description;
+\item HOM for virtual mixture description;
+\item VSE for areas and volumes results;
+\item VSI for local surfaces and regions identification;
+\item RID for final region numbering;
+\item SID for final surface numbering
+\item PNT for pin contents description;
+\item PIN for pins location.
+\end{itemize}
+
+In Table~\ref{tabl:NXTCell}, a description of the  additional \namedir{NXTRecords} records associated with cell $i=1$ can be found.
+
+\begin{DescriptionEnregistrement}{Cell $i=1$ records in \namedir{NXTRecords}}{8.0cm}\label{tabl:NXTCell}
+\IntEnr
+  {C00000001DIM}{$40$}
+  {array $N^{\text{GC}}_{j}$ containing the dimensioning information required to rebuilt the cell} 
+\DbleEnr
+  {C00000001SMR}{$N^{\text{GC}}_{2}$}{cm}
+  {array $r^{\text{GC}}_{j}$ containing the cell radial mesh description}
+\DbleEnr
+  {C00000001SMX}{$N^{\text{GC}}_{3}$}{cm}
+  {array $x^{\text{GC}}_{j}$ containing the cell $x$-directed mesh description}
+\DbleEnr
+  {C00000001SMY}{$N^{\text{GC}}_{4}$}{cm}
+  {array $y^{\text{GC}}_{j}$ containing the cell $y$-directed mesh description}
+\DbleEnr
+  {C00000001SMZ}{$N^{\text{GC}}_{5}$}{cm}
+  {array $z^{\text{GC}}_{j}$ containing the cell $z$-directed mesh description}
+\IntEnr
+  {C00000001MIX}{$N^{\text{GC}}_{6}$}
+  {array $M^{\text{GC}}_{j}$ containing the cell physical mixture for each region}
+\IntEnr
+  {C00000001HOM}{$N^{\text{GC}}_{6}$}
+  {array $H^{\text{GC}}_{j}$ containing the cell virtual mixture for each region}
+\DbleEnr
+  {C00000001VSE}{$-N^{\text{GC}}_{9}:N^{\text{GC}}_{8}$}{}
+  {array $\text{SV}^{\text{GC}}_{j}$ containing surface area $j$ ($\text{SV}^{\text{GC}}_{-j}=S^{\text{GC}}_{j}$ in cm for 2-D and cm$^{2}$ for 3-D problems) and
+regional volumes $j$ ($\text{SV}^{\text{GC}}_{j}=V^{\text{GC}}_{j}$ in cm$^{2}$ for 2-D and cm$^{3}$ for 3-D problems)}
+\IntEnr
+  {C00000001VSI}{$4,-N^{\text{GC}}_{9}:N^{\text{GC}}_{8}$}
+  {array $\text{VSI}^{\text{GC}}_{k,j}$ containing the location of a surface ($j<0$) and a region ($j>0$)}
+\IntEnr
+  {C00000001RID}{$N^{\text{GC}}_{8}$}
+  {index array $\text{RID}^{\text{GC}}_{j}$ associating local and global region numbering}
+\IntEnr
+  {C00000001SID}{$N^{\text{GC}}_{9}$}
+  {index array $\text{SID}^{\text{GC}}_{j,i}$ associating local and global outer surface numbering}
+\IntEnr
+  {C00000001PNT}{$3,N^{\text{GC}}_{18}$}
+  {array $\text{PC}^{\text{GC}}_{k,j}$ containing the cell pin contents}
+\DbleEnr
+  {C00000001PIN}{$-1:4,N^{\text{GC}}_{18}$}{}
+  {array $\text{p}^{\text{GC}}_{k,j}$ containing the location of the pins in cell}
+\end{DescriptionEnregistrement}
+Note that the record names above are built using the following FORTRAN instructions:
+  \begin{quote}
+    \verb|WRITE(NAMREC,'(A1,I8.8,A3)') 'C',|$i$\verb|,NAMEXT|
+  \end{quote}
+
+The cell dimensioning array $N^{\text{GC}}_{i}$ for cell $i$ contains the following information:
+\begin{itemize}
+\item $N^{\text{GC}}_{1}$ cell geometry type (see the definition of $\mathcal{S}^{G}_{1}$ in \Sect{geometrydirmain});
+\item $N^{\text{GC}}_{2}$ dimensions of the radial mesh array; 
+\item $N^{\text{GC}}_{3}$ dimensions of the $x$-directed mesh array;  
+\item $N^{\text{GC}}_{4}$ dimensions of the $y$-directed mesh array;  
+\item $N^{\text{GC}}_{5}$ dimensions of the $z$-directed mesh array;  
+\item $N^{\text{GC}}_{6}$ dimensions of the mixture record;
+\item $N^{\text{GC}}_{7}$ geometry level (1 for cell);
+\item $N^{\text{GC}}_{8}$ number of regions in the cell before symmetry considerations;
+\item $N^{\text{GC}}_{9}$ number of surfaces in the cell before symmetry considerations;
+\item $N^{\text{GC}}_{10}$ number of regions in the cell after symmetry considerations;
+\item $N^{\text{GC}}_{11}$ number of surfaces in the cell after symmetry considerations;
+\item $N^{\text{GC}}_{12}$ first global region number for cell;
+\item $N^{\text{GC}}_{13}$ last global region number for cell;
+\item $N^{\text{GC}}_{14}$ first global surface number for cell;
+\item $N^{\text{GC}}_{15}$ last global surface number for cell;
+\item $N^{\text{GC}}_{16}$ number of pin cluster geometries in cell;
+\item $N^{\text{GC}}_{17}$ first pin cluster geometry associated with cell;
+\item $N^{\text{GC}}_{18}$ total number of pins in cell;
+\item $N^{\text{GC}}_{19}$ number of times this cell is used in the global cell.
+\end{itemize}
+while the remaining elements are not used. 
+
+The array $x^{\text{GC}}_{j}$ contains the following information:
+\begin{itemize}
+\item $x^{\text{GC}}_{-1}$ contains the displacement of the center of the cylindrical region with respect to the center of the Cartesian mesh in direction $x$. This
+center is located at:
+$$
+x_{c}=\frac{x^{\text{Gc}}_{n^{x}}+x^{\text{GC}}_{0}}{2}
+$$
+where we have used $n^{x}=N^{\text{GC}}_{3}$.
+\item $x^{\text{GC}}_{j-1}$ and $x^{\text{GC}}_{j}$ are the lower and upper $x$ limits of mesh element $j$ ($j=1,n^{x}$).
+\end{itemize}
+The array $y^{\text{GC}}_{j}$ contains the following information:
+\begin{itemize}
+\item $y^{\text{GC}}_{-1}$ contains the displacement of the center of the cylindrical region with respect to the center of the Cartesian mesh in direction $y$. This
+center is located at:
+$$
+y_{c}=\frac{y^{\text{GC}}_{n^{y}}+y^{\text{GC}}_{0}}{2}
+$$
+where we have used $n^{y}=N^{\text{GC}}_{4}$.
+\item $y^{\text{GC}}_{j-1}$ and $y^{\text{GC}}_{j}$ are the lower and upper $y$ limits of mesh element $j$ ($j=1,n^{y}$).
+\end{itemize}
+The array $z^{\text{GC}}_{j}$ contains the following information:
+\begin{itemize}
+\item $z^{\text{GC}}_{-1}$ contains the displacement of the center of the cylindrical region with respect to the center of the Cartesian mesh in direction $z$. This
+center is located at:
+$$
+z_{c}=\frac{z^{\text{GC}}_{n^{z}}+z^{\text{GC}}_{0}}{2}
+$$
+where we have used $n^{z}=N^{\text{GC}}_{5}$.
+\item $z^{\text{GC}}_{j-1}$ and $z^{\text{GC}}_{j}$ are the lower and upper $z$ limits of mesh element $j$ ($j=1,n^{z}$).
+\end{itemize}
+  The array $r^{\text{GC}}_{j}$ contains the following information:
+\begin{itemize}
+\item $r^{\text{GC}}_{-1}=r^{\text{GC}}_{0}=0$. 
+\item $r^{\text{GC}}_{j-1}\le r\le r^{\text{GC}}_{j}$ describes the position in $r$ of mesh element $j$ ($j=1,N^{\text{GC}}_{2}$). 
+\end{itemize}
+The array $p^{\text{GC}}_{j}$ contains the following information:
+\begin{itemize}
+\item $p^{\text{GC}}_{-1}$ is the angular position of $z$-, $x$- or $y$-directed pin with respect to the $x$, $y$ or $z$ axis. 
+\item $p^{\text{GC}}_{0}$ is the radial position of $z$-, $x$- or $y$-directed pin with respect to the $x-y$, $y-z$ or $z-x$ center of the cell where the pin is located. 
+\item $p^{\text{GC}}_{1}$ is the height of a $x$-directed pin.
+\item $p^{\text{GC}}_{2}$ is the height of a $y$-directed pin.
+\item $p^{\text{GC}}_{3}$ is the height of a $z$-directed pin.
+\item $p^{\text{GC}}_{4}$ is the outer radius of the pin.
+\end{itemize}
+
+In Table~\ref{tabl:NXTCluster}, a description of the additional \namedir{NXTRecords} records associated with pin $i=1$ can be found. These records are identified using a procedure similar to that used for cell records, namely
+\begin{quote}
+\begin{verbatim}
+INTEGER      I
+CHARACTER*12 NAMREC
+CHARACTER*3  NREC
+WRITE(NAMREC,'(A1,I8.8,A3)') 'P',I,NREC
+\end{verbatim}
+\end{quote}
+where the variable \verb|NREC| can take the same values as for cell records, except for \moc{NREC}=\moc{PNT} and \moc{NREC}=\moc{PIN} which are now forbidden.
+
+\begin{DescriptionEnregistrement}{Pin $i=1$ records in \namedir{NXTRecords}}{8.0cm}\label{tabl:NXTCluster}
+\IntEnr
+  {P00000001DIM}{$40$}
+  {array $N^{\text{GP}}_{j}$ containing the dimensioning information required to rebuilt the pin} 
+\DbleEnr
+  {P00000001SMR}{$N^{\text{GP}}_{2}$}{cm}
+  {array $r^{\text{GP}}_{j}$ containing the pin radial mesh description}
+\DbleEnr
+  {P00000001SMX}{$N^{\text{GP}}_{3}$}{cm}
+  {array $x^{\text{GP}}_{j}$ containing the pin $x$-directed mesh description}
+\DbleEnr
+  {P00000001SMY}{$N^{\text{GP}}_{4}$}{cm}
+  {array $y^{\text{GP}}_{j}$ containing the pin $y$-directed mesh description}
+\DbleEnr
+  {P00000001SMZ}{$N^{\text{GP}}_{5}$}{cm}
+  {array $z^{\text{GP}}_{j}$ containing the pin $z$-directed mesh description}
+\IntEnr
+  {P00000001MIX}{$N^{\text{GP}}_{6}$}
+  {array $M^{\text{GP}}_{j}$ containing the pin physical mixture for each region}
+\IntEnr
+  {P00000001HOM}{$N^{\text{GP}}_{6}$}
+  {array $H^{\text{GP}}_{j}$ containing the pin virtual mixture for each region}
+\DbleEnr
+  {P00000001VSE}{$-N^{\text{GP}}_{9}:N^{\text{GP}}_{8}$}{}
+  {array $\text{SV}^{\text{GP}}_{j}$ containing surface area $j$ ($\text{SV}^{\text{GP}}_{-j}=S^{\text{GP}}_{j}$ in cm for 2-D and cm$^{2}$ for 3-D problems) and
+regional volumes $j$ ($\text{SV}^{\text{GP}}_{j}=V^{\text{GP}}_{j}$ in cm$^{2}$ for 2-D and cm$^{3}$ for 3-D problems)}
+\IntEnr
+  {P00000001VSI}{$4,-N^{\text{GP}}_{9}:N^{\text{GP}}_{8}$}
+  {array $\text{VSI}^{\text{GP}}_{k,j}$ containing the location of a surface ($j<0$) and a region ($j>0$)}
+\IntEnr
+  {P00000001RID}{$N^{\text{GP}}_{8}$}
+  {index array $\text{RID}^{\text{GP}}_{j}$ associating local and global region numbering}
+\IntEnr
+  {P00000001SID}{$N^{\text{GP}}_{9}$}
+  {index array $\text{SID}^{\text{GP}}_{j,i}$ associating local and global outer surface numbering}
+\end{DescriptionEnregistrement}
+
+The pin dimensioning array $N^{\text{GP}}$ contains the following information:
+\begin{itemize}
+\item $N^{\text{GP}}_{1,}$ pin geometry type (see the definition of $\mathcal{S}^{G}_{1}$ in \Sect{geometrydirmain});
+\item $N^{\text{GP}}_{2,}$ dimensions of the radial mesh array;
+\item $N^{\text{GP}}_{3,}$ dimensions of the $x$-directed mesh array;  
+\item $N^{\text{GP}}_{4,}$ dimensions of the $y$-directed mesh array;  
+\item $N^{\text{GP}}_{5,}$ dimensions of the $z$-directed mesh array;  
+\item $N^{\text{GP}}_{6,}$ dimensions of the mixture record;
+\item $N^{\text{GP}}_{7,}$ geometry level (2 for pins);
+\item $N^{\text{GP}}_{8,}$ number of regions in the pin before symmetry considerations;
+\item $N^{\text{GP}}_{9,}$ number of surfaces in the pin before symmetry considerations;
+\item $N^{\text{GP}}_{10}$ number of regions in the pin after symmetry considerations;
+\item $N^{\text{GP}}_{11}$ number of surfaces in the pin after symmetry considerations;
+\item $N^{\text{GP}}_{12}$ first global region number for pins in cluster;
+\item $N^{\text{GP}}_{13}$ last global region number for pins in cluster;
+\item $N^{\text{GP}}_{14}$ first global surface number for pins in cluster;
+\item $N^{\text{GP}}_{15}$ last global surface number for pins in cluster;
+\item $N^{\text{GP}}_{16}$ first pin cluster geometry for pins in cluster.
+\item $N^{\text{GP}}_{17}$ total number of pins in cluster.
+\end{itemize}
+while the remaining elements are not used. 
+The array $x^{\text{GP}}_{j}$ contains the following information:
+\begin{itemize}
+\item $x^{\text{GP}}_{-1}$ contains the displacement of the center of the cylindrical region with respect to the center of the Cartesian mesh in direction $x$. This
+center is located at:
+$$
+x_{c}=\frac{x^{\text{GP}}_{n^{x}}+x^{\text{GP}}_{0}}{2}
+$$
+where we have used $n^{x}=N^{\text{GP}}_{3}$.
+\item $x^{\text{GP}}_{j-1}$ and $x^{\text{GP}}_{j}$ re the lower and upper $x$ limits of mesh element $j$ ($j=1,n^{x}$).
+\end{itemize}
+The array $y^{\text{GP}}_{j}$ contains the following information:
+\begin{itemize}
+\item $y^{\text{GP}}_{-1}$ contains the displacement of the center of the cylindrical region with respect to the center of the Cartesian mesh in direction $y$. This
+center is located at:
+$$
+y_{c}=\frac{y^{\text{GP}}_{n^{y}}+y^{\text{GP}}_{0}}{2}
+$$
+where we have used $n^{y}=N^{\text{GP}}_{4}$.
+\item $y^{\text{GP}}_{j-1}$ and $y^{\text{GP}}_{j}$ are the lower and upper $y$ limits of mesh element $j$ ($j=1,n^{y}$).
+\end{itemize}
+The array $z^{\text{GP}}_{j}$ contains the following information:
+\begin{itemize}
+\item $z^{\text{GP}}_{-1}$ contains the displacement of the center of the cylindrical region with respect to the center of the Cartesian mesh in direction $z$. This
+center is located at:
+$$
+z_{c}=\frac{z^{\text{GP}}_{n^{z}}+z^{\text{GP}}_{0}}{2}
+$$
+where we have used $n^{z}=N^{\text{GP}}_{5}$.
+\item $z^{\text{GP}}_{j-1}$ and $z^{\text{GP}}_{j}$ are the lower and upper $z$ limits of mesh element $j$ ($j=1,n^{z}$).
+\end{itemize}
+  The array $r^{\text{GP}}_{j}$ contains the following information:
+\begin{itemize}
+\item $r^{\text{GP}}_{-1}=r^{\text{GP}}_{0}=0$. 
+\item $r^{\text{GP}}_{j-1}\le r\le r^{\text{GP}}_{j}$ describes the position in $r$ of mesh element $j$ with $j=1,N^{\text{GP}}_{2}$. 
+\end{itemize}
+
+Finally the \namedir{NXTRecords} directory also contains records associated with global identification of the surfaces and volumes as illustrated in Table~\ref{tabl:NXTGlobal}.
+
+\begin{DescriptionEnregistrement}{Global geometry records in \namedir{NXTRecords}}{8.0cm}\label{tabl:NXTGlobal}
+\DbleEnr
+  {TrackingDnsA}{$\mathcal{S}^{t}_{20}$}{cm}
+  {array $D_{i}$ containing the spatial spacing for each track direction}
+\DbleEnr
+  {TrackingDirc}{$N^{\text{GG}}_{1},\mathcal{S}^{t}_{20}$}{}
+  {array $\alpha_{j,i}$ containing the director cosine for axis $j$ for each track direction}
+\DbleEnr
+  {TrackingOrig}{$N^{\text{GG}}_{1},N_{p},\mathcal{S}^{t}_{20}$}{cm}
+  {array $L_{k,j,i}$ containing the origin in space ($k=1,N^{\text{GG}}_{1}$) and the direction of the normal plan for each plane $j$ and track direction $i$}
+\DbleEnr
+  {TrackingWgtD}{$\mathcal{S}^{t}_{20}$}{}
+  {array $W_{i}$ containing the integration weight for each track direction}
+\IntEnr
+  {NumMerge\blank{4}}{$N^{\text{GG}}_{23}$}
+  {merging indices used for normalizing the tracks}
+\DbleEnr
+  {VolMerge\blank{4}}{$N^{\text{GG}}_{23}$}{}
+  {volumes used for normalizing the tracks}
+\DbleEnr
+  {VTNormalize\blank{1}}{$N^{\text{GG}}_{23}$}{}
+  {array $R_{i}$ containing the ratio of the analytical and numerical volume for each region}
+\OptDbleEnr
+  {VTNormalizeD}{$N^{\text{GG}}_{23},\mathcal{S}^{t}_{20}$}{$\mathcal{S}^{t}_{8}=-1$}{}
+  {array $R_{i}$ containing the ratio of the analytical and numerical volume for region $i$ for each track direction}
+\IntEnr
+  {KEYMRG\blank{6}}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}
+  {array $\text{MRG}_{i}$ containing the global merging index}
+\IntEnr
+  {MATALB\blank{6}}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}
+  {array containing the albedo number associated with each surface and the physical mixture number associated with each region}
+\IntEnr
+  {HOMMATALB\blank{3}}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}
+  {array containing the albedo number associated with each surface and the virtual mixture number associated with each region}
+\DbleEnr
+  {SAreaRvolume}{$-N^{\text{GG}}_{22}:N^{\text{GG}}_{23}$}{}
+  {array containing the area ($S_{\alpha}$ in cm for 2-D and cm$^{2}$ for 3-D problems) and volumes ($V_{i}$ cm$^{2}$ for 2-D and
+cm$^{3}$ for 3-D problems) of each external surface and region in the geometry}
+\end{DescriptionEnregistrement}
+
+\clearpage