do_residual - barotropic residual currents. (Source File: residual.F90)

INTERFACE:

    subroutine do_residual(finish)
DESCRIPTION:

Here, the residual transports and depths are integrated up every time step. At the end of the simulation, the Eulerian residual currents are calculated from:

$\displaystyle u_{res} =
\frac{\displaystyle\int_{t_0}^{t_1}U\,\mbox{d}\tau}
{\d...
...t_{t_0}^{t_1}V\,\mbox{d}\tau}
{\displaystyle\int_{t_0}^{t_1}D^v\,\mbox{d}\tau},$ (94)

where $t_0$ is the time when the residual calculation begins (to be chosen from namelist) and $t_1$ is the finishing time of the model simulation.

USES:

    use variables_2d, only: u,v,res_du,res_u,res_dv,res_v,du,dv
    IMPLICIT NONE
INPUT PARAMETERS:
    integer, intent(in)                 :: finish
REVISION HISTORY:
    Original author(s): Karsten Bolding & Hans Burchard