Index: /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_DOMAINcfg.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_DOMAINcfg.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_DOMAINcfg.tex	(revision 11598)
@@ -2,18 +2,28 @@
 
 \begin{document}
+
 \chapter{A brief guide to the DOMAINcfg tool}
 \label{apdx:DOMCFG}
 
+%    {\em 4.0} & {\em Andrew Coward} & {\em Created at v4.0 from materials removed from chap\_DOM that are still relevant to the \forcode{DOMAINcfg} tool and which illustrate and explain the choices to be made by the user when setting up new domains }  \\
+
+\thispagestyle{plain}
+
 \chaptertoc
-\vfill
-\begin{figure}[b]
-%% =================================================================================================
-\subsubsection*{Changes record}
-\begin{tabular}{m{0.08\linewidth}||m{0.32\linewidth}|m{0.6\linewidth}}
-    Release   & Author(s)     & Modifications \\
-\hline
-    {\em 4.0} & {\em Andrew Coward} & {\em Created at v4.0 from materials removed from chap\_DOM that are still relevant to the \forcode{DOMAINcfg} tool and which illustrate and explain the choices to be made by the user when setting up new domains }  \\
-\end{tabular}
-\end{figure}
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 This appendix briefly describes some of the options available in the
@@ -34,5 +44,4 @@
 \section{Choice of horizontal grid}
 \label{sec:DOMCFG_hor}
-
 
 \begin{listing}
@@ -298,13 +307,10 @@
 \np{nn_bathy}{nn\_bathy} (found in \nam{dom}{dom} namelist (\texttt{DOMAINCFG} variant) ):
 \begin{description}
-\item [{\np[=0]{nn_bathy}{nn\_bathy}}]:
-  a flat-bottom domain is defined.
+\item [{\np[=0]{nn_bathy}{nn\_bathy}}]: a flat-bottom domain is defined.
   The total depth $z_w (jpk)$ is given by the coordinate transformation.
   The domain can either be a closed basin or a periodic channel depending on the parameter \np{jperio}{jperio}.
-\item [{\np[=-1]{nn_bathy}{nn\_bathy}}]:
-  a domain with a bump of topography one third of the domain width at the central latitude.
+\item [{\np[=-1]{nn_bathy}{nn\_bathy}}]: a domain with a bump of topography one third of the domain width at the central latitude.
   This is meant for the "EEL-R5" configuration, a periodic or open boundary channel with a seamount.
-\item [{\np[=1]{nn_bathy}{nn\_bathy}}]:
-  read a bathymetry and ice shelf draft (if needed).
+\item [{\np[=1]{nn_bathy}{nn\_bathy}}]: read a bathymetry and ice shelf draft (if needed).
   The \ifile{bathy\_meter} file (Netcdf format) provides the ocean depth (positive, in meters) at
   each grid point of the model grid.
@@ -326,14 +332,10 @@
 After reading the bathymetry, the algorithm for vertical grid definition differs between the different options:
 \begin{description}
-\item [\forcode{ln_zco = .true.}]
-  set a reference coordinate transformation $z_0(k)$, and set $z(i,j,k,t) = z_0(k)$ where $z_0(k)$ is the closest match to the depth at $(i,j)$.
-\item [\forcode{ln_zps = .true.}]
-  set a reference coordinate transformation $z_0(k)$, and calculate the thickness of the deepest level at
+\item [\forcode{ln_zco = .true.}] set a reference coordinate transformation $z_0(k)$, and set $z(i,j,k,t) = z_0(k)$ where $z_0(k)$ is the closest match to the depth at $(i,j)$.
+\item [\forcode{ln_zps = .true.}] set a reference coordinate transformation $z_0(k)$, and calculate the thickness of the deepest level at
   each $(i,j)$ point using the bathymetry, to obtain the final three-dimensional depth and scale factor arrays.
-\item [\forcode{ln_sco = .true.}]
-  smooth the bathymetry to fulfill the hydrostatic consistency criteria and
+\item [\forcode{ln_sco = .true.}] smooth the bathymetry to fulfill the hydrostatic consistency criteria and
   set the three-dimensional transformation.
-\item [\forcode{s-z and s-zps}]
-  smooth the bathymetry to fulfill the hydrostatic consistency criteria and
+\item [\forcode{s-z and s-zps}] smooth the bathymetry to fulfill the hydrostatic consistency criteria and
   set the three-dimensional transformation $z(i,j,k)$,
   and possibly introduce masking of extra land points to better fit the original bathymetry file.
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_algos.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_algos.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_algos.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Note on some algorithms}
 \label{apdx:ALGOS}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 This appendix some on going consideration on algorithms used or planned to be used in \NEMO.
@@ -356,7 +374,5 @@
 This expression of the iso-neutral diffusion has been chosen in order to satisfy the following six properties:
 \begin{description}
-\item [$\bullet$ horizontal diffusion]
-  The discretization of the diffusion operator recovers the traditional five-point Laplacian in
-  the limit of flat iso-neutral direction:
+\item [Horizontal diffusion] The discretization of the diffusion operator recovers the traditional five-point Laplacian in the limit of flat iso-neutral direction:
   \[
     % \label{eq:ALGOS_Gf_property1a}
@@ -366,7 +382,5 @@
     { _i^k \mathbb{R}_{i_p}^{k_p} }=0
   \]
-
-\item [$\bullet$ implicit treatment in the vertical]
-  In the diagonal term associated with the vertical divergence of the iso-neutral fluxes
+\item [Implicit treatment in the vertical] In the diagonal term associated with the vertical divergence of the iso-neutral fluxes
   \ie\ the term associated with a second order vertical derivative)
   appears only tracer values associated with a single water column.
@@ -380,7 +394,5 @@
   \]
   can be quite large.
-
-\item [$\bullet$ pure iso-neutral operator]
-  The iso-neutral flux of locally referenced potential density is zero, \ie
+\item [Pure iso-neutral operator] The iso-neutral flux of locally referenced potential density is zero, \ie
   \begin{align*}
     % \label{eq:ALGOS_Gf_property2}
@@ -396,7 +408,5 @@
   This result is trivially obtained using the \autoref{eq:ALGOS_Gf_triads} applied to $T$ and $S$ and
   the definition of the triads' slopes \autoref{eq:ALGOS_Gf_slopes}.
-
-\item [$\bullet$ conservation of tracer]
-  The iso-neutral diffusion term conserve the total tracer content, \ie
+\item [Conservation of tracer] The iso-neutral diffusion term conserve the total tracer content, \ie
   \[
     % \label{eq:ALGOS_Gf_property1}
@@ -404,7 +414,5 @@
   \]
 This property is trivially satisfied since the iso-neutral diffusive operator is written in flux form.
-
-\item [$\bullet$ decrease of tracer variance]
-  The iso-neutral diffusion term does not increase the total tracer variance, \ie
+\item [Decrease of tracer variance] The iso-neutral diffusion term does not increase the total tracer variance, \ie
   \[
     % \label{eq:ALGOS_Gf_property1}
@@ -417,7 +425,5 @@
 It therfore ensures that, when the diffusivity coefficient is large enough,
 the field on which it is applied become free of grid-point noise.
-
-\item [$\bullet$ self-adjoint operator]
-  The iso-neutral diffusion operator is self-adjoint, \ie
+\item [Self-adjoint operator] The iso-neutral diffusion operator is self-adjoint, \ie
   \[
     % \label{eq:ALGOS_Gf_property1}
@@ -583,5 +589,4 @@
 \ie\ it does not include a diffusive component but is a "pure" advection term.
 
-$\ $\newpage      %force an empty line
 %% =================================================================================================
 \subsection{Discrete invariants of the iso-neutral diffrusion}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_diff_opers.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_diff_opers.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_diff_opers.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Diffusive Operators}
 \label{apdx:DIFFOPERS}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 %% =================================================================================================
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_invariants.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_invariants.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_invariants.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Discrete Invariants of the Equations}
 \label{apdx:INVARIANTS}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 %%%  Appendix put in gmcomment as it has not been updated for \zstar and s coordinate
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_s_coord.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_s_coord.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_s_coord.tex	(revision 11598)
@@ -6,18 +6,26 @@
 \label{apdx:SCOORD}
 
+%    {\em 4.0} & {\em Mike Bell} & {\em review}  \\
+%    {\em 3.x} & {\em Gurvan Madec} & {\em original}  \\
+
+\thispagestyle{plain}
+
 \chaptertoc
 
-\vfill
-\begin{figure}[b]
-%% =================================================================================================
-\subsubsection*{Changes record}
-\begin{tabular}{l||l|m{0.65\linewidth}}
-    Release   & Author        & Modifications \\
-    {\em 4.0} & {\em Mike Bell} & {\em review}  \\
-    {\em 3.x} & {\em Gurvan Madec} & {\em original}  \\
-\end{tabular}
-\end{figure}
-
-%% =================================================================================================
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
+
 \section{Chain rule for $s-$coordinates}
 \label{sec:SCOORD_chain}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_triads.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_triads.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/apdx_triads.tex	(revision 11598)
@@ -12,12 +12,29 @@
 
 \begin{document}
+
 \chapter{Iso-Neutral Diffusion and Eddy Advection using Triads}
 \label{apdx:TRIADS}
 
+\thispagestyle{plain}
+
 \chaptertoc
 
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
+
 %% =================================================================================================
 \section[Choice of \forcode{namtra\_ldf} namelist parameters]{Choice of \protect\nam{tra_ldf}{tra\_ldf} namelist parameters}
-
 
 Two scheme are available to perform the iso-neutral diffusion.
@@ -36,6 +53,5 @@
 The options specific to the Griffies scheme include:
 \begin{description}
-\item [{\np{ln_triad_iso}{ln\_triad\_iso}}]
-  See \autoref{sec:TRIADS_taper}.
+\item [{\np{ln_triad_iso}{ln\_triad\_iso}}] See \autoref{sec:TRIADS_taper}.
   If this is set false (the default),
   then `iso-neutral' mixing is accomplished within the surface mixed-layer along slopes linearly decreasing with
@@ -47,17 +63,15 @@
   giving an almost pure horizontal diffusive tracer flux within the mixed layer.
   This is similar to the tapering suggested by \citet{gerdes.koberle.ea_CD91}. See \autoref{subsec:TRIADS_Gerdes-taper}
-\item [{\np{ln_botmix_triad}{ln\_botmix\_triad}}]
-  See \autoref{sec:TRIADS_iso_bdry}.
+\item [{\np{ln_botmix_triad}{ln\_botmix\_triad}}] See \autoref{sec:TRIADS_iso_bdry}.
   If this is set false (the default) then the lateral diffusive fluxes
   associated with triads partly masked by topography are neglected.
   If it is set true, however, then these lateral diffusive fluxes are applied,
   giving smoother bottom tracer fields at the cost of introducing diapycnal mixing.
-\item [{\np{rn_sw_triad}{rn\_sw\_triad}}]
-  blah blah to be added....
+\item [{\np{rn_sw_triad}{rn\_sw\_triad}}] blah blah to be added....
 \end{description}
 The options shared with the Standard scheme include:
 \begin{description}
-\item [{\np{ln_traldf_msc}{ln\_traldf\_msc}}]   blah blah to be added
-\item [{\np{rn_slpmax}{rn\_slpmax}}]  blah blah to be added
+\item [{\np{ln_traldf_msc}{ln\_traldf\_msc}}] blah blah to be added
+\item [{\np{rn_slpmax}{rn\_slpmax}}]          blah blah to be added
 \end{description}
 
@@ -196,5 +210,4 @@
 % Instead of multiplying the mean slope calculated at the $u$-point by
 % the mean vertical gradient at the $u$-point,
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 \begin{figure}[tb]
   \centering
@@ -207,5 +220,4 @@
   \label{fig:TRIADS_ISO_triad}
 \end{figure}
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 They get the skew flux from the products of the vertical gradients at each $w$-point surrounding the $u$-point with
 the corresponding `triad' slope calculated from the lateral density gradient across the $u$-point divided by
@@ -258,5 +270,4 @@
 while the metrics are calculated at the $u$- and $w$-points on the arms.
 
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 \begin{figure}[tb]
   \centering
@@ -269,5 +280,4 @@
   \label{fig:TRIADS_qcells}
 \end{figure}
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 
 Each triad $\{_i^{k}\:_{i_p}^{k_p}\}$ is associated (\autoref{fig:TRIADS_qcells}) with the quarter cell that is
@@ -549,7 +559,7 @@
 where $b_T= e_{1T}\,e_{2T}\,e_{3T}$ is the volume of $T$-cells.
 The diffusion scheme satisfies the following six properties:
+
 \begin{description}
-\item [$\bullet$ horizontal diffusion]
-  The discretization of the diffusion operator recovers the traditional five-point Laplacian
+\item [Horizontal diffusion] The discretization of the diffusion operator recovers the traditional five-point Laplacian
   \autoref{eq:TRIADS_lat-normal} in the limit of flat iso-neutral direction:
   \[
@@ -560,7 +570,5 @@
     \text{when} \quad { _i^k \mathbb{R}_{i_p}^{k_p} }=0
   \]
-
-\item [$\bullet$ implicit treatment in the vertical]
-  Only tracer values associated with a single water column appear in the expression \autoref{eq:TRIADS_i33} for
+\item [Implicit treatment in the vertical] Only tracer values associated with a single water column appear in the expression \autoref{eq:TRIADS_i33} for
   the $_{33}$ fluxes, vertical fluxes driven by vertical gradients.
   This is of paramount importance since it means that a time-implicit algorithm can be used to
@@ -576,11 +584,7 @@
   \]
   (where $b_w= e_{1w}\,e_{2w}\,e_{3w}$ is the volume of $w$-cells) can be quite large.
-
-\item [$\bullet$ pure iso-neutral operator]
-  The iso-neutral flux of locally referenced potential density is zero.
+\item [Pure iso-neutral operator] The iso-neutral flux of locally referenced potential density is zero.
   See \autoref{eq:TRIADS_latflux-rho} and \autoref{eq:TRIADS_vertflux-triad2}.
-
-\item [$\bullet$ conservation of tracer]
-  The iso-neutral diffusion conserves tracer content, \ie
+\item [Conservation of tracer] The iso-neutral diffusion conserves tracer content, \ie
   \[
     % \label{eq:TRIADS_iso_property1}
@@ -588,7 +592,5 @@
   \]
   This property is trivially satisfied since the iso-neutral diffusive operator is written in flux form.
-
-\item [$\bullet$ no increase of tracer variance]
-  The iso-neutral diffusion does not increase the tracer variance, \ie
+\item [No increase of tracer variance] The iso-neutral diffusion does not increase the tracer variance, \ie
   \[
     % \label{eq:TRIADS_iso_property2}
@@ -601,7 +603,5 @@
   It therefore ensures that, when the diffusivity coefficient is large enough,
   the field on which it is applied becomes free of grid-point noise.
-
-\item [$\bullet$ self-adjoint operator]
-  The iso-neutral diffusion operator is self-adjoint, \ie
+\item [Self-adjoint operator] The iso-neutral diffusion operator is self-adjoint, \ie
   \begin{equation}
     \label{eq:TRIADS_iso_property3}
@@ -655,5 +655,4 @@
 (\np[=.true.]{ln_trabbl}{ln\_trabbl}, with \np[=1]{nn_bbl_ldf}{nn\_bbl\_ldf}), or for simple idealized problems.
 For setups with topography without bbl mixing, \np[=.true.]{ln_botmix_triad}{ln\_botmix\_triad} may be necessary.
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 \begin{figure}[h]
   \centering
@@ -679,5 +678,4 @@
   \label{fig:TRIADS_bdry_triads}
 \end{figure}
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 
 %% =================================================================================================
@@ -808,5 +806,4 @@
 \end{enumerate}
 
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 \begin{figure}[h]
   \centering
@@ -831,5 +828,4 @@
   \label{fig:TRIADS_MLB_triad}
 \end{figure}
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 
 %% =================================================================================================
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ASM.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ASM.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ASM.tex	(revision 11598)
@@ -2,19 +2,29 @@
 
 \begin{document}
+
 \chapter{Apply Assimilation Increments (ASM)}
 \label{chap:ASM}
 
+%    {\em 4.0} & {\em D. J. Lea} & {\em \NEMO\ 4.0 updates}  \\
+%    {\em 3.4} & {\em D. J. Lea, M. Martin, K. Mogensen, A. Weaver} & {\em Initial version}  \\
+
+\thispagestyle{plain}
+
 \chaptertoc
 
-\vfill
-\begin{figure}[b]
-%% =================================================================================================
-\subsubsection*{Changes record}
-\begin{tabular}{l||l|m{0.65\linewidth}}
-    Release   & Author        & Modifications \\
-    {\em 4.0} & {\em D. J. Lea} & {\em \NEMO\ 4.0 updates}  \\
-    {\em 3.4} & {\em D. J. Lea, M. Martin, K. Mogensen, A. Weaver} & {\em Initial version}  \\
-\end{tabular}
-\end{figure}
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 The ASM code adds the functionality to apply increments to the model variables: temperature, salinity,
@@ -25,6 +35,4 @@
 There is a brief description of all the namelist options provided.
 To build the ASM code \key{asminc} must be set.
-
-%===============================================================
 
 %% =================================================================================================
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIA.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIA.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIA.tex	(revision 11598)
@@ -2,22 +2,32 @@
 
 \begin{document}
+
 \chapter{Output and Diagnostics (IOM, DIA, TRD, FLO)}
 \label{chap:DIA}
 
+%    {\em 4.0} & {\em Mirek Andrejczuk, Massimiliano Drudi} & {\em }  \\
+%    {\em }      & {\em Dorotea Iovino, Nicolas Martin} & {\em }  \\
+%    {\em 3.6} & {\em Gurvan Madec, Sebastien Masson } & {\em }  \\
+%    {\em 3.4} & {\em Gurvan Madec, Rachid Benshila, Andrew Coward } & {\em }  \\
+%    {\em }      & {\em Christian Ethe, Sebastien Masson } & {\em }  \\
+
+\thispagestyle{plain}
+
 \chaptertoc
 
-\vfill
-\begin{figure}[b]
-%% =================================================================================================
-\subsubsection*{Changes record}
-\begin{tabular}{l||l|m{0.65\linewidth}}
-    Release   & Author        & Modifications \\
-    {\em 4.0} & {\em Mirek Andrejczuk, Massimiliano Drudi} & {\em }  \\
-    {\em }      & {\em Dorotea Iovino, Nicolas Martin} & {\em }  \\
-    {\em 3.6} & {\em Gurvan Madec, Sebastien Masson } & {\em }  \\
-    {\em 3.4} & {\em Gurvan Madec, Rachid Benshila, Andrew Coward } & {\em }  \\
-    {\em }      & {\em Christian Ethe, Sebastien Masson } & {\em }  \\
-\end{tabular}
-\end{figure}
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 %% =================================================================================================
@@ -134,13 +144,13 @@
 
 If an additional variable must be written to a restart file, the following steps are needed:
-\begin{description}
-   \item [step 1:] add variable name to a list of restart variables (in subroutine \rou{iom\_set\_rst\_vars,} \mdl{iom}) and
+\begin{enumerate}
+\item Add variable name to a list of restart variables (in subroutine \rou{iom\_set\_rst\_vars,} \mdl{iom}) and
 define correct grid for the variable (\forcode{grid_N_3D} - 3D variable, \forcode{grid_N} - 2D variable, \forcode{grid_vector} -
 1D variable, \forcode{grid_scalar} - scalar),
-   \item [step 2:] add variable to the list of fields written by restart.  This can be done either in subroutine
+\item Add variable to the list of fields written by restart.  This can be done either in subroutine
 \rou{iom\_set\_rstw\_core} (\mdl{iom}) or by calling  \rou{iom\_set\_rstw\_active} (\mdl{iom}) with the name of a variable
 as an argument. This convention follows approach for writing restart using iom, where variables are
 written either by \rou{rst\_write} or by calling \rou{iom\_rstput} from individual routines.
-\end{description}
+\end{enumerate}
 
 An older versions of XIOS do not support reading functionality. It's recommended to use at least XIOS2@1451.
@@ -266,14 +276,10 @@
 
 \begin{enumerate}
-\item [1.]
-  in \NEMO\ code, add a \forcode{CALL iom_put( 'identifier', array )} where you want to output a 2D or 3D array.
-\item [2.]
-  If necessary, add \forcode{USE iom ! I/O manager library} to the list of used modules in
+\item in \NEMO\ code, add a \forcode{CALL iom_put( 'identifier', array )} where you want to output a 2D or 3D array.
+\item If necessary, add \forcode{USE iom ! I/O manager library} to the list of used modules in
   the upper part of your module.
-\item [3.]
-  in the field\_def.xml file, add the definition of your variable using the same identifier you used in the f90 code
+\item in the field\_def.xml file, add the definition of your variable using the same identifier you used in the f90 code
   (see subsequent sections for a details of the XML syntax and rules).
   For example:
-
 \begin{xmllines}
 <field_definition>
@@ -284,5 +290,4 @@
 </field_definition>
 \end{xmllines}
-
 Note your definition must be added to the field\_group whose reference grid is consistent with the size of
 the array passed to iomput.
@@ -291,17 +296,13 @@
 (iom\_set\_domain\_attr and iom\_set\_axis\_attr in \mdl{iom}) or defined in the domain\_def.xml file.
 \eg:
-
 \begin{xmllines}
 <grid id="grid_T_3D" domain_ref="grid_T" axis_ref="deptht"/>
 \end{xmllines}
-
 Note, if your array is computed within the surface module each \np{nn_fsbc}{nn\_fsbc} time\_step,
 add the field definition within the field\_group defined with the id "SBC":
 \xmlcode{<field_group id="SBC" ...>} which has been defined with the correct frequency of operations
 (iom\_set\_field\_attr in \mdl{iom})
-\item [4.]
-  add your field in one of the output files defined in iodef.xml
+\item add your field in one of the output files defined in iodef.xml
   (again see subsequent sections for syntax and rules)
-
 \begin{xmllines}
 <file id="file1" .../>
@@ -311,5 +312,4 @@
 </file>
 \end{xmllines}
-
 \end{enumerate}
 
@@ -1342,5 +1342,4 @@
 \NEMO\ executables linked with NetCDF4 libraries can be made to produce NetCDF3 files by
 setting the \np{ln_nc4zip}{ln\_nc4zip} logical to false in the \nam{nc4}{nc4} namelist:
-
 
 \begin{listing}
@@ -1442,5 +1441,4 @@
 \label{sec:DIA_trd}
 
-
 \begin{listing}
   \nlst{namtrd}
@@ -1458,22 +1456,15 @@
 
 \begin{description}
-\item [{\np{ln_glo_trd}{ln\_glo\_trd}}]:
-  at each \np{nn_trd}{nn\_trd} time-step a check of the basin averaged properties of
+\item [{\np{ln_glo_trd}{ln\_glo\_trd}}]: at each \np{nn_trd}{nn\_trd} time-step a check of the basin averaged properties of
   the momentum and tracer equations is performed.
   This also includes a check of $T^2$, $S^2$, $\tfrac{1}{2} (u^2+v2)$,
   and potential energy time evolution equations properties;
-\item [{\np{ln_dyn_trd}{ln\_dyn\_trd}}]:
-  each 3D trend of the evolution of the two momentum components is output;
-\item [{\np{ln_dyn_mxl}{ln\_dyn\_mxl}}]:
-  each 3D trend of the evolution of the two momentum components averaged over the mixed layer is output;
-\item [{\np{ln_vor_trd}{ln\_vor\_trd}}]:
-  a vertical summation of the moment tendencies is performed,
+\item [{\np{ln_dyn_trd}{ln\_dyn\_trd}}]: each 3D trend of the evolution of the two momentum components is output;
+\item [{\np{ln_dyn_mxl}{ln\_dyn\_mxl}}]: each 3D trend of the evolution of the two momentum components averaged over the mixed layer is output;
+\item [{\np{ln_vor_trd}{ln\_vor\_trd}}]: a vertical summation of the moment tendencies is performed,
   then the curl is computed to obtain the barotropic vorticity tendencies which are output;
-\item [{\np{ln_KE_trd}{ln\_KE\_trd}}] :
-  each 3D trend of the Kinetic Energy equation is output;
-\item [{\np{ln_tra_trd}{ln\_tra\_trd}}]:
-  each 3D trend of the evolution of temperature and salinity is output;
-\item [{\np{ln_tra_mxl}{ln\_tra\_mxl}}]:
-  each 2D trend of the evolution of temperature and salinity averaged over the mixed layer is output;
+\item [{\np{ln_KE_trd}{ln\_KE\_trd}}]  : each 3D trend of the Kinetic Energy equation is output;
+\item [{\np{ln_tra_trd}{ln\_tra\_trd}}]: each 3D trend of the evolution of temperature and salinity is output;
+\item [{\np{ln_tra_mxl}{ln\_tra\_mxl}}]: each 2D trend of the evolution of temperature and salinity averaged over the mixed layer is output;
 \end{description}
 
@@ -1638,5 +1629,4 @@
 \section[Transports across sections (\texttt{\textbf{key\_diadct}})]{Transports across sections (\protect\key{diadct})}
 \label{sec:DIA_diag_dct}
-
 
 \begin{listing}
@@ -1993,7 +1983,4 @@
 \end{figure}
 
-% -----------------------------------------------------------
-%       CMIP specific diagnostics
-% -----------------------------------------------------------
 %% =================================================================================================
 \subsection[CMIP specific diagnostics (\textit{diaar5.F90}, \textit{diaptr.F90})]{CMIP specific diagnostics (\protect\mdl{diaar5})}
@@ -2015,5 +2002,4 @@
 the Indo-Pacific mask been deduced from the sum of the Indian and Pacific mask (\autoref{fig:DIA_mask_subasins}).
 
-
 \begin{listing}
   \nlst{namptr}
@@ -2022,10 +2008,6 @@
 \end{listing}
 
-% -----------------------------------------------------------
-%       25 hour mean and hourly Surface, Mid and Bed
-% -----------------------------------------------------------
 %% =================================================================================================
 \subsection{25 hour mean output for tidal models}
-
 
 \begin{listing}
@@ -2040,10 +2022,6 @@
 This diagnostic is actived with the logical $ln\_dia25h$.
 
-% -----------------------------------------------------------
-%     Top Middle and Bed hourly output
-% -----------------------------------------------------------
 %% =================================================================================================
 \subsection{Top middle and bed hourly output}
-
 
 \begin{listing}
@@ -2059,7 +2037,4 @@
 This diagnostic is actived with the logical $ln\_diatmb$.
 
-% -----------------------------------------------------------
-%     Courant numbers
-% -----------------------------------------------------------
 %% =================================================================================================
 \subsection{Courant numbers}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIU.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIU.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIU.tex	(revision 11598)
@@ -2,10 +2,26 @@
 
 \begin{document}
+
 \chapter{Diurnal SST Models (DIU)}
 \label{chap:DIU}
 
+\thispagestyle{plain}
+
 \chaptertoc
 
-$\ $\newline % force a new line
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 Code to produce an estimate of the diurnal warming and cooling of the sea surface skin
@@ -37,8 +53,6 @@
 This namelist contains only two variables:
 \begin{description}
-\item [{\np{ln_diurnal}{ln\_diurnal}}]
-  A logical switch for turning on/off both the cool skin and warm layer.
-\item [{\np{ln_diurnal_only}{ln\_diurnal\_only}}]
-  A logical switch which if \forcode{.true.} will run the diurnal model without the other dynamical parts of \NEMO.
+\item [{\np{ln_diurnal}{ln\_diurnal}}] A logical switch for turning on/off both the cool skin and warm layer.
+\item [{\np{ln_diurnal_only}{ln\_diurnal\_only}}] A logical switch which if \forcode{.true.} will run the diurnal model without the other dynamical parts of \NEMO.
   \np{ln_diurnal_only}{ln\_diurnal\_only} must be \forcode{.false.} if \np{ln_diurnal}{ln\_diurnal} is \forcode{.false.}.
 \end{description}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex	(revision 11598)
@@ -2,8 +2,7 @@
 
 \begin{document}
+
 \chapter{Space Domain (DOM)}
 \label{chap:DOM}
-
-\chaptertoc
 
 % Missing things:
@@ -15,21 +14,30 @@
 %     - domclo:  closed sea and lakes.... management of closea sea area : specific to global configuration, both forced and coupled
 
-\vfill
-
-\begin{table}[b]
-  \footnotesize
-  \caption*{Changes record}
+%    {\em 4.0} & {\em Simon M\"{u}ller \& Andrew Coward} &
+%    {\em
+%      Compatibility changes Major simplification has moved many of the options to external domain configuration tools.
+%      (see \autoref{apdx:DOMCFG})
+%    }                                                                                            \\
+%    {\em 3.x} & {\em Rachid Benshila, Gurvan Madec \& S\'{e}bastien Masson} &
+%    {\em First version}                                                                          \\
+
+\thispagestyle{plain}
+
+\chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
   \begin{tabularx}{\textwidth}{l||X|X}
-    Release & Author(s) & Modifications                                                          \\
-    \hline
-    {\em 4.0} & {\em Simon M\"{u}ller \& Andrew Coward} &
-    {\em
-      Compatibility changes Major simplification has moved many of the options to external domain configuration tools.
-      (see \autoref{apdx:DOMCFG})
-    }                                                                                            \\
-    {\em 3.x} & {\em Rachid Benshila, Gurvan Madec \& S\'{e}bastien Masson} &
-    {\em First version}                                                                          \\
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
   \end{tabularx}
-\end{table}
+}
+
+\clearpage
 
 Having defined the continuous equations in \autoref{chap:MB} and chosen a time discretisation \autoref{chap:TD},
@@ -259,7 +267,4 @@
 Furthermore, the direction of the vertical indexing has been reversed and the surface level set at $k = 1$.
 
-% -----------------------------------
-%        Horizontal Indexing
-% -----------------------------------
 %% =================================================================================================
 \subsubsection{Horizontal indexing}
@@ -272,7 +277,4 @@
 A $t$-point and its nearest north-east $f$-point have the same $i$-and $j$-indices.
 
-% -----------------------------------
-%        Vertical indexing
-% -----------------------------------
 %% =================================================================================================
 \subsubsection{Vertical indexing}
@@ -342,7 +344,4 @@
 and the vertical grid (\autoref{subsec:DOM_zgr}).
 
-% -----------------------------------
-%        Domain Size
-% -----------------------------------
 %% =================================================================================================
 \subsection{Domain size}
@@ -639,12 +638,10 @@
 
 \begin{description}
-\item [{\np[=.true.]{ln_tsd_init}{ln\_tsd\_init}}]
-  Use T and S input files that can be given on the model grid itself or on their native input data grids.
+\item [{\np[=.true.]{ln_tsd_init}{ln\_tsd\_init}}] Use T and S input files that can be given on the model grid itself or on their native input data grids.
   In the latter case, the data will be interpolated on-the-fly both in the horizontal and the vertical to the model grid
   (see \autoref{subsec:SBC_iof}).
   The information relating to the input files are specified in the \np{sn_tem}{sn\_tem} and \np{sn_sal}{sn\_sal} structures.
   The computation is done in the \mdl{dtatsd} module.
-\item [{\np[=.false.]{ln_tsd_init}{ln\_tsd\_init}}]
-  Initial values for T and S are set via a user supplied \rou{usr\_def\_istate} routine contained in \mdl{userdef\_istate}.
+\item [{\np[=.false.]{ln_tsd_init}{ln\_tsd\_init}}] Initial values for T and S are set via a user supplied \rou{usr\_def\_istate} routine contained in \mdl{userdef\_istate}.
   The default version sets horizontally uniform T and profiles as used in the GYRE configuration
   (see \autoref{sec:CFGS_gyre}).
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Ocean Dynamics (DYN)}
 \label{chap:DYN}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 Using the representation described in \autoref{chap:DOM},
@@ -149,5 +167,4 @@
 (see \autoref{subsec:DOM_Num_Index_vertical}).
 
-
 %% =================================================================================================
 \section{Coriolis and advection: vector invariant form}
@@ -301,5 +318,4 @@
   \label{fig:DYN_een_triad}
 \end{figure}
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
 
 A key point in \autoref{eq:DYN_een_e3f} is how the averaging in the \textbf{i}- and \textbf{j}- directions is made.
@@ -387,5 +403,4 @@
 a similar split-explicit time stepping should be used on vertical advection of tracer to ensure a better stability,
 an option which is only available with a TVD scheme (see \np{ln_traadv_tvd_zts}{ln\_traadv\_tvd\_zts} in \autoref{subsec:TRA_adv_tvd}).
-
 
 %% =================================================================================================
@@ -834,5 +849,4 @@
 (see their figure 12, lower left).
 
-%>   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >
 \begin{figure}[!t]
   \centering
@@ -859,5 +873,4 @@
   \label{fig:DYN_spg_ts}
 \end{figure}
-%>   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >
 
 In the default case (\np[=.true.]{ln_bt_fw}{ln\_bt\_fw}),
@@ -909,5 +922,4 @@
 it is still significant as shown by \citet{levier.treguier.ea_rpt07} in the case of an analytical barotropic Kelvin wave.
 
-%>>>>>===============
 \gmcomment{               %%% copy from griffies Book
 
@@ -1028,6 +1040,4 @@
 
 }            %%end gm comment (copy of griffies book)
-
-%>>>>>===============
 
 %% =================================================================================================
@@ -1251,5 +1261,4 @@
 the snow-ice mass is taken into account when computing the surface pressure gradient.
 
-
 \gmcomment{ missing : the lateral boundary condition !!!   another external forcing
  }
@@ -1307,5 +1316,4 @@
 The final sub-section covers some additional considerations that are relevant to both schemes.
 
-
 %   Iterative limiters
 %% =================================================================================================
@@ -1336,5 +1344,4 @@
 to be its depth times its velocity. This depth is considered to be zero at ``dry'' $u$-points consistent with its
 treatment in the calculation of the flux of mass across the cell face.
-
 
 \cite{warner.defne.ea_CG13} state that in their scheme the velocity masks at the cell faces for the baroclinic
@@ -1462,5 +1469,4 @@
 directional limiter does.
 
-
 %      Surface pressure gradients
 %% =================================================================================================
@@ -1481,5 +1487,4 @@
 column.  The three possible combinations are illustrated in \autoref{fig:DYN_WAD_dynhpg}.
 
-%>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
 \begin{figure}[!ht]
   \centering
@@ -1491,5 +1496,4 @@
   \label{fig:DYN_WAD_dynhpg}
 \end{figure}
-%>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
 
 The first logical, $\mathrm{ll\_tmp1}$, is set to true if and only if the water depth at
@@ -1554,10 +1558,7 @@
 See the WAD tests MY\_DOC documention for details of the WAD test cases.
 
-
-
 %% =================================================================================================
 \section[Time evolution term (\textit{dynnxt.F90})]{Time evolution term (\protect\mdl{dynnxt})}
 \label{sec:DYN_nxt}
-
 
 Options are defined through the \nam{dom}{dom} namelist variables.
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LBC.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LBC.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LBC.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Lateral Boundary Condition (LBC)}
 \label{chap:LBC}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 %gm% add here introduction to this chapter
@@ -164,16 +182,13 @@
 \begin{description}
 
-\item [For closed boundary (\jp{jperio}\forcode{=0})],
-  solid walls are imposed at all model boundaries:
+\item [For closed boundary (\jp{jperio}\forcode{=0})], solid walls are imposed at all model boundaries:
   first and last rows and columns are set to zero.
 
-\item [For cyclic east-west boundary (\jp{jperio}\forcode{=1})],
-  first and last rows are set to zero (closed) whilst the first column is set to
+\item [For cyclic east-west boundary (\jp{jperio}\forcode{=1})], first and last rows are set to zero (closed) whilst the first column is set to
   the value of the last-but-one column and the last column to the value of the second one
   (\autoref{fig:LBC_jperio}-a).
   Whatever flows out of the eastern (western) end of the basin enters the western (eastern) end.
 
-\item [For cyclic north-south boundary (\jp{jperio}\forcode{=2})],
-  first and last columns are set to zero (closed) whilst the first row is set to
+\item [For cyclic north-south boundary (\jp{jperio}\forcode{=2})], first and last columns are set to zero (closed) whilst the first row is set to
   the value of the last-but-one row and the last row to the value of the second one
   (\autoref{fig:LBC_jperio}-a).
@@ -215,5 +230,4 @@
 \section[Exchange with neighbouring processors (\textit{lbclnk.F90}, \textit{lib\_mpp.F90})]{Exchange with neighbouring processors (\protect\mdl{lbclnk}, \protect\mdl{lib\_mpp})}
 \label{sec:LBC_mpp}
-
 
 \begin{listing}
@@ -325,5 +339,4 @@
 \section{Unstructured open boundary conditions (BDY)}
 \label{sec:LBC_bdy}
-
 
 \begin{listing}
@@ -645,5 +658,4 @@
 \label{subsec:LBC_bdy_tides}
 
-
 \begin{listing}
   \nlst{nambdy_tide}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex	(revision 11598)
@@ -6,5 +6,22 @@
 \label{chap:LDF}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 The lateral physics terms in the momentum and tracer equations have been described in \autoref{eq:MB_zdf} and
@@ -22,6 +39,4 @@
 is described in \autoref{apdx:TRIADS}
 
-
-
 %% =================================================================================================
 \section[Lateral mixing operators]{Lateral mixing operators}
@@ -144,18 +159,12 @@
 
 \begin{description}
-
-\item [$z$-coordinate with full step: ]
-  in \autoref{eq:LDF_slp_iso} the densities appearing in the $i$ and $j$ derivatives  are taken at the same depth,
+\item [$z$-coordinate with full step:] in \autoref{eq:LDF_slp_iso} the densities appearing in the $i$ and $j$ derivatives  are taken at the same depth,
   thus the $in situ$ density can be used.
   This is not the case for the vertical derivatives: $\delta_{k+1/2}[\rho]$ is replaced by $-\rho N^2/g$,
   where $N^2$ is the local Brunt-Vais\"{a}l\"{a} frequency evaluated following \citet{mcdougall_JPO87}
   (see \autoref{subsec:TRA_bn2}).
-
-\item [$z$-coordinate with partial step: ]
-  this case is identical to the full step case except that at partial step level,
+\item [$z$-coordinate with partial step:] this case is identical to the full step case except that at partial step level,
   the \emph{horizontal} density gradient is evaluated as described in \autoref{sec:TRA_zpshde}.
-
-\item [$s$- or hybrid $s$-$z$- coordinate: ]
-  in the current release of \NEMO, iso-neutral mixing is only employed for $s$-coordinates if
+\item [$s$- or hybrid $s$-$z$- coordinate:] in the current release of \NEMO, iso-neutral mixing is only employed for $s$-coordinates if
   the Griffies scheme is used (\np[=.true.]{ln_traldf_triad}{ln\_traldf\_triad};
   see \autoref{apdx:TRIADS}).
@@ -207,5 +216,4 @@
 
 Note that such a formulation could be also used in the $z$-coordinate and $z$-coordinate with partial steps cases.
-
 \end{description}
 
@@ -464,5 +472,4 @@
 \label{sec:LDF_eiv}
 
-
 \begin{listing}
   \nlst{namtra_eiv}
@@ -470,5 +477,4 @@
   \label{lst:namtra_eiv}
 \end{listing}
-
 
 %%gm  from Triad appendix  : to be incorporated....
@@ -529,5 +535,4 @@
 \label{sec:LDF_mle}
 
-
 \begin{listing}
   \nlst{namtra_mle}
@@ -536,5 +541,4 @@
 \end{listing}
 
-
 If  \np[=.true.]{ln_mle}{ln\_mle} in \nam{tra_mle}{tra\_mle} namelist, a parameterization of the mixing due to unresolved mixed layer instabilities is activated (\citet{foxkemper.ferrari_JPO08}). Additional transport is computed in \rou{ldf\_mle\_trp} and added to the eulerian transport in \rou{tra\_adv} as done for eddy induced advection.
 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_OBS.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_OBS.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_OBS.tex	(revision 11598)
@@ -2,21 +2,35 @@
 
 \begin{document}
+
 \chapter{Observation and Model Comparison (OBS)}
 \label{chap:OBS}
 
+%\subsubsection*{Changes record}
+%\begin{tabular}{l||l|m{0.65\linewidth}}
+%    Release   & Author        & Modifications \\
+%    {\em 4.0} & {\em D. J. Lea} & {\em \NEMO\ 4.0 updates}  \\
+%    {\em 3.6} & {\em M. Martin, A. Ryan} & {\em Add averaging operator, standalone obs oper} \\
+%    {\em 3.4} & {\em D. J. Lea, M. Martin, ...} & {\em Initial version}  \\
+%    {\em --\texttt{"}--} & {\em ... K. Mogensen, A. Vidard, A. Weaver} & {\em ---\texttt{"}---}  \\
+%\end{tabular}
+
+\thispagestyle{plain}
+
 \chaptertoc
 
-\vfill
-\begin{figure}[b]
-%% =================================================================================================
-\subsubsection*{Changes record}
-\begin{tabular}{l||l|m{0.65\linewidth}}
-    Release   & Author        & Modifications \\
-    {\em 4.0} & {\em D. J. Lea} & {\em \NEMO\ 4.0 updates}  \\
-    {\em 3.6} & {\em M. Martin, A. Ryan} & {\em Add averaging operator, standalone obs oper} \\
-    {\em 3.4} & {\em D. J. Lea, M. Martin, ...} & {\em Initial version}  \\
-    {\em --\texttt{"}--} & {\em ... K. Mogensen, A. Vidard, A. Weaver} & {\em ---\texttt{"}---}  \\
-\end{tabular}
-\end{figure}
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 The observation and model comparison code, the observation operator (OBS), reads in observation files
@@ -111,5 +125,4 @@
 Here we show a more complete example namelist \nam{obs}{obs} and also show the NetCDF headers of
 the observation files that may be used with the observation operator.
-
 
 \begin{listing}
@@ -608,6 +621,5 @@
 
 \begin{enumerate}
-
-\item [1.] {\bfseries Great-Circle distance-weighted interpolation.}
+\item {\bfseries Great-Circle distance-weighted interpolation.}
   The weights are computed as a function of the great-circle distance $s(P, \cdot)$ between $P$ and
   the model grid points $A$, $B$ etc.
@@ -654,5 +666,5 @@
    \end{alignat*}
 
-\item [2.] {\bfseries Great-Circle distance-weighted interpolation with small angle approximation.}
+\item {\bfseries Great-Circle distance-weighted interpolation with small angle approximation.}
   Similar to the previous interpolation but with the distance $s$ computed as
   \begin{alignat*}{2}
@@ -664,8 +676,8 @@
   where $M$ corresponds to $A$, $B$, $C$ or $D$.
 
-\item [3.] {\bfseries Bilinear interpolation for a regular spaced grid.}
+\item {\bfseries Bilinear interpolation for a regular spaced grid.}
   The interpolation is split into two 1D interpolations in the longitude and latitude directions, respectively.
 
-\item [4.] {\bfseries Bilinear remapping interpolation for a general grid.}
+\item {\bfseries Bilinear remapping interpolation for a general grid.}
   An iterative scheme that involves first mapping a quadrilateral cell into
   a cell with coordinates (0,0), (1,0), (0,1) and (1,1).
@@ -705,5 +717,4 @@
   \label{fig:OBS_avgrec}
 \end{figure}
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 
 \begin{figure}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_SBC.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_SBC.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_SBC.tex	(revision 11598)
@@ -6,6 +6,22 @@
 \label{chap:SBC}
 
+\thispagestyle{plain}
+
 \chaptertoc
 
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 \begin{listing}
@@ -212,11 +228,9 @@
 where
 \begin{description}
-\item [File name]:
-  the stem name of the NetCDF file to be opened.
+\item [File name]: the stem name of the NetCDF file to be opened.
   This stem will be completed automatically by the model, with the addition of a '.nc' at its end and
   by date information and possibly a prefix (when using AGRIF).
   \autoref{tab:SBC_fldread} provides the resulting file name in all possible cases according to
   whether it is a climatological file or not, and to the open/close frequency (see below for definition).
-
   \begin{table}[htbp]
     \centering
@@ -245,17 +259,11 @@
     \label{tab:SBC_fldread}
   \end{table}
-
-\item [Record frequency]:
-  the frequency of the records contained in the input file.
+\item [Record frequency]: the frequency of the records contained in the input file.
   Its unit is in hours if it is positive (for example 24 for daily forcing) or in months if negative
   (for example -1 for monthly forcing or -12 for annual forcing).
   Note that this frequency must REALLY be an integer and not a real.
   On some computers, setting it to '24.' can be interpreted as 240!
-
-\item [Variable name]:
-  the name of the variable to be read in the input NetCDF file.
-
-\item [Time interpolation]:
-  a logical to activate, or not, the time interpolation.
+\item [Variable name]: the name of the variable to be read in the input NetCDF file.
+\item [Time interpolation]: a logical to activate, or not, the time interpolation.
   If set to 'false', the forcing will have a steplike shape remaining constant during each forcing period.
   For example, when using a daily forcing without time interpolation, the forcing remaining constant from
@@ -265,13 +273,9 @@
   For example, when using a daily forcing with time interpolation,
   linear interpolation will be performed between mid-day of two consecutive days.
-
-\item [Climatological forcing]:
-  a logical to specify if a input file contains climatological forcing which can be cycle in time,
+\item [Climatological forcing]: a logical to specify if a input file contains climatological forcing which can be cycle in time,
   or an interannual forcing which will requires additional files if
   the period covered by the simulation exceeds the one of the file.
   See the above file naming strategy which impacts the expected name of the file to be opened.
-
-\item [Open/close frequency]:
-  the frequency at which forcing files must be opened/closed.
+\item [Open/close frequency]: the frequency at which forcing files must be opened/closed.
   Four cases are coded:
   'daily', 'weekLLL' (with 'LLL' the first 3 letters of the first day of the week), 'monthly' and 'yearly' which
@@ -280,9 +284,6 @@
   For example, the first record of a yearly file containing daily data is Jan 1st even if
   the experiment is not starting at the beginning of the year.
-
-\item [Others]:
-  'weights filename', 'pairing rotation' and 'land/sea mask' are associated with
+\item [Others]:  'weights filename', 'pairing rotation' and 'land/sea mask' are associated with
   on-the-fly interpolation which is described in \autoref{subsec:SBC_iof}.
-
 \end{description}
 
@@ -449,5 +450,4 @@
 \label{subsec:SBC_SAS}
 
-
 \begin{listing}
   \nlst{namsbc_sas}
@@ -477,22 +477,16 @@
 
 \begin{itemize}
-\item \mdl{nemogcm}:
-  This routine initialises the rest of the model and repeatedly calls the stp time stepping routine (\mdl{step}).
+\item \mdl{nemogcm}: This routine initialises the rest of the model and repeatedly calls the stp time stepping routine (\mdl{step}).
   Since the ocean state is not calculated all associated initialisations have been removed.
-\item \mdl{step}:
-  The main time stepping routine now only needs to call the sbc routine (and a few utility functions).
-\item \mdl{sbcmod}:
-  This has been cut down and now only calculates surface forcing and the ice model required.
+\item \mdl{step}: The main time stepping routine now only needs to call the sbc routine (and a few utility functions).
+\item \mdl{sbcmod}: This has been cut down and now only calculates surface forcing and the ice model required.
   New surface modules that can function when only the surface level of the ocean state is defined can also be added
   (\eg\ icebergs).
-\item \mdl{daymod}:
-  No ocean restarts are read or written (though the ice model restarts are retained),
+\item \mdl{daymod}: No ocean restarts are read or written (though the ice model restarts are retained),
   so calls to restart functions have been removed.
   This also means that the calendar cannot be controlled by time in a restart file,
   so the user must check that nn\_date0 in the model namelist is correct for his or her purposes.
-\item \mdl{stpctl}:
-  Since there is no free surface solver, references to it have been removed from \rou{stp\_ctl} module.
-\item \mdl{diawri}:
-  All 3D data have been removed from the output.
+\item \mdl{stpctl}: Since there is no free surface solver, references to it have been removed from \rou{stp\_ctl} module.
+\item \mdl{diawri}: All 3D data have been removed from the output.
   The surface temperature, salinity and velocity components (which have been read in) are written along with
   relevant forcing and ice data.
@@ -502,6 +496,5 @@
 
 \begin{itemize}
-\item \mdl{sbcsas}:
-  This module initialises the input files needed for reading temperature, salinity and
+\item \mdl{sbcsas}: This module initialises the input files needed for reading temperature, salinity and
   velocity arrays at the surface.
   These filenames are supplied in namelist namsbc\_sas.
@@ -621,6 +614,5 @@
 their neutral transfer coefficients relationships with neutral wind.
 \begin{itemize}
-\item NCAR (\np[=.true.]{ln_NCAR}{ln\_NCAR}):
-  The NCAR bulk formulae have been developed by \citet{large.yeager_rpt04}.
+\item NCAR (\np[=.true.]{ln_NCAR}{ln\_NCAR}): The NCAR bulk formulae have been developed by \citet{large.yeager_rpt04}.
   They have been designed to handle the NCAR forcing, a mixture of NCEP reanalysis and satellite data.
   They use an inertial dissipative method to compute the turbulent transfer coefficients
@@ -630,10 +622,7 @@
   Note that substituting ERA40 to NCEP reanalysis fields does not require changes in the bulk formulea themself.
   This is the so-called DRAKKAR Forcing Set (DFS) \citep{brodeau.barnier.ea_OM10}.
-\item COARE 3.0 (\np[=.true.]{ln_COARE_3p0}{ln\_COARE\_3p0}):
-  See \citet{fairall.bradley.ea_JC03} for more details
-\item COARE 3.5 (\np[=.true.]{ln_COARE_3p5}{ln\_COARE\_3p5}):
-  See \citet{edson.jampana.ea_JPO13} for more details
-\item ECMWF (\np[=.true.]{ln_ECMWF}{ln\_ECMWF}):
-  Based on \href{https://www.ecmwf.int/node/9221}{IFS (Cy31)} implementation and documentation.
+\item COARE 3.0 (\np[=.true.]{ln_COARE_3p0}{ln\_COARE\_3p0}): See \citet{fairall.bradley.ea_JC03} for more details
+\item COARE 3.5 (\np[=.true.]{ln_COARE_3p5}{ln\_COARE\_3p5}): See \citet{edson.jampana.ea_JPO13} for more details
+\item ECMWF (\np[=.true.]{ln_ECMWF}{ln\_ECMWF}): Based on \href{https://www.ecmwf.int/node/9221}{IFS (Cy31)} implementation and documentation.
   Surface roughness lengths needed for the Obukhov length are computed following \citet{beljaars_QJRMS95}.
 \end{itemize}
@@ -741,5 +730,4 @@
 \label{sec:SBC_tide}
 
-
 \begin{listing}
   \nlst{nam_tide}
@@ -924,7 +912,5 @@
 
 \begin{description}
-
-  \item [{\np[=1]{nn_isf}{nn\_isf}}]:
-  The ice shelf cavity is represented (\np[=.true.]{ln_isfcav}{ln\_isfcav} needed).
+  \item [{\np[=1]{nn_isf}{nn\_isf}}]: The ice shelf cavity is represented (\np[=.true.]{ln_isfcav}{ln\_isfcav} needed).
   The fwf and heat flux are depending of the local water properties.
 
@@ -932,9 +918,7 @@
 
 	\begin{description}
-	\item [{\np[=1]{nn_isfblk}{nn\_isfblk}}]:
-     The melt rate is based on a balance between the upward ocean heat flux and
+	\item [{\np[=1]{nn_isfblk}{nn\_isfblk}}]: The melt rate is based on a balance between the upward ocean heat flux and
      the latent heat flux at the ice shelf base. A complete description is available in \citet{hunter_rpt06}.
-	\item [{\np[=2]{nn_isfblk}{nn\_isfblk}}]:
-     The melt rate and the heat flux are based on a 3 equations formulation
+	\item [{\np[=2]{nn_isfblk}{nn\_isfblk}}]: The melt rate and the heat flux are based on a 3 equations formulation
      (a heat flux budget at the ice base, a salt flux budget at the ice base and a linearised freezing point temperature equation).
      A complete description is available in \citet{jenkins_JGR91}.
@@ -952,6 +936,5 @@
      There are 3 different ways to compute the exchange coeficient:
 	\begin{description}
-        \item [{\np[=0]{nn_gammablk}{nn\_gammablk}}]:
-     The salt and heat exchange coefficients are constant and defined by \np{rn_gammas0}{rn\_gammas0} and \np{rn_gammat0}{rn\_gammat0}.
+        \item [{\np[=0]{nn_gammablk}{nn\_gammablk}}]: The salt and heat exchange coefficients are constant and defined by \np{rn_gammas0}{rn\_gammas0} and \np{rn_gammat0}{rn\_gammat0}.
      \begin{gather*}
        % \label{eq:SBC_isf_gamma_iso}
@@ -960,6 +943,5 @@
      \end{gather*}
      This is the recommended formulation for ISOMIP.
-	\item [{\np[=1]{nn_gammablk}{nn\_gammablk}}]:
-     The salt and heat exchange coefficients are velocity dependent and defined as
+	\item [{\np[=1]{nn_gammablk}{nn\_gammablk}}]: The salt and heat exchange coefficients are velocity dependent and defined as
      \begin{gather*}
        \gamma^{T} = rn\_gammat0 \times u_{*} \\
@@ -968,6 +950,5 @@
      where $u_{*}$ is the friction velocity in the top boundary layer (ie first \np{rn_hisf_tbl}{rn\_hisf\_tbl} meters).
      See \citet{jenkins.nicholls.ea_JPO10} for all the details on this formulation. It is the recommended formulation for realistic application.
-	\item [{\np[=2]{nn_gammablk}{nn\_gammablk}}]:
-     The salt and heat exchange coefficients are velocity and stability dependent and defined as:
+	\item [{\np[=2]{nn_gammablk}{nn\_gammablk}}]: The salt and heat exchange coefficients are velocity and stability dependent and defined as:
 \[
 \gamma^{T,S} = \frac{u_{*}}{\Gamma_{Turb} + \Gamma^{T,S}_{Mole}}
@@ -979,6 +960,5 @@
      This formulation has not been extensively tested in \NEMO\ (not recommended).
    \end{description}
-  \item [{\np[=2]{nn_isf}{nn\_isf}}]:
-   The ice shelf cavity is not represented.
+  \item [{\np[=2]{nn_isf}{nn\_isf}}]: The ice shelf cavity is not represented.
    The fwf and heat flux are computed using the \citet{beckmann.goosse_OM03} parameterisation of isf melting.
    The fluxes are distributed along the ice shelf edge between the depth of the average grounding line (GL)
@@ -986,12 +966,10 @@
    (\np{sn_depmin_isf}{sn\_depmin\_isf}) as in (\np[=3]{nn_isf}{nn\_isf}).
    The effective melting length (\np{sn_Leff_isf}{sn\_Leff\_isf}) is read from a file.
-  \item [{\np[=3]{nn_isf}{nn\_isf}}]:
-   The ice shelf cavity is not represented.
+  \item [{\np[=3]{nn_isf}{nn\_isf}}]: The ice shelf cavity is not represented.
    The fwf (\np{sn_rnfisf}{sn\_rnfisf}) is prescribed and distributed along the ice shelf edge between
    the depth of the average grounding line (GL) (\np{sn_depmax_isf}{sn\_depmax\_isf}) and
    the base of the ice shelf along the calving front (\np{sn_depmin_isf}{sn\_depmin\_isf}).
    The heat flux ($Q_h$) is computed as $Q_h = fwf \times L_f$.
-  \item [{\np[=4]{nn_isf}{nn\_isf}}]:
-   The ice shelf cavity is opened (\np[=.true.]{ln_isfcav}{ln\_isfcav} needed).
+  \item [{\np[=4]{nn_isf}{nn\_isf}}]: The ice shelf cavity is opened (\np[=.true.]{ln_isfcav}{ln\_isfcav} needed).
    However, the fwf is not computed but specified from file \np{sn_fwfisf}{sn\_fwfisf}).
    The heat flux ($Q_h$) is computed as $Q_h = fwf \times L_f$.
@@ -1037,11 +1015,11 @@
 At each restart step:
 
-\begin{description}
-\item [Step 1]: the ice sheet model send a new bathymetry and ice shelf draft netcdf file.
-\item [Step 2]: a new domcfg.nc file is built using the DOMAINcfg tools.
-\item [Step 3]: \NEMO\ run for a specific period and output the average melt rate over the period.
-\item [Step 4]: the ice sheet model run using the melt rate outputed in step 4.
-\item [Step 5]: go back to 1.
-\end{description}
+\begin{enumerate}
+\item the ice sheet model send a new bathymetry and ice shelf draft netcdf file.
+\item a new domcfg.nc file is built using the DOMAINcfg tools.
+\item \NEMO\ run for a specific period and output the average melt rate over the period.
+\item the ice sheet model run using the melt rate outputed in step 4.
+\item go back to 1.
+\end{enumerate}
 
 If \np[=.true.]{ln_iscpl}{ln\_iscpl}, the isf draft is assume to be different at each restart step with
@@ -1050,20 +1028,14 @@
 
 \begin{description}
-\item [Thin a cell down]:
-  T/S/ssh are unchanged and U/V in the top cell are corrected to keep the barotropic transport (bt) constant
+\item [Thin a cell down]: T/S/ssh are unchanged and U/V in the top cell are corrected to keep the barotropic transport (bt) constant
   ($bt_b=bt_n$).
-\item [Enlarge  a cell]:
-  See case "Thin a cell down"
-\item [Dry a cell]:
-  mask, T/S, U/V and ssh are set to 0.
+\item [Enlarge  a cell]: See case "Thin a cell down"
+\item [Dry a cell]: mask, T/S, U/V and ssh are set to 0.
   Furthermore, U/V into the water column are modified to satisfy ($bt_b=bt_n$).
-\item [Wet a cell]:
-  mask is set to 1, T/S is extrapolated from neighbours, $ssh_n = ssh_b$ and U/V set to 0.
+\item [Wet a cell]: mask is set to 1, T/S is extrapolated from neighbours, $ssh_n = ssh_b$ and U/V set to 0.
   If no neighbours, T/S is extrapolated from old top cell value.
   If no neighbours along i,j and k (both previous test failed), T/S/U/V/ssh and mask are set to 0.
-\item [Dry a column]:
-	mask, T/S, U/V are set to 0 everywhere in the column and ssh set to 0.
-\item [Wet a column]:
-  set mask to 1, T/S is extrapolated from neighbours, ssh is extrapolated from neighbours and U/V set to 0.
+\item [Dry a column]: mask, T/S, U/V are set to 0 everywhere in the column and ssh set to 0.
+\item [Wet a column]: set mask to 1, T/S is extrapolated from neighbours, ssh is extrapolated from neighbours and U/V set to 0.
   If no neighbour, T/S/U/V and mask set to 0.
 \end{description}
@@ -1109,6 +1081,5 @@
 Two initialisation schemes are possible.
 \begin{description}
-\item [{\np{nn_test_icebergs}{nn\_test\_icebergs}~$>$~0}]
-  In this scheme, the value of \np{nn_test_icebergs}{nn\_test\_icebergs} represents the class of iceberg to generate
+\item [{\np{nn_test_icebergs}{nn\_test\_icebergs}~$>$~0}] In this scheme, the value of \np{nn_test_icebergs}{nn\_test\_icebergs} represents the class of iceberg to generate
   (so between 1 and 10), and \np{nn_test_icebergs}{nn\_test\_icebergs} provides a lon/lat box in the domain at each grid point of
   which an iceberg is generated at the beginning of the run.
@@ -1116,6 +1087,5 @@
   \np{nn_test_icebergs}{nn\_test\_icebergs} is defined by four numbers in \np{nn_test_box}{nn\_test\_box} representing the corners of
   the geographical box: lonmin,lonmax,latmin,latmax
-\item [{\np[=-1]{nn_test_icebergs}{nn\_test\_icebergs}}]
-  In this scheme, the model reads a calving file supplied in the \np{sn_icb}{sn\_icb} parameter.
+\item [{\np[=-1]{nn_test_icebergs}{nn\_test\_icebergs}}] In this scheme, the model reads a calving file supplied in the \np{sn_icb}{sn\_icb} parameter.
   This should be a file with a field on the configuration grid (typically ORCA)
   representing ice accumulation rate at each model point.
@@ -1184,8 +1154,4 @@
 \end{description}
 
-% ----------------------------------------------------------------
-% Neutral drag coefficient from wave model (ln_cdgw)
-
-% ----------------------------------------------------------------
 %% =================================================================================================
 \subsection[Neutral drag coefficient from wave model (\forcode{ln_cdgw})]{Neutral drag coefficient from wave model (\protect\np{ln_cdgw}{ln\_cdgw})}
@@ -1198,7 +1164,4 @@
 air-sea interface following \citet{large.yeager_rpt04}.
 
-% ----------------------------------------------------------------
-% 3D Stokes Drift (ln_sdw, nn_sdrift)
-% ----------------------------------------------------------------
 %% =================================================================================================
 \subsection[3D Stokes Drift (\forcode{ln_sdw} \& \forcode{nn_sdrift})]{3D Stokes Drift (\protect\np{ln_sdw}{ln\_sdw} \& \np{nn_sdrift}{nn\_sdrift})}
@@ -1294,7 +1257,4 @@
 \]
 
-% ----------------------------------------------------------------
-% Stokes-Coriolis term (ln_stcor)
-% ----------------------------------------------------------------
 %% =================================================================================================
 \subsection[Stokes-Coriolis term (\forcode{ln_stcor})]{Stokes-Coriolis term (\protect\np{ln_stcor}{ln\_stcor})}
@@ -1308,7 +1268,4 @@
 \np[=.true.]{ln_stcor}{ln\_stcor} has to be set.
 
-% ----------------------------------------------------------------
-% Waves modified stress (ln_tauwoc, ln_tauw)
-% ----------------------------------------------------------------
 %% =================================================================================================
 \subsection[Wave modified stress (\forcode{ln_tauwoc} \& \forcode{ln_tauw})]{Wave modified sress (\protect\np{ln_tauwoc}{ln\_tauwoc} \& \np{ln_tauw}{ln\_tauw})}
@@ -1357,5 +1314,4 @@
 \label{subsec:SBC_dcy}
 %
-
 
 \begin{figure}[!t]
@@ -1475,11 +1431,9 @@
 the value of the \np{nn_ice}{nn\_ice} namelist parameter found in \nam{sbc}{sbc} namelist.
 \begin{description}
-\item [nn\_ice = 0]
-  there will never be sea-ice in the computational domain.
+\item [nn\_ice = 0] there will never be sea-ice in the computational domain.
   This is a typical namelist value used for tropical ocean domain.
   The surface fluxes are simply specified for an ice-free ocean.
   No specific things is done for sea-ice.
-\item [nn\_ice = 1]
-  sea-ice can exist in the computational domain, but no sea-ice model is used.
+\item [nn\_ice = 1] sea-ice can exist in the computational domain, but no sea-ice model is used.
   An observed ice covered area is read in a file.
   Below this area, the SST is restored to the freezing point and
@@ -1492,6 +1446,5 @@
   is usually referred as the \textit{ice-if} model.
   It can be found in the \mdl{sbcice\_if} module.
-\item [nn\_ice = 2 or more]
-  A full sea ice model is used.
+\item [nn\_ice = 2 or more] A full sea ice model is used.
   This model computes the ice-ocean fluxes,
   that are combined with the air-sea fluxes using the ice fraction of each model cell to
@@ -1545,13 +1498,10 @@
 
 \begin{description}
-\item [{\np[=0]{nn_fwb}{nn\_fwb}}]
-  no control at all.
+\item [{\np[=0]{nn_fwb}{nn\_fwb}}] no control at all.
   The mean sea level is free to drift, and will certainly do so.
-\item [{\np[=1]{nn_fwb}{nn\_fwb}}]
-  global mean \textit{emp} set to zero at each model time step.
+\item [{\np[=1]{nn_fwb}{nn\_fwb}}] global mean \textit{emp} set to zero at each model time step.
   %GS: comment below still relevant ?
   %Note that with a sea-ice model, this technique only controls the mean sea level with linear free surface and no mass flux between ocean and ice (as it is implemented in the current ice-ocean coupling).
-\item [{\np[=2]{nn_fwb}{nn\_fwb}}]
-  freshwater budget is adjusted from the previous year annual mean budget which
+\item [{\np[=2]{nn_fwb}{nn\_fwb}}] freshwater budget is adjusted from the previous year annual mean budget which
   is read in the \textit{EMPave\_old.dat} file.
   As the model uses the Boussinesq approximation, the annual mean fresh water budget is simply evaluated from
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_STO.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_STO.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_STO.tex	(revision 11598)
@@ -2,8 +2,24 @@
 
 \begin{document}
+
 \chapter{Stochastic Parametrization of EOS (STO)}
 \label{chap:STO}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
 
 % \vfill
@@ -18,7 +34,5 @@
 % \end{figure}
 
-Authors: \\
-C. Levy release 4.0.1 update \\
-P.-A. Bouttier release 3.6 inital version
+\clearpage
 
 As a result of the nonlinearity of the seawater equation of state, unresolved scales represent a major source of uncertainties in the computation of the large-scale horizontal density gradient from the large-scale temperature and salinity fields. Following  \cite{brankart_OM13}, the impact of these uncertainties can be simulated by random processes representing unresolved T/S fluctuations. The Stochastic Parametrization of EOS (STO) module implements this parametrization.
@@ -178,11 +192,13 @@
 
 \begin{description}
-\item [{\np{nn_sto_eos}{nn\_sto\_eos}:}]   number of independent random walks
-\item [{\np{rn_eos_stdxy}{rn\_eos\_stdxy}:}] random walk horizontal standard deviation (in grid points)
-\item [{\np{rn_eos_stdz}{rn\_eos\_stdz}:}]  random walk vertical standard deviation (in grid points)
-\item [{\np{rn_eos_tcor}{rn\_eos\_tcor}:}]  random walk time correlation (in timesteps)
-\item [{\np{nn_eos_ord}{nn\_eos\_ord}:}]   order of autoregressive processes
-\item [{\np{nn_eos_flt}{nn\_eos\_flt}:}]   passes of Laplacian filter
-\item [{\np{rn_eos_lim}{rn\_eos\_lim}:}]   limitation factor (default = 3.0)
+\item [{\np{nn_sto_eos}{nn\_sto\_eos}:}]     number of independent random walks
+\item [{\np{rn_eos_stdxy}{rn\_eos\_stdxy}:}] random walk horizontal standard deviation
+  (in grid points)
+\item [{\np{rn_eos_stdz}{rn\_eos\_stdz}:}]   random walk vertical standard deviation
+  (in grid points)
+\item [{\np{rn_eos_tcor}{rn\_eos\_tcor}:}]   random walk time correlation (in timesteps)
+\item [{\np{nn_eos_ord}{nn\_eos\_ord}:}]     order of autoregressive processes
+\item [{\np{nn_eos_flt}{nn\_eos\_flt}:}]     passes of Laplacian filter
+\item [{\np{rn_eos_lim}{rn\_eos\_lim}:}]     limitation factor (default = 3.0)
 \end{description}
 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_TRA.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_TRA.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_TRA.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Ocean Tracers (TRA)}
 \label{chap:TRA}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 % missing/update
@@ -111,6 +129,5 @@
 
 \begin{description}
-\item [linear free surface:]
-  (\np[=.true.]{ln_linssh}{ln\_linssh})
+\item [linear free surface:] (\np[=.true.]{ln_linssh}{ln\_linssh})
   the first level thickness is constant in time:
   the vertical boundary condition is applied at the fixed surface $z = 0$ rather than on
@@ -119,6 +136,5 @@
   $\tau_w|_{k = 1/2} = T_{k = 1}$, \ie\ the product of surface velocity (at $z = 0$) by
   the first level tracer value.
-\item [non-linear free surface:]
-  (\np[=.false.]{ln_linssh}{ln\_linssh})
+\item [non-linear free surface:] (\np[=.false.]{ln_linssh}{ln\_linssh})
   convergence/divergence in the first ocean level moves the free surface up/down.
   There is no tracer advection through it so that the advective fluxes through the surface are also zero.
@@ -435,14 +451,11 @@
 
 \begin{description}
-\item [{\np[=.true.]{ln_traldf_OFF}{ln\_traldf\_OFF}}]
-  no operator selected, the lateral diffusive tendency will not be applied to the tracer equation.
+\item [{\np[=.true.]{ln_traldf_OFF}{ln\_traldf\_OFF}}] no operator selected, the lateral diffusive tendency will not be applied to the tracer equation.
   This option can be used when the selected advection scheme is diffusive enough (MUSCL scheme for example).
-\item [{\np[=.true.]{ln_traldf_lap}{ln\_traldf\_lap}}]
-  a laplacian operator is selected.
+\item [{\np[=.true.]{ln_traldf_lap}{ln\_traldf\_lap}}] a laplacian operator is selected.
   This harmonic operator takes the following expression:  $\mathcal{L}(T) = \nabla \cdot A_{ht} \; \nabla T $,
   where the gradient operates along the selected direction (see \autoref{subsec:TRA_ldf_dir}),
   and $A_{ht}$ is the eddy diffusivity coefficient expressed in $m^2/s$ (see \autoref{chap:LDF}).
-\item [{\np[=.true.]{ln_traldf_blp}{ln\_traldf\_blp}}]:
-  a bilaplacian operator is selected.
+\item [{\np[=.true.]{ln_traldf_blp}{ln\_traldf\_blp}}] a bilaplacian operator is selected.
   This biharmonic operator takes the following expression:
   $\mathcal{B} = - \mathcal{L}(\mathcal{L}(T)) = - \nabla \cdot b \nabla (\nabla \cdot b \nabla T)$
@@ -597,5 +610,4 @@
 \section[Tracer vertical diffusion (\textit{trazdf.F90})]{Tracer vertical diffusion (\protect\mdl{trazdf})}
 \label{sec:TRA_zdf}
-
 
 Options are defined through the \nam{zdf}{zdf} namelist variables.
@@ -774,14 +786,9 @@
 
 \begin{description}
-\item [{\np[=0]{nn_chldta}{nn\_chldta}}]
-  a constant 0.05 g.Chl/L value everywhere ;
-\item [{\np[=1]{nn_chldta}{nn\_chldta}}]
-  an observed time varying chlorophyll deduced from satellite surface ocean color measurement spread uniformly in
-  the vertical direction;
-\item [{\np[=2]{nn_chldta}{nn\_chldta}}]
-  same as previous case except that a vertical profile of chlorophyl is used.
+\item [{\np[=0]{nn_chldta}{nn\_chldta}}] a constant 0.05 g.Chl/L value everywhere ;
+\item [{\np[=1]{nn_chldta}{nn\_chldta}}] an observed time varying chlorophyll deduced from satellite surface ocean color measurement spread uniformly in the vertical direction;
+\item [{\np[=2]{nn_chldta}{nn\_chldta}}] same as previous case except that a vertical profile of chlorophyl is used.
   Following \cite{morel.berthon_LO89}, the profile is computed from the local surface chlorophyll value;
-\item [{\np[=.true.]{ln_qsr_bio}{ln\_qsr\_bio}}]
-  simulated time varying chlorophyll by TOP biogeochemical model.
+\item [{\np[=.true.]{ln_qsr_bio}{ln\_qsr\_bio}}] simulated time varying chlorophyll by TOP biogeochemical model.
   In this case, the RGB formulation is used to calculate both the phytoplankton light limitation in
   PISCES and the oceanic heating rate.
@@ -1134,6 +1141,5 @@
 
 \begin{description}
-\item [{\np[=.true.]{ln_teos10}{ln\_teos10}}]
-  the polyTEOS10-bsq equation of seawater \citep{roquet.madec.ea_OM15} is used.
+\item [{\np[=.true.]{ln_teos10}{ln\_teos10}}] the polyTEOS10-bsq equation of seawater \citep{roquet.madec.ea_OM15} is used.
   The accuracy of this approximation is comparable to the TEOS-10 rational function approximation,
   but it is optimized for a boussinesq fluid and the polynomial expressions have simpler and
@@ -1153,6 +1159,5 @@
   either computing the air-sea and ice-sea fluxes (forced mode) or
   sending the SST field to the atmosphere (coupled mode).
-\item [{\np[=.true.]{ln_eos80}{ln\_eos80}}]
-  the polyEOS80-bsq equation of seawater is used.
+\item [{\np[=.true.]{ln_eos80}{ln\_eos80}}] the polyEOS80-bsq equation of seawater is used.
   It takes the same polynomial form as the polyTEOS10, but the coefficients have been optimized to
   accurately fit EOS80 (Roquet, personal comm.).
@@ -1165,6 +1170,5 @@
   Nevertheless, a severe assumption is made in order to have a heat content ($C_p T_p$) which
   is conserved by the model: $C_p$ is set to a constant value, the TEOS10 value.
-\item [{\np[=.true.]{ln_seos}{ln\_seos}}]
-  a simplified EOS (S-EOS) inspired by \citet{vallis_bk06} is chosen,
+\item [{\np[=.true.]{ln_seos}{ln\_seos}}] a simplified EOS (S-EOS) inspired by \citet{vallis_bk06} is chosen,
   the coefficients of which has been optimized to fit the behavior of TEOS10
   (Roquet, personal comm.) (see also \citet{roquet.madec.ea_JPO15}).
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex	(revision 11598)
@@ -5,8 +5,26 @@
 
 \begin{document}
+
 \chapter{Vertical Ocean Physics (ZDF)}
 \label{chap:ZDF}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 %gm% Add here a small introduction to ZDF and naming of the different physics (similar to what have been written for TRA and DYN.
@@ -39,5 +57,4 @@
 %and thus of the formulation used (see \autoref{chap:TD}).
 
-
 \begin{listing}
   \nlst{namzdf}
@@ -69,5 +86,4 @@
 \subsection[Richardson number dependent (\forcode{ln_zdfric})]{Richardson number dependent (\protect\np{ln_zdfric}{ln\_zdfric})}
 \label{subsec:ZDF_ric}
-
 
 \begin{listing}
@@ -405,5 +421,4 @@
 \label{subsec:ZDF_gls}
 
-
 \begin{listing}
   \nlst{namzdf_gls}
@@ -859,5 +874,4 @@
   \label{lst:namdrg_bot}
 \end{listing}
-
 
 Options to define the top and bottom friction are defined through the \nam{drg}{drg} namelist variables.
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_cfgs.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_cfgs.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_cfgs.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Configurations}
 \label{chap:CFGS}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 %% =================================================================================================
@@ -18,5 +36,4 @@
 though by no means are all options exercised in the reference configurations.
 Configuration is defined manually through the \nam{cfg}{cfg} namelist variables.
-
 
 \begin{listing}
@@ -48,15 +65,12 @@
 and second, the two components of the velocity are moved on a $T$-point.
 Therefore, defining \key{c1d} changes some things in the code behaviour:
-\begin{description}
-\item [(1)]
-  a simplified \rou{stp} routine is used (\rou{stp\_c1d}, see \mdl{step\_c1d} module) in which
+\begin{enumerate}
+\item a simplified \rou{stp} routine is used (\rou{stp\_c1d}, see \mdl{step\_c1d} module) in which
   both lateral tendancy terms and lateral physics are not called;
-\item [(2)]
-  the vertical velocity is zero
+\item the vertical velocity is zero
   (so far, no attempt at introducing a Ekman pumping velocity has been made);
-\item [(3)]
-  a simplified treatment of the Coriolis term is performed as $U$- and $V$-points are the same
+\item a simplified treatment of the Coriolis term is performed as $U$- and $V$-points are the same
   (see \mdl{dyncor\_c1d}).
-\end{description}
+\end{enumerate}
 All the relevant \textit{\_c1d} modules can be found in the src/OCE/C1D directory of
 the \NEMO\ distribution.
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_conservation.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_conservation.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_conservation.tex	(revision 11598)
@@ -6,5 +6,22 @@
 \label{chap:CONS}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 The continuous equations of motion have many analytic properties.
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_misc.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_misc.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_misc.tex	(revision 11598)
@@ -2,8 +2,26 @@
 
 \begin{document}
+
 \chapter{Miscellaneous Topics}
 \label{chap:MISC}
 
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 %% =================================================================================================
@@ -397,2 +415,5 @@
 increment also applies to the time.step file which is otherwise updated every timestep.
 
+\onlyinsubfile{\input{../../global/epilogue}}
+
+\end{document}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics.tex	(revision 11598)
@@ -6,12 +6,27 @@
 \label{chap:MB}
 
+\thispagestyle{plain}
+
 \chaptertoc
 
-%% =================================================================================================
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
+
 %% =================================================================================================
 \section{Primitive equations}
 \label{sec:MB_PE}
 
-%% =================================================================================================
 %% =================================================================================================
 \subsection{Vector invariant formulation}
@@ -92,5 +107,4 @@
 Their nature and formulation are discussed in \autoref{sec:MB_zdf_ldf} and \autoref{subsec:MB_boundary_condition}.
 
-%% =================================================================================================
 %% =================================================================================================
 \subsection{Boundary conditions}
@@ -166,9 +180,7 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \section{Horizontal pressure gradient}
 \label{sec:MB_hor_pg}
 
-%% =================================================================================================
 %% =================================================================================================
 \subsection{Pressure formulation}
@@ -204,5 +216,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsection{Free surface formulation}
 \label{subsec:MB_free_surface}
@@ -255,9 +266,7 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \section{Curvilinear \textit{z-}coordinate system}
 \label{sec:MB_zco}
 
-%% =================================================================================================
 %% =================================================================================================
 \subsection{Tensorial formalism}
@@ -301,5 +310,4 @@
 \end{equation}
 
-% >>>>>>>>>>>>>>>>>>>>>>>>>>>>
 \begin{figure}[!tb]
   \centering
@@ -344,5 +352,4 @@
 where $q$ is a scalar quantity and $\vect A = (a_1,a_2,a_3)$ a vector in the $(i,j,k)$ coordinates system.
 
-%% =================================================================================================
 %% =================================================================================================
 \subsection{Continuous model equations}
@@ -530,5 +537,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \section{Curvilinear generalised vertical coordinate system}
 \label{sec:MB_gco}
@@ -611,5 +617,4 @@
 %}
 
-%% =================================================================================================
 %% =================================================================================================
 \subsection{\textit{S}-coordinate formulation}
@@ -701,5 +706,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsection{Curvilinear \zstar-coordinate system}
 \label{subsec:MB_zco_star}
@@ -788,9 +792,7 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsection{Curvilinear terrain-following \textit{s}--coordinate}
 \label{subsec:MB_sco}
 
-%% =================================================================================================
 %% =================================================================================================
 \subsubsection{Introduction}
@@ -875,5 +877,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsection{\texorpdfstring{Curvilinear \ztilde-coordinate}{}}
 \label{subsec:MB_zco_tilde}
@@ -884,5 +885,4 @@
 Its use is therefore not recommended.
 
-%% =================================================================================================
 %% =================================================================================================
 \section{Subgrid scale physics}
@@ -909,5 +909,4 @@
 The formulation of these terms and their underlying physics are briefly discussed in the next two subsections.
 
-%% =================================================================================================
 %% =================================================================================================
 \subsection{Vertical subgrid scale physics}
@@ -943,5 +942,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsection{Formulation of the lateral diffusive and viscous operators}
 \label{subsec:MB_ldf}
@@ -998,5 +996,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsubsection{Lateral laplacian tracer diffusive operator}
 
@@ -1040,5 +1037,4 @@
 while in $s$-coordinates $\pd[]{k}$ is replaced by $\pd[]{s}$.
 
-%% =================================================================================================
 %% =================================================================================================
 \subsubsection{Eddy induced velocity}
@@ -1079,5 +1075,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsubsection{Lateral bilaplacian tracer diffusive operator}
 
@@ -1091,5 +1086,4 @@
 the harmonic eddy diffusion coefficient set to the square root of the biharmonic one.
 
-%% =================================================================================================
 %% =================================================================================================
 \subsubsection{Lateral Laplacian momentum diffusive operator}
@@ -1126,5 +1120,4 @@
 
 %% =================================================================================================
-%% =================================================================================================
 \subsubsection{Lateral bilaplacian momentum diffusive operator}
 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics_zstar.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics_zstar.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics_zstar.tex	(revision 11598)
@@ -2,5 +2,26 @@
 
 \begin{document}
+
 \chapter{ essai \zstar \sstar}
+
+\thispagestyle{plain}
+
+\chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
+
 %% =================================================================================================
 \section{Curvilinear \zstar- or \sstar coordinate system}
@@ -126,5 +147,4 @@
 is a multiple of \np{rdtbt}{rdtbt} in the \nam{dom}{dom} namelist (Figure III.3).
 
-%>   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >
 \begin{figure}[!t]
   \centering
@@ -150,5 +170,4 @@
   \label{fig:MBZ_dyn_dynspg_ts}
 \end{figure}
-%>   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >
 
 The split-explicit formulation has a damping effect on external gravity waves,
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_time_domain.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_time_domain.tex	(revision 11597)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_time_domain.tex	(revision 11598)
@@ -5,5 +5,23 @@
 \chapter{Time Domain}
 \label{chap:TD}
+
+\thispagestyle{plain}
+
 \chaptertoc
+
+\paragraph{Changes record} ~\\
+
+{\footnotesize
+  \begin{tabularx}{\textwidth}{l||X|X}
+    Release & Author(s) & Modifications \\
+    \hline
+    {\em   4.0} & {\em ...} & {\em ...} \\
+    {\em   3.6} & {\em ...} & {\em ...} \\
+    {\em   3.4} & {\em ...} & {\em ...} \\
+    {\em <=3.4} & {\em ...} & {\em ...}
+  \end{tabularx}
+}
+
+\clearpage
 
 % Missing things:
@@ -13,5 +31,4 @@
   would help  ==> to be added}
 %%%%
-
 
 Having defined the continuous equations in \autoref{chap:MB}, we need now to choose a time discretization,
@@ -291,6 +308,4 @@
 \gmcomment{       % add a subsection here
 
-%        Time Domain
-% ------------------------------------------------------------------------------------------------------------
 %% =================================================================================================
 \subsection{Time domain}
@@ -305,6 +320,4 @@
 
 }     	%% end add
-
-
 
 %%
Index: /NEMO/trunk/doc/tools/shr_func.sh
===================================================================
--- /NEMO/trunk/doc/tools/shr_func.sh	(revision 11597)
+++ /NEMO/trunk/doc/tools/shr_func.sh	(revision 11598)
@@ -2,11 +2,6 @@
 
 clean() {
-    ## Not sure if this step is needed, guess latexmk should be able to detect a change
     printf "\t¤ Clean previous build"
-    find latex/$1/build -mindepth 1 -prune -not -name $1_manual.pyg -exec rm -rf {} \;
-
-    ## HTML exports
-    #printf '   - possible HTML export'
-    #find latex/$1 -type d -name 'html*'    -exec rm -r {} \;
+    find latex/$1/build -mindepth 1 -delete
 
     echo
@@ -15,9 +10,7 @@
 build() {
     printf "\t¤ Generation of the PDF format\n"
-    latexmk -r  ./latex/global/latexmkrc  \
-	    -cd ./latex/$1/main/$1_manual \
-	    1> /dev/null
+    latexmk -r  ./latex/global/latexmk.pl -pdfxe ./latex/$1/main/$1_manual \
+#	1> /dev/null
     [ -f ./latex/$1/build/$1_manual.pdf ] && mv ./latex/$1/build/$1_manual.pdf .
     echo
 }
-
