Index: /NEMO/trunk/doc/latex/NEMO/subfiles/annex_A.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/annex_A.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/annex_A.tex	(revision 10406)
@@ -259,10 +259,10 @@
 Applying the time derivative chain rule (first equation of (\autoref{apdx:A_s_chain_rule})) to $u$ and
 using (\autoref{apdx:A_w_in_s}) provides the expression of the last term of the right hand side,
-\begin{equation*} {\begin{array}{*{20}l} 
+\[ {\begin{array}{*{20}l} 
 w_s  \;\frac{\partial u}{\partial s} 
 	= \frac{\partial s}{\partial t} \;  \frac{\partial u }{\partial s}
 	= \left. {\frac{\partial u }{\partial t}} \right|_s  - \left. {\frac{\partial u }{\partial t}} \right|_z \quad , 
 \end{array} } 		
-\end{equation*}
+\]
 leads to the $s-$coordinate formulation of the total $z-$coordinate time derivative,
 $i.e.$ the total $s-$coordinate time derivative :
@@ -370,22 +370,22 @@
 
 The horizontal pressure gradient term can be transformed as follows:
-\begin{equation*}
+\[
 \begin{split}
- -\frac{1}{\rho _o \, e_1 }\left. {\frac{\partial p}{\partial i}} \right|_z
- & =-\frac{1}{\rho _o e_1 }\left[ {\left. {\frac{\partial p}{\partial i}} \right|_s -\frac{e_1 }{e_3 }\sigma _1 \frac{\partial p}{\partial s}} \right] \\
-& =-\frac{1}{\rho _o \,e_1 }\left. {\frac{\partial p}{\partial i}} \right|_s +\frac{\sigma _1 }{\rho _o \,e_3 }\left( {-g\;\rho \;e_3 } \right) \\
-&=-\frac{1}{\rho _o \,e_1 }\left. {\frac{\partial p}{\partial i}} \right|_s -\frac{g\;\rho }{\rho _o }\sigma _1
+ -\frac{1}{\rho_o \, e_1 }\left. {\frac{\partial p}{\partial i}} \right|_z
+ & =-\frac{1}{\rho_o e_1 }\left[ {\left. {\frac{\partial p}{\partial i}} \right|_s -\frac{e_1 }{e_3 }\sigma _1 \frac{\partial p}{\partial s}} \right] \\
+& =-\frac{1}{\rho_o \,e_1 }\left. {\frac{\partial p}{\partial i}} \right|_s +\frac{\sigma _1 }{\rho_o \,e_3 }\left( {-g\;\rho \;e_3 } \right) \\
+&=-\frac{1}{\rho_o \,e_1 }\left. {\frac{\partial p}{\partial i}} \right|_s -\frac{g\;\rho }{\rho_o }\sigma _1
 \end{split}
-\end{equation*}
+\]
 Applying similar manipulation to the second component and
 replacing $\sigma _1$ and $\sigma _2$ by their expression \autoref{apdx:A_s_slope}, it comes:
 \begin{equation} \label{apdx:A_grad_p_1}
 \begin{split}
- -\frac{1}{\rho _o \, e_1 } \left. {\frac{\partial p}{\partial i}} \right|_z
-&=-\frac{1}{\rho _o \,e_1 } \left(     \left.              {\frac{\partial p}{\partial i}} \right|_s 
+ -\frac{1}{\rho_o \, e_1 } \left. {\frac{\partial p}{\partial i}} \right|_z
+&=-\frac{1}{\rho_o \,e_1 } \left(     \left.              {\frac{\partial p}{\partial i}} \right|_s 
                                                   + g\;\rho  \;\left. {\frac{\partial z}{\partial i}} \right|_s    \right) \\
 %
- -\frac{1}{\rho _o \, e_2 }\left. {\frac{\partial p}{\partial j}} \right|_z
-&=-\frac{1}{\rho _o \,e_2 } \left(    \left.               {\frac{\partial p}{\partial j}} \right|_s 
+ -\frac{1}{\rho_o \, e_2 }\left. {\frac{\partial p}{\partial j}} \right|_z
+&=-\frac{1}{\rho_o \,e_2 } \left(    \left.               {\frac{\partial p}{\partial j}} \right|_s 
                                                    + g\;\rho \;\left. {\frac{\partial z}{\partial j}} \right|_s   \right) \\
 \end{split}
@@ -400,10 +400,10 @@
 and a hydrostatic pressure anomaly, $p_h'$, by $p_h'= g \; \int_z^\eta d \; e_3 \; dk$.
 The pressure is then given by:
-\begin{equation*} 
+\[ 
 \begin{split}
 p &= g\; \int_z^\eta \rho \; e_3 \; dk = g\; \int_z^\eta \left(  \rho_o \, d + 1 \right) \; e_3 \; dk   \\
    &= g \, \rho_o \; \int_z^\eta d \; e_3 \; dk + g \, \int_z^\eta e_3 \; dk    
 \end{split}
-\end{equation*}
+\]
 Therefore, $p$ and $p_h'$ are linked through:
 \begin{equation} \label{apdx:A_pressure}
@@ -411,7 +411,7 @@
 \end{equation}
 and the hydrostatic pressure balance expressed in terms of $p_h'$ and $d$ is:
-\begin{equation*} 
+\[ 
 \frac{\partial p_h'}{\partial k} = - d \, g \, e_3
-\end{equation*}
+\]
 
 Substituing \autoref{apdx:A_pressure} in \autoref{apdx:A_grad_p_1} and
@@ -419,9 +419,9 @@
 \begin{equation} \label{apdx:A_grad_p_2}
 \begin{split}
- -\frac{1}{\rho _o \, e_1 } \left. {\frac{\partial p}{\partial i}} \right|_z
+ -\frac{1}{\rho_o \, e_1 } \left. {\frac{\partial p}{\partial i}} \right|_z
 &=-\frac{1}{e_1 } \left(     \left.              {\frac{\partial p_h'}{\partial i}} \right|_s 
                                        + g\; d  \;\left. {\frac{\partial z}{\partial i}} \right|_s    \right)  - \frac{g}{e_1 } \frac{\partial \eta}{\partial i} \\
 %
- -\frac{1}{\rho _o \, e_2 }\left. {\frac{\partial p}{\partial j}} \right|_z
+ -\frac{1}{\rho_o \, e_2 }\left. {\frac{\partial p}{\partial j}} \right|_z
 &=-\frac{1}{e_2 } \left(    \left.               {\frac{\partial p_h'}{\partial j}} \right|_s 
                                         + g\; d \;\left. {\frac{\partial z}{\partial j}} \right|_s   \right)  - \frac{g}{e_2 } \frac{\partial \eta}{\partial j}\\
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/annex_B.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/annex_B.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/annex_B.tex	(revision 10406)
@@ -20,5 +20,5 @@
 \subsubsection*{In z-coordinates}
 In $z$-coordinates, the horizontal/vertical second order tracer diffusion operator is given by:
-\begin{eqnarray} \label{apdx:B1}
+\begin{align} \label{apdx:B1}
  &D^T = \frac{1}{e_1 \, e_2}      \left[
   \left. \frac{\partial}{\partial i} \left( 	\frac{e_2}{e_1}A^{lT} \;\left. \frac{\partial T}{\partial i} \right|_z   \right)   \right|_z      \right.
@@ -26,5 +26,5 @@
 + \left. \frac{\partial}{\partial j} \left(  \frac{e_1}{e_2}A^{lT} \;\left. \frac{\partial T}{\partial j} \right|_z   \right)   \right|_z      \right]
 + \frac{\partial }{\partial z}\left( {A^{vT} \;\frac{\partial T}{\partial z}} \right)
-\end{eqnarray}
+\end{align}
 
 \subsubsection*{In generalized vertical coordinates}
@@ -156,10 +156,10 @@
 \end{equation}
 where ($a_1$, $a_2$) are the isopycnal slopes in ($\textbf{i}$, $\textbf{j}$) directions, relative to geopotentials:
-\begin{equation*}
+\[
 a_1 =\frac{e_3 }{e_1 }\left( {\frac{\partial \rho }{\partial i}} \right)\left( {\frac{\partial \rho }{\partial k}} \right)^{-1}
 \qquad , \qquad
 a_2 =\frac{e_3 }{e_2 }\left( {\frac{\partial \rho }{\partial j}}
 \right)\left( {\frac{\partial \rho }{\partial k}} \right)^{-1}
-\end{equation*}
+\]
 
 In practice, isopycnal slopes are generally less than $10^{-2}$ in the ocean,
@@ -204,8 +204,8 @@
 The demonstration of the first property is trivial as \autoref{apdx:B4} is the divergence of fluxes.
 Let us demonstrate the second one:
-\begin{equation*}
+\[
 \iiint\limits_D T\;\nabla .\left( {\textbf{A}}_{\textbf{I}} \nabla T \right)\,dv
           = -\iiint\limits_D \nabla T\;.\left( {\textbf{A}}_{\textbf{I}} \nabla T \right)\,dv,
-\end{equation*}
+\]
 and since
 \begin{subequations}
@@ -249,10 +249,10 @@
 where ($r_1$, $r_2$) are the isopycnal slopes in ($\textbf{i}$, $\textbf{j}$) directions,
 relative to $s$-coordinate surfaces:
-\begin{equation*}
+\[
 r_1 =\frac{e_3 }{e_1 }\left( {\frac{\partial \rho }{\partial i}} \right)\left( {\frac{\partial \rho }{\partial s}} \right)^{-1}
 \qquad , \qquad
 r_2 =\frac{e_3 }{e_2 }\left( {\frac{\partial \rho }{\partial j}}
 \right)\left( {\frac{\partial \rho }{\partial s}} \right)^{-1}.
-\end{equation*}
+\]
 
 To prove \autoref{apdx:B5} by direct re-expression of \autoref{apdx:B_ldfiso} is straightforward, but laborious.
@@ -325,7 +325,7 @@
 Using \autoref{eq:PE_div}, the definition of the horizontal divergence,
 the third componant of the second vector is obviously zero and thus :
-\begin{equation*}
+\[
 \Delta {\textbf{U}}_h = \nabla _h \left( \chi \right) - \nabla _h \times \left( \zeta \right) + \frac {1}{e_3 } \frac {\partial }{\partial k} \left( {\frac {1}{e_3 } \frac{\partial {\textbf{ U}}_h }{\partial k}} \right)
-\end{equation*}
+\]
 
 Note that this operator ensures a full separation between
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/annex_C.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/annex_C.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/annex_C.tex	(revision 10406)
@@ -25,12 +25,12 @@
 
 fluxes at the faces of a $T$-box:
-\begin{equation*}
+\[
 U = e_{2u}\,e_{3u}\; u  \qquad  V = e_{1v}\,e_{3v}\; v  \qquad W = e_{1w}\,e_{2w}\; \omega     \\
-\end{equation*}
+\]
 
 volume of cells at $u$-, $v$-, and $T$-points:
-\begin{equation*}
+\[
 b_u = e_{1u}\,e_{2u}\,e_{3u}  \qquad  b_v = e_{1v}\,e_{2v}\,e_{3v}  \qquad b_t = e_{1t}\,e_{2t}\,e_{3t}     \\
-\end{equation*}
+\]
 
 partial derivative notation: $\partial_\bullet = \frac{\partial}{\partial \bullet}$
@@ -47,17 +47,17 @@
 
 Continuity equation with the above notation:
-\begin{equation*}
+\[
 \frac{1}{e_{3t}} \partial_t (e_{3t})+ \frac{1}{b_t}  \biggl\{  \delta_i [U] + \delta_j [V] + \delta_k [W] \biggr\} = 0
-\end{equation*}
+\]
 
 A quantity, $Q$ is conserved when its domain averaged time change is zero, that is when:
-\begin{equation*}
+\[
 \partial_t \left( \int_D{ Q\;dv } \right) =0
-\end{equation*}
+\]
 Noting that the coordinate system used ....  blah blah
-\begin{equation*}
+\[
 \partial_t \left( \int_D {Q\;dv} \right) =  \int_D { \partial_t \left( e_3 \, Q \right) e_1e_2\;di\,dj\,dk }
                                                        =  \int_D { \frac{1}{e_3} \partial_t \left( e_3 \, Q \right) dv } =0
-\end{equation*}
+\]
 equation of evolution of $Q$ written as
 the time evolution of the vertical content of $Q$ like for tracers, or momentum in flux form,
@@ -162,9 +162,9 @@
 $\ $\newline    % force a new ligne
 The prognostic ocean dynamics equation can be summarized as follows:
-\begin{equation*}
+\[
 \text{NXT} = \dbinom	{\text{VOR} + \text{KEG} + \text {ZAD} }
 						{\text{COR} + \text{ADV}                       }
 			+ \text{HPG} + \text{SPG} + \text{LDF} + \text{ZDF}
-\end{equation*}
+\]
 $\ $\newline    % force a new ligne
 
@@ -365,5 +365,5 @@
 
 For the ENE scheme, the two components of the vorticity term are given by:
-\begin{equation*}
+\[
 - e_3 \, q \;{\textbf{k}}\times {\textbf {U}}_h    \equiv 
 	\left( {{  \begin{array} {*{20}c}
@@ -373,5 +373,5 @@
 		\overline {\, q \ \overline {\left( e_{2u}\,e_{3u}\,u \right)}^{\,j+1/2}} ^{\,i}       \hfill \\
 	\end{array}} }    \right)
-\end{equation*}
+\]
 
 This formulation does not conserve the enstrophy but it does conserve the total kinetic energy.
@@ -471,9 +471,9 @@
 
 The change of Kinetic Energy (KE) due to the vertical advection is exactly balanced by the change of KE due to the horizontal gradient of KE~:
-\begin{equation*}
+\[
     \int_D \textbf{U}_h \cdot \frac{1}{e_3 } \omega \partial_k \textbf{U}_h \;dv
 =  -   \int_D \textbf{U}_h \cdot \nabla_h \left( \frac{1}{2}\;{\textbf{U}_h}^2 \right)\;dv
    +   \frac{1}{2} \int_D {  \frac{{\textbf{U}_h}^2}{e_3} \partial_t ( e_3) \;dv }  \\
-\end{equation*}
+\]
 Indeed, using successively \autoref{eq:DOM_di_adj} ($i.e.$ the skew symmetry property of the $\delta$ operator)
 and the continuity equation, then \autoref{eq:DOM_di_adj} again,
@@ -544,5 +544,5 @@
 For example KE can also be discretized as $1/2\,({\overline u^{\,i}}^2 + {\overline v^{\,j}}^2)$.
 This leads to the following expression for the vertical advection:
-\begin{equation*}
+\[
 \frac{1} {e_3 }\; \omega\; \partial_k \textbf{U}_h
 \equiv \left( {{\begin{array} {*{20}c}
@@ -552,5 +552,5 @@
 \left[ \overline v^{\,j+1/2} \right]}}^{\,j+1/2,k} \hfill \\
 \end{array}} } \right)
-\end{equation*}
+\]
 a formulation that requires an additional horizontal mean in contrast with the one used in NEMO.
 Nine velocity points have to be used instead of 3.
@@ -588,9 +588,9 @@
 the change of KE due to the work of pressure forces is balanced by
 the change of potential energy due to buoyancy forces: 
-\begin{equation*}
+\[
 - \int_D  \left. \nabla p \right|_z \cdot \textbf{U}_h \;dv 
 = - \int_D \nabla \cdot \left( \rho \,\textbf {U} \right) \,g\,z \;dv
   + \int_D g\, \rho \; \partial_t (z)  \;dv
-\end{equation*}
+\]
 
 This property can be satisfied in a discrete sense for both $z$- and $s$-coordinates.
@@ -771,10 +771,10 @@
 This altered Coriolis parameter is discretised at an f-point.
 It is given by:
-\begin{equation*}
+\[
 f+\frac{1} {e_1 e_2 } \left( v \frac{\partial e_2 } {\partial i} - u \frac{\partial e_1 } {\partial j}\right)\;
 \equiv \;
 f+\frac{1} {e_{1f}\,e_{2f}} \left( \overline v^{\,i+1/2} \delta_{i+1/2} \left[ e_{2u} \right] 
                                                -\overline u^{\,j+1/2} \delta_{j+1/2} \left[ e_{1u}  \right] \right)
-\end{equation*}
+\]
 
 Either the ENE or EEN scheme is then applied to obtain the vorticity term in flux form.
@@ -842,5 +842,5 @@
 \end{flalign*}
 Applying similar manipulation applied to the second term of the scalar product leads to:
-\begin{equation*}
+\[
  -  \int_D \textbf{U}_h \cdot     \left(                 {{\begin{array} {*{20}c}
 \nabla \cdot \left( \textbf{U}\,u \right) \hfill \\
@@ -848,5 +848,5 @@
 \equiv + \sum\limits_{i,j,k}  \frac{1}{2}  \left( \overline {u^2}^{\,i} + \overline {v^2}^{\,j} \right)
 	\biggl\{     \left(   \frac{1}{e_{3t}} \frac{\partial e_{3t}}{\partial t}   \right) \; b_t     \biggr\}    
-\end{equation*}
+\]
 which is the discrete form of $ \frac{1}{2} \int_D u \cdot \nabla \cdot \left(   \textbf{U}\,u   \right) \; dv $.
 \autoref{eq:C_ADV_KE_flux} is thus satisfied.
@@ -1032,8 +1032,8 @@
 
 conservation of a tracer, $T$:
-\begin{equation*}
+\[
 \frac{\partial }{\partial t} \left(   \int_D {T\;dv}   \right) 
 =  \int_D { \frac{1}{e_3}\frac{\partial \left( e_3 \, T \right)}{\partial t} \;dv }=0
-\end{equation*}
+\]
 
 conservation of its variance:
@@ -1156,5 +1156,5 @@
 The lateral momentum diffusion term dissipates the horizontal kinetic energy:
 %\begin{flalign*}
-\begin{equation*}
+\[
 \begin{split}
 \int_D \textbf{U}_h \cdot 
@@ -1198,5 +1198,5 @@
 \quad \leq 0       \\
 \end{split}
-\end{equation*}
+\]
 
 % -------------------------------------------------------------------------------------------------------------
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/annex_E.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/annex_E.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/annex_E.tex	(revision 10406)
@@ -26,22 +26,22 @@
 For example, in the $i$-direction:
 \begin{equation} \label{eq:tra_adv_ubs2}
-\tau _u^{ubs} = \left\{	 \begin{aligned}
-  & \tau _u^{cen4} + \frac{1}{12} \,\tau"_i	   & \quad \text{if }\ u_{i+1/2} \geqslant 0 \\
-  & \tau _u^{cen4} - \frac{1}{12} \,\tau"_{i+1} & \quad \text{if }\ u_{i+1/2}       <       0
+\tau_u^{ubs} = \left\{	 \begin{aligned}
+  & \tau_u^{cen4} + \frac{1}{12} \,\tau"_i	   & \quad \text{if }\ u_{i+1/2} \geqslant 0 \\
+  & \tau_u^{cen4} - \frac{1}{12} \,\tau"_{i+1} & \quad \text{if }\ u_{i+1/2}       <       0
   						 \end{aligned}    \right.
 \end{equation}
 or equivalently, the advective flux is
 \begin{equation} \label{eq:tra_adv_ubs2}
-U_{i+1/2} \ \tau _u^{ubs} 
+U_{i+1/2} \ \tau_u^{ubs} 
 =U_{i+1/2} \ \overline{ T_i - \frac{1}{6}\,\tau"_i }^{\,i+1/2}
 - \frac{1}{2}\, |U|_{i+1/2} \;\frac{1}{6} \;\delta_{i+1/2}[\tau"_i]
 \end{equation}
 where $U_{i+1/2} = e_{1u}\,e_{3u}\,u_{i+1/2}$ and
-$\tau "_i =\delta _i \left[ {\delta _{i+1/2} \left[ \tau \right]} \right]$.
+$\tau "_i =\delta_i \left[ {\delta_{i+1/2} \left[ \tau \right]} \right]$.
 By choosing this expression for $\tau "$ we consider a fourth order approximation of $\partial_i^2$ with
 a constant i-grid spacing ($\Delta i=1$).
 
 Alternative choice: introduce the scale factors:  
-$\tau "_i =\frac{e_{1T}}{e_{2T}\,e_{3T}}\delta _i \left[ \frac{e_{2u} e_{3u} }{e_{1u} }\delta _{i+1/2}[\tau] \right]$.
+$\tau "_i =\frac{e_{1T}}{e_{2T}\,e_{3T}}\delta_i \left[ \frac{e_{2u} e_{3u} }{e_{1u} }\delta_{i+1/2}[\tau] \right]$.
 
 
@@ -76,5 +76,5 @@
 
 NB 2: In a forthcoming release four options will be proposed for the vertical component used in the UBS scheme.
-$\tau _w^{ubs}$ will be evaluated using either \textit{(a)} a centered $2^{nd}$ order scheme,
+$\tau_w^{ubs}$ will be evaluated using either \textit{(a)} a centered $2^{nd}$ order scheme,
 or \textit{(b)} a TVD scheme, or \textit{(c)} an interpolation based on conservative parabolic splines following
 \citet{Shchepetkin_McWilliams_OM05} implementation of UBS in ROMS, or \textit{(d)} an UBS.
@@ -83,7 +83,7 @@
 NB 3: It is straight forward to rewrite \autoref{eq:tra_adv_ubs} as follows:
 \begin{equation} \label{eq:tra_adv_ubs2}
-\tau _u^{ubs} = \left\{	 \begin{aligned}
-	& \tau _u^{cen4} + \frac{1}{12} \tau"_i		& \quad \text{if }\ u_{i+1/2} \geqslant 0 \\
-	& \tau _u^{cen4} - \frac{1}{12} \tau"_{i+1}	& \quad \text{if }\ u_{i+1/2}       <       0
+\tau_u^{ubs} = \left\{	 \begin{aligned}
+	& \tau_u^{cen4} + \frac{1}{12} \tau"_i		& \quad \text{if }\ u_{i+1/2} \geqslant 0 \\
+	& \tau_u^{cen4} - \frac{1}{12} \tau"_{i+1}	& \quad \text{if }\ u_{i+1/2}       <       0
   						 \end{aligned}    \right.
 \end{equation}
@@ -91,5 +91,5 @@
 \begin{equation} \label{eq:tra_adv_ubs2}
 \begin{split}
-e_{2u} e_{3u}\,u_{i+1/2} \ \tau _u^{ubs} 
+e_{2u} e_{3u}\,u_{i+1/2} \ \tau_u^{ubs} 
 &= e_{2u} e_{3u}\,u_{i+1/2} \ \overline{ T - \frac{1}{6}\,\tau"_i }^{\,i+1/2} \\
 & - \frac{1}{2} e_{2u} e_{3u}\,|u|_{i+1/2} \;\frac{1}{6} \;\delta_{i+1/2}[\tau"_i]
@@ -107,10 +107,10 @@
 \begin{equation} \label{eq:tra_ldf_lap}
 \begin{split}
-D_T^{lT} =\frac{1}{e_{1T} \; e_{2T}\;  e_{3T} } &\left[ {\quad \delta _i 
-\left[ {A_u^{lT} \frac{e_{2u} e_{3u} }{e_{1u} }\;\delta _{i+1/2} 
+D_T^{lT} =\frac{1}{e_{1T} \; e_{2T}\;  e_{3T} } &\left[ {\quad \delta_i 
+\left[ {A_u^{lT} \frac{e_{2u} e_{3u} }{e_{1u} }\;\delta_{i+1/2} 
 \left[ T \right]} \right]} \right.
 \\
-&\ \left. {+\; \delta _j \left[ 
-{A_v^{lT} \left( {\frac{e_{1v} e_{3v} }{e_{2v} }\;\delta _{j+1/2} \left[ T 
+&\ \left. {+\; \delta_j \left[ 
+{A_v^{lT} \left( {\frac{e_{1v} e_{3v} }{e_{2v} }\;\delta_{j+1/2} \left[ T 
 \right]} \right)} \right]\quad } \right]
 \end{split}
@@ -121,7 +121,7 @@
 \begin{split}
 D_T^{lT} =&-\frac{1}{e_{1T} \; e_{2T}\;  e_{3T}} \\
-& \delta _i \left[  \sqrt{A_u^{lT}}\ \frac{e_{2u}\,e_{3u}}{e_{1u}}\;\delta _{i+1/2} 
+& \delta_i \left[  \sqrt{A_u^{lT}}\ \frac{e_{2u}\,e_{3u}}{e_{1u}}\;\delta_{i+1/2} 
 		  \left[ \frac{1}{e_{1T}\,e_{2T}\, e_{3T}}
-    \delta _i \left[ \sqrt{A_u^{lT}}\ \frac{e_{2u}\,e_{3u}}{e_{1u}}\;\delta _{i+1/2} 
+    \delta_i \left[ \sqrt{A_u^{lT}}\ \frac{e_{2u}\,e_{3u}}{e_{1u}}\;\delta_{i+1/2} 
     		  [T] \right] \right] \right]
 \end{split}
@@ -133,7 +133,7 @@
 \begin{split}
 D_T^{lT} =&-\frac{1}{12}\,\frac{1}{e_{1T} \; e_{2T}\;  e_{3T}} \\
-& \delta _i \left[ e_{2u}\,e_{3u}\,\sqrt{ e_{1u}\,|u|\,}\;\delta _{i+1/2} 
+& \delta_i \left[ e_{2u}\,e_{3u}\,\sqrt{ e_{1u}\,|u|\,}\;\delta_{i+1/2} 
 		 \left[ \frac{1}{e_{1T}\,e_{2T}\, e_{3T}} 
-    \delta _i \left[ e_{2u}\,e_{3u}\,\sqrt{ e_{1u}\,|u|\,}\;\delta _{i+1/2} 
+    \delta_i \left[ e_{2u}\,e_{3u}\,\sqrt{ e_{1u}\,|u|\,}\;\delta_{i+1/2} 
     		[T] \right] \right] \right]
 \end{split}
@@ -143,7 +143,7 @@
 \begin{split}
 F_u^{lT} = - \frac{1}{12}
- e_{2u}\,e_{3u}\,|u| \;\sqrt{ e_{1u}}\,\delta _{i+1/2} 
+ e_{2u}\,e_{3u}\,|u| \;\sqrt{ e_{1u}}\,\delta_{i+1/2} 
 		 \left[ \frac{1}{e_{1T}\,e_{2T}\, e_{3T}} 
-    \delta _i \left[ e_{2u}\,e_{3u}\,\sqrt{ e_{1u}}\:\delta _{i+1/2} 
+    \delta_i \left[ e_{2u}\,e_{3u}\,\sqrt{ e_{1u}}\:\delta_{i+1/2} 
     		[T] \right] \right]
 \end{split}
@@ -158,5 +158,5 @@
 \end{equation}
 if the velocity is uniform ($i.e.$ $|u|=cst$) and
-choosing $\tau "_i =\frac{e_{1T}}{e_{2T}\,e_{3T}}\delta _i \left[ \frac{e_{2u} e_{3u} }{e_{1u} } \delta _{i+1/2}[\tau] \right]$
+choosing $\tau "_i =\frac{e_{1T}}{e_{2T}\,e_{3T}}\delta_i \left[ \frac{e_{2u} e_{3u} }{e_{1u} } \delta_{i+1/2}[\tau] \right]$
 
 sol 1 coefficient at T-point ( add $e_{1u}$ and $e_{1T}$ on both side of first $\delta$):
@@ -164,5 +164,5 @@
 \begin{split}
 F_u^{lT}
-&= - \frac{1}{12} \frac{e_{2u} e_{3u}}{e_{1u}}\;\delta_{i+1/2}\left[ \frac{e_{1T}^3\,|u|}{e_{1T}e_{2T}\,e_{3T}}\,\delta _i \left[ \frac{e_{2u} e_{3u} }{e_{1u} } \delta _{i+1/2}[\tau] \right] \right]
+&= - \frac{1}{12} \frac{e_{2u} e_{3u}}{e_{1u}}\;\delta_{i+1/2}\left[ \frac{e_{1T}^3\,|u|}{e_{1T}e_{2T}\,e_{3T}}\,\delta_i \left[ \frac{e_{2u} e_{3u} }{e_{1u} } \delta_{i+1/2}[\tau] \right] \right]
 \end{split}
 \end{equation}
@@ -173,6 +173,6 @@
 \begin{split}
 F_u^{lT}
-&= - \frac{1}{12} {e_{1u}}^1 \sqrt{e_{1u}|u|} \frac{e_{2u} e_{3u}}{e_{1u}}\;\delta_{i+1/2}\left[ \frac{1}{e_{2T}\,e_{3T}}\,\delta _i \left[ \sqrt{e_{1u}|u|} \frac{e_{2u} e_{3u} }{e_{1u} } \delta _{i+1/2}[\tau] \right] \right] \\
-&= - \frac{1}{12} e_{1u} \sqrt{e_{1u}|u|\,} \frac{e_{2u} e_{3u}}{e_{1u}}\;\delta_{i+1/2}\left[ \frac{1}{e_{1T}\,e_{2T}\,e_{3T}}\,\delta _i \left[ e_{1u} \sqrt{e_{1u}|u|\,} \frac{e_{2u} e_{3u} }{e_{1u}} \delta _{i+1/2}[\tau] \right] \right]
+&= - \frac{1}{12} {e_{1u}}^1 \sqrt{e_{1u}|u|} \frac{e_{2u} e_{3u}}{e_{1u}}\;\delta_{i+1/2}\left[ \frac{1}{e_{2T}\,e_{3T}}\,\delta_i \left[ \sqrt{e_{1u}|u|} \frac{e_{2u} e_{3u} }{e_{1u} } \delta_{i+1/2}[\tau] \right] \right] \\
+&= - \frac{1}{12} e_{1u} \sqrt{e_{1u}|u|\,} \frac{e_{2u} e_{3u}}{e_{1u}}\;\delta_{i+1/2}\left[ \frac{1}{e_{1T}\,e_{2T}\,e_{3T}}\,\delta_i \left[ e_{1u} \sqrt{e_{1u}|u|\,} \frac{e_{2u} e_{3u} }{e_{1u}} \delta_{i+1/2}[\tau] \right] \right]
 \end{split}
 \end{equation}
@@ -191,5 +191,5 @@
 \begin{subequations} \label{eq:dt_mt}
 \begin{align}
- \delta _{t+\rdt/2} [q]     &=  \  \ \,   q^{t+\rdt}  - q^{t}		\\
+ \delta_{t+\rdt/2} [q]     &=  \  \ \,   q^{t+\rdt}  - q^{t}		\\
  \overline q^{\,t+\rdt/2} &= \left\{ q^{t+\rdt} + q^{t} \right\} \; / \; 2
 \end{align}
@@ -202,5 +202,5 @@
 \begin{equation} \label{eq:LF}
    \frac{\partial q}{\partial t} 
-   		\equiv \frac{1}{\rdt} \overline{ \delta _{t+\rdt/2}[q]}^{\,t} 
+   		\equiv \frac{1}{\rdt} \overline{ \delta_{t+\rdt/2}[q]}^{\,t} 
 		=         \frac{q^{t+\rdt}-q^{t-\rdt}}{2\rdt}
 \end{equation} 
@@ -219,9 +219,9 @@
 \int_{t_0}^{t_1} {q\, \frac{\partial q}{\partial t} \;dt} 
 	&\equiv \sum\limits_{0}^{N} 
-			{\frac{1}{\rdt} q^t \ \overline{ \delta _{t+\rdt/2}[q]}^{\,t} \ \rdt} 
-	   \equiv \sum\limits_{0}^{N}  { q^t \ \overline{ \delta _{t+\rdt/2}[q]}^{\,t} } \\
-	&\equiv \sum\limits_{0}^{N}  { \overline{q}^{\,t+\Delta/2}{ \delta _{t+\rdt/2}[q]}}
-	   \equiv \sum\limits_{0}^{N}  { \frac{1}{2} \delta _{t+\rdt/2}[q^2] }\\
-	&\equiv \sum\limits_{0}^{N}  { \frac{1}{2} \delta _{t+\rdt/2}[q^2] }
+			{\frac{1}{\rdt} q^t \ \overline{ \delta_{t+\rdt/2}[q]}^{\,t} \ \rdt} 
+	   \equiv \sum\limits_{0}^{N}  { q^t \ \overline{ \delta_{t+\rdt/2}[q]}^{\,t} } \\
+	&\equiv \sum\limits_{0}^{N}  { \overline{q}^{\,t+\Delta/2}{ \delta_{t+\rdt/2}[q]}}
+	   \equiv \sum\limits_{0}^{N}  { \frac{1}{2} \delta_{t+\rdt/2}[q^2] }\\
+	&\equiv \sum\limits_{0}^{N}  { \frac{1}{2} \delta_{t+\rdt/2}[q^2] }
 	   \equiv \frac{1}{2} \left( {q_{t_1}}^2 - {q_{t_0}}^2 \right)
 \end{split} \end{equation}
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/annex_iso.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/annex_iso.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/annex_iso.tex	(revision 10406)
@@ -155,15 +155,15 @@
 noting that the $e_{{3}_{i+1/2}^k}$ in the area $e{_{3}}_{i+1/2}^k{e_{2}}_{i+1/2}i^k$ at the $u$-point cancels out with
 the $1/{e_{3}}_{i+1/2}^k$ associated with the vertical tracer gradient, is then \autoref{eq:tra_ldf_iso}
-\begin{equation*}
+\[
   \left(F_u^{13} \right)_{i+\hhalf}^k = \Alts_{i+\hhalf}^k
   {e_{2}}_{i+1/2}^k \overline{\overline
     r_1} ^{\,i,k}\,\overline{\overline{\delta_k T}}^{\,i,k},
-\end{equation*}
+\]
 where
-\begin{equation*}
+\[
   \overline{\overline
    r_1} ^{\,i,k} = -\frac{{e_{3u}}_{i+1/2}^k}{{e_{1u}}_{i+1/2}^k}
   \frac{\delta_{i+1/2} [\rho]}{\overline{\overline{\delta_k \rho}}^{\,i,k}},
-\end{equation*}
+\]
 and here and in the following we drop the $^{lT}$ superscript from $\Alt$ for simplicity.
 Unfortunately the resulting combination $\overline{\overline{\delta_k\bullet}}^{\,i,k}$ of a $k$ average and
@@ -200,9 +200,9 @@
   \label{eq:i13}
   \left( F_u^{13}  \right)_{i+\frac{1}{2}}^k = \Alts_{i+1}^k a_1 s_1
-  \delta _{k+\frac{1}{2}} \left[ T^{i+1}
+  \delta_{k+\frac{1}{2}} \left[ T^{i+1}
   \right]/e_{{3w}_{i+1}}^{k+\frac{1}{2}}  + \Alts _i^k a_2 s_2 \delta
   _{k+\frac{1}{2}} \left[ T^i
   \right]/e_{{3w}_{i+1}}^{k+\frac{1}{2}} \\
-   +\Alts _{i+1}^k a_3 s_3 \delta _{k-\frac{1}{2}} \left[ T^{i+1}
+   +\Alts _{i+1}^k a_3 s_3 \delta_{k-\frac{1}{2}} \left[ T^{i+1}
   \right]/e_{{3w}_{i+1}}^{k+\frac{1}{2}}  +\Alts _i^k a_4 s_4 \delta
   _{k-\frac{1}{2}} \left[ T^i \right]/e_{{3w}_{i+1}}^{k+\frac{1}{2}},
@@ -218,8 +218,8 @@
   \label{eq:i31}
   \left( F_w^{31} \right) _i ^{k+\frac{1}{2}} =  \Alts_i^{k+1} a_{1}'
-  s_{1}' \delta _{i-\frac{1}{2}} \left[ T^{k+1} \right]/{e_{3u}}_{i-\frac{1}{2}}^{k+1}
-   +\Alts_i^{k+1} a_{2}' s_{2}' \delta _{i+\frac{1}{2}} \left[ T^{k+1} \right]/{e_{3u}}_{i+\frac{1}{2}}^{k+1}\\
-  + \Alts_i^k a_{3}' s_{3}' \delta _{i-\frac{1}{2}} \left[ T^k\right]/{e_{3u}}_{i-\frac{1}{2}}^k
-  +\Alts_i^k a_{4}' s_{4}' \delta _{i+\frac{1}{2}} \left[ T^k \right]/{e_{3u}}_{i+\frac{1}{2}}^k.
+  s_{1}' \delta_{i-\frac{1}{2}} \left[ T^{k+1} \right]/{e_{3u}}_{i-\frac{1}{2}}^{k+1}
+   +\Alts_i^{k+1} a_{2}' s_{2}' \delta_{i+\frac{1}{2}} \left[ T^{k+1} \right]/{e_{3u}}_{i+\frac{1}{2}}^{k+1}\\
+  + \Alts_i^k a_{3}' s_{3}' \delta_{i-\frac{1}{2}} \left[ T^k\right]/{e_{3u}}_{i-\frac{1}{2}}^k
+  +\Alts_i^k a_{4}' s_{4}' \delta_{i+\frac{1}{2}} \left[ T^k \right]/{e_{3u}}_{i+\frac{1}{2}}^k.
 \end{multline}
 
@@ -275,5 +275,5 @@
   - \left( \Alts_i^{k+1} a_{1} + \Alts_i^{k+1} a_{2} + \Alts_i^k
     a_{3} + \Alts_i^k a_{4} \right)
-  \frac{\delta _{i+1/2} \left[ T^k\right]}{{e_{1u}}_{\,i+1/2}^{\,k}},
+  \frac{\delta_{i+1/2} \left[ T^k\right]}{{e_{1u}}_{\,i+1/2}^{\,k}},
 \end{equation}
 where the areas $a_i$ are as in \autoref{eq:i13}.
@@ -647,7 +647,7 @@
 (here the $\sigma_i$ are the slopes of the coordinate surfaces relative to geopotentials)
 \autoref{eq:PE_slopes_eiv} rather than the slope $r_i$ relative to coordinate surfaces, so we require
-\begin{equation*}
+\[
   |\tilde{r}_i|\leq \tilde{r}_\mathrm{max}=0.01.
-\end{equation*}
+\]
 and then recalculate the slopes $r_i$ relative to coordinates.
 Each individual triad slope
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ASM.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ASM.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ASM.tex	(revision 10406)
@@ -45,8 +45,8 @@
 is corrected by adding the analysis increments for temperature, salinity, horizontal velocity and SSH as
 additional tendency terms to the prognostic equations:
-\begin{eqnarray}     \label{eq:wa_traj_iau}
+\begin{align}     \label{eq:wa_traj_iau}
 {\bf x}^{a}(t_{i}) = M(t_{i}, t_{0})[{\bf x}^{b}(t_{0})] 
 \; + \; F_{i} \delta \tilde{\bf x}^{a} 
-\end{eqnarray}
+\end{align}
 where $F_{i}$ is a weighting function for applying the increments $\delta\tilde{\bf x}^{a}$ defined such that
 $\sum_{i=1}^{N} F_{i}=1$.
@@ -58,5 +58,5 @@
 In addition, two different weighting functions have been implemented.
 The first function employs constant weights, 
-\begin{eqnarray}    \label{eq:F1_i}
+\begin{align}    \label{eq:F1_i}
 F^{(1)}_{i}
 =\left\{ \begin{array}{ll}
@@ -65,9 +65,9 @@
    0     &    {\rm if} \; \; \; t_{i} > t_{n}
   \end{array} \right. 
-\end{eqnarray}
+\end{align}
 where $M = m-n$.
 The second function employs peaked hat-like weights in order to give maximum weight in the centre of the sub-window,
 with the weighting reduced linearly to a small value at the window end-points:
-\begin{eqnarray}     \label{eq:F2_i}
+\begin{align}     \label{eq:F2_i}
 F^{(2)}_{i}
 =\left\{ \begin{array}{ll}
@@ -77,5 +77,5 @@
    0                            &   {\rm if} \; \; \; t_{i}        >    t_{n}
   \end{array} \right.
-\end{eqnarray}
+\end{align}
 where $\alpha^{-1} = \sum_{i=1}^{M/2} 2i$ and $M$ is assumed to be even. 
 The weights described by \autoref{eq:F2_i} provide a smoother transition of the analysis trajectory from
@@ -91,8 +91,8 @@
 \begin{equation} \label{eq:asm_dmp}
 \left\{ \begin{aligned}
- u^{n}_I = u^{n-1}_I + \frac{1}{e_{1u} } \delta _{i+1/2} \left( {A_D
+ u^{n}_I = u^{n-1}_I + \frac{1}{e_{1u} } \delta_{i+1/2} \left( {A_D
 \;\chi^{n-1}_I } \right) \\
 \\
- v^{n}_I = v^{n-1}_I + \frac{1}{e_{2v} } \delta _{j+1/2} \left( {A_D
+ v^{n}_I = v^{n-1}_I + \frac{1}{e_{2v} } \delta_{j+1/2} \left( {A_D
 \;\chi^{n-1}_I } \right) \\
 \end{aligned} \right.,
@@ -101,6 +101,6 @@
 \begin{equation} \label{eq:asm_div}
 \chi^{n-1}_I = \frac{1}{e_{1t}\,e_{2t}\,e_{3t} }
-                \left( {\delta _i \left[ {e_{2u}\,e_{3u}\,u^{n-1}_I} \right]
-                       +\delta _j \left[ {e_{1v}\,e_{3v}\,v^{n-1}_I} \right]} \right).
+                \left( {\delta_i \left[ {e_{2u}\,e_{3u}\,u^{n-1}_I} \right]
+                       +\delta_j \left[ {e_{1v}\,e_{3v}\,v^{n-1}_I} \right]} \right).
 \end{equation}
 By the application of \autoref{eq:asm_dmp} and \autoref{eq:asm_dmp} the divergence is filtered in each iteration,
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIU.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIU.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DIU.tex	(revision 10406)
@@ -61,10 +61,10 @@
 The warm layer is calculated using the model of \citet{Takaya_al_JGR10} (TAKAYA10 model hereafter).
 This is a simple flux based model that is defined by the equations
-\begin{eqnarray}
+\begin{align}
 \frac{\partial{\Delta T_{\rm{wl}}}}{\partial{t}}&=&\frac{Q(\nu+1)}{D_T\rho_w c_p
 \nu}-\frac{(\nu+1)ku^*_{w}f(L_a)\Delta T}{D_T\Phi\!\left(\frac{D_T}{L}\right)} \mbox{,}
 \label{eq:ecmwf1} \\
 L&=&\frac{\rho_w c_p u^{*^3}_{w}}{\kappa g \alpha_w Q }\mbox{,}\label{eq:ecmwf2}
-\end{eqnarray}
+\end{align}
 where $\Delta T_{\rm{wl}}$ is the temperature difference between the top of the warm layer and the depth $D_T=3$\,m at which there is assumed to be no diurnal signal.
 In equation (\autoref{eq:ecmwf1}) $\alpha_w=2\times10^{-4}$ is the thermal expansion coefficient of water,
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex	(revision 10406)
@@ -113,5 +113,5 @@
 \begin{subequations} \label{eq:di_mi}
 \begin{align}
- \delta _i [q]       &=  \  \    q(i+1/2)  - q(i-1/2)		\\
+ \delta_i [q]       &=  \  \    q(i+1/2)  - q(i-1/2)		\\
  \overline q^{\,i} &= \left\{ q(i+1/2) + q(i-1/2) \right\} \; / \; 2
 \end{align}
@@ -124,7 +124,7 @@
 These operators have the following discrete forms in the curvilinear $s$-coordinate system:
 \begin{equation} \label{eq:DOM_grad}
-\nabla q\equiv 	\frac{1}{e_{1u} } \delta _{i+1/2 } [q] \;\,\mathbf{i}
-		+	\frac{1}{e_{2v} } \delta _{j+1/2 } [q] \;\,\mathbf{j}
-		+	\frac{1}{e_{3w}} \delta _{k+1/2} [q] \;\,\mathbf{k}
+\nabla q\equiv 	\frac{1}{e_{1u} } \delta_{i+1/2 } [q] \;\,\mathbf{i}
+		+	\frac{1}{e_{2v} } \delta_{j+1/2 } [q] \;\,\mathbf{j}
+		+	\frac{1}{e_{3w}} \delta_{k+1/2} [q] \;\,\mathbf{k}
 \end{equation}
 \begin{multline} \label{eq:DOM_lap}
@@ -138,15 +138,15 @@
 defined at vector points $(u,v,w)$ has its three curl components defined at $vw$-, $uw$, and $f$-points,
 and its divergence defined at $t$-points:
-\begin{eqnarray}  \label{eq:DOM_curl}
+\begin{align}  \label{eq:DOM_curl}
  \nabla \times {\rm{\bf A}}\equiv &
       \frac{1}{e_{2v}\,e_{3vw} } \ \left( \delta_{j +1/2} \left[e_{3w}\,a_3 \right] -\delta_{k+1/2} \left[e_{2v} \,a_2 \right] \right)  &\ \mathbf{i} \\ 
  +& \frac{1}{e_{2u}\,e_{3uw}} \ \left( \delta_{k+1/2} \left[e_{1u}\,a_1  \right] -\delta_{i +1/2} \left[e_{3w}\,a_3 \right] \right)  &\ \mathbf{j} \\
  +& \frac{1}{e_{1f} \,e_{2f}    } \ \left( \delta_{i +1/2} \left[e_{2v}\,a_2  \right] -\delta_{j +1/2} \left[e_{1u}\,a_1 \right] \right)  &\ \mathbf{k}
- \end{eqnarray}
-\begin{eqnarray} \label{eq:DOM_div}
+ \end{align}
+\begin{align} \label{eq:DOM_div}
 \nabla \cdot \rm{\bf A} \equiv 
     \frac{1}{e_{1t}\,e_{2t}\,e_{3t}} \left( \delta_i \left[e_{2u}\,e_{3u}\,a_1 \right]
                                            +\delta_j \left[e_{1v}\,e_{3v}\,a_2 \right] \right)+\frac{1}{e_{3t} }\delta_k \left[a_3 \right]
-\end{eqnarray}
+\end{align}
 
 The vertical average over the whole water column denoted by an overbar becomes for a quantity $q$ which
@@ -181,6 +181,6 @@
 \begin{align} 
 \label{eq:DOM_di_adj}
-\sum\limits_i { a_i \;\delta _i \left[ b \right]} 
-	&\equiv -\sum\limits_i {\delta _{i+1/2} \left[ a \right]\;b_{i+1/2} }      \\
+\sum\limits_i { a_i \;\delta_i \left[ b \right]} 
+	&\equiv -\sum\limits_i {\delta_{i+1/2} \left[ a \right]\;b_{i+1/2} }      \\
 \label{eq:DOM_mi_adj}
 \sum\limits_i { a_i \;\overline b^{\,i}} 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex	(revision 10406)
@@ -19,9 +19,9 @@
 
 The prognostic ocean dynamics equation can be summarized as follows:
-\begin{equation*}
+\[
 \text{NXT} = \dbinom	{\text{VOR} + \text{KEG} + \text {ZAD} }
 						{\text{COR} + \text{ADV}                       }
 			+ \text{HPG} + \text{SPG} + \text{LDF} + \text{ZDF}
-\end{equation*}
+\]
 NXT stands for next, referring to the time-stepping.
 The first group of terms on the rhs of this equation corresponds to the Coriolis and advection terms that
@@ -73,6 +73,6 @@
 The vorticity is defined at an $f$-point ($i.e.$ corner point) as follows:
 \begin{equation} \label{eq:divcur_cur}
-\zeta =\frac{1}{e_{1f}\,e_{2f} }\left( {\;\delta _{i+1/2} \left[ {e_{2v}\;v} \right]
-						        -\delta _{j+1/2} \left[ {e_{1u}\;u} \right]\;} \right)
+\zeta =\frac{1}{e_{1f}\,e_{2f} }\left( {\;\delta_{i+1/2} \left[ {e_{2v}\;v} \right]
+						        -\delta_{j+1/2} \left[ {e_{1u}\;u} \right]\;} \right)
 \end{equation} 
 
@@ -81,6 +81,6 @@
 \begin{equation} \label{eq:divcur_div}
 \chi =\frac{1}{e_{1t}\,e_{2t}\,e_{3t} }
-		\left( {\delta _i \left[ {e_{2u}\,e_{3u}\,u} \right]
-		       +\delta _j \left[ {e_{1v}\,e_{3v}\,v} \right]} \right)
+		\left( {\delta_i \left[ {e_{2u}\,e_{3u}\,u} \right]
+		       +\delta_j \left[ {e_{1v}\,e_{3v}\,v} \right]} \right)
 \end{equation} 
 
@@ -109,13 +109,13 @@
 \begin{aligned}
 \frac{\partial \eta }{\partial t}
-&\equiv    \frac{1}{e_{1t} e_{2t} }\sum\limits_k { \left\{  \delta _i \left[ {e_{2u}\,e_{3u}\;u} \right]
-                                                                                  +\delta _j \left[ {e_{1v}\,e_{3v}\;v} \right]  \right\} } 
-           -    \frac{\textit{emp}}{\rho _w }   \\
-&\equiv    \sum\limits_k {\chi \ e_{3t}}  -  \frac{\textit{emp}}{\rho _w }
+&\equiv    \frac{1}{e_{1t} e_{2t} }\sum\limits_k { \left\{  \delta_i \left[ {e_{2u}\,e_{3u}\;u} \right]
+                                                                                  +\delta_j \left[ {e_{1v}\,e_{3v}\;v} \right]  \right\} } 
+           -    \frac{\textit{emp}}{\rho_w }   \\
+&\equiv    \sum\limits_k {\chi \ e_{3t}}  -  \frac{\textit{emp}}{\rho_w }
 \end{aligned}
 \end{equation}
 where \textit{emp} is the surface freshwater budget (evaporation minus precipitation), 
 expressed in Kg/m$^2$/s (which is equal to mm/s),
-and $\rho _w$=1,035~Kg/m$^3$ is the reference density of sea water (Boussinesq approximation).
+and $\rho_w$=1,035~Kg/m$^3$ is the reference density of sea water (Boussinesq approximation).
 If river runoff is expressed as a surface freshwater flux (see \autoref{chap:SBC}) then
 \textit{emp} can be written as the evaporation minus precipitation, minus the river runoff. 
@@ -355,6 +355,6 @@
 \begin{equation} \label{eq:dynkeg}
 \left\{ \begin{aligned}
- -\frac{1}{2 \; e_{1u} }  & \ \delta _{i+1/2} \left[ {\overline {u^2}^{\,i} + \overline{v^2}^{\,j}} \right]   \\
- -\frac{1}{2 \; e_{2v} }  & \ \delta _{j+1/2} \left[ {\overline {u^2}^{\,i} + \overline{v^2}^{\,j}} \right]    
+ -\frac{1}{2 \; e_{1u} }  & \ \delta_{i+1/2} \left[ {\overline {u^2}^{\,i} + \overline{v^2}^{\,j}} \right]   \\
+ -\frac{1}{2 \; e_{2v} }  & \ \delta_{j+1/2} \left[ {\overline {u^2}^{\,i} + \overline{v^2}^{\,j}} \right]    
 \end{aligned} \right.
 \end{equation} 
@@ -373,6 +373,6 @@
 \begin{equation} \label{eq:dynzad}
 \left\{ 		\begin{aligned}
--\frac{1} {e_{1u}\,e_{2u}\,e_{3u}} &\ \overline{\ \overline{ e_{1t}\,e_{2t}\;w } ^{\,i+1/2}  \;\delta _{k+1/2} \left[ u \right]\  }^{\,k}  \\
--\frac{1} {e_{1v}\,e_{2v}\,e_{3v}}  &\ \overline{\ \overline{ e_{1t}\,e_{2t}\;w } ^{\,j+1/2}  \;\delta _{k+1/2} \left[ u \right]\  }^{\,k} 
+-\frac{1} {e_{1u}\,e_{2u}\,e_{3u}} &\ \overline{\ \overline{ e_{1t}\,e_{2t}\;w } ^{\,i+1/2}  \;\delta_{k+1/2} \left[ u \right]\  }^{\,k}  \\
+-\frac{1} {e_{1v}\,e_{2v}\,e_{3v}}  &\ \overline{\ \overline{ e_{1t}\,e_{2t}\;w } ^{\,j+1/2}  \;\delta_{k+1/2} \left[ u \right]\  }^{\,k} 
 \end{aligned}         \right.
 \end{equation} 
@@ -414,6 +414,6 @@
 \begin{multline} \label{eq:dyncor_metric}
 f+\frac{1}{e_1 e_2 }\left( {v\frac{\partial e_2 }{\partial i}  -  u\frac{\partial e_1 }{\partial j}} \right)  \\
-   \equiv   f + \frac{1}{e_{1f} e_{2f} } \left( { \ \overline v ^{i+1/2}\delta _{i+1/2} \left[ {e_{2u} } \right]  
-                                                                 -  \overline u ^{j+1/2}\delta _{j+1/2} \left[ {e_{1u} } \right]  }  \ \right)
+   \equiv   f + \frac{1}{e_{1f} e_{2f} } \left( { \ \overline v ^{i+1/2}\delta_{i+1/2} \left[ {e_{2u} } \right]  
+                                                                 -  \overline u ^{j+1/2}\delta_{j+1/2} \left[ {e_{1u} } \right]  }  \ \right)
 \end{multline} 
 
@@ -434,12 +434,12 @@
 \begin{aligned}
 \frac{1}{e_{1u}\,e_{2u}\,e_{3u}} 
-\left(      \delta _{i+1/2} \left[ \overline{e_{2u}\,e_{3u}\;u }^{i       }  \ u_t      \right]    
-          + \delta _{j       } \left[ \overline{e_{1u}\,e_{3u}\;v }^{i+1/2}  \ u_f      \right] \right.  \ \;   \\
-\left.   + \delta _{k      } \left[ \overline{e_{1w}\,e_{2w}\;w}^{i+1/2}  \ u_{uw} \right] \right)   \\
+\left(      \delta_{i+1/2} \left[ \overline{e_{2u}\,e_{3u}\;u }^{i       }  \ u_t      \right]    
+          + \delta_{j       } \left[ \overline{e_{1u}\,e_{3u}\;v }^{i+1/2}  \ u_f      \right] \right.  \ \;   \\
+\left.   + \delta_{k      } \left[ \overline{e_{1w}\,e_{2w}\;w}^{i+1/2}  \ u_{uw} \right] \right)   \\
 \\
 \frac{1}{e_{1v}\,e_{2v}\,e_{3v}} 
-\left(     \delta _{i       } \left[ \overline{e_{2u}\,e_{3u }\;u }^{j+1/2} \ v_f       \right] 
-         + \delta _{j+1/2} \left[ \overline{e_{1u}\,e_{3u }\;v }^{i       } \ v_t       \right] \right.  \ \, \, \\
-\left.  + \delta _{k      } \left[ \overline{e_{1w}\,e_{2w}\;w}^{j+1/2} \ v_{vw}  \right] \right) \\
+\left(     \delta_{i       } \left[ \overline{e_{2u}\,e_{3u }\;u }^{j+1/2} \ v_f       \right] 
+         + \delta_{j+1/2} \left[ \overline{e_{1u}\,e_{3u }\;v }^{i       } \ v_t       \right] \right.  \ \, \, \\
+\left.  + \delta_{k      } \left[ \overline{e_{1w}\,e_{2w}\;w}^{j+1/2} \ v_{vw}  \right] \right) \\
 \end{aligned}
 \right.
@@ -490,5 +490,5 @@
 \end{cases}
 \end{equation}
-where $u"_{i+1/2} =\delta _{i+1/2} \left[ {\delta _i \left[ u \right]} \right]$.
+where $u"_{i+1/2} =\delta_{i+1/2} \left[ {\delta_i \left[ u \right]} \right]$.
 This results in a dissipatively dominant ($i.e.$ hyper-diffusive) truncation error
 \citep{Shchepetkin_McWilliams_OM05}.
@@ -562,8 +562,8 @@
 \begin{equation} \label{eq:dynhpg_zco_surf}
 \left\{ \begin{aligned}
-					\left. \delta _{i+1/2} \left[  p^h 			 \right] \right|_{k=km} 
-&= \frac{1}{2} g \ 	\left. \delta _{i+1/2} \left[  e_{3w} \ \rho \right] \right|_{k=km}   \\
-     					\left. \delta _{j+1/2} \left[  p^h  			 \right] \right|_{k=km} 
-&= \frac{1}{2} g \ 	\left. \delta _{j+1/2} \left[  e_{3w} \ \rho \right] \right|_{k=km}   \\
+					\left. \delta_{i+1/2} \left[  p^h 			 \right] \right|_{k=km} 
+&= \frac{1}{2} g \ 	\left. \delta_{i+1/2} \left[  e_{3w} \ \rho \right] \right|_{k=km}   \\
+     					\left. \delta_{j+1/2} \left[  p^h  			 \right] \right|_{k=km} 
+&= \frac{1}{2} g \ 	\left. \delta_{j+1/2} \left[  e_{3w} \ \rho \right] \right|_{k=km}   \\
 \end{aligned} \right.
 \end{equation} 
@@ -572,10 +572,10 @@
 \begin{equation} \label{eq:dynhpg_zco}
 \left\{ \begin{aligned}
-					\left. \delta _{i+1/2} \left[  p^h 			 \right] \right|_{k} 
-&=					\left. \delta _{i+1/2} \left[  p^h 			 \right] \right|_{k-1} 
-+    \frac{1}{2}\;g\;	\left. \delta _{i+1/2} \left[  e_{3w} \ \overline {\rho}^{k+1/2} \right] \right|_{k}   \\
-     					\left. \delta _{j+1/2} \left[  p^h  			 \right] \right|_{k} 
-&=     				\left. \delta _{j+1/2} \left[  p^h  			 \right] \right|_{k-1} 
-+    \frac{1}{2}\;g\;	\left. \delta _{j+1/2} \left[  e_{3w} \ \overline {\rho}^{k+1/2} \right] \right|_{k}   \\
+					\left. \delta_{i+1/2} \left[  p^h 			 \right] \right|_{k} 
+&=					\left. \delta_{i+1/2} \left[  p^h 			 \right] \right|_{k-1} 
++    \frac{1}{2}\;g\;	\left. \delta_{i+1/2} \left[  e_{3w} \ \overline {\rho}^{k+1/2} \right] \right|_{k}   \\
+     					\left. \delta_{j+1/2} \left[  p^h  			 \right] \right|_{k} 
+&=     				\left. \delta_{j+1/2} \left[  p^h  			 \right] \right|_{k-1} 
++    \frac{1}{2}\;g\;	\left. \delta_{j+1/2} \left[  e_{3w} \ \overline {\rho}^{k+1/2} \right] \right|_{k}   \\
 \end{aligned} \right.
 \end{equation} 
@@ -622,8 +622,8 @@
 \begin{equation} \label{eq:dynhpg_sco}
 \left\{ \begin{aligned}
- - \frac{1}    					{\rho_o \, e_{1u}} \;	\delta _{i+1/2} \left[  p^h  \right] 
-+ \frac{g\; \overline {\rho}^{i+1/2}}	{\rho_o \, e_{1u}} \;	\delta _{i+1/2} \left[  z_t   \right]    \\
- - \frac{1}    					{\rho_o \, e_{2v}} \;	\delta _{j+1/2} \left[  p^h  \right]  
-+ \frac{g\; \overline {\rho}^{j+1/2}}	{\rho_o \, e_{2v}} \;	\delta _{j+1/2} \left[  z_t   \right]    \\
+ - \frac{1}    					{\rho_o \, e_{1u}} \;	\delta_{i+1/2} \left[  p^h  \right] 
++ \frac{g\; \overline {\rho}^{i+1/2}}	{\rho_o \, e_{1u}} \;	\delta_{i+1/2} \left[  z_t   \right]    \\
+ - \frac{1}    					{\rho_o \, e_{2v}} \;	\delta_{j+1/2} \left[  p^h  \right]  
++ \frac{g\; \overline {\rho}^{j+1/2}}	{\rho_o \, e_{2v}} \;	\delta_{j+1/2} \left[  z_t   \right]    \\
 \end{aligned} \right.
 \end{equation} 
@@ -695,5 +695,5 @@
 \begin{equation} \label{eq:dynhpg_lf}
 \frac{u^{t+\rdt}-u^{t-\rdt}}{2\rdt} = \;\cdots \;
-	-\frac{1}{\rho _o \,e_{1u} }\delta _{i+1/2} \left[ {p_h^t } \right]
+	-\frac{1}{\rho_o \,e_{1u} }\delta_{i+1/2} \left[ {p_h^t } \right]
 \end{equation}
 
@@ -701,5 +701,5 @@
 \begin{equation} \label{eq:dynhpg_imp}
 \frac{u^{t+\rdt}-u^{t-\rdt}}{2\rdt} = \;\cdots \;
-	-\frac{1}{4\,\rho _o \,e_{1u} } \delta_{i+1/2} \left[ p_h^{t+\rdt} +2\,p_h^t +p_h^{t-\rdt}  \right]
+	-\frac{1}{4\,\rho_o \,e_{1u} } \delta_{i+1/2} \left[ p_h^{t+\rdt} +2\,p_h^t +p_h^{t-\rdt}  \right]
 \end{equation}
 
@@ -783,6 +783,6 @@
 \begin{equation} \label{eq:dynspg_exp}
 \left\{ \begin{aligned}
- - \frac{1}{e_{1u}\,\rho_o} \;	\delta _{i+1/2} \left[  \,\rho \,\eta\,  \right] 	\\
- - \frac{1}{e_{2v}\,\rho_o} \;	\delta _{j+1/2} \left[  \,\rho \,\eta\,  \right]  
+ - \frac{1}{e_{1u}\,\rho_o} \;	\delta_{i+1/2} \left[  \,\rho \,\eta\,  \right] 	\\
+ - \frac{1}{e_{2v}\,\rho_o} \;	\delta_{j+1/2} \left[  \,\rho \,\eta\,  \right]  
 \end{aligned} \right.
 \end{equation} 
@@ -1093,10 +1093,10 @@
 \begin{equation} \label{eq:dynldf_lap}
 \left\{ \begin{aligned}
- D_u^{l{\rm {\bf U}}} =\frac{1}{e_{1u} }\delta _{i+1/2} \left[ {A_T^{lm} 
-\;\chi } \right]-\frac{1}{e_{2u} {\kern 1pt}e_{3u} }\delta _j \left[ 
+ D_u^{l{\rm {\bf U}}} =\frac{1}{e_{1u} }\delta_{i+1/2} \left[ {A_T^{lm} 
+\;\chi } \right]-\frac{1}{e_{2u} {\kern 1pt}e_{3u} }\delta_j \left[ 
 {A_f^{lm} \;e_{3f} \zeta } \right] \\ 
 \\
- D_v^{l{\rm {\bf U}}} =\frac{1}{e_{2v} }\delta _{j+1/2} \left[ {A_T^{lm} 
-\;\chi } \right]+\frac{1}{e_{1v} {\kern 1pt}e_{3v} }\delta _i \left[ 
+ D_v^{l{\rm {\bf U}}} =\frac{1}{e_{2v} }\delta_{j+1/2} \left[ {A_T^{lm} 
+\;\chi } \right]+\frac{1}{e_{1v} {\kern 1pt}e_{3v} }\delta_i \left[ 
 {A_f^{lm} \;e_{3f} \zeta } \right] \\ 
 \end{aligned} \right.
@@ -1127,12 +1127,12 @@
 \begin{split}
  D_u^{l\textbf{U}} &= \frac{1}{e_{1u} \, e_{2u} \, e_{3u} }	\\
-&  \left\{\quad  {\delta _{i+1/2} \left[ {A_T^{lm}  \left( 
-	 {\frac{e_{2t} \; e_{3t} }{e_{1t} } \,\delta _{i}[u]
-	-e_{2t} \; r_{1t} \,\overline{\overline {\delta _{k+1/2}[u]}}^{\,i,\,k}}
+&  \left\{\quad  {\delta_{i+1/2} \left[ {A_T^{lm}  \left( 
+	 {\frac{e_{2t} \; e_{3t} }{e_{1t} } \,\delta_{i}[u]
+	-e_{2t} \; r_{1t} \,\overline{\overline {\delta_{k+1/2}[u]}}^{\,i,\,k}}
  \right)} \right]} 	\right.
 \\ 
 & \qquad +\ \delta_j \left[ {A_f^{lm} \left( {\frac{e_{1f}\,e_{3f} }{e_{2f} 
-}\,\delta _{j+1/2} [u] - e_{1f}\, r_{2f} 
-\,\overline{\overline {\delta _{k+1/2} [u]}} ^{\,j+1/2,\,k}} 
+}\,\delta_{j+1/2} [u] - e_{1f}\, r_{2f} 
+\,\overline{\overline {\delta_{k+1/2} [u]}} ^{\,j+1/2,\,k}} 
 \right)} \right] 
 \\ 
@@ -1150,12 +1150,12 @@
 \\
  D_v^{l\textbf{V}} &= \frac{1}{e_{1v} \, e_{2v} \, e_{3v} }    \\
-&  \left\{\quad  {\delta _{i+1/2} \left[ {A_f^{lm}  \left( 
-	 {\frac{e_{2f} \; e_{3f} }{e_{1f} } \,\delta _{i+1/2}[v]
-	-e_{2f} \; r_{1f} \,\overline{\overline {\delta _{k+1/2}[v]}}^{\,i+1/2,\,k}}
+&  \left\{\quad  {\delta_{i+1/2} \left[ {A_f^{lm}  \left( 
+	 {\frac{e_{2f} \; e_{3f} }{e_{1f} } \,\delta_{i+1/2}[v]
+	-e_{2f} \; r_{1f} \,\overline{\overline {\delta_{k+1/2}[v]}}^{\,i+1/2,\,k}}
  \right)} \right]} 	\right.
 \\ 
 & \qquad +\ \delta_j \left[ {A_T^{lm} \left( {\frac{e_{1t}\,e_{3t} }{e_{2t} 
-}\,\delta _{j} [v] - e_{1t}\, r_{2t} 
-\,\overline{\overline {\delta _{k+1/2} [v]}} ^{\,j,\,k}} 
+}\,\delta_{j} [v] - e_{1t}\, r_{2t} 
+\,\overline{\overline {\delta_{k+1/2} [v]}} ^{\,j,\,k}} 
 \right)} \right] 
 \\ 
@@ -1213,9 +1213,9 @@
 \begin{equation} \label{eq:dynzdf}
 \left\{   \begin{aligned}
-D_u^{vm} &\equiv \frac{1}{e_{3u}} \ \delta _k \left[ \frac{A_{uw}^{vm} }{e_{3uw} }
-                              \ \delta _{k+1/2} [\,u\,]         \right]     \\
+D_u^{vm} &\equiv \frac{1}{e_{3u}} \ \delta_k \left[ \frac{A_{uw}^{vm} }{e_{3uw} }
+                              \ \delta_{k+1/2} [\,u\,]         \right]     \\
 \\
-D_v^{vm} &\equiv \frac{1}{e_{3v}} \ \delta _k \left[ \frac{A_{vw}^{vm} }{e_{3vw} }
-                              \ \delta _{k+1/2} [\,v\,]         \right]
+D_v^{vm} &\equiv \frac{1}{e_{3v}} \ \delta_k \left[ \frac{A_{vw}^{vm} }{e_{3vw} }
+                              \ \delta_{k+1/2} [\,v\,]         \right]
 \end{aligned}   \right.
 \end{equation} 
@@ -1228,7 +1228,7 @@
 \begin{equation} \label{eq:dynzdf_sbc}
 \left.{\left( {\frac{A^{vm} }{e_3 }\ \frac{\partial \textbf{U}_h}{\partial k}} \right)} \right|_{z=1}
-	 = \frac{1}{\rho _o} \binom{\tau _u}{\tau _v }
-\end{equation}
-where $\left( \tau _u ,\tau _v \right)$ are the two components of the wind stress vector in
+	 = \frac{1}{\rho_o} \binom{\tau_u}{\tau_v }
+\end{equation}
+where $\left( \tau_u ,\tau_v \right)$ are the two components of the wind stress vector in
 the (\textbf{i},\textbf{j}) coordinate system.
 The high mixing coefficients in the surface mixed layer ensure that the surface wind stress is distributed in 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LBC.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LBC.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_LBC.tex	(revision 10406)
@@ -46,5 +46,5 @@
 \begin{equation} \label{eq:lbc_aaaa}
 \frac{A^{lT} }{e_1 }\frac{\partial T}{\partial i}\equiv \frac{A_u^{lT} 
-}{e_{1u} } \; \delta _{i+1 / 2} \left[ T \right]\;\;mask_u 
+}{e_{1u} } \; \delta_{i+1 / 2} \left[ T \right]\;\;mask_u 
 \end{equation}
 (where mask$_{u}$ is the mask array at a $u$-point) ensures that the heat flux is zero inside land and
@@ -106,13 +106,13 @@
   Therefore, the vorticity along the coastlines is given by: 
 
-\begin{equation*}
+\[
 \zeta \equiv 2 \left(\delta_{i+1/2} \left[e_{2v} v \right] - \delta_{j+1/2} \left[e_{1u} u \right] \right) / \left(e_{1f} e_{2f} \right) \ ,
-\end{equation*}
+\]
 where $u$ and $v$ are masked fields.
 Setting the mask$_{f}$ array to $2$ along the coastline provides a vorticity field computed with
 the no-slip boundary condition, simply by multiplying it by the mask$_{f}$ :
 \begin{equation} \label{eq:lbc_bbbb}
-\zeta \equiv \frac{1}{e_{1f} {\kern 1pt}e_{2f} }\left( {\delta _{i+1/2} 
-\left[ {e_{2v} \,v} \right]-\delta _{j+1/2} \left[ {e_{1u} \,u} \right]} 
+\zeta \equiv \frac{1}{e_{1f} {\kern 1pt}e_{2f} }\left( {\delta_{i+1/2} 
+\left[ {e_{2v} \,v} \right]-\delta_{j+1/2} \left[ {e_{1u} \,u} \right]} 
 \right)\;\mbox{mask}_f 
 \end{equation}
@@ -279,8 +279,8 @@
 The whole domain dimensions are named \np{jpiglo}, \np{jpjglo} and \jp{jpk}.
 The relationship between the whole domain and a sub-domain is:
-\begin{eqnarray} 
+\begin{align} 
       jpi & = & ( jpiglo-2*jpreci + (jpni-1) ) / jpni + 2*jpreci  \nonumber \\
       jpj & = & ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj  \label{eq:lbc_jpi}
-\end{eqnarray}
+\end{align}
 where \jp{jpni}, \jp{jpnj} are the number of processors following the i- and j-axis.
 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_OBS.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_OBS.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_OBS.tex	(revision 10406)
@@ -578,5 +578,5 @@
 All horizontal interpolation methods implemented in NEMO estimate the value of a model variable $x$ at point $P$ as
 a weighted linear combination of the values of the model variables at the grid points ${\rm A}$, ${\rm B}$ etc.:
-\begin{eqnarray}
+\begin{align}
 {x_{}}_{\rm P} & \hspace{-2mm} = \hspace{-2mm} & 
 \frac{1}{w} \left( {w_{}}_{\rm A} {x_{}}_{\rm A} + 
@@ -584,5 +584,5 @@
                    {w_{}}_{\rm C} {x_{}}_{\rm C} + 
                    {w_{}}_{\rm D} {x_{}}_{\rm D} \right)
-\end{eqnarray}
+\end{align}
 where ${w_{}}_{\rm A}$, ${w_{}}_{\rm B}$ etc. are the respective weights for the model field at
 points ${\rm A}$, ${\rm B}$ etc., and $w = {w_{}}_{\rm A} + {w_{}}_{\rm B} + {w_{}}_{\rm C} + {w_{}}_{\rm D}$.
@@ -597,10 +597,10 @@
   For example, the weight given to the field ${x_{}}_{\rm A}$ is specified as the product of the distances
   from ${\rm P}$ to the other points:
-  \begin{eqnarray}
+  \begin{align}
   {w_{}}_{\rm A} = s({\rm P}, {\rm B}) \, s({\rm P}, {\rm C}) \, s({\rm P}, {\rm D})
   \nonumber
-  \end{eqnarray}
+  \end{align}
   where 
-  \begin{eqnarray}
+  \begin{align}
    s\left ({\rm P}, {\rm M} \right ) 
      & \hspace{-2mm} = \hspace{-2mm} & 
@@ -610,22 +610,22 @@
                \cos ({\lambda_{}}_{\rm M} - {\lambda_{}}_{\rm P}) 
                    \right\}
-   \end{eqnarray}
+   \end{align}
    and $M$ corresponds to $B$, $C$ or $D$.
    A more stable form of the great-circle distance formula for small distances ($x$ near 1)
    involves the arcsine function ($e.g.$ see p.~101 of \citet{Daley_Barker_Bk01}:
-   \begin{eqnarray}
+   \begin{align}
    s\left( {\rm P}, {\rm M} \right) 
      & \hspace{-2mm} = \hspace{-2mm} & 
       \sin^{-1} \! \left\{ \sqrt{ 1 - x^2 } \right\}
    \nonumber
-   \end{eqnarray}
+   \end{align}
    where
-   \begin{eqnarray}
+   \begin{align}
     x & \hspace{-2mm} = \hspace{-2mm} & 
       {a_{}}_{\rm M} {a_{}}_{\rm P} + {b_{}}_{\rm M} {b_{}}_{\rm P} + {c_{}}_{\rm M} {c_{}}_{\rm P}
    \nonumber
-   \end{eqnarray}
+   \end{align}
    and 
-   \begin{eqnarray}
+   \begin{align}
       {a_{}}_{\rm M} & \hspace{-2mm} = \hspace{-2mm} & \sin {\phi_{}}_{\rm M}, 
       \nonumber \\
@@ -641,9 +641,9 @@
       \nonumber
    \nonumber
-  \end{eqnarray}
+  \end{align}
 
 \item[2.] {\bf Great-Circle distance-weighted interpolation with small angle approximation.}
   Similar to the previous interpolation but with the distance $s$ computed as
-  \begin{eqnarray}
+  \begin{align}
     s\left( {\rm P}, {\rm M} \right) 
      & \hspace{-2mm} = \hspace{-2mm} & 
@@ -651,5 +651,5 @@
       + \left( {\lambda_{}}_{\rm M} - {\lambda_{}}_{\rm P} \right)^{2}
         \cos^{2} {\phi_{}}_{\rm M} }
-  \end{eqnarray}
+  \end{align}
   where $M$ corresponds to $A$, $B$, $C$ or $D$.
 
@@ -719,5 +719,5 @@
 denote the bottom left, bottom right, top left and top right corner points of the cell, respectively. 
 To determine if P is inside the cell, we verify that the cross-products 
-\begin{eqnarray}
+\begin{align}
 \begin{array}{lllll}
 {{\bf r}_{}}_{\rm PA} \times {{\bf r}_{}}_{\rm PC}
@@ -743,5 +743,5 @@
 \end{array}
 \label{eq:cross}
-\end{eqnarray}
+\end{align}
 point in the opposite direction to the unit normal $\widehat{\bf k}$
 (i.e., that the coefficients of $\widehat{\bf k}$ are negative),
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_SBC.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_SBC.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_SBC.tex	(revision 10406)
@@ -19,5 +19,5 @@
 \begin{itemize}
 \item
-  the two components of the surface ocean stress $\left( {\tau _u \;,\;\tau _v} \right)$
+  the two components of the surface ocean stress $\left( {\tau_u \;,\;\tau_v} \right)$
 \item
   the incoming solar and non solar heat fluxes $\left( {Q_{ns} \;,\;Q_{sr} } \right)$
@@ -391,5 +391,5 @@
 The symbolic algorithm used to calculate values on the model grid is now:
 
-\begin{equation*} \begin{split}
+\[ \begin{split}
 f_{m}(i,j) =  f_{m}(i,j) +& \sum_{k=1}^{4} {wgt(k)f(idx(src(k)))}     
               +   \sum_{k=5}^{8} {wgt(k)\left.\frac{\partial f}{\partial i}\right| _{idx(src(k))} }    \\
@@ -397,5 +397,5 @@
               +   \sum_{k=13}^{16} {wgt(k)\left.\frac{\partial ^2 f}{\partial i \partial j}\right| _{idx(src(k))} }
 \end{split}
-\end{equation*}
+\]
 The gradients here are taken with respect to the horizontal indices and not distances since
 the spatial dependency has been absorbed into the weights.
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_TRA.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_TRA.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_TRA.tex	(revision 10406)
@@ -27,8 +27,8 @@
 The two active tracers are potential temperature and salinity.
 Their prognostic equations can be summarized as follows:
-\begin{equation*}
+\[
 \text{NXT} = \text{ADV}+\text{LDF}+\text{ZDF}+\text{SBC}
                    \ (+\text{QSR})\ (+\text{BBC})\ (+\text{BBL})\ (+\text{DMP})
-\end{equation*}
+\]
 
 NXT stands for next, referring to the time-stepping.
@@ -76,7 +76,7 @@
 \begin{equation} \label{eq:tra_adv}
 ADV_\tau =-\frac{1}{b_t} \left( 
-\;\delta _i \left[ e_{2u}\,e_{3u} \;  u\; \tau _u  \right]
-+\delta _j \left[ e_{1v}\,e_{3v}  \;  v\; \tau _v  \right] \; \right)
--\frac{1}{e_{3t}} \;\delta _k \left[ w\; \tau _w \right]
+\;\delta_i \left[ e_{2u}\,e_{3u} \;  u\; \tau_u  \right]
++\delta_j \left[ e_{1v}\,e_{3v}  \;  v\; \tau_v  \right] \; \right)
+-\frac{1}{e_{3t}} \;\delta_k \left[ w\; \tau_w \right]
 \end{equation}
 where $\tau$ is either T or S, and $b_t= e_{1t}\,e_{2t}\,e_{3t}$ is the volume of $T$-cells.
@@ -125,5 +125,5 @@
   the vertical boundary condition is applied at the fixed surface $z=0$ rather than on the moving surface $z=\eta$.
   There is a non-zero advective flux which is set for all advection schemes as
-  $\left. {\tau _w } \right|_{k=1/2} =T_{k=1} $,
+  $\left. {\tau_w } \right|_{k=1/2} =T_{k=1} $,
   $i.e.$ the product of surface velocity (at $z=0$) by the first level tracer value.
 \item[non-linear free surface:]
@@ -194,5 +194,5 @@
 For example, in the $i$-direction :
 \begin{equation} \label{eq:tra_adv_cen2}
-\tau _u^{cen2} =\overline T ^{i+1/2}
+\tau_u^{cen2} =\overline T ^{i+1/2}
 \end{equation}
 
@@ -213,6 +213,6 @@
 For example, in the $i$-direction:
 \begin{equation} \label{eq:tra_adv_cen4}
-\tau _u^{cen4} 
-=\overline{   T - \frac{1}{6}\,\delta _i \left[ \delta_{i+1/2}[T] \,\right]   }^{\,i+1/2}
+\tau_u^{cen4} 
+=\overline{   T - \frac{1}{6}\,\delta_i \left[ \delta_{i+1/2}[T] \,\right]   }^{\,i+1/2}
 \end{equation}
 In the vertical direction (\np{nn\_cen\_v}\forcode{ = 4}),
@@ -238,5 +238,5 @@
 
 At a $T$-grid cell adjacent to a boundary (coastline, bottom and surface),
-an additional hypothesis must be made to evaluate $\tau _u^{cen4}$.
+an additional hypothesis must be made to evaluate $\tau_u^{cen4}$.
 This hypothesis usually reduces the order of the scheme.
 Here we choose to set the gradient of $T$ across the boundary to zero.
@@ -260,10 +260,10 @@
 \begin{equation} \label{eq:tra_adv_fct}
 \begin{split}
-\tau _u^{ups}&= \begin{cases}
+\tau_u^{ups}&= \begin{cases}
  					T_{i+1} 	& \text{if $\ u_{i+1/2} <     0$} \hfill \\
  					T_i   		& \text{if $\ u_{i+1/2} \geq 0$} \hfill \\
 				  \end{cases}     \\
 \\
-\tau _u^{fct}&=\tau _u^{ups} +c_u \;\left( {\tau _u^{cen} -\tau _u^{ups} } \right)
+\tau_u^{fct}&=\tau_u^{ups} +c_u \;\left( {\tau_u^{cen} -\tau_u^{ups} } \right)
 \end{split}
 \end{equation}
@@ -287,6 +287,6 @@
 
 For stability reasons (see \autoref{chap:STP}),
-$\tau _u^{cen}$ is evaluated in (\autoref{eq:tra_adv_fct}) using the \textit{now} tracer while
-$\tau _u^{ups}$ is evaluated using the \textit{before} tracer.
+$\tau_u^{cen}$ is evaluated in (\autoref{eq:tra_adv_fct}) using the \textit{now} tracer while
+$\tau_u^{ups}$ is evaluated using the \textit{before} tracer.
 In other words, the advective part of the scheme is time stepped with a leap-frog scheme
 while a forward scheme is used for the diffusive part. 
@@ -306,8 +306,8 @@
 For example, in the $i$-direction :
 \begin{equation} \label{eq:tra_adv_mus}
-   \tau _u^{mus} = \left\{      \begin{aligned}
-         &\tau _i  &+ \frac{1}{2} \;\left( 1-\frac{u_{i+1/2} \;\rdt}{e_{1u}} \right)
+   \tau_u^{mus} = \left\{      \begin{aligned}
+         &\tau_i  &+ \frac{1}{2} \;\left( 1-\frac{u_{i+1/2} \;\rdt}{e_{1u}} \right)
          &\ \widetilde{\partial _i \tau}  & \quad \text{if }\;u_{i+1/2} \geqslant 0      \\
-         &\tau _{i+1/2} &+\frac{1}{2}\;\left( 1+\frac{u_{i+1/2} \;\rdt}{e_{1u} } \right)
+         &\tau_{i+1/2} &+\frac{1}{2}\;\left( 1+\frac{u_{i+1/2} \;\rdt}{e_{1u} } \right)
          &\ \widetilde{\partial_{i+1/2} \tau } & \text{if }\;u_{i+1/2} <0
    \end{aligned}    \right.
@@ -317,5 +317,5 @@
 
 The time stepping is performed using a forward scheme,
-that is the \textit{before} tracer field is used to evaluate $\tau _u^{mus}$.
+that is the \textit{before} tracer field is used to evaluate $\tau_u^{mus}$.
 
 For an ocean grid point adjacent to land and where the ocean velocity is directed toward land,
@@ -339,5 +339,5 @@
 For example, in the $i$-direction:
 \begin{equation} \label{eq:tra_adv_ubs}
-   \tau _u^{ubs} =\overline T ^{i+1/2}-\;\frac{1}{6} \left\{      
+   \tau_u^{ubs} =\overline T ^{i+1/2}-\;\frac{1}{6} \left\{      
    \begin{aligned}
          &\tau"_i        	& \quad \text{if }\ u_{i+1/2} \geqslant 0      \\
@@ -345,5 +345,5 @@
    \end{aligned}    \right.
 \end{equation}
-where $\tau "_i =\delta _i \left[ {\delta _{i+1/2} \left[ \tau \right]} \right]$.
+where $\tau "_i =\delta_i \left[ {\delta_{i+1/2} \left[ \tau \right]} \right]$.
 
 This results in a dissipatively dominant (i.e. hyper-diffusive) truncation error
@@ -374,5 +374,5 @@
 Note that it is straightforward to rewrite \autoref{eq:tra_adv_ubs} as follows:
 \begin{equation} \label{eq:traadv_ubs2}
-\tau _u^{ubs} = \tau _u^{cen4} + \frac{1}{12} \left\{	 
+\tau_u^{ubs} = \tau_u^{cen4} + \frac{1}{12} \left\{	 
    \begin{aligned}
 	& + \tau"_i			& \quad \text{if }\ u_{i+1/2} \geqslant 0 \\
@@ -382,6 +382,6 @@
 or equivalently 
 \begin{equation} \label{eq:traadv_ubs2b}
-u_{i+1/2} \ \tau _u^{ubs} 
-=u_{i+1/2} \ \overline{ T - \frac{1}{6}\,\delta _i\left[ \delta_{i+1/2}[T] \,\right] }^{\,i+1/2}
+u_{i+1/2} \ \tau_u^{ubs} 
+=u_{i+1/2} \ \overline{ T - \frac{1}{6}\,\delta_i\left[ \delta_{i+1/2}[T] \,\right] }^{\,i+1/2}
 - \frac{1}{2} |u|_{i+1/2} \;\frac{1}{6} \;\delta_{i+1/2}[\tau"_i]
 \end{equation}
@@ -521,6 +521,6 @@
 \begin{equation} \label{eq:tra_ldf_lap}
 D_t^{lT} =\frac{1}{b_t} \left( \;
-   \delta _{i}\left[ A_u^{lT} \; \frac{e_{2u}\,e_{3u}}{e_{1u}} \;\delta _{i+1/2} [T] \right] 
-+ \delta _{j}\left[ A_v^{lT} \;  \frac{e_{1v}\,e_{3v}}{e_{2v}} \;\delta _{j+1/2} [T] \right]  \;\right)
+   \delta_{i}\left[ A_u^{lT} \; \frac{e_{2u}\,e_{3u}}{e_{1u}} \;\delta_{i+1/2} [T] \right] 
++ \delta_{j}\left[ A_v^{lT} \;  \frac{e_{1v}\,e_{3v}}{e_{2v}} \;\delta_{j+1/2} [T] \right]  \;\right)
 \end{equation}
 where  $b_t$=$e_{1t}\,e_{2t}\,e_{3t}$  is the volume of $T$-cells and
@@ -728,6 +728,6 @@
 \begin{equation} \label{eq:tra_sbc}
 \begin{aligned}
- &F^T = \frac{ 1 }{\rho _o \;C_p \,\left. e_{3t} \right|_{k=1} }  &\overline{ Q_{ns}       }^t  & \\ 
-& F^S =\frac{ 1 }{\rho _o  \,      \left. e_{3t} \right|_{k=1} }  &\overline{ \textit{sfx} }^t   & \\   
+ &F^T = \frac{ 1 }{\rho_o \;C_p \,\left. e_{3t} \right|_{k=1} }  &\overline{ Q_{ns}       }^t  & \\ 
+& F^S =\frac{ 1 }{\rho_o  \,      \left. e_{3t} \right|_{k=1} }  &\overline{ \textit{sfx} }^t   & \\   
  \end{aligned}
 \end{equation} 
@@ -743,8 +743,8 @@
 \begin{equation} \label{eq:tra_sbc_lin}
 \begin{aligned}
- &F^T = \frac{ 1 }{\rho _o \;C_p \,\left. e_{3t} \right|_{k=1} }   
+ &F^T = \frac{ 1 }{\rho_o \;C_p \,\left. e_{3t} \right|_{k=1} }   
            &\overline{ \left( Q_{ns} - \textit{emp}\;C_p\,\left. T \right|_{k=1} \right) }^t  & \\ 
 %
-& F^S =\frac{ 1 }{\rho _o \,\left. e_{3t} \right|_{k=1} } 
+& F^S =\frac{ 1 }{\rho_o \,\left. e_{3t} \right|_{k=1} } 
            &\overline{ \left( \;\textit{sfx} - \textit{emp} \;\left. S \right|_{k=1}  \right) }^t   & \\   
  \end{aligned}
@@ -1410,5 +1410,5 @@
       A linear interpolation is used to estimate $\widetilde{T}_k^{i+1}$,
       the tracer value at the depth of the shallower tracer point of the two adjacent bottom $T$-points.
-      The horizontal difference is then given by: $\delta _{i+1/2} T_k=  \widetilde{T}_k^{\,i+1} -T_k^{\,i}$ and
+      The horizontal difference is then given by: $\delta_{i+1/2} T_k=  \widetilde{T}_k^{\,i+1} -T_k^{\,i}$ and
       the average by: $\overline{T}_k^{\,i+1/2}= ( \widetilde{T}_k^{\,i+1/2} - T_k^{\,i} ) / 2$.
     }
@@ -1416,17 +1416,17 @@
 \end{figure}
 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-\begin{equation*}
+\[
 \widetilde{T}= \left\{  \begin{aligned}  
-&T^{\,i+1}      -\frac{ \left( e_{3w}^{i+1} -e_{3w}^i \right)}{ e_{3w}^{i+1} }\;\delta _k T^{i+1}	
+&T^{\,i+1}      -\frac{ \left( e_{3w}^{i+1} -e_{3w}^i \right)}{ e_{3w}^{i+1} }\;\delta_k T^{i+1}	
  								&& \quad\text{if  $\ e_{3w}^{i+1} \geq e_{3w}^i$   } 	\\
 										\\
-&T^{\,i} \ \ \ \,+\frac{ \left( e_{3w}^{i+1} -e_{3w}^i \right) }{e_{3w}^i       }\;\delta _k T^{i+1}
+&T^{\,i} \ \ \ \,+\frac{ \left( e_{3w}^{i+1} -e_{3w}^i \right) }{e_{3w}^i       }\;\delta_k T^{i+1}
 			 					&& \quad\text{if  $\ e_{3w}^{i+1}    <   e_{3w}^i$   } 
             \end{aligned}   \right.
-\end{equation*}
+\]
 and the resulting forms for the horizontal difference and the horizontal average value of $T$ at a $U$-point are: 
 \begin{equation} \label{eq:zps_hde}
 \begin{aligned}
- \delta _{i+1/2} T= 	\begin{cases}
+ \delta_{i+1/2} T= 	\begin{cases}
 \ \ \ \widetilde {T}\quad\ -T^i	 	& \ \ \quad\quad\text{if  $\ e_{3w}^{i+1} \geq e_{3w}^i$ } \\
 										\\
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex	(revision 10406)
@@ -1329,7 +1329,7 @@
 In the above formula, $h_{ab}$ denotes the height above bottom,
 $h_{wkb}$ denotes the WKB-stretched height above bottom, defined by
-\begin{equation*}
+\[
 h_{wkb} = H \, \frac{ \int_{-H}^{z} N \, dz' } { \int_{-H}^{\eta} N \, dz'  } \; ,
-\end{equation*}
+\]
 The $n_p$ parameter (given by \np{nn\_zpyc} in \ngn{namzdf\_tmx\_new} namelist)
 controls the stratification-dependence of the pycnocline-intensified dissipation.
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_conservation.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_conservation.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_conservation.tex	(revision 10406)
@@ -123,5 +123,5 @@
 
 \begin{equation} \label{eq:hpg_pe}
-\int_D {-\frac{1}{\rho _o }\left. {\nabla p^h} \right|_z \cdot {\textbf {U}}_h \;dv} \;=\;\int_D {\nabla .\left( {\rho \,{\textbf{U}}} \right)\;g\;z\;\;dv}
+\int_D {-\frac{1}{\rho_o }\left. {\nabla p^h} \right|_z \cdot {\textbf {U}}_h \;dv} \;=\;\int_D {\nabla .\left( {\rho \,{\textbf{U}}} \right)\;g\;z\;\;dv}
 \end{equation}
 
@@ -143,5 +143,5 @@
 surface pressure forces is exactly zero:
 \begin{equation} \label{eq:spg}
-\int_D {-\frac{1}{\rho _o }\nabla _h } \left( {p_s } \right)\cdot {\textbf{U}}_h \;dv=0
+\int_D {-\frac{1}{\rho_o }\nabla _h } \left( {p_s } \right)\cdot {\textbf{U}}_h \;dv=0
 \end{equation}
 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics.tex	(revision 10406)
@@ -68,5 +68,5 @@
             +\frac{1}{2}\nabla \left( {{\rm {\bf U}}^2} \right)}    \right]_h
  -f\;{\rm {\bf k}}\times {\rm {\bf U}}_h 
--\frac{1}{\rho _o }\nabla _h p + {\rm {\bf D}}^{\rm {\bf U}} + {\rm {\bf F}}^{\rm {\bf U}}
+-\frac{1}{\rho_o }\nabla _h p + {\rm {\bf D}}^{\rm {\bf U}} + {\rm {\bf F}}^{\rm {\bf U}}
   \end{equation}
   \begin{equation}     \label{eq:PE_hydrostatic}
@@ -570,5 +570,5 @@
    -   \frac{1}{2\,e_1}           \frac{\partial}{\partial i} \left(  u^2+v^2   \right) 
    -   \frac{1}{e_3    }  w     \frac{\partial u}{\partial k}      &      \\
-   -   \frac{1}{e_1    }            \frac{\partial}{\partial i} \left( \frac{p_s+p_h }{\rho _o}    \right)    
+   -   \frac{1}{e_1    }            \frac{\partial}{\partial i} \left( \frac{p_s+p_h }{\rho_o}    \right)    
    &+   D_u^{\vect{U}}  +   F_u^{\vect{U}}      \\
 \\
@@ -577,5 +577,5 @@
        -   \frac{1}{2\,e_2 }        \frac{\partial }{\partial j}\left(  u^2+v^2  \right)   
        -   \frac{1}{e_3     }   w  \frac{\partial v}{\partial k}     &      \\
-       -   \frac{1}{e_2     }        \frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho _o}  \right)    
+       -   \frac{1}{e_2     }        \frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho_o}  \right)    
     &+  D_v^{\vect{U}}  +   F_v^{\vect{U}}
 \end{split} \end{equation}
@@ -595,5 +595,5 @@
 		+			\frac{\partial \left( {e_1 \,v\,u} \right)}{\partial j}	\right)
                  - \frac{1}{e_3 }\frac{\partial \left( {         w\,u} \right)}{\partial k}    \\
--   \frac{1}{e_1 }\frac{\partial}{\partial i}\left( \frac{p_s+p_h }{\rho _o}   \right)
+-   \frac{1}{e_1 }\frac{\partial}{\partial i}\left( \frac{p_s+p_h }{\rho_o}   \right)
 +   D_u^{\vect{U}} +   F_u^{\vect{U}}
 \end{multline}
@@ -607,5 +607,5 @@
 		+			\frac{\partial \left( {e_1 \,v\,v} \right)}{\partial j}	\right)
                  - \frac{1}{e_3 } \frac{\partial \left( {        w\,v} \right)}{\partial k}    \\
--   \frac{1}{e_2 }\frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho _o}    \right)
+-   \frac{1}{e_2 }\frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho_o}    \right)
 +  D_v^{\vect{U}} +  F_v^{\vect{U}} 
 \end{multline}
@@ -771,6 +771,6 @@
 	-   \frac{1}{2\,e_1} \frac{\partial}{\partial i} \left(  u^2+v^2   \right) 
 	-   \frac{1}{e_3} \omega \frac{\partial u}{\partial k}       \\
-	-   \frac{1}{e_1} \frac{\partial}{\partial i} \left( \frac{p_s + p_h}{\rho _o}    \right)    
-	+  g\frac{\rho }{\rho _o}\sigma _1 
+	-   \frac{1}{e_1} \frac{\partial}{\partial i} \left( \frac{p_s + p_h}{\rho_o}    \right)    
+	+  g\frac{\rho }{\rho_o}\sigma _1 
 	+   D_u^{\vect{U}}  +   F_u^{\vect{U}} \quad
 \end{multline}
@@ -780,6 +780,6 @@
 	-   \frac{1}{2\,e_2 }\frac{\partial }{\partial j}\left(  u^2+v^2  \right)     	
 	-   \frac{1}{e_3 } \omega \frac{\partial v}{\partial k}         \\
-	-   \frac{1}{e_2 }\frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho _o}  \right) 
-	 +  g\frac{\rho }{\rho _o }\sigma _2   
+	-   \frac{1}{e_2 }\frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho_o}  \right) 
+	 +  g\frac{\rho }{\rho_o }\sigma _2   
 	+  D_v^{\vect{U}}  +   F_v^{\vect{U}} \quad
 \end{multline}
@@ -796,6 +796,6 @@
 		+			\frac{\partial \left( {e_1 \, e_3 \, v\,u} \right)}{\partial j}	\right)
    - \frac{1}{e_3 }\frac{\partial \left( { \omega\,u} \right)}{\partial k}    \\
-	- \frac{1}{e_1} \frac{\partial}{\partial i} \left( \frac{p_s + p_h}{\rho _o}    \right)    
-	+  g\frac{\rho }{\rho _o}\sigma _1 
+	- \frac{1}{e_1} \frac{\partial}{\partial i} \left( \frac{p_s + p_h}{\rho_o}    \right)    
+	+  g\frac{\rho }{\rho_o}\sigma _1 
 	+   D_u^{\vect{U}}  +   F_u^{\vect{U}} \quad
 \end{multline}
@@ -809,6 +809,6 @@
 		+			\frac{\partial \left( {e_1 \; e_3  \,v\,v} \right)}{\partial j}	\right)
                  - \frac{1}{e_3 } \frac{\partial \left( { \omega\,v} \right)}{\partial k}    \\
-	-   \frac{1}{e_2 }\frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho _o}  \right) 
-	 +  g\frac{\rho }{\rho _o }\sigma _2   
+	-   \frac{1}{e_2 }\frac{\partial }{\partial j}\left( \frac{p_s+p_h }{\rho_o}  \right) 
+	 +  g\frac{\rho }{\rho_o }\sigma _2   
 	+  D_v^{\vect{U}}  +   F_v^{\vect{U}} \quad
 \end{multline}
@@ -896,8 +896,8 @@
 $\textit{z*} = 0$ and  $\textit{z*} = -H$ respectively.
 Also the divergence of the flow field is no longer zero as shown by the continuity equation:
-\begin{equation*} 
+\[ 
 \frac{\partial r}{\partial t} = \nabla_{\textit{z*}} \cdot \left( r \; \rm{\bf U}_h \right)
 		\left( r \; w\textit{*} \right) = 0 
-\end{equation*} 
+\] 
 %}
 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics_zstar.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics_zstar.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics_zstar.tex	(revision 10406)
@@ -106,12 +106,12 @@
 The sea surface height is given by:
 \begin{equation} \label{eq:dynspg_ssh}
-\frac{\partial \eta }{\partial t}\equiv \frac{\text{EMP}}{\rho _w }+\frac{1}{e_{1T} 
-e_{2T} }\sum\limits_k {\left( {\delta _i \left[ {e_{2u} e_{3u} u} 
-\right]+\delta _j \left[ {e_{1v} e_{3v} v} \right]} \right)} 
+\frac{\partial \eta }{\partial t}\equiv \frac{\text{EMP}}{\rho_w }+\frac{1}{e_{1T} 
+e_{2T} }\sum\limits_k {\left( {\delta_i \left[ {e_{2u} e_{3u} u} 
+\right]+\delta_j \left[ {e_{1v} e_{3v} v} \right]} \right)} 
 \end{equation}
 
 where EMP is the surface freshwater budget (evaporation minus precipitation, and minus river runoffs
 (if the later are introduced as a surface freshwater flux, see \autoref{chap:SBC}) expressed in $Kg.m^{-2}.s^{-1}$,
-and $\rho _w =1,000\,Kg.m^{-3}$ is the volumic mass of pure water.
+and $\rho_w =1,000\,Kg.m^{-3}$ is the volumic mass of pure water.
 The sea-surface height is evaluated using a leapfrog scheme in combination with an Asselin time filter,
 i.e. the velocity appearing in (\autoref{eq:dynspg_ssh}) is centred in time (\textit{now} velocity). 
@@ -120,11 +120,11 @@
 \begin{equation} \label{eq:dynspg_exp}
 \left\{ \begin{aligned}
- - \frac{1}    						{e_{1u}} \;	\delta _{i+1/2} \left[  \,\eta\,  \right] 	\\
+ - \frac{1}    						{e_{1u}} \;	\delta_{i+1/2} \left[  \,\eta\,  \right] 	\\
  \\
- - \frac{1}    						{e_{2v}} \;	\delta _{j+1/2} \left[  \,\eta\,  \right]  
+ - \frac{1}    						{e_{2v}} \;	\delta_{j+1/2} \left[  \,\eta\,  \right]  
 \end{aligned} \right.
 \end{equation} 
 
-Consistent with the linearization, a $\left. \rho \right|_{k=1} / \rho _o$ factor is omitted in
+Consistent with the linearization, a $\left. \rho \right|_{k=1} / \rho_o$ factor is omitted in
 (\autoref{eq:dynspg_exp}). 
 
Index: /NEMO/trunk/doc/latex/NEMO/subfiles/chap_time_domain.tex
===================================================================
--- /NEMO/trunk/doc/latex/NEMO/subfiles/chap_time_domain.tex	(revision 10405)
+++ /NEMO/trunk/doc/latex/NEMO/subfiles/chap_time_domain.tex	(revision 10406)
@@ -163,5 +163,5 @@
 \begin{equation} \label{eq:STP_imp_zdf}
 \frac{T(k)^{t+1}-T(k)^{t-1}}{2\;\rdt}\equiv \text{RHS}+\frac{1}{e_{3t} }\delta 
-_k \left[ {\frac{A_w^{vT} }{e_{3w} }\delta _{k+1/2} \left[ {T^{t+1}} \right]} 
+_k \left[ {\frac{A_w^{vT} }{e_{3w} }\delta_{k+1/2} \left[ {T^{t+1}} \right]} 
 \right]
 \end{equation}
@@ -352,8 +352,8 @@
 \begin{flalign*}
 &\frac{\left( e_{3t}\,T \right)_k^{t+1}-\left( e_{3t}\,T \right)_k^{t-1}}{2\rdt}
-\equiv \text{RHS}+ \delta _k \left[ {\frac{A_w^{vt} }{e_{3w}^{t+1} }\delta _{k+1/2} \left[ {T^{t+1}} \right]} 
+\equiv \text{RHS}+ \delta_k \left[ {\frac{A_w^{vt} }{e_{3w}^{t+1} }\delta_{k+1/2} \left[ {T^{t+1}} \right]} 
 \right]      \\
 &\left( e_{3t}\,T \right)_k^{t+1}-\left( e_{3t}\,T \right)_k^{t-1}
-\equiv {2\rdt} \ \text{RHS}+ {2\rdt} \ \delta _k \left[ {\frac{A_w^{vt} }{e_{3w}^{t+1} }\delta _{k+1/2} \left[ {T^{t+1}} \right]} 
+\equiv {2\rdt} \ \text{RHS}+ {2\rdt} \ \delta_k \left[ {\frac{A_w^{vt} }{e_{3w}^{t+1} }\delta_{k+1/2} \left[ {T^{t+1}} \right]} 
 \right]      \\
 &\left( e_{3t}\,T \right)_k^{t+1}-\left( e_{3t}\,T \right)_k^{t-1}
