advect_vel.F

subroutine advect_vel

Description

The subroutine calls advection routines htmlrefadvpuuadvpuu and htmlrefadvpuvadvpuv to do momentum advection. The order of advection changes on alternate time steps, such that on odd timesteps u advection is followed by v advection and on even timesteps v advection is followed by u advection.

Subroutine Arguments

oddstep logical variable telling whether timestep is odd

Local variables

i,j,k grid indices

u2 temporary storage for intermediate velocities

Logical units

none

Global variables changed

u, v

Order of Things

1. Initialise depths
2. Call momentum advection routines in order determined by odd or even timestep

Calls

advpuu, advpuv

b3ddbg (with debug option)

gsynchwrk

Called By

b3drun.F

Options - Logical

advect_v turns on northward advection.

Options - Compiler

NOPHYSADV turns off all physics advection

NOVELADV turns off velocity advection

DEBUG outputs debugging information

Known Issues

Limitations of directional splitting not resolved by alternating order call to gsynchwrk is likely to ineffient on some systems - not sure whether it is still neded.