133 |
\begin{eqnarray} |
\begin{eqnarray} |
134 |
\label{EQ:model_equations} |
\label{EQ:model_equations} |
135 |
\frac{Du}{Dt} - fv + |
\frac{Du}{Dt} - fv + |
136 |
\frac{1}{\rho}\frac{\partial p^{'}}{\partial \lambda} - |
\frac{1}{\rho}\frac{\partial p^{\prime}}{\partial \lambda} - |
137 |
A_{h}\nabla_{h}^2u - A_{z}\frac{\partial^{2}u}{\partial z^{2}} |
A_{h}\nabla_{h}^2u - A_{z}\frac{\partial^{2}u}{\partial z^{2}} |
138 |
& = & |
& = & |
139 |
\cal{F} |
\cal{F} |
140 |
\\ |
\\ |
141 |
\frac{Dv}{Dt} + fu + |
\frac{Dv}{Dt} + fu + |
142 |
\frac{1}{\rho}\frac{\partial p^{'}}{\partial \varphi} - |
\frac{1}{\rho}\frac{\partial p^{\prime}}{\partial \varphi} - |
143 |
A_{h}\nabla_{h}^2v - A_{z}\frac{\partial^{2}v}{\partial z^{2}} |
A_{h}\nabla_{h}^2v - A_{z}\frac{\partial^{2}v}{\partial z^{2}} |
144 |
& = & |
& = & |
145 |
0 |
0 |
154 |
& = & |
& = & |
155 |
0 |
0 |
156 |
\\ |
\\ |
157 |
p^{'} & = & |
p^{\prime} & = & |
158 |
g\rho_{0} \eta + \int^{0}_{-z}\rho^{'} dz |
g\rho_{0} \eta + \int^{0}_{-z}\rho^{\prime} dz |
159 |
\\ |
\\ |
160 |
\rho^{'} & = &- \alpha_{\theta}\rho_{0}\theta^{'} |
\rho^{\prime} & = &- \alpha_{\theta}\rho_{0}\theta^{\prime} |
161 |
\\ |
\\ |
162 |
{\cal F} |_{s} & = & \frac{\tau_{x}}{\rho_{0}\Delta z_{s}} |
{\cal F} |_{s} & = & \frac{\tau_{x}}{\rho_{0}\Delta z_{s}} |
163 |
\\ |
\\ |
166 |
|
|
167 |
\noindent where $u$ and $v$ are the components of the horizontal |
\noindent where $u$ and $v$ are the components of the horizontal |
168 |
flow vector $\vec{u}$ on the sphere ($u=\dot{\lambda},v=\dot{\varphi}$). |
flow vector $\vec{u}$ on the sphere ($u=\dot{\lambda},v=\dot{\varphi}$). |
169 |
|
The terms $H\hat{u}$ and $H\hat{v}$ are the components of the term |
170 |
|
integrated in eqaution \ref{eq:free-surface}, as descirbed in section |
171 |
|
|
172 |
The suffices ${s},{i}$ indicate surface and interior of the domain. |
The suffices ${s},{i}$ indicate surface and interior of the domain. |
173 |
The forcing $\cal F$ is only applied at the surface. |
The forcing $\cal F$ is only applied at the surface. |
174 |
The pressure field $p^{'}$ is separated into a barotropic part |
The pressure field, $p^{\prime}$, is separated into a barotropic part |
175 |
due to variations in sea-surface height, $\eta$, and a hydrostatic |
due to variations in sea-surface height, $\eta$, and a hydrostatic |
176 |
part due to variations in density, $\rho^{'}$, over the water column. |
part due to variations in density, $\rho^{\prime}$, over the water column. |
177 |
|
|
178 |
\subsection{Discrete Numerical Configuration} |
\subsection{Discrete Numerical Configuration} |
179 |
|
|
182 |
$\Delta \lambda=\Delta \varphi=1^{\circ}$, so |
$\Delta \lambda=\Delta \varphi=1^{\circ}$, so |
183 |
that there are sixty grid cells in the zonal and meridional directions. |
that there are sixty grid cells in the zonal and meridional directions. |
184 |
Vertically the |
Vertically the |
185 |
model is configured with a four layers with constant depth, |
model is configured with four layers with constant depth, |
186 |
$\Delta z$, of $500$~m. The internal, locally orthogonal, model coordinate |
$\Delta z$, of $500$~m. The internal, locally orthogonal, model coordinate |
187 |
variables $x$ and $y$ are initialised from the values of |
variables $x$ and $y$ are initialised from the values of |
188 |
$\lambda$, $\varphi$, $\Delta \lambda$ and $\Delta \varphi$ in |
$\lambda$, $\varphi$, $\Delta \lambda$ and $\Delta \varphi$ in |