\begin{figure}
 %\textcolor{magenta}{*** Figure 3 ***}
 %\textcolor{magenta}{*** Left side of Figure 2 of the CFD4NRS-7 paper ***}
%   \begin{tabular}{cc}
%      \adjincludegraphics[width=60mm,valign=c]{figures/figure03/figure3-320by320-Code_Saturne-LiquidFraction.pdf}\label{figure:03:a} 
%    &
%      \adjincludegraphics[width=87.5mm,valign=c]{figures/figure03/contoursatepsiloncomparison.pdf}\label{figure:03:b}
%  \end{tabular}
   \begin{center}
      \includegraphics[width=0.6\textwidth]{figures/figure03/figure3-320by320-Code_Saturne-LiquidFraction.pdf}\label{figure:03:a} \\
      \includegraphics[width=0.8\textwidth]{figures/figure03/contoursatepsiloncomparison.pdf}\label{figure:03:b}
  \end{center}


 \caption[
    Comparison between between results obtained by \CSAT~ with the study of \citep{Voller:1987}. 
    The field at the top is the liquid fraction.  
    The plot below gives the contours at $T = \pm\epsilon=\pm$\num{0.1}, which indicate the boundaries of the mushy area. Note that one in every twenty data points is plotted for the 320$\times$320 mesh.
   ]{}\label{fig:03} % first figure in \S\,\ref{casestud:02}
\end{figure}
