next up previous
Next: advect_vel.F Up: The Proudman Oceanographic Laboratory Previous: advdif_spm.F

Subsections


advect_sca.F


subroutine advect_sca

Description

The subroutine calls advection routines advpbu and advpbv to advect the scalars temperature and/or salinity. 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.

Depth conservation can be checked in the debugging option to see that the new depth given by advection routines agrees with that given by the barotropic time steps.

Subroutine Arguments

oddstep logical variable telling whether timestep is odd

Local variables

i,j grid indices
idif,jdif position of maximum depth differences in conservation check
dif,diffg,diffl,hdiff used in depth conservation check

Logical units

Only used in debugging option.

Global variables changed

tmp, sal

Order of Things

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

Calls

advpbu
advpbv
b3ddbg (with debug option)
gmax1d
MPI_Gather, copy1di , subvol3d

Called By

b3drun

Options - Logical

no_tmp turns off temperature integration
no_sal turns off salinity integration
advect_v turns on northward advection.

Options - Compiler

NOPHYSADV turns off all physics advection
NOSCAADV turns off scalar advection
DEBUG outputs debugging information

Known Issues

Limitations of directional splitting not resolved by alternating order


next up previous
Next: advect_vel.F Up: The Proudman Oceanographic Laboratory Previous: advdif_spm.F
The AMMP Project 2005-04-20