June 1973
Users' Guides
to the
Interactive Versions
of
Three Point Source Dispersion Programs:
PTMAX, PTDIS, and PTMTP
by '
D. Bruce Turner*
and
Adrian D. Busse*
Meteorology Laboratory
Program Element 21ADN
Notice
This document is a preliminary draft. It has not been formally
released by EPA and should not at this stage be construed to represent
Agency policy. It is being circulated for comment on its technical
accuracy and policy implications.
NATIONAL ENVIRONMENTAL RESEARCH CENTER
OFFICE OF RESEARCH AND MONITORING
U.S. ENVIRONMENTAL PROTECTION AGENCY
RESEARCH TRIANGLE PARK, N. C. 27711
*0n Assignment from the National Oceanic & Atmospheric Administration,
U.S. Department of Commerce
-------
Three steady—state Gaussian plume point source models have
recently been added to EPA ’s UNANAP (Users’ Network for Applied
Modeling of Air Pollution) and to the UNAMAP system available to
non—EPA users.
EPA users access UNAMAP models by entering “UNAMAP” as the
procedure name at log on. Non—EPA users may obtain information
as to access to UNANAP models by contacting Mr. Peter Loux,
Computer Sciences Corporation, 1701 North Port Meyer Drive,
Arlington, Va. 22209. Phone (703) 527—6080. On both systems
the command “TJNAMAP” after log on is accomplished, produces a
brief description of each model currently available.
Each of the three programs is briefly described including
a program abstract for each program.
-------
READY
unama p
‘TSO.UNAMAP(CATAL)’
UNAMAP- Users Network for Applied Modeling of Air Pollution
Contact: Dr. Ron Ruff, Chief, Computer Techniques Group
Division of Meteorology
Environmental Protection Agency
Research Triangle Park, North Carolina 27711
Phone 919/5’ 9—t 566
CATALOG OF PROGRAMS as of 06/01/73
APRAC — A short-term Urban fliffusion Model that calculates the
automotive contribution to Carbon Monoxide. The model
was developed by Stanford Research Institute (SRI).
A 120 page manual is available on the model.
HIWAY — A model that calculates a pollutant concentration
in the vicinity of a roadway. The model is
self-documenting.
PTMAX - An interactive program which performs an analysis of
the maximum, short—term concentration from a point
source as a function of stability and wind speed.
PTDIS — An interactive program which computes short-term
concentrations downwind from a point source at
distances specified by the user.
PThTP - An interactive program which computes, at multiple
receptors, short term concentrations resultin from
multiple point sources.
-------
Users’ Guide to PTMAX (the Interactive Version of DBT 52)
Program Abstract
PTMAX produces an analysis of maximum concentration
as the function of wind speed and stability. A separate
analysis is made for each individual stack. Input to
the program consists of ambient air temperature, and
characteristics of the source, such as emission rate,
physical stack height, and stack gas temperature.
Either the stack gas volume flow or both the stack
gas velocity and inside diameter at the top are also
required. Outputs of the program consist of effective
height of emission, maximum ground level concentration,
and distance of maximum concentration for each condition
of stability and wind speed.
This program determines for each wind speed and stability the
final plume rise using methods suggested by Briggs. This plume rise
is added to the physical stack height to determine the effective
height of emission. The effective height is used to determine both
the maximum concentration and the distance to maximum concentration.
The following assumptions are made: a steady—state Gaussian
plume model is applicable to determine ground level concentrations.
Computations can be performed according to the “Workbook of Atmospheric
Dispersion Estimates.” The dispersion parameter values used for the
horizontal dispersion coefficient, sigma y, and the vertical
dispersion coefficient, sigma z, are those given in Figures 3—2 and
3-3 of the workbook. The stated wind speed occurs at the stack top
for dilution of the plume and through the layer that the plume rise
occurs. The stated stability occurs from ground level to well above
the top of the plume. If there is a limit to vertical mixing, it
-------
2
occurs far enough above the top of the plume so that it has no influence
upon the maximum concentration. There are no topographic obstructions
in the vicinity of the source. The source exists in either flat or
gently rolling terrain.
Use of this program is applicable where single sources exist in
relatively uniform terrain. It is not applicable if aerodynamic
downwash around buildings in the vicliaity of the source affects the
plume emitted from the stack. The calculated concentrations are for
the single source considered. Where multiple stacks exist for a given
single plant this program can be applied to each individual stack. It
cannot give the maximum concentrations of the combination of the
stacks however. This program is useful in determining what combinations
of wind speed and stability produce maximum concentrations. For a
given stability the critical wind velocity, that is, the wind speed
that causes the maximum concentration, can be determined. This can
be done by seeing which wind speed produces the highest concentration
for that stability.
The use of the interactive version of the program is quite
simple. An alpha—numeric title is used to put a heading on the
output. The ambient air temperature in degrees—K is asked for.
If a zero is entered a value of 293 will be used. The stability
class is asked for next. If a zero is entered, all stabilities
1 through 6 are considered. Numbers 1 through 6 correspond to
stabilities A through F. Source strength of the pollutant considered
-------
3
in grams per second is asked for next. The physical stack height in
meters is then required. The stack gas temperature in degrees—IC is
next entered. If the volume flow in cubic meters per second is known,
it is next entered. If it Is not known, a zero is entered, and both
stack gas velocity in meters per second, and stack diameter at the
top of the stack in meters are required.
The above inputs are then used to perform the calculations and
are repeated on the output. The output table is in the form of 5
columns of information: stability, wind speed, maximum concentration,
the distance of the maximum concentration, and plume height. The
maximum concentration is in grams per cubic meter. The distance of
the maximum concentration is in kilometers and the plume height which
is final plume height according to Briggs method is in meters. When
the table is completed for each stability or the stated stability the
program asks for input for a second run, starting with the title. This
is repeated until the user wishes to terminate the run. Termination
is by entering END when the title is asked for. For second and
subsequent runs, parameter values which are to remain unchanged from
the preceding run may be entered by using a comma. Concentrations of
9.90 E+01 and distances of 999.000 are indicators that no concentration
estimates were attempted. Also numbers in parentheses refer to
footnotes. These are printed at the bottom of the first output If
they are needed. The user is offered the option of another run
starting with a title or concluding the run by typing END.
-------
READY
ptrnax
ENTER ALPHANUMERIC TITLE OF UP TO 6 1 e CHARACTERS.. OR “END”.
test of ptmax 7/5/73
ENTER AMBIENT AI TEMPERATURE (DEG K) OR ZERO TO USE DEFAULT VALUE
OF 293.
0
ENTER SELECTED STABILITY CLASS OR ZERO (0) FOR ALL STABILITIES
0
ENTER SOURCE STRENGTH (( /SEC)
287
ENTER PHYSICAL STACK HEIGHT CM)
7
30
ENTER STACK GAS TEMPERATURE (DEG K)
I
350
ENTER VOLUME FLOW CM**3/SEC) IF KNOWN, OR ZERO (0) IF NOT KNOtIN
0
ENTER STACK GAS VELOCITY (M/SEC)
I
20
ENTER STACK DIAMETER CM)
0.6
TEST OF PTt4AX 7/5/73
ANALYSIS OF CONCENTRATION AS A FUNCTION OF STABILITY AND WIN SPEED.
1971 VERSION, 0. 9. TURNER.
EMISSION RATE (G/SEC) = 287.00
PHYSICAL STACK HEI ’HT CM) = 30.00
STACK GAS TE 1P (DEG K) = 350.00
AMBIENT AI TEMPEPATURE (OF.G K) = 293.
STACK GAS VEL (M/SEC) = 20.00
STACK OIAMETER CM) = 0.60
VOLUME FLOW (CU M/SEC) = 5.65
-------
STAB I LIlY
WIND SPEED
(M/SEC)
MAX CONC
(G/CU M)
01ST OF MAX
(KM)
PLUME HEIGHT
CM)
1
0.5
7.7219E—03
0.486
124.6
1
0.8
7.9789E—03
0.397
89.1
1
1.0
8.0089E—03
0.357
77.3
1
1.5
7.761 14E—03
0.297
61.5
1
2.0
7.3243E—03
0.264
53.6
1
2.5
6.8569E—03
0.243
148.9
1
3.0
6. 4141 49E—03
0.229
45.8
2
0.5
5,9313E—03
0.858
124.6
2
0.8
6.8110E—03
0.632
80.1
2
1.0
7.0593E—03
0.555
77.3
2
1.5
7.1534E—03
0.445
61.5
2
2.0
6.9109E—03
0.387
53.6
2
2.5
6.5671E—03
0.352
48.0
2
3.0
6.2193E—03
0.329
45.8
2
4.0
5.5448E—03
0.300
41.8
2
5.0
li.9602E—03
0.233
39.5
3
2.0
6.9162E—03
0.593
53.5
3
2.5
6.65L&5E-03
0.536
48.9
3
3.0
6.3366E—03
0.498
45.8
3
14.0
5.6913E-03
0.452
41.8
3
5.0
5.].156E—03
0.424
39.5
3
7.0
Ie.2116E— 03
0.392
36.8
3
10.0
3.3027E—03
0.368
34.7
3
12.0
2.8817E—03
0.359
33.9
3
15.0
2.4164E—03
0.350
33.2
-------
‘4 0.5 2.7076E—03 Lj.218 121i.5
4 0.8 3.8535E—03 2.469 89.1
‘ 4 1.0 1 ,.3 419E—03 1.978 71.3
4 1.5 5.J153E—03 1.387 61.5
‘ 4 2.0 5.2362E—03 1.120 53.6
‘4 2.5 5.20 ’ 49E—03 0.980 48.9
‘4 3.0 5.00L 2E—03 0.912 45.8
‘4 4.0 ‘ 4 .5535E—03 0.827 41.8
4 5.0 4.1275E—03 0.777 39•5
‘4 7.0 3 . 4331E— 03 0.719 36.8
4 10.0 2.7143E—03 0.676 314 7
4 12.0 2.3762E—03 0.660 33•()
‘4 15.0 1.qgq3E—o3 0.643 33.2
4 20.0 1.5790E—03 0.626 32.4
5 2.0 2.8291E—03 2.547 61.0
5 2.5 2.4885E—03 2.403 58.7
5 3.0 2.2350E—03 2.294 57.0
5 4.0 1.8775E—03 2.139 5’4.6
5 5.0 1.6334E—03 2.031 52.8
6 2.0 2.4535E—03 4.410 55.7
6 2.5 2.1665E—03 4.120 53.9
6 3.0 1.9514E—03 3.906 52.4
6 4.0 1.6456E—03 3.603 50.
6 5.0 1.4270E—03 3.434 48.9
ENTER ALPHA;JUMERIC TITLE OF U TO 64 CHP RACTE’ S, q lt r) t
end
-------
Users’ Guide to PTDIS (the Interactive Version of DBT 43)
Program Abstract
PTDIS calculates downwind ground—level concentrations
for various downwind distances for the input meteorological
conditions. Only an individual source can be considered.
Inputs to the program consist of information on the
source and Information on the meteorological conditions
to be considered. Primary output of the program consists
of a table with height of emission and concentration
given for each downwind distance. Also included in
this table but not frequently needed are the values of
the dispersion parameters sigma y and sigma z for each
distance and also a relative concentration normalized
for wind speed and source strength commonly called
xu/Q. An optional feature of the program allows the
user to enter a value of concentration to be used for
the determination of half—width of isopleths. For
each distance, if the concentration exceeds the stated
isopleth value the half—width of an isopleth will be
determined. Also this half—width will be compared in
the form of a ratio to the half—width of a sector of
given angular size in terms of degrees. The user is
also given the option of either specifying effective
height of emission or having it calculated using
Briggs’ plume rise methods.
This program determines the concentration at ground level from
a single point source using a steady—state Gaussian model. The
computations used are similar to those shown in the Workbook of
Atmospheric Dispersion Estimates. The dispersion parameter values
are also those given in Figures 3—2 and 3—3 of this Workbook. The
concentrations are for a single meteorological condition defined
by a stability class using the numbers 1 through 6 to represent the
Pasquill stability types A through F. The single wind speed used
is assumed to be representative of the top of the stack, as well
as through the layer that plume rise occurs. The effect of a
-------
2
definite limit to vertical dispersion or mixing height is included
in the computations. It is assumed that complete eddy reflection
occurs at this barrier. It is assumed that the given stability occurs
from ground level to the mixing height. The concept of a mixing
height is not employed for stabilities 5 or 6. It is assumed that
there are no topographic obstructions in the vicinity of the
source and that the source is in an area of either flat or gently
rolling terrain. No consideration of the possibility of aerodynamic
downwash is included.
Given a stability and wind speed condition, this program is
useful in obtaining the variation of ground level concentration
with distance. Concentratons derived from these calculations can
be considered to be valid for averaging times from ten minutes to
an hour. The meteorological input must also be valid for these
averaging times. Although the computational system will yield
numerical values for effective plume heights of above 1,000 meters,
any computations performed for plume heights of greater than three
or four hundred meters should be considered with a certain amount
of skepticism.
The use of the interactive version of this program is quite
simple. Some cautionary statements are printed out first
indicating where the computations are most valid. The number of
distances for which calculations are to be made is entered first.
-------
3
There is a maximum of 50 distances. After this number of distances
has been entered, the distances, in kilometers, are entered separated
by commas or spaces. Although the distances can be entered in any
order, it is handy in interpreting output if the smallest distances
are first and the distances ordered with increasing distance. The
user Is asked if he wants to use the isopleth option. He can answer
yes or no. If he says yes he is asked how many isopleths he wants to
consider. There is a maximum of eight. Next, the isopleth values
in grams per cubic meter are entered separated by coas and he is
also asked to enter a wind segment size in degrees. This is
commonly 100 or 22.5. This is used to compare the isopleth
widths at each distance. Whether the Isopleth option is used or
not, the next variable entered is the stability class, a number
from 1 to 6. Next the wind speed in meters per second is entered,
and finally the mixing height in meters. This concludes the entries
of the meteorological variables. The source strength In grams per
second of the pollutant being considered is entered next. If a
specific effective height of emission Is desired this is entered,
in meters. A zero is entered here if the plume rise is to be
calculated. The physical stack height in meters is entered. The
stack gas temperature in degrees Kelvin is entered. The ambient
air temperature in degrees Kelvin is entered. If a zero is used
here a value of 293 is used. If the stack gas volume flow, in
cubic meters per second, is known it is entered next; if it is
-------
4
unknown, a zero is entered. If the volume flow was not entered,
the stack gas velocity, in meters per second, is next entered and
the stack diameter in meters. Outputs consist of a listing of
the source conditions, the meteorological conditions, and a table
giving distance, effective height, and concentration for each
distance. If the isopleth option is used an additional table
giving the half—width, and a ratio of the isopleth half—width to
the sector half—width for each isopleth value, for each distance.
This concludes one run of the program. The user has options as
to where he can reenter the program or to conclude the run.
-------
• READY
ptdis
** NOTICE ***
USE OF THIS MODEL PRIOR TO JUNE 14, 1973 MAY HAVE PPOTJCED Er n EO1JS
RESULTS IF MORE THAN ONE PASS WAS MADE AND A COMMA WAS U.SEP TO O EFIrJE
MIXING HEIGHT.
DO YOU WANT THE PRECAUTIONARY MESSAGE PRINTED? EUTE YES OR NO
yes
CARE SHOUL’) BE EXERCISED IN THE INTERPRETATION OF THESE C4LCULATEP
CONCENTRATIONS. CONCENTRATION ESTI!IATES PlAY BE EXPECTE ) TO BE
WITHIN A FACTOR OF THREE FOR: 1) ALL STABILITIES OR DISTA 1CES
OF TRAVEL OUT TO A FEW HUNDRED METERS. 2) NEUTRAL TO MODERATELY
UNSTABLE CONDITIONS FOR DISTANCES OUT TO A FEW KILOIETERS.
3) UNSTABLE CONDITIONS IN THE LOWER 1000 METERS OF THE ATIIOSP’IERE
WITH A MARKED INVERSION ABOVE FOR DISTANCES OUT TO TEN
KILOMETERS OR MORE. FOR OTHER CONDITIONS THESE ESTIM4TES BECOME LESS
RELIABLE FOR EXTREMES OF STABILITY AND AS TRAVEL DISTANCE
INCREASES.
ENTER ALPHANUMERIC TITLE (UP TO 64 CHARACTERS)
test of ptdis 7/5/73
ENTER NUMBER OF DISTANCES FOR WHICH CALCULATIONS AR’E TO BE
MADE. MAXIMUM 50
22
ENTER DISTANCES (KM) SEPARATED BY COMMAS OR SPACES
0.1,0.2,0.3,0.4,0.452,0.5,0.6,0.7,0.8,0.9,1.0,2.,3.,5.,7.,10.,
15., 20., 30. ,50. , 70., 100.
DO YOU WANT THE ISOPLETH OPTION? ENTER YES OR NO
yes
ENTER NUMBER OF ISOPLETHS TO BE CONSIDERED. MAX!MU ’1 S
3
ENTER ISOPLETH VALUES (G/’ 1**3) SEPARATED BY COMMAS OR SPACES
1. Oe—03,2 .Oe—4, 1. Oe—04
ENTER WIND SEGMENT SIZE (DEG)
22.
-------
ENTER SOURCE STRENGTH (G/SEC)
287
ENTER EFFECTIVE HEIGHT OF EMISSION CM) IF YOU WISH OR ENTER ZERO
(0) TO HAVE PLUME RISE CALCULATED
0
ENTER PHYSICAL STACK HEIGHT (M)
30
ENTER STACK GAS TEMPERATURE (DEG K)
‘7
350
ENTER VOLUME FLOW (M**3/SF C) IF KNOWN, OR ZERO (0) IF NOT KNOWN
‘7
0
ENTER STACK GAS VELOCITY (M/SEC)
‘7
20
ENTER STACK DIAMETER CM)
‘7
0.6
ENTER AMBIENT AIR TEMPERATURE (DEG K), OR ZERO (0) TO USE DEFAULT
VALUE OF 293
‘7
0
ENTER STABILITY CLASS (1-6)
‘7
3
ENTER WIND SPEED (M/SEC)
‘7
14
ENTER MIXING HEIGHT CM)
‘7
700
-------
DOWNWIND CONCENTRATIONS FOP SPECIFIC DISTANCES
DBTI43 - JUNE 1973 VERSION, 0. B. TUR 4ER
TEST OF PTDIS 7/5/73
* * * SOURCE CONDITIONS * * *
SOURCE STRENGTH (G/SEC) = 287.0
PHYSICAL STACK HEICHT (Ii) = 30.0
STACK GAS TEMPERATURE (DEG K) = 350.0
STACK GAS VELOCITY (M/SEC) = 20.0
STACK DIAMETER CM) = 0.6
VOLUME FLOW (M**3/SEC) = 5.7
* * * METEOROLOGICAL CONDITIONS * * *
AMBIENT AIR TEMPERATURE (DEG K) = 293.0
STABILITY CLASS = 3
WIND SPEED (M/SEC) = 11.0
HEIGHT OF MIXING LAYER (M) = 700.0
FINAL EFFECTIVE HEIGHT OF EMISSION CM) =
DISTANCE TO FINAL EFFECTIVE HEIGHT (KM) =
lel. 8
0.095
DISTANCE
HEIGHT
CONCENTRATION
SIGY
SIGZ
C II*U/Q
(KM)
CM)
(G/CU M)
CM)
CM)
(SEC/ .1**3)
0.100
41.8
3.t 1 1E—08
12.46
7.44
4.75E—10
0.200
41.8
8.1OE—04
23.62
14.03
1.13E—05
0.300
41.8
3.95E—03
34.29
20.33
5.50E—05
0.400
41.8
5.54E—03
44.65
26. 45
7.72E—05
0.452
41.8
5.69E—03
49.94
29.57
7.93E—05
0.500
41.8
5.60E—03
54.77
32.43
7.80E—05
0.600
41.8
5.08E—03
64.71
38.32
7.08E—05
0.700
41.8
4 . 43E—03
74.49
44.12
6.18E—05
0.800
41.8
3.83E—03
84.14
49.85
5.34E—05
0.900
41.8
3.31E—03
93.68
55.52
Ii.61—05
1.000
41.8
2.87E—03
103.11
61.14
4. OOE—05
2.000
41.8
q.59E—04
193.45
115.26
1.34E—05
3.000
41.8
4.75E—04
279.00
167.01
6.62E—06
5.000
41.8
1 g2 —oL
441.64
266.’i7
2.67E—0
7.000
41.8
1. O SE—04
5g6.81
362.4)
1.46E—06
10.000
41.8
5.76E—05
820.13
502.32
8.03E—07
15.000
41.8
3.51E—05
1175.00
727.85
!i.90E—07
20.000
41.8
2.70E—05
1514.57
9t 1 5 ’J3
3.76E—07
30.000
41.3
1.80E—05
2161.94
1372.00
2. IiE —07
5i).000
41.8
1.21E—05
3373.06
2189.25
1.69E—07
70.000
41.8
9.07E—06
4510.52
2978.13
1.26E—07
100.000
41.8
6.68E—06
6123.50
4126.98
9 .31E—08
-------
RATIO IS THE HALF-WIDTH OF THE ISOPLETH COMPARED TO THE HALF-WIDTH OF
A SECTOR OF 22.5 DEGREES AT THIS DISTANCE.
0. 10000E—02
ISOPLETH VALUES C(RAMS PER CU3IC METER)
0. 20000E—03
HA IF —
WI t )TH
CM)
RATIO
0.
10000E—03
HALF-
HALF-
HALF-
DISTANCE WIDTH
RATIO WIDTH
RATIO WIDTH
(KM) CM)
CM)
CM)
RATIO
0.100 0.
0.0 0.
0.0 0.
0.200 0.
0.0 40.
1.006 149•
0.300 57.
0.961 84.
1.412 93.
1.227
0.400 83.
1.045 116.
1.453 127.
1.567
1.597
0. 1452 9 1 4•
1.042 130.
1.443 142.
0.500 102.
1.028 142.
1.427 156.
1.585
0.600 117.
0.982 165.
1.384 182.
0.700 129.
0.927 186.
1.336 206.
0.800 138.
0.870 205.
1.288 228.
1.477
0.900 145.
0.813 222.
1.243 248.
1.431
1.000 150.
0.755 238.
1.199 268.
2.000 0.
0.0 3143.
0.863 412.
1.346
3.000 0.
0.0 367.
0.616 493.
1.036
5.000 0.
0.0 0.
0.0 501i.
7.000 0.
0.0 0.
0.0 187.
10.000 0.
0.0 0.
0.0 0.
0.0
15.000 0.
0.0 0.
0.0 0.
0.0
20.000 0.
0.0 0.
0.0 0.
30.000 0.
0.0 0.
0.0 0.
0.0
0.0
50.000 0.
0.0 0.
0.0 0.
0.0
70.000 0.
0.0 0.
0.0 0.
100.000 0.
0.0 0.
0.0 0.
0.0
0.0
ENTER “DI STANCE”
OR “SOURCE” OR
“METEOROLOGY”
OR
“END”
?
end
r r nv
-------
Users’ Guide to PThTP (The Interactive Version of DBT 51)
Program Abstract
PThTP produces hourly concentrations at up to 30
receptors whose locations are specified from up to 25
point sources. A Gaussian plume model is used. Inputs
to the program consist of the number of sources to be
considered, and for each source the emission rate,
physical height, stack gas temperature, volume flow,
or stack gas velocity and diameter, the location, in
coordinates. The number of receptors, the coordinates
of each and the height above ground of each receptor
are also required. Concentrations for a number of
hours up to 24 can be estimated, and an average
concentration over this time period is calculated.
For each hour the meteorological information required
is: wind direction, wind speed, stability class,
mixing height, and ambient air temperature.
The assumptions that are made in this model follow:
Meteorological conditions are steady—state for each hour and a
Gaussian plume model is applicable to determine ground level
concentrations. Computations can be performed according to the
“Workbook of Atmospheric Dispersion Estimates.” The dispersion
parameter values used for the horizontal dispersion coefficient,
sigma y, and the vertical dispersion coefficient, sigma z, are those
given in Figures 3—2 and 3—3 of the Workbook. The sources and
receptors exist in either flat or gently rolling terrain, and the
stacks are tall enough to be free from building turbulence so that
no aerodynamic downwash occurs. The wind speed and wind direction
apply from the shortest to the tallest plume height. No wind
direction shear or wind speed shear occurs. The given stability
exists from ground—level to well above the top of the plume.
-------
2
Calculations for each hour are made by considering each source—
receptor pair. Plume rise is calculated according to Briggs’ plume
rise estimates. For each source—receptor pair, the downwind and
crosswind distances are determined. If the downwind distance is
closer than the distance to final rise, the plume rise for this
distance is calculated. The concentration from this source upon
this receptor is determined using these distances by the Gaussian
model.
The use of the Interactive version of the program is relatively
straightforward. First, an alphanumeric title to identify the output
is entered. Next, the number of sources to be considered is given.
The source strength, physical height, stack gas temperature, and
volume flow is entered for each stack. If the volume flow is not
known the stack gas velocity and diameter are required. The coordinates
based on a coordinate system having units of one kilometer are
required for each source. Next, the number of receptors to be
processed, the coordinates of each and the height above ground for
each are entered. The meteorological information includes the number
of hours to be averaged up to 24, the wind direction, wind speed,
stability class, mixing height, and ambient air temperature are
entered for each hour. An option exists to print the partial
concentrations, that is, the concentration from each source at each
receptor. Also, an option exists to print the hourly concentrations.
-------
3
The output is quite simple, consisting of title followed by
input information on the sources, receptors, and meteorology. This
is followed by hour by hour partial concentrations if desired and
total concentrations. If partial concentrations are printed the
final plume height for that hour for each source is also printed.
Then average concentrations for the time period are printed including
partial concentrations if desired. When the output is complete,
the user is offered the option of ending the run or entering at
3 different points. He may go back to enter new sources or he
may keep the same sources and enter new receptors or he may keep
both the same sources and receptors and enter only different
meteorological conditions.
-------
ptmtp
* * * NOT ICE * * *
USE OF THIS MODEL PRIOR TO JUNE 7, 1973 lAY !IAVE PRO )UCE!
ERRONEOUS RESULTS IF ANY RECEPTORS WERE CLOSER Tfl SOURCES
THAN THE OISTANCE TO FINAL PLUME RISE.
ENTER ALPHANUMERIC TITLE CUP TO 611 CHARACTERS)
test of ptmtp 7/5/73
ENTER NUMBER OF SOURCES TO BE CONSIDERED. MAX 25
1
ENTER SOURCE STRENGTH (G/SEC) FOR EACH STACK
287,287,287,287
ENTER PHYSICAL HEIGHT (Ii) OF EACH STACK
30, 30, 30, 30
ENTER GAS TEMPERATURE (DEG K) OF EACH STACK
Le*35 0
IS VOLU 1E FLOW KNOWN FOR EACH STACK? YES OR NO
no
ENTER GAS VELOCITY (t 1/SEC) FOR EACH STACK
*20
ENTER DIAMETER (M) OF EACH STACK
4*0.5
ENTER COORDINATES (KM) OF EACH STACK. OROERE PAIRS
1.10., 1.05,0., 1.10,0., 1.15,0.
-------
ENTER NUMBER OF RECEPTORS TO BE PROCESSED. MAX 30
7
‘I ’
ENTER COORDINATES CKM) OF EACH RECEPTOR. ORDERED PAIRS
7
0.8,0., 1.02,0., 1.07,0., 1.12,0., 1.17,0., 1.2,0.,
9
1.3,0., 1.14,0., 1.5,0., 1.6,0., 1.7,0., 1.8,0.,
9
1.9,0., 2.0,0.,
ENTER HEIGHT CM) ABOVE GROUND FOR EACH RECEPTOR
7
114*0.
ENTER NUMBER OF HOURS TO BE AVERAGED. t.IAX 214
9
3
ENTER WIND DIRECTION (DEG) FOR EACH HOUR
7
265,270,275
ENTER WIND SPEED (M/SEC) FOR EACH HOUR
7
‘4,14,14
ENTER STABILITY CLASS FOR EACH HOUR
7
3*3
ENTER MIXING HEIGHT CM) FOP EACH HOUR
7
3*7 00
ENTER AMBIENT AIR TEMPERATURE (DEG K) FOR EACH HOUR
7
3*2 93
DO YOU WANT PARTIAL CONCENTRATIONS PRINTED? YES OR NO
7
yes
DO YOU IANT HOURLY CONCENTRATIONS PRINTED? YES OR NO
7
yes
-------
CEPTOR NUMBER***
1 2 3
0.0 0.0 1i.650E—13 1.554E—05 1.733E—0 [ i 6.0114E—0 14
0.0 0.0 0.0 Ii.650E—13 1.5514E—06 4.653E—05
0.0 0.0 0.0 0.0 4.650E—13 2.LiO5E—08
0.0 0.0 0.0 0.0 0.0 8.15)E—21
TOTAL CONCENT ATIOFJ (G/I1**3)
0.0 14.650E—13 1.55LiE—06 1.7149E—0 1 4 6.Z480E—04
HOUR # 1
***RECEPTOR NUMBER***
7 8 9 10
S EFFHT PARTIAL CONCENTRATIONS (G/M**3)
2.925E—03 ti.072E—03 4.080E—03 3.669E—03 3.179E—03 2.726E—03
1.772E—03 3.695E—03 4.164E—03 3.898E—03 3.42 4E—03 2.9Li5E—03
6.014E—04 2.925E—03 k.072E—03 Le.080E—03 3.669E—03 3.179E—03
Li.653E—05 1.772E—03 3.695E—03 Ii.164E—03 3.898E—03 3.14214E—03
TOTAL CONCENTRATION (G/M**3)
5.31e 1 4E—03 1.2le6E—02 1.60].E—02 1.581E—02 1.Lel7E—02 1.227E—02
2.337E—03 2.0114E—03
2.523E—03 2.168E—03
2.726E—03 2.337E—03
2.945E—03 2.523E—03
TOTAL CONCENTRATION (G/M**3)
hOUR # 1
S EFFHT PARTIAL CONCENTRATIONS (G/M**3)
1 42.
2 142.
3 42.
le 42.
4
5
6
0.0
11
12
1 42.
2 42.
3 42.
14 142.
HOUR #
S EFFHT
1 42.
2 42.
3 42.
14 42.
CE PT() R
13
l ie
N U I’•l B E R * * *
PARTIAL CONCENTRATIONS ((/M**3)
l.053E—02 9.OLe3E—03
-------
TEST OF PTMTP 7/5/73
MULTIPLE SOURCE MODEL DBT51, JUNE 1973 VERSION
*
I
(DEG K)
CES *
HP
(M)
30.0
30.0
30.0
30.0
* * *SOUR
NO Q
(GISEC)
1 287.0
2 287.0
3 287.0
4 287.0
***REC
N’) RREC
(KM)
* *
TS
VS
0
VF
(DEG
K)
(M/SEC)
(N)
(M**3/sEc)
350.0
350.0
350.0
350.0
20.0
20.0
20.0
20.0
0.6
0.6
0.6
0.6
E PTO RS *
SREC
(KM)
R
(KM)
1.000
1.050
1.100
1.150
5.7
5.7
5.7
5.7
S
(KM)
0.0
0.0
0.0
0.0
1
0.800
0.0
2
1.020
0.0
3
1.070
0.0
4
1.120
0.0
5
1.170
0.0
6
1.200
0.0
7
1.300
0.0
8
1.400
0.0
9
1.500
0.0
10
1.600
0.0
11
1.700
0.0
12
1.800
0.0
13
1.900
0.0
14
2.000
0.0
* *
z
CM)
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
HL
(H)
***METE
NO THETA
(DEG)
OR’) L
U
(M/SEC)
OG
KST
1
265.0
4.0
3
700.
293.
2
270.0
4.0
3
700.
293.
3
275.0
4.0
3
700.
293.
-------
HOUR # 2
***RECEPTOR NUMBER***
1 2 3
S EFFHT PARTIAL CONCENTRATIONS (G/M**3)
0.0 0.0 6.813E—13
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
TOTAL CONCENTRATION (G/M**3)
0.0 0.0
HOUR # 2
***RECEPTOR NUMBER***
7 8
S EFFHT PARTIAL COflCENTRATIONS (G/M**3)
11 12
1 42.
2 ti2.
3 42.
4 42.
3.946E—03 5.538E—03 5.598E—03 5.077E—03
2.384E—03 5.004E—03 5.68gE—03 5.371E—03
8.096E—Ole 3.946E—03 5.538E—03 5.598E—03
.6.323E—05 2.384E—03 5.OOliE—03 5.689E—03
TOTAL CONCENTRATION (G/M**3)
4.431&E—03 3.829E—03
4.756E—03 4.123E—03
5.077E—03 Le.434E—03
5.371E—03 Ii.756E—03
7.203E—03 1.687E—02 2.183E—02 2.173E—02 1.964E—02 1.714E—02
HOUR # 2
S EFFHT
1 42.
2 42.
3 42.
‘4 42.
CEPTOR NUM3ER***
13 14
PAPTIAL CONCENTRATIONS (G/M**3)
3.306E—03 2.867E—03
3.557E—03 3.077E—03
3.829E—03 3.306E—03
Li.123E—03 3.557E—03
TOTAL CONCENTRATION (G/M**3)
1 42.
2 42.
3 42.
4 42.
‘4
2. 148E—06
6. 813E—13
0.0
0.0
5 6
2.342E—0 ’ 4 8.Og6E—o ’ 4
2.lZe8E—06 6.322E—05
6.813E—].3 3. ’eO7E—08
0.0 1.32LiE—20
6.813E—13 2.1 ’ 48E—06 2.364E-04 8.728E—Ole
9
10
1.432E—02 1.281E—02
-------
0.0 0.0 14.637E—13
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
TOTAL CONCENTRATION (G/M**3)
0.0 0.0
CEPTOR NUMBER***
7 8 9 10
PARTIAL CONCENTRATIONS (G/M**3)
1 42.
2 42.
3 42.
4 42.
2.918E—03 k.062E—03 4.070E—03 3.660E—03
1.768E—03 3.686E—03 Li.154E—03 3.888E—03
6.OO1E—04 2.918E—03 Li.062E—03 k.070E—03
4.642E—05 1.768E—03 3.686E—03 le.l5teE—03
TOTAL CONCENTRATION (G/M**3)
3.171E—03 2.719E—03
3.415E—03 2.937E—03
3.660E—03 3.171E—03
3.888E—03 3.415E—03
5.332E—03 1.2fi3E—02 1.597E—02 1.577E-02 1.’413E—02 1.224E—02
HOUR # 3
S EFFHT
CEPTOR NUt.IBER***
13 14
PARTIAL CONCENTRATIONS (GIM**3)
2.331E—03 2.009E—03
2.517E—03 2.162E—03
2.719E—03 2.331E—03
2.937E—03 2.517E—03
TOTAL CONCENTRATION (( /M**3)
HOUR # 3
**RECEPTOR
1 2 3
S EFFHT PARTIAL CONCENTRATIONS (G/M**3)
N U M B E R * * *
II
1 42.
2 42.
3 42.
l 42.
HOUR 3
S EFFHT
1. 550E—06
£i.637E—13
0.0
0.0
5 6
1.729E—04 6.OO1E—04
1.550E—06 4.642E—05
4.637E—13 2.398 —08
0.0 8.121E—21
4.637E—13 1.550E—06 1.745E—OL1 6.465E—04
11 12
1 42.
2 42.
3 42.
4 42.
1.050E—02 9.019E—03
-------
AVERAGE CONCENTRATIONS FOR 3 HOURS.
***RECEPTOR NUMflER***
1 2 3 14
S PARTIAL CONCENTRATIONS (G/M**3)
0.0 0.0 5.367E—13
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
TOTAL CONCENTRATION CG/M**3)
5.367E—13 1.751E—06 1.952E—Oli
NUMBER * * *
S
1 3.263E—03 !4.557E—03 k.582E—03 1 4.135E—03 3.595E—03
2 ]..975E—03 Le.128E—03 Li.669E—03 i..386E—03 3.865E—03
3 6.7014E—014 3.263E—03 1 .557E—03 Ii.582E—03 L .135E—O3
5.206E—05 1.975E—03 1s.128E—03 14.669E—03 L*.386E—03
TOTAL CONCENTRATION (G/M**3)
5.960E—03 1.392E—02 1.794E—02 1.777E—02 1.598E—02
CEPTOR NUMBER***
13 114
PARTIAL CONCENTRATIONS (G/M**3)
2.658E—03 2.297E—03
2.866E—03 2.1469E—03
3.091E—03 2.658E—03
3.335E—03 2.866E—03
TOTAL CONCENTRATION (G/M**3)
1.195E—02 1.O2gE—02
1
2
3
14
5
1.751E—06 1.935E—014
5.367E-13 1.751E—06
0.0 5.367E—13
0.0 0.0
0.0
CE PTO R
7
0.0
8
9
PARTIAL CONCENTRATIONS (G/M**3)
10
1].
6
6. 7014E—014
5.206E—05
2. 737E—08
9. 836E—21
7. 2214E—014
12
3. 091E—03
3. 33 5E—03
. “‘ _
_, . a _J S L_
3. 865 E—03
1. 389E—02
S
1
2
3
ENTER “SOURCES” OR “RECEPTORS” OR “METEOROLOGY” OR “END”
end
READY
‘ogoff
-------
REFERENCES
Briggs, Gary A., 1971: Some recent analyses of plume rise observation.
1029—1032 in Proceedings of the Second International Clean Air
Congress , edited by H. N. Englund and W. T. Berry. Academic
Press, New York.
Turner, D. Bruce, 1970: Workbook of Atmospheric Dispersion Estimates .
Office of Air Programs Publication No. P2-26. Environmental
Protection Agency, Research Triangle Park, N. C., 84 p.
------- |