s_phys_pkgs/text/gridalt.tex

\section{Gridalt - Alternate Grid Package}

\subsection {Introduction} 
To take advantage of a `high end' turbulence parameterization
(and convection parameterization), the vertical resolution near the surface
must be increased substantially as compared to the vertical resolution needed
aloft. This cannot be accomplished if the high end physics is computed using 
the $p^*$ coordinate currently in use in the MIT gcm. 

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