drydepo_mod Module

Arguments

Return Value real

Arguments

Arguments

Belward, A.S., Estes, J.E., and Kline, K.D., 1999, * The IGBP-DIS 1-Km Land-Cover Data Set DISCover: * A Project Overview: Photogrammetric Engineering and Remote Sensing , * v. 65, no. 9, p. 1013-1020 * * if there are no data in the inventory * the ECMWF land/sea mask is used to distinguish * between sea (-> ocean) and land (-> grasslands). * * Author: A. Stohl * * 5 December 1996 * 8 February 1999 Additional use of nests, A. Stohl * 29 December 2006 new landuse inventory, S. Eckhardt *

Arguments

Arguments

Variables: * * nunc uncertainty class of the respective particle * nage age class of the respective particle * deposit amount (kg) to be deposited * *

Arguments

Variables: * * nunc uncertainty class of the respective particle * nage age class of the respective particle * deposit amount (kg) to be deposited * *

Arguments

Variables: * alpha help variable * cun 'slip-flow' correction after Cunningham * d01 [um] upper diameter * d02 [um] lower diameter * dc [m2/s] coefficient of Brownian diffusion * delta distance given in standard deviation units * density [kg/m3] density of the particle * dmean geometric mean diameter of interval * dquer [um] geometric mass mean particle diameter * dsigma e.g. dsigma=10 or dsigma=0.1 means that 68% of the mass * are between 0.1dquer and 10dquer * fract(ni) mass fraction of each diameter interval * kn Knudsen number * ni number of diameter intervals, for which deposition * is calculated * schmidt Schmidt number * schmi schmidt**2/3 * vsh [m/s] gravitational settling velocity of the particle * x01 normalized upper diameter * x02 normalized lower diameter * * Constants: * g [m/s2] Acceleration of gravity * kb [J/K] Stefan-Boltzmann constant * lam [m] mean free path of air molecules * myl [kg/m/s] dynamical viscosity of air * nyl [m2/s] kinematic viscosity of air * tr reference temperature * * Function: * erf calculates the integral of the Gauss function * *

Arguments

Calculation of the probability for dry deposition * * Particle positions are read in - prob returned * *

Arguments

Arguments

This routine calculates the dry deposition velocities. * * Author: A. Stohl * * 20 December 1996 * Sabine Eckhardt, Jan 07 * if the latitude is negative: add half a year to the julian day * *

Arguments

This routine calculates the dry deposition velocities. * * Author: A. Stohl * * 20 December 1996 * Sabine Eckhardt, Jan 07 * if the latitude is negative: add half a year to the julian day * *

Arguments

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

Arguments

Calculation of the quasilaminar sublayer resistance to dry deposition. * * AUTHOR: Andreas Stohl, 20 May 1995 * *

Arguments

Calculation of the surface resistance according to the procedure given * in: * Wesely (1989): Parameterization of surface resistances to gaseous * dry deposition in regional-scale numerical models. * Atmos. Environ. 23, 1293-1304. * * * AUTHOR: Andreas Stohl, 19 May 1995 * *

Arguments