*
Calculation of the dry deposition velocities of particles. *
This routine is based on Stokes' law for considering settling and *
assumes constant dynamic viscosity of the air. *
*
AUTHOR: Andreas Stohl, 12 November 1993 *
Update: 20 December 1996 *
*
Literature: *
[1] Hicks/Baldocchi/Meyers/Hosker/Matt (1987), A Preliminary *
Multiple Resistance Routine for Deriving Dry Deposition *
Velocities from Measured Quantities. *
Water, Air and Soil Pollution 36 (1987), pp.311-330. *
[2] Slinn (1982), Predictions for Particle Deposition to *
Vegetative Canopies. Atm.Env.16-7 (1982), pp.1785-1794. *
[3] Slinn/Slinn (1980), Predictions for Particle Deposition on *
Natural Waters. Atm.Env.14 (1980), pp.1013-1016. *
[4] Scire/Yamartino/Carmichael/Chang (1989), *
CALGRID: A Mesoscale Photochemical Grid Model. *
Vol II: User's Guide. (Report No.A049-1, June, 1989) *
[5] Langer M. (1992): Ein einfaches Modell zur Abschaetzung der *
Depositionsgeschwindigkeit von Teilchen und Gasen. *
Internal report. *
*
*
Variables: * alpha help variable * fract(nc,ni) mass fraction of each diameter interval * lpdep(nc) 1 for particle deposition, 0 else * nc actual number of chemical components * ni number of diameter intervals, for which vdepj is calc. rdp [s/m] deposition layer resistance * ra [s/m] aerodynamical resistance * schmi(nc,ni) Schmidt number*2/3 of each diameter interval * stokes Stokes number * ustar [m/s] friction velocity * vdep_tmp(nc) [m/s] deposition velocities of all components * vdepj [m/s] help, deposition velocity of 1 interval * vset(nc,ni) gravitational settling velocity of each interval * * Constants: * nc number of chemical species * ni number of diameter intervals, for which deposition * is calculated * *
Use this stokes number for different shapes
Use this stokes number for different shapes
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer | :: | nc | ||||
| real, | intent(in) | :: | density(maxspec) | |||
| real, | intent(in) | :: | fract(maxspec,maxndia) | |||
| real | :: | schmi(maxspec,maxndia) | ||||
| real, | intent(in) | :: | vset(maxspec,maxndia) | |||
| real, | intent(in) | :: | ra | |||
| real, | intent(in) | :: | ustar | |||
| real, | intent(in) | :: | nyl | |||
| real, | intent(in) | :: | rhoa | |||
| real, | intent(inout) | :: | vdep_tmp(maxspec) |