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:

$$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