\def\legsz{\footnotesize}
\def\legendes{\legsz $\theta K$ &\legsz $w'\theta'$ &\legsz B &\legsz C &\legsz u (m~s$^{-1}$) &\legsz v (m~s$^{-1}$)\normalsize}
\def\incun#1{\includegraphics[width=3.4cm,angle=-90]{#1.epsi}}
\def\incn#1#2{#2 &
\incun{#1/#2/temp} &
\incun{#1/#2/flux} &
\incun{#1/#2/B} &
\incun{#1/#2/C} &
\incun{#1/#2/u} &
\incun{#1/#2/v} \\
& \legendes}

\begin{figure*}
\begin{center}
\begin{tabular}{ccccccc}
\incn{\thermdir/Kref}{05WC}\\
\incn{\thermdir/Kref}{05SC}\\
\incn{\thermdir/Kref}{24SC}\\
\incn{\thermdir/Kref}{24B}
\end{tabular}
\end{center}
\caption{Comparison of LES profiles with the results of the
thermals model with nominal values of the parameters
($\lambda=20$, $\mu=2$, $r=2$ and $\gamma=0$)
for 4 cases (05WC, 05SC, 24SC and 24B).
For each case, the initial profile (dashed) and the average
between times $t_1$ and $t_2$ for the LES (thin solid curve) and for the
parameterization (heavy solid curve) are shown.
For the heat flux, there is no initial value.
\label{fg:profref}}
\end{figure*}