| 364 |
$\chi_{ob}$ (note: passive tracers are currently not implemented and |
$\chi_{ob}$ (note: passive tracers are currently not implemented and |
| 365 |
the code stops when package \code{ptracers} is used together with this |
the code stops when package \code{ptracers} is used together with this |
| 366 |
option). Currently, the code vertically averages the normal velocity |
option). Currently, the code vertically averages the normal velocity |
| 367 |
as specified. From these prescribed values the code computes the |
as specified in \code{OB[E,W]u} or \code{OB[N,S]v}. From these |
| 368 |
boundary values for the next timestep $n+1$ as follows (as an |
prescribed values the code computes the boundary values for the next |
| 369 |
example, we use the notation for an eastern or western boundary): |
timestep $n+1$ as follows (as an example, we use the notation for an |
| 370 |
|
eastern or western boundary): |
| 371 |
\begin{itemize} |
\begin{itemize} |
| 372 |
\item $u^{n+1}(y,z) = \bar{u}_{ob}(y) + u'(y,z)$, where $u_{n}'$ is the |
\item $u^{n+1}(y,z) = \bar{u}_{ob}(y) + (u')^{n}(y,z)$, where $(u')^{n}$ |
| 373 |
deviation from the vertically averaged velocity one grid point |
is the deviation from the vertically averaged velocity at timestep |
| 374 |
inward from the boundary. |
$n$ one grid point inward from the boundary. |
| 375 |
\item If $u^{n+1}$ is directed into the model domain, the boudary |
\item If $u^{n+1}$ is directed into the model domain, the boudary |
| 376 |
value for tracer $\chi$ is restored to the prescribed values: |
value for tracer $\chi$ is restored to the prescribed values: |
| 377 |
\[\chi^{n+1} = \chi^{n} + \frac{\Delta{t}}{\tau_\chi} (\chi_{ob} - |
\[\chi^{n+1} = \chi^{n} + \frac{\Delta{t}}{\tau_\chi} (\chi_{ob} - |
| 378 |
\chi^{n}),\] where $\tau_\chi$ is the relaxation time |
\chi^{n}),\] where $\tau_\chi$ is the relaxation time |
| 379 |
scale \texttt{T/SrelaxStevens}. |
scale \texttt{T/SrelaxStevens}. The new $\chi^{n+1}$ is then subject |
| 380 |
\item If $u^{n+1}$ is directed out of the model domain, the tracer is |
to the advection by $u^{n+1}$. |
| 381 |
advected out of the domain with $u^{n+1}+c$, where $c$ is a phase |
\item If $u^{n+1}$ is directed out of the model domain, the tracer |
| 382 |
velocity estimated as |
$\chi^{n+1}$ on the boundary at timestep $n+1$ is estimated from |
| 383 |
$\frac{1}{2}\frac{\partial\chi}{\partial{t}}/\frac{\partial\chi}{\partial{x}}$. |
advection advected out of the domain with $u^{n+1}+c$, where $c$ is |
| 384 |
|
a phase velocity estimated as |
| 385 |
|
$\frac{1}{2}\frac{\partial\chi}{\partial{t}}/\frac{\partial\chi}{\partial{x}}$. The |
| 386 |
|
numerical scheme is (as an example for an eastern boundary): |
| 387 |
|
\[\chi_{i,j,k}^{n+1} = \chi_{i,j,k}^{n} + \Delta{t} |
| 388 |
|
(u^{n+1}+c)_{i_{b},j,k}\frac{\chi_{i_{b},j,k}^{n} |
| 389 |
|
- \chi_{i_{b}-1,j,k}^{n}}{\Delta{x}_{i_{b},j}^{C}}\mbox{, if }u_{i_{b},j,k}^{n+1}>0, |
| 390 |
|
\] where $i_{b}$ is the boundary index. |
| 391 |
|
|
| 392 |
For test purposes, the phase velocity contribution or the entire |
For test purposes, the phase velocity contribution or the entire |
| 393 |
advection can |
advection can be turned off by setting the corresponding parameters |
|
be turned off by setting the corresponding parameters |
|
| 394 |
\texttt{useStevensPhaseVel} and \texttt{useStevensAdvection} to |
\texttt{useStevensPhaseVel} and \texttt{useStevensAdvection} to |
| 395 |
\texttt{.FALSE.}.\end{itemize} See \citet{stevens:90} for details. |
\texttt{.FALSE.}.\end{itemize} See \citet{stevens:90} for details. |
| 396 |
|
|
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch