s_phys_pkgs/text/gridalt.tex

\section{Gridalt - Alternate Grid Package}
\label{sec:pkg:gridalt}
\begin{rawhtml}
<!-- CMIREDIR:package_gridalt: -->
\end{rawhtml}

\subsection {Introduction} 

The gridalt package was developed to allow the high end atmospheric physics 
(fizhi) physics to be run on a separate grid from the hydrodynamics. The package
could (with some user modification) be used in conjunction with other packages 
or for other calculations within the GCM. For the case of the atmospheric 
physics, a modified $p^*$ coordinate, which adds additional levels between
the lower levels of the existing $p^*$ grid (and perhaps between the levels near 
the tropopause as well), is implemented. The vertical discretization is
different for each grid point, although it consist of the same number of 
levels. This is illustrated as follows:
\begin{figure}[htbp]
\vspace*{-0.4in}
\begin{center}
\includegraphics[height=2.4in]{part6/vertical.eps}
\end{center}
\end{figure}

\vspace*{-0.5in}
In addition to computing the physical forcing terms of the momentum,
thermodynamic and humidity equations on the modified (higher resolution)
grid, the higher resolution structure of the atmosphere (the boundary
layer) is retained between calculations. This neccessitates a second 
set of evolution equations for the atmospheric state variables on the 
modified grid. If the equations for the evolution of the state
on $p^*$ can be expressed as:
\[
\left . {\partial U \over {\partial t}} \right |_{p^*}^{total} = 
\left . {\partial U \over {\partial t}} \right |_{p^*}^{dynamics} + 
\left . {\partial U \over {\partial t}} \right |_{p^*}^{physics}
\]
where the physics forcing terms on $p^*$ have been computed from a
mapping from the modified grid, then an additional set of equations
to govern the evolution of $U$ on the modified grid are written:
\[
\left . {\partial U \over {\partial t}} \right |_{p^{*m}}^{total} = 
\left . {\partial U \over {\partial t}} \right |_{p^{*m}}^{dynamics} + 
\left . {\partial U \over {\partial t}} \right |_{p^{*m}}^{physics} +
\gamma ({\left . U \right |_{p^*}} - {\left . U \right |_{p^{*m}}})
\]
where $p^{*m}$ refers to the modified higher resolution grid, and
the dynamics forcing terms have been mapped from the $p^*$ space.
The last term on the RHS is a relaxation term, meant to constrain
the state variables on the modified vertical grid to `track' the
state variables on the $p^*$ grid on some time scale, $\gamma$.

\subsection {Key subroutines, parameters and files } 

\subsection {Dos and donts}

In the context of a Held-Suarez type of model experiment (located
in the fizhi-gridalt-hs verification experiment) with
topography, the forcing terms which represent the physics are computed on 
the modified grid. The forcing terms are computed as functions of the
state variables on the modified grid. The tendencies are then interpolated
to the standard grid

\subsection {Gridalt Reference} 
1	molod	1.1	\section{Gridalt - Alternate Grid Package}
2	edhill	1.2	\label{sec:pkg:gridalt}
3			\begin{rawhtml}
4			<!-- CMIREDIR:package_gridalt: -->
5			\end{rawhtml}
6	molod	1.1
7			\subsection {Introduction}
8
9			The gridalt package was developed to allow the high end atmospheric physics
10			(fizhi) physics to be run on a separate grid from the hydrodynamics. The package
11			could (with some user modification) be used in conjunction with other packages
12			or for other calculations within the GCM. For the case of the atmospheric
13			physics, a modified $p^*$ coordinate, which adds additional levels between
14			the lower levels of the existing $p^*$ grid (and perhaps between the levels near
15			the tropopause as well), is implemented. The vertical discretization is
16			different for each grid point, although it consist of the same number of
17			levels. This is illustrated as follows:
18			\begin{figure}[htbp]
19			\vspace*{-0.4in}
20			\begin{center}
21	molod	1.3	\includegraphics[height=2.4in]{part6/vertical.eps}
22	molod	1.1	\end{center}
23			\end{figure}
24
25			\vspace*{-0.5in}
26			In addition to computing the physical forcing terms of the momentum,
27			thermodynamic and humidity equations on the modified (higher resolution)
28			grid, the higher resolution structure of the atmosphere (the boundary
29			layer) is retained between calculations. This neccessitates a second
30			set of evolution equations for the atmospheric state variables on the
31			modified grid. If the equations for the evolution of the state
32			on $p^*$ can be expressed as:
33			\[
34			\left . {\partial U \over {\partial t}} \right \|_{p^*}^{total} =
35			\left . {\partial U \over {\partial t}} \right \|_{p^*}^{dynamics} +
36			\left . {\partial U \over {\partial t}} \right \|_{p^*}^{physics}
37			\]
38			where the physics forcing terms on $p^*$ have been computed from a
39			mapping from the modified grid, then an additional set of equations
40			to govern the evolution of $U$ on the modified grid are written:
41			\[
42			\left . {\partial U \over {\partial t}} \right \|_{p^{*m}}^{total} =
43			\left . {\partial U \over {\partial t}} \right \|_{p^{*m}}^{dynamics} +
44			\left . {\partial U \over {\partial t}} \right \|_{p^{*m}}^{physics} +
45			\gamma ({\left . U \right \|_{p^}} - {\left . U \right \|_{p^{m}}})
46			\]
47			where $p^{*m}$ refers to the modified higher resolution grid, and
48			the dynamics forcing terms have been mapped from the $p^*$ space.
49			The last term on the RHS is a relaxation term, meant to constrain
50			the state variables on the modified vertical grid to `track' the
51			state variables on the $p^*$ grid on some time scale, $\gamma$.
52
53			\subsection {Key subroutines, parameters and files }
54
55			\subsection {Dos and donts}
56
57			In the context of a Held-Suarez type of model experiment (located
58	molod	1.3	in the fizhi-gridalt-hs verification experiment) with
59	molod	1.1	topography, the forcing terms which represent the physics are computed on
60			the modified grid. The forcing terms are computed as functions of the
61			state variables on the modified grid. The tendencies are then interpolated
62			to the standard grid
63
64			\subsection {Gridalt Reference}