!********************************************************************** ! 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 . * !********************************************************************** 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