\documentclass[12pt]{article}
\usepackage{jas}
\usepackage{a4}
\usepackage{epsf}
\usepackage{graphicx}
\usepackage[round]{natbib}
\usepackage{equations}

\textheight 29cm
\pagestyle{empty}
\topmargin -3.cm
\textwidth 18cm
\oddsidemargin -1cm

\def\der#1#2{\frac{\partial #1}{\partial #2}}
\def\dep#1{\left(#1\right)}
\def\depb#1{\left[#1\right]}
\def\dem{1/2}
\def\eq#1{Eq.~\ref{eq:#1}}
\def\fg#1{Fig.~\ref{fg:#1}}
\def\dq{{\dep{\delta q}}}
\def\dt{\delta t}
\def\eg{e.~g.}


\begin{document}

\bibliographystyle{jas}


\large
\thispagestyle{empty}

\clearpage
\pagenumbering{arabic}

\begin{figure}
% Viewport original de uvref 44 80 586 239
\centerline{\includegraphics[width=17cm,trim=0 10 0 0]{\EPS/uvref.eps}}
%viewport original de uv.eps 42 63 586 539
\centerline{\includegraphics[width=17cm]{\EPS/uv.eps}}
\caption{Test d'advection horizontal avec un champ de vent
calculé à partir d'une fonction de courant analytique:\label{fg:w}}
$\Psi= a U_0 \cos^2\frac{\lambda}{2}\cos^2\phi$
($a$ rayon planétaire, $U_0$ une vitesse de référence, $\lambda$ la
longitude et $\phi$ la latitude).
Les vents correspondants sont représentés sur la figure en haut à gauche.
La distribution initiale de traceurs (également en haut à gauche) est
une gaussienne en fonction de la longitude sans variations latitudinales.
Les autres graphiques du haut montrent la même distribution advectée
analytiquement aux temps $T/2$ et $T$ ($T$ étant la période de rotation
au centre du graphe).
Les graphiques du bas montrent les résultats de l'advection avec les
différents schémas de transport pour trois résolutions spatiales~:
$120\times 60$ points, 1/2 avec $60 \times 30$
et 1/3 avec $40\times 20$ en longiutde-latitude.
\end{figure}

\newpage

\begin{figure}
\centerline{\includegraphics[width=18.cm]{\EPS/new1.eps}}
\caption{Advection d'un traceur idéalisé avec source purement
continentale sur toute la colonne d'atmosphère et une décroissance
radioactive avec une constante de temps de 10 jours.
Simulation LMD-Z avec une grille de 3 degrés par 2 degrés
et 20 niveaux verticaux.}
%``Res 1$/$2'' signifie que les mailles ont
été regroupées 2 par 2 dans les deux directions horizontales pour
calculer le transport du traceur.
Le champ de traceur est tracé au jour 20 à 500 mbar.
\end{figure}


\newpage

\begin{figure}
\centerline{\includegraphics[width=20cm]{\EPS/evol.eps}}
\caption{Evolution temporelle sur 40 jours de l'erreur
quadratique moyenne calculée avec comme référence le résultat de l'advection
avec le schéma de Prather sur la grille la plus fine.
Les simulations sont les mêmes que pour la Fig.~2 et les
cartes correspondent au niveau 500 mbar.}

La figure de gauche montre les erreurs calculées sur la grille la plus fine.\\
La figure de droite montre les erreurs calculées sur la grille
%%``Res. 1$/$2''.
}
\end{figure}

\newpage

\begin{figure}
\vspace{-1.5cm}
\centerline{\includegraphics[width=14.cm]{\EPS/alex.eps}}
\caption{Distribution moyenne sur 2 ans de la distribution de $^{222}Rn$
dans des simulations effectuées avec
le modèle LMD-cycle~5 (A: Schéma de pentes) et le modèle LMDZ
(B: Van Leer I et C: Godunov) avec contrôle de l'émission par la température
du sol et une grille identique de 64 par 50 point en sinus de la latitude.}
\end{figure}

\end{document}
