Next: advdif_spm.F
Up: The Proudman Oceanographic Laboratory
Previous: Subroutines
Subsections
aaksetr.F
subroutine aaksetr
Sets the eddy viscosity and eddy diffusivity from a a Richardson number
fromulation.
It has a mixing length scale l to be assumed to be a fraction of the water depth
in this case l = 0.05H and the richardson number is:
A basic value of aa or ak resulting from surface or bed stress is assumed to be given
by , where is the appropiate friction velocity. Then this is reduced
liniearly until a cutoff value (Ri = 0.23) is reached. Hence for the bed stress
contribution to aa

(36) 
with similar formulation for the surface stress contribution, which decreases monotonically
from the surface downwards.The contributions from surface and bed stress are then added
together. A minimum value 0.002 ms is assumed for due to surface stress,
and a minimum value of 10 is used for aa and ak.
none
 ff Fraction of
 i,icg,ipn,j,jcg,k,ki Grid indeces.
 tagb,tagu,tagv,tagfb,tagfs,taggb,taggs Message passing taggs.
 alb,als Length scale for surface and bed stress.
 ustb,usts,ustl for bottom and surface and
 wsp3 Wind speed cubed.
 ri,rif Richardson number and Richardson number function.
 ab,as bottom and surface contributions for aa or ak
 aa Eddy viscosity
 ak Eddy diffusivity
 none
 Exhanges variables needed between processors.
 Calculates Ri at u points and makes sure that is positive.
 Calculates rif.
 Calculates for bottom and surface.
 Calculates the mixing length for surface and bottom.
 Calculates aa due to bed stress decreasing upwards.
 Calculates aa due to surface stress decreasing downwards.
 Combine the surface stress and bottom stress contribution of aa.
 Calculates Ri at b points and makes sure that is positive.
 Calculates rif.
 Calculates for bottom and surface.
 Calculates the mixing length for surface and bottom.
 Calculates ak due to bed stress decreasing upwards.
 Calculates ak due to surface stress decreasing downwards.
 Combine the surface stress and bottom stress contribution of ak.
 Exch3DS,Exch3DR
 b3drun.F
 none
 none
 Not complitely parallelised
Next: advdif_spm.F
Up: The Proudman Oceanographic Laboratory
Previous: Subroutines
The AMMP Project
20050420