\begin{figure}
\footnotesize
\begin{tabular}{c@{\extracolsep{1mm}}c@{\extracolsep{1mm}}c}
\textit{Benchmark 1:} &\textit{Benchmark 2:} &\textit{Benchmark 3:} \\
%
\epsfig{file=../plots/01bench/ap1t2r351806m1r129l128j257/12.um.11.epsi,width=0.332\columnwidth,clip=}&
\epsfig{file=../plots/01bench/ap1t1r80000m1p50Nrho3D/12.um.19.epsi,width=0.332\columnwidth,clip=}&
\epsfig{file=../plots/01bench/ae035p2t40r25000000m2p2n2Nrho3/12.um.47.epsi,width=0.332\columnwidth,clip=}\\
%
\epsfig{file=../plots/01bench/ap1t2r351806m1r129l128j257/plot_pol.parms.plot.11.epsi,width=0.332\columnwidth,clip=}& 
\epsfig{file=../plots/01bench/ap1t1r80000m1p50Nrho3D/plot_pol.parms.plot.19.epsi,width=0.332\columnwidth,clip=}&
\epsfig{file=../plots/01bench/ae035p2t40r25000000m2p2n2Nrho3/plot_pol.parms.plot.47.epsi,width=0.332\columnwidth,clip=}\\
%
\raisebox{0.15\columnwidth}{--}&
\epsfig{file=../plots/01bench/ap1t1r80000m1p50Nrho3D/plot_4.parms.plot.19.epsi,width=0.332\columnwidth,clip=} &
\epsfig{file=../plots/01bench/ae035p2t40r25000000m2p2n2Nrho3/plot_4.parms.plot.47.epsi,width=0.332\columnwidth,clip=} \\
%
\normalsize
\end{tabular}
\caption{Solutions structures of benchmark cases 1, 2 and 3 (left to
right). The first plot in each column shows azimuthally-averaged 
isocontours of $\overline{u}_\varphi$ (left half) and of the
streamlines $r\sin\theta(\partial_\theta \overline{\pol})$ (right
half) in the meridional plane. The second plot in each column shows
isocontours of $u_r$ in the equatorial plane. The third plot in each
column shows isocontours of $B_r$ at $r=r_o$. }
\label{fig01}
\end{figure}



%$\eta=$, $n=$, $N_\rho=$, $k=$, $\R=$, $\Pr=$ $\Pm=$, $\tau$,
%and stress-free velocity BCs,
%fixed-value entropy BC,
%and insulating magnetic field BCs.
