1 files changed, 137 insertions, 0 deletions
diff --git a/doc/IGE335/Section3.05_shi.tex b/doc/IGE335/Section3.05_shi.tex
new file mode 100644
index 0000000..92f7f22
--- /dev/null
+++ b/doc/IGE335/Section3.05_shi.tex
@@ -0,0 +1,137 @@
+\subsection{The {\tt SHI:} module}\label{sect:SHIData}
+
+The {\tt SHI:} module perform self-shielding calculations in DRAGON, using
+the generalized Stamm'ler method.\cite{SHIBA}  This approach is based on an heterogeneous-homogeneous equivalence principle. In this case, an {\sl equivalent dilution parameter} $\sigma_{{\rm e},g}$ is computed for each resonant isotope, in each resonant region and
+each resonant energy group $g$. This dilution parameter is used to interpolate pretabulated effective cross sections for the infinite homogeneous medium, previously obtained with the {\sl flux calculator} of the {\tt GROUPR} module in code NJOY.\cite{njoy2010}
+Each resonant isotope, identified as such by the \dusa{inrs}
+parameter defined in \Sect{LIBData}, is to be recalculated. The general format of
+the data for this module is:
+
+\begin{DataStructure}{Structure \dstr{SHI:}}
+\dusa{MICLIB} \moc{:=} \moc{SHI:} $\{$ \dusa{MICLIB} $|$ \dusa{OLDLIB} $\}$ 
+\dusa{TRKNAM} $[$ \dusa{TRKFIL} $]$ \moc{::} \dstr{descshi}
+\end{DataStructure}
+
+\noindent
+where
+
+\begin{ListeDeDescription}{mmmmmmmm}
+
+\item[\dusa{MICLIB}] {\tt character*12} name of the \dds{microlib} that will
+contain the microscopic and macroscopic cross sections updated by the
+self-shielding module. If
+\dusa{MICLIB} appears on both LHS and RHS, it is updated; otherwise, the
+internal library \dusa{OLDLIB} is copied into
+\dusa{MICLIB} and \dusa{MICLIB} is updated.
+
+\item[\dusa{OLDLIB}] {\tt character*12} name of a read-only \dds{microlib} 
+that is copied into \dusa{MICLIB}.
+
+\item[\dusa{TRKNAM}] {\tt character*12} name of the required \dds{tracking}
+data structure.
+
+\item[\dusa{TRKFIL}] {\tt character*12} name of the sequential binary tracking
+file used to store the tracks lengths. This file is given if and only if it was
+required in the previous tracking module call (see \Sect{TRKData}).
+
+\item[\dstr{descshi}] structure describing the self-shielding options.
+
+\end{ListeDeDescription}
+
+Each time the \moc{SHI:} module is called, a sub-directory is updated in the
+\dds{microlib} data structure to hold the last values defined in the
+\dstr{descshi} structure. The next time this module is called,
+these values will be used as floating defaults.
+
+\subsubsection{Data input for module {\tt SHI:}}\label{sect:descshi}
+
+\begin{DataStructure}{Structure \dstr{descshi}}
+$[$ \moc{EDIT} \dusa{iprint} $]$ \\
+$[$ \moc{GRMIN}  \dusa{lgrmin} $]~~[$ \moc{GRMAX}  \dusa{lgrmax} $]$ \\
+$[$ \moc{MXIT} \dusa{imxit} $]~~[$ \moc{EPS}  \dusa{valeps} $]$  \\
+$[~\{$ \moc{LJ} $|$ \moc{NOLJ} $\}~]$ $[~\{$ \moc{GC} $|$ \moc{NOGC} $\}$
+$[~\{$ \moc{TRAN} $|$ \moc{NOTR} $\}~]$
+$[$ \moc{LEVEL} \dusa{ilev} $]$ \\
+$[~\{$ \moc{PIJ} $|$ \moc{ARM} $\}~]$ \\
+{\tt ;}
+\end{DataStructure}
+
+\noindent  where
+
+\begin{ListeDeDescription}{mmmmmmmm}
+
+\item[\moc{EDIT}] keyword used to modify the print level \dusa{iprint}.
+
+\item[\dusa{iprint}] index used to control the printing of this module. The
+amount of output produced by this tracking module will vary substantially
+depending on the print level specified. 
+
+\item[\moc{GRMIN}] keyword to specify the minimum group number considered
+during the self-shielding process.
+
+\item[\dusa{lgrmin}] first group number considered during the
+self-shielding process. By default, \dusa{lgrmin} is set to the first group
+number containing self-shielding data in the library.
+
+\item[\moc{GRMAX}]  keyword to specify the maximum group number considered
+during the self-shielding process.
+
+\item[\dusa{lgrmax}] last group number considered during the self-shielding
+process. By default, \dusa{lgrmax} is set is set to the last group
+number containing self-shielding data in the library.
+
+\item[\moc{MXIT}]  keyword to specify the maximum number of iterations during
+the self-shielding process.
+
+\item[\dusa{imxit}] the maximum number of iterations. The default is
+\dusa{imxit}=20.
+
+\item[\moc{EPS}] keyword to specify the convergence criterion for the
+self-shielding iteration.
+
+\item[\dusa{valeps}] the convergence criterion for the self-shielding iteration.
+By default, \dusa{valeps}=$1.0\times 10^{-4}$.
+
+\item[\moc{LJ}] keyword to activate the Livolant and Jeanpierre
+normalization scheme which modifies the self-shielded averaged neutron fluxes in
+heterogeneous geometries. By default the Livolant and Jeanpierre
+normalization scheme is not activated.
+
+\item[\moc{NOLJ}] keyword to deactivate the Livolant and Jeanpierre
+normalization scheme which modifies the self-shielded averaged neutron fluxes in
+heterogeneous geometries. This is the default option.
+
+\item[\moc{GC}] keyword to activate the Goldstein-Cohen approximation in
+cases where Goldstein-Cohen parameters are stored on the internal library. These
+parameters may not be available with some libraries (e.g., {\tt APLIB1}, {\tt
+APLIB2} or MATXS-type libraries). The Goldstein-Cohen parameters can always be
+imposed using the \moc{IRSET} keyword of the \moc{LIB:} module (see
+\Sect{LIBData}). This is the default option.
+
+\item[\moc{NOGC}] keyword to deactivate the Goldstein-Cohen approximation in
+cases where Goldstein-Cohen parameters are stored on the internal library.
+
+\item[\moc{TRAN}] keyword to activate the transport correction option for
+self-shielding calculations (see \moc{CTRA} in \Sectand{MACData}{LIBData}). This is the default option.
+
+\item[\moc{NOTR}] keyword to deactivate the transport correction option for
+self-shielding calculations (see \moc{CTRA} in \Sectand{MACData}{LIBData}).
+
+\item[\moc{LEVEL}] keyword to specify the self-shielding model.
+
+\item[\dusa{ilev}] $=0$: original Stamm'ler model (without distributed effects);
+ $=1$: use the Nordheim (PIC) distributed self-shielding model\cite{toronto04};
+ $=2$: use both Nordheim (PIC) distributed self-shielding model and Riemann integration
+ method\cite{hasan}. By default, \dusa{ilev}$\,=0$.
+
+\item[\moc{PIJ}] keyword to specify the use of complete collision
+probabilities in the self-shielding calculations of {\tt SHI:}.
+This is the default option for \moc{EXCELT:} and \moc{SYBILT:} trackings.
+This option is not available for \moc{MCCGT:} trackings.
+
+\item[\moc{ARM}] keyword to specify the use of iterative flux techniques
+in the self-shielding calculations of {\tt SHI:}.
+This is the default option for \moc{MCCGT:} trackings.
+
+\end{ListeDeDescription}
+\eject