1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
!**********************************************************************
! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010 *
! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, *
! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann *
! *
! This file is part of FLEXPART. *
! *
! FLEXPART is free software: you can redistribute it and/or modify *
! it under the terms of the GNU General Public License as published by*
! the Free Software Foundation, either version 3 of the License, or *
! (at your option) any later version. *
! *
! FLEXPART is distributed in the hope that it will be useful, *
! but WITHOUT ANY WARRANTY; without even the implied warranty of *
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
! GNU General Public License for more details. *
! *
! You should have received a copy of the GNU General Public License *
! along with FLEXPART. If not, see <http://www.gnu.org/licenses/>. *
!**********************************************************************
subroutine partdep(nc,density,fract,schmi,vset,ra,ustar,nyl,vdep)
! i i i i i i i i i/o
!*****************************************************************************
! *
! 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(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 par_mod
implicit none
real :: density(maxspec),schmi(maxspec,ni),fract(maxspec,ni)
real :: vset(maxspec,ni)
real :: vdep(maxspec),stokes,vdepj,rdp,ustar,alpha,ra,nyl
real,parameter :: eps=1.e-5
integer :: ic,j,nc
do ic=1,nc ! loop over all species
if (density(ic).gt.0.) then
do j=1,ni ! loop over all diameter intervals
if (ustar.gt.eps) then
! Stokes number for each diameter interval
!*****************************************
stokes=vset(ic,j)/ga*ustar*ustar/nyl
alpha=-3./stokes
! Deposition layer resistance
!****************************
if (alpha.le.log10(eps)) then
rdp=1./(schmi(ic,j)*ustar)
else
rdp=1./((schmi(ic,j)+10.**alpha)*ustar)
endif
vdepj=vset(ic,j)+1./(ra+rdp+ra*rdp*vset(ic,j))
else
vdepj=vset(ic,j)
endif
! deposition velocities of each interval are weighted with mass fraction
!***********************************************************************
vdep(ic)=vdep(ic)+vdepj*fract(ic,j)
end do
endif
end do
end subroutine partdep
|
