source: trunk/documentation/disvert.tex @ 109

Last change on this file since 109 was 109, checked in by slebonnois, 14 years ago

SLebonnois: discretisation verticale: cohabitation entre
la methode Terre et les autres.

File size: 2.9 KB
Line 
1\documentclass[a4paper,10pt]{article}
2%\usepackage{graphicx}
3\usepackage{natbib}  % si appel à bibtex
4%\usepackage[francais]{babel}
5%\usepackage[latin1]{inputenc}  % accents directs (é...), avec babel
6%\usepackage{rotating}
7
8\setlength{\hoffset}{-1.in}
9\setlength{\oddsidemargin}{3.cm}
10\setlength{\textwidth}{15.cm}
11\setlength{\marginparsep}{0.mm}
12\setlength{\marginparwidth}{0.mm}
13
14\setlength{\voffset}{-1.in}
15\setlength{\topmargin}{0.mm}
16\setlength{\headheight}{0.mm}
17\setlength{\headsep}{30.mm}
18\setlength{\textheight}{24.cm}
19\setlength{\footskip}{1.cm}
20
21\setlength{\parindent}{0.mm}
22\setlength{\parskip}{1 em}
23\newcommand{\ten}[1]{$\times 10^{#1}$~} 
24\renewcommand{\baselinestretch}{1.}
25
26\begin{document}
27\pagestyle{plain}
28
29\begin{center}
30{\bf \LARGE 
31Documentation for LMDZ, Planets version
32
33\vspace{1cm}
34\Large
35The vertical discretization
36}
37
38\vspace{1cm}
39S\'ebastien Lebonnois
40
41\vspace{1cm}
42Latest version: \today
43\end{center}
44
45
46\section{Theoretical aspects}
47
48The position of the layers:
49\begin{itemize}
50\item pressure limit between two layers,
51\item pressure within the layers
52\end{itemize}
53
54The Exner function: definition.
55It corresponds to the pressure levels within the layers.
56Used for the computation of the potential temperature.
57For the Earth, we use a specific scheme that computes these positions so that
58it maintains a condition of proportionality between total,
59internal and potential energy (cf. a note from F. Hourdin).
60
61\section{Pratical aspects in the code}
62
63\begin{itemize}
64\item \textsf{disvert\_[no]terre.F[90]}:
65position of the interface pressure levels from an input file
66(several possibilities).
67Definition of ap, bp and presnivs.
68In the planetary version, definition of aps and bps.
69
70This is done only once, called at the beginning from \textsf{iniconst.F}.
71
72\item Interface pressures:
73computed in \textsf{caldyn0.F, caldyn.F, integrd.F, leapfrog.F}
74through the \textsf{pression.F} routine.
75
76\item Exner function (and therefore pressure within the layers):
77computed at three different places in \textsf{leapfrog.F} through the
78\textsf{exner\_[hyb/milieu].F} routine.
79For the Earth, we use \textsf{exner\_hyb.F}, that computes the positions in a
80specific way to maintain a condition of proportionality between total,
81internal and potential energy (cf. a note from F. Hourdin).
82For other planets, we use \textsf{exner\_milieu.F}, that computes the positions
83of these pressure levels exactly in the middle of each layer.
84Though this fails to maintain the previous condition, there is no evidence of
85any significant influence on the results, and it makes it a lot easier to
86define correctly the level positions with the input file.
87\end{itemize}
88
89%\begin{thebibliography}{2}
90%\providecommand{\natexlab}[1]{#1}
91%\expandafter\ifx\csname urlstyle\endcsname\relax
92%  \providecommand{\doi}[1]{doi:\discretionary{}{}{}#1}\else
93%  \providecommand{\doi}{doi:\discretionary{}{}{}\begingroup
94%  \urlstyle{rm}\Url}\fi
95
96%\end{thebibliography}
97
98\end{document}
Note: See TracBrowser for help on using the repository browser.