do_rivers - updating river points

INTERFACE:

    subroutine do_rivers(do_3d)
DESCRIPTION:

Here, the temperature, salinity, sea surface elevation and layer heights are updated in the river inflow grid boxes. Temperature and salinity are mixed with riverine values proportional to the old volume and the river inflow volume at that time step, sea surface elevation is simply increased by the inflow volume divided by the grid box area, and the layer heights are increased proportionally. USES:

    IMPLICIT NONE
INPUT PARAMETERS:
    logical, intent(in)                 :: do_3d
LOCAL VARIABLES:
    integer                   :: i,j,k,m,n
    integer                   :: kl,kh
    integer, save             :: nn=0
    REALTYPE                  :: ramp=_ONE_
    REALTYPE                  :: rvol,height
    REALTYPE                  :: river_depth,x