-------
IPERD). After all hourly meteorological data have been processed by the
program, the "N"-day tables, highest, second-highest and third-highest tables
and the maximum 50 tables are alternately printed for each source group for
each specified time interval. The number of tables is governed by the number
of source groups (NGROUP) and time periods (ISW(7) through ISW(14)) specified.
b. Output File. The ISCST program is capable of generating an output
file containing the calculated average concentration or total deposition
values based on the selected time periods and source groups. If ISW(5) equals
"1", this output file is generated. The user must assign an output file and
associate the logical unit number specified in parameter ITAP to the output
file (see Section 3.2.3.a).
The output file is written with a FORTRAN unformatted (binary) WRITE
statement and consists of constant length records whose lengths equal the
total number of receptor points (NPNTS) plus 3 words. Word 1 of each record
contains the hour at which the corresponding values were calculated in words 4
to NPNTS +3. Word 2 contains the Julian Day and word 3 contains the source
group number. Words 4 through NPNTS + 3 contain the calculated average
concentration or total deposition values for all receptors. The values
calculated for the receptor grid (if any) are written first followed by the
values calculated at the discrete receptors (if any). Starting with the first
Y point (direction radial) of the Y-axis (radial) grid, the calculated values
are written for the X-axis (ranges) in the same order that receptor locations
were entered in parameter GRIDX (see Section 3.2.3.a). For each successive
Y-axis (radial), the values are written for the X-axis (ranges). After the
calculated values have been written for the receptor grid, the calculated
values are written for the discrete points in the order the discrete points
were entered in parameters XDIS and YDIS (see Section 3.2:3.a).
3-60
-------
The content and number of records produced is governed by the number of
source groups (specified in parameter NGROUP) and time periods (specified in
parameters ISW(7) through ISW(14)}. For each day of meteorological data
processed by the program and for each hour, the program generates records of
calculated values for all applicable time period intervals for all source
groups. For hour one, a 1-hour record of calculated values for source group
1, followed by 1-hour records of calculated values for each remaining source
groups are written to the output file. For hour two, a 1-hour and a 2-hour
record are written to the output file for each source group. For hour three,
a 1-hour and 3-hour record are written to the output file for each source
group. For hour four, a 1-hour, 2-hour, and 4-hour record of calculated
values are written to the output file for each source group. This format is
continued for each hour of the day. For example, if there is one source group
and only 24-hour average concentrations are calculated, only one record per
day is written to the output file. If ISW(15) equals "1", records of the
"N"-day average concentration or total deposition values are additionally
written to the output file for all source groups after the program has
processed all "N"-days of meteorological data.
3.2.5 Program Run Time, Page and Tape Output Estimates
This section provides the user with equations which estimate the amount of
run time required and program output generated for a given problem run. The
.equations describing the amount of printed output data (in pages) and tape
output data (in words) can be quite accurately estimated. The run time
estimate is less accurate because of unknowns such as the nature of the hourly
meteorology and wake effects. These unknowns may affect the run time estimate
significantly for a large problem run.
3-61
-------
a. Run Time. The amount of time a problem takes to execute is. primarily
governed by si:-: factors. These factors are: (1) the number of hours in a day
of meteorological data (NHOURS); (2) the number of days of meteorological data
processed (NDAYS); (3) the number of sources (NSOURC); (4) the number of
source groups (NGROUP); (5) the number of receptor points (NPNTS); and (6) the
number of time periods (NAVG). Using these factors, the following equation
estimates the run time in minutes:
No. of Minutes * C • (NDAYS +!)•(!+ NHOURS • (1 + 0.8 NSOURC
(3-2)
• (1 + 0.6 • NPNTS + 0.1 • NGROUP • NAVG)))
where
C = 2.1 • 10~5
The constant, C, is derived from problem runs made on a UN I VAC 1103 computer
and is different for other computers.
b. Page Output. The number of pages of printer output produced by a
problem run is primarily controlled by which categories of output are desired
by the user. The content of these categories of program print output are
discussed in Section 3.2.4.a. Input parameters ISW(6), ISW(15), ISW(16),
ISW(17), and ISW(18), discussed in Section 3.2.3.a., control which categories
of program print output are produced. Other factors which determine the
amount of print output are the number of receptor points, number of source
groups, and the number of time periods for which average concentration or
total deposition values are computed.
If ISW(6) equals "1", all input data are printed, producing about 5 pages
of print output. For source with gravitational settling categories (NVS
greater than zero) or variational emission rates (QFLG greater than zero), add
one-third of a page per source. For sources with direction specific building
dimensions (HB less than zero), add one-fifth of a page per source. If ISW(6)
3-62 12/87
-------
equals "2", all meteorological data processed by the program are printed. Add
one page for every day of meteorological data processed.
If ISW(15) equals "1", tables of the "N"-day average concentration or
total deposition values are printed. The number of tables printed equals the
number of source groups desired by the user (NGROUP). If parameter NGROUP is
specified as "0", one table will be printed. The number of pages produced for
each "N"-day table is given the following equation:
Number of Pages = (NXPNTS/9) (NYPNTS/38) + (NXWYPT/114) ' (3-3)
where
NXPNTS = the number of X points on the' X-axis grid or
the number of grid ranges
NYPNTS = the number of Y points on the Y-axis grid or
the number of grid direction radials
NXWYPT = the number of discrete receptor points
Round up any fractional number in each term to the nearest whole number.
If ISW(16) equals "1", tables of average concentration or total deposition
for user-defined combinations of source groups and time periods for each day
of meteorological data processed by the program are printed. The number of
tables produced by this output category for each day is given by the following
equation:
No. of Tables = NGROUP • (24/IPERD) • ISW(7)
+ (12/IPERD) • ISW(8) + (8/IPERD) • ISW{9)
+ (6/IPERD) • ISW(IO) + (4/IPERD) • ISW(ll) (3-4)
+ (3/IPERD) • ISW{12) + (2/IPERD) • ISW(13)
+ (1/IPERD) • ISW(14)
where
NGROUP = number of source groups as specified by input parameters
NGROUP. If NGROUP is specified as "0", assume a value of
"1" for this equation.
3-63 12/87
-------
IPERD = "N"th time interval for all time periods as specified by
input parameter IPERD. Note that if IPERD is not set to
"0", the term (3/IPERD) • ISW(i) equals (3) • ISW(i). If
IPERD is set greater than "0", 'the term (j /IPERD) •
ISW(i) equals (1) • ISW(i) if (j/IPERD) is greater than
or equal to "1"; otherwise, it equals (0) • ISW(i) if
(j/IPERD) is less than "1".
ISW(7)- = the corresponding 1-, 2-, 3-, 4-, 6-, 8-, 12-, and
ISW(14) 24-hour time periods as specified by input parameters
ISW(7) through ISW(14). The "1" or "0" values specified
by the user in these parameters are the numeric values
used in the equation
The number of pages produced by each table is given in Equation (3-3). Hence,
the total number of pages generated by the print output option ISW(16) equals
the product of the number of days processed by the program for a probler, run,
the number of tables printed according to Equation (3-4) and the number of
pages produced per table according to Equation (3-3).
If ISW(17) equals "1", tables of the highest and second-highest average
concentration or total deposition values found at each receptor are printed
for all user-defined combinations of source groups and time periods. If
ISW(17) equals "2" tables of highest, second-highest, and third-highest are
printed. The number of tables printed equals two or three (depending on
ISW{17)) times the number of time periods specified (the number of input
parameters ISW(7) through ISW(14) set to "1") multiplied by the number of
source groups desired. If no source groups are specified (input parameter
NGROUP equals "0"), assume one source group for the purpose of computing the
number of tables printed by this option (ISW(17». The number of pages each
table produces is given by the following equation:
Number of Pages = (NXPNTS/5) (NYPNTS/38) + (NXWYPT/76) (3-5)
where NXPNTS, NYPNTS, and NXWYPT are defined following Equation (3-3). Round
up any fractional number in each term to the nearest whole nvunber. Hence, the
3-64 . 12/87
-------
number of pages printed by this output category equals two or three, times the
product of the number of time periods., the number of source groups, and the
number of pages produced per table according to Equation (3-5).
If ISW(18) equals "1", tables of the maximum 50 average concentration or
total deposition values calculated are printed for all user-defined
combinations of source groups and time periods. Because each table printed
produces only one page of output, the total number of pages printed by this
output category equals the number of time periods specified (the number of
input parameters (ISW(7) through ISW(14) set to "1") multiplied by the number
of source groups specified. Again, if no source groups are specified (input
parameter NGROUP equal to zero), assume one source group.
Thus, the total number of pages of output produced by the program equals
the sum of the number of pages produced by each optional print output category
desired by the user for a problem run.
c. Output to Disc File. Values of average concentration or total
deposition are written by a FORTRAN unformatted WRITE statement to an output
file only if parameter ISW(5) equals "1". Otherwise, the program does not
generate an output file. It is not practical to discuss the physical amount
(length of magnetic tape or number of tracks or sectors of mass storage)
generated since this introduces factors which depend on the computer
installation. Instead, the number of computer words generated by a problem
run is discussed. The user may then equate this number to a physical amount
for the particular storage device being used.
The output file is written in records, where the length of each record
equals the number of receptor points (NPNTS) plus 3 for a total of NPNTS + 3
computer words for a given problem run. For each day of meteorological data
processed, the number of records written to the file is governed by the number
3-65
-------
of source groups and time.periods specified by the user. If we substitute the
term "Tables" used in Equation (3-4) with the word, "Records" and set IPERD
equal to "0", Equation (3-4) gives the number of records written to the file
for each day of meteorological data processed. All variables used to
formulate Equation (3-4) maintain the same definition. Hence, the number of
records equals the value computed from Equation (3-4) multiplied by the number
of days of meteorological data processed by the program for a problem run.
Also, if input parameter ISW(15) equals "1", .additional records containing
"N"-day average concentration or total deposition values are written to the
file depending on the number of source groups specified by the input parameter
NGROUP. If NGROUP equals "0", assume one source group. Hence, the total
number of computer words written to the file equals the number of records
generated, multiplied by (NPNTS + 3) computer words per record for a problem
run.
3.2.6 Program Diagnostic Messages
The ISCST program prints diagnostic messages when certain conditions occur
during a problem run. The diagnostic messages consist of two types. The
first type is a table format that informs the user of the conditions found,
but does not terminate program execution. The second type is an error message
which informs the user of the condition. The run is terminated after the
error message is printed.
The diagnostic message in a table format informs the user when a receptor
is located within one meter or three building heights (or three effective
building widths) of a source. As shown in Figure C-4 in Appendix C, the table
lists all source-receptor combinations for which this condition has occurred.
The table lists the source number, receptor location, and calculated distance
3-66
-------
between the corresponding source and receptor. A negative distance value
implies that the receptor is located within the dimensions of a volume or area
source.
Four types of diagnostic error messages may be printed by the program. If
the allocated data storage is not sufficient for the data required by a
problem run, an error message is printed (Figure 3-2(a)). An error message is
printed if the station numbers or years read from the meteorological data
input tape do not match the 'corresponding station numbers or years specified
by the user in parameters ISS, ISY, IUS, IUY (Figure 3-2(b». If the number
of input sources equals "0", an error message is printed (Figure 3-2(c)).
Finally, if there are no gravitational settling categories to calculate
deposition for any source, an error message is printed as shown in
Figure 3-2(d).
3.2.7 Program Modification for Computers Other than UKIVAC 1100 Series
Computers
The ISCST program, which is written in FORTRAN 77, provides easy transport
and adaption for use on other computers. The program design requires that:
(1) at least four Hollerith characters can be stored in one computer word;
(2) the computer word lengths of integer and real type variables are the same;
and, (3) at least 132 characters per line can be printed on a page with 57
lines per page. The program requires about 70,000 words of executable
storage, 26,500'of which consist of the program itself compiled on a UNIVAC
1100 Computer. The size of the compiled program will vary depending on the
FORTRAN compiler and computer installation. The remaining 43,500 words
consist of data storage used by the program for storing the input data values,
intermediate values, and output results of a given problem run.
3-67 12/87
-------
TABLE 3-8
I SCSI ERROR MESSAGES
**ERROR**CALCULATED STORAGE ALLOCATION LIMIT EQUALS n AND EXCEEDS THE MAXIMUM
STORAGE ALLOCATION LIMIT OF m
RUN TERMINATED. The program has determined that n locations of core storage
are required for the run, but only m locations are available. See equation
(3-1) in Section 3.2.3
(a)
***ERROR***MET DATA -REQUESTED DOES NOT MATCH MET DATA READ.
'REQUESTED/READ' VALUES ARE:
SURFACE STATION NO. = isisis/jsjsjs YEAR OF SURFACE DATA = iys/jys
UPPER AIR STATION NO. = iuiuiu/jujuju YEAR OF UPPER AIR DATA = iuy/juy
RUN TERMINATED. The surface or upper air station identifiers or the years
selected on the card input deck do not match the identifiers or years in the
preprocessed meteorological tape file. Correct the identifiers and years for
the proper values and rerun.
(b)
***ERROR***NUMBER OF SOURCES TO BE READ EQUALS ZERO. RUN TERMINATED. The
parameter NSOURC on the input card deck has been set to zero. The program
requires at least one source to execute properly.
(c)
***ERROR***SOURCE NUMBER HAS NO GRAVITATIONAL SETTLING CATEGORIES WITH WHICH
TO CALCULATE DEPOSITION. RUN TERMINATED. For deposition calculations, the
program requires particulate settling parameters to be entered for each
source. Check the input card deck and rerun.
(d)
***ERROR***DIRECTION SPECIFIC BUILDING HEIGHT OR WIDTH IS GREATER THAN 9999
FOR SOURCE NO: n, RUN TERMINATEJ. When the stack height is less than or
equal to the building height plus one-half the lesser of the building height
or width, the program expects the building height to be set as a negative
value in order to read direction specific building dimensions. The program
attempted to read direction specific building dimensions and thus encountered
values which were out of range. Reset HB to a negative value, insert 36
direction specific building heights and widths and rerun.
(e)
FIGURE 3-2. (a) through (e) show the five types of error messages printed by
the ISCST Program. The run is terminated after an error message
is printed.
3-68 12/87
-------
If it is necessary to adjust the current allotment of 43,500 words of data
storage, only two FORTRAN statements in the ISCST program need to be
modified. The FORTRAN statement with sequence number ISC07540 (in columns
73-80) in the main program allocates the data storage in array QF. Also, the
value assigned to the variable LIMIT at sequence number ISC07620 must agree
with the value used in array QF.
The program assumes FORTRAN logical unit 5 for the card reader and logical
unit 6 for the printer. These logical unit numbers may be modified on
sequence numbers ISC07690 and ISC07700 in the main program.
3-69 12/87
-------
SECTION 4
USER'S INSTRUCTION FOR THE ISC LONG-TERM
(ISCLT) MODEL PROGRAM
4.1 Summary of Program Options, Data Requirements and Output
4.1.1 Summary of ISCLT Program Options
The program options of the ISC Dispersion Model long-term computer program
ISCLT consist of three general categories:
• Meteorological data input options
• ( Dispersion-model options
• Output options
Each category is discussed separately below.
a. Meteorological Data Input Options. Table 4-1 lists the meteorological
data input options for the ISCLT computer program. All meteorological data
may be input by card deck or by a magnetic tape inventory previously generated
by ISCLT (see Section 4.1.1.C below). ISCLT accepts STAR summaries with six
Pasquill stability categories (A through F) or five Pasguill stability
categories (A through E with the E and F categories combined). It does not
accept STAR summaries with separate day and night neutral categories.
Site-specific mixing heights and ambient air temperature are ISCLT input
requirements rather than options. Suggested procedures for developing these
inputs are given in Section 2.2.1.2. The remaining meteorological data input
options listed in Table 4-1 are identical to the ISCST meteorological data
input options discussed in Section 3.1.1.a.
b. Dispersion Model Options. Table 4-2 lists the dispersion model options
for the ISCLT computer program. In general, these options correspond to the
4-1
-------
TABLE 4-1
METEOROLOGICAL DATA INPUT OPTIONS FOR ISCLT
Input of all meteorological data by card deck or by magnetic tape inventory
previously generated by ISCLT
STAR summaries with five or six Pasguill stability categories
Site-specific mixing heights
Site-specific ambient air temperatures
Site-specific wind-profile exponents
Site-specific vertical potential temperature gradients
Rural Mode or Urban Mode 1, 2 or 3
Final or distance dependent plume rise
Wind system measurement height if other than 10 meters
TABLE 4-2
DISPERSION-MODEL OPTIONS FOR ISCLT
Concentration or dry deposition calculations
Inclusion of the effects of gravitational settling and/or deposition in
concentration calculations
Inclusion of terrain effects (concentration calculations only)
Grid or discrete receptrrs (Cartesian or polar system), with capability to
model receptor heights above ground.
Stack, volume and area sources
Pollutant emission rates held constant or varied by season or by wind speed
and stability
Time-dependent exponential decay of pollutants
Inclusion of building wake, stack-tip downwash and buoyancy-induced dispersion
effects
Time periods for which concentration or deposition calculations are to be made
(seasonal and/or annual)
4-2 12/87
-------
ISCST dispersion-model options discussed in Section S.l.l.b. Pollutant
emission rates may be held constant or varied by season or by wind speed and
stability in ISCLT calculations. The program uses seasonal STAR summaries to
calculate seasonal and/or annual concentration or deposition values.
Additionally, monthly STAR summaries may be used to calculate monthly
concentration or deposition values.
c. Output Options. Table 4-3 lists the ISCLT program output options. A
more detailed discussion of the ISCLT output information is given in Section
4.1.3.
The ISCLT program has the capability to generate a master file inventory
containing all meteorological and source inputs and the results of all
concentration or deposition calculations. This file can then be used as input
to future update runs. For example, assume that the user wishes to add a new
source and modify an existing source at a previously modeled industrial source
complex. Concentration or deposition calculations are made for these or
modified sources alone and the results of these calculations in combination
with select sources from the original file inventory are used to generate an
updated inventory. That is, it is not necessary to repeat the concentration
or deposition calculations for the unaffected sources in the industrial source
complex in order to obtain an updated estimate of the concentration or
deposition values for the combined emissions. The optional master file
inventory is discussed in detail in Section 4.2.4.b.
The ISCLT user may elect to print one or more of the following tables:
• The program control parameters, meteorological input data and
receptor data
• The source input data
• The seasonal and/or annual average concentration or total
deposition values calculated at each receptor for each source
or for the combined emissions from select groups or all sources
4-3
-------
TABLE 4-3
ISCLT OUTPUT OPTIONS
Master file inventory of meteorological and source inputs and the results of
the concentration or deposition calculations
Printout of program control parameters, meteorological data and receptor data
Printout of tables of source input data
Printout of seasonal and/or annual average concentrations or total seasonal
and/or annual deposition values calculated at each receptor for each source or
for the combined emissions from a select group or all sources
Printout of the contributions of the individual sources to the 10 highest
concentration or deposition values calculated for the combined emissions from
a select group of all sources or the contributions of the individual sources
to the total concentration or deposition values calculated for the combined
emissions from a select group of all sources at 10 user-specified receptors
4-4
-------
• The contributions of the individual sources to the 10 receptors
with highest concentration (or deposition) values obtained from
the combined emissions of select groups of sources; or the
contributions of each individual source, as well as the
combined sources, to a select -group of user specified receptor
points; or the maximum 10 concentration (or deposition) values
for each source and for the combined sources, determined
independently of each other
4.1.2 Data Input Requirements
This section provides a description of all input data parameters required
by the ISCLT program. The user should note that some input parameters are not
read or are ignored by the program, depending on the values assigned to the
control parameters (options) by the user.
a.Program control Parameter Data. These data contain parameters which
provide user-control Parameter Data.
Parameter
Name
ISW(l)
ISW(2)
ISW(3)
Concentration/Deposition Option—Directs the program to
calculate either average concentration or total deposition. A
value of "1" indicates average concentration is to be
calculated and a value of "2" indicates total deposition is to
be calculated. If this parameter is not punched, the program
defaults to "1" or concentration.
Receptor Reference Grid System Option—Specifies whether a
right-handed rectangular Cartesian coordinate system or a
polar system is to be input to the program to form the
receptor reference grid system. A value of "1" indicates a
Cartesian- reference grid system is being input and a value of
"2" indicates a polar reference grid system is being input.
If this parameter is not punched, the program will default to
a value of "1." '
Discrete Receptor Option—Specifies whether a right-handed
rectangular Cartesian reference system or polar reference
system is used to reference the input discrete receptor
points. A value of "1" indicates that the Cartesian reference
system is used and a value of "2" indicates that a polar
reference system is used. If this parameter is not punched,
the program will default to a value of "1."
4-5
-------
Parameter
Name
ISW(4) Receptor Terrain Elevation Option—Specifies whether the user
desires to input the terrain elevations for each receptor
point or to use the program as a flat terrain model. A value
of "0" indicates terrain elevations are not to be input and a
value of "1" indicates terrain elevations for each receptor
point are to be input. Note that terrain elevations cannot be
used with the deposition model. The default for this
parameter is no terrain or "0." If equal to "-1," the program
assumes input elevations are in meters rather than feet.
ISW(5) Input/Output File Option—Specifies whether disc file input
and/or output is to be used. A value of "0" indicates no file
input or output. A value of "1" indicates an output file is
to be produced on the output unit specified by ISW(15). A
value of "2" indicates an input file is required on the input
unit specified by ISW(14). A value of "3" indicates both
input and output files are being used. Default for this
parameter is "0". It is the user's responsibility to ensure
that the correct tapes or files are mounted on the correct
units.
ISW(6) Print Input Data Option—Specifies what input data are to be
printed. A value of "0" indicates no input data are to be
printed. A value of "1" indicates only the control
parameters, receptor points and meteorological data are to be
printed. A value of "2" indicates only the source input data
are to be printed and a value of "3" indicates all input data
are to be printed. The default for this parameter is "0."
ISW{7) Seasonal/Annual Print Option—Specifies whether seasonal
concentration (or deposition) values are to be printed, or
annual values only, or both seasonal and annual values. An
ISW(7) value of "1" indicates only seasonal output is to be
printed, a value of "2" indicates only annual output is to be
printed, and a value of "3" indicates both seasonal and annual
output are to be printed. If this parameter is not punched or
is "0," the ptogram defaults to "3."
ISW(8) Individual/Combined Sources Print Option—Specifies whether
output for individual sources or the combined sources (sum of
sources) or both is to be.printed. An ISW(8) value of "1"
indicates output for individual sources only is to be printed,
a value of "2" indicates output for the combined sources only
is to be printed, and a value of "3" indicates output for both
individual and combined sources is to be printed. The default
for this parameter is "3." This parameter is used in
conjunction with the parameter NGROUP below. If NGROUP equals
"0,", all sources input to the program are considered for
output under ISW(8). However, if NGROUP is greater than "0,"
only those sources explicitly or implicitly defined under
NGROUP are considered for output unddr ISW(8). Also, a single
source defined under NGROUP is logically treated as combined
source output when ISW(8) equals "2" or "3."
4-6
-------
Parameter
Name
ISW(9>
ISW(IO)
ISW(ll)
Rural/Urban Option—Specifies whether rural or urban modes are
to be used (see Table 2-3). A value of "1" specifies Urban
Mode 1 and the E and F stability categories are redefined as
D. A value of "2" specifies Urban Mode 2 and stability
categories A and B are redefined as A, C becomes B, D becomes
C, and E and F become D. A value of "3" specifies the Rural
Mode and does not redefine the stability categories. The
rural Pasquill-Gifford dispersion curves are used with values
of 1 through 3. A value of "4" specifies Urban Mode 3, with
no.stability category adjustment and use of the urban Briggs
dispersion curves. If this parameter is not punched or is
"0," the program defaults to "3." If file input is used, the
program defaults to the value saved on file. The parameter
ISW(9) is only used for card input sources and/or tape input
sources when ISW(12) equals "I." It should be noted that the
use of Urban Modes 1 and 2 are not recommended for regulatory
purposes.
Maximum 10 Print Option—Specifies whether the maximum 10
values of concentration or deposition only are to be printed,
or the results of the calculations for all receptors only, or
both are to be printed. A value of "1" directs the program to
calculate and print only the maximum 10 values and receptors
according to ISW(ll) or ISW(12) below. Values at receptors
other than the maximum 10 are not printed if this option
equals "1." A value of "0" directs the program to print the
results of the calculations at all receptors; the maximum 10
values are not produced. A value of "2" directs the program
to print the results of. the calculations at all receptor
locations as well as the maximum 10. The default for this
parameter is "0." The ISCLT program will print less than 10
values in cases where there are less than 10 concentration
(deposition) values greater than zero calculated.
Maximum 10 Calculation Option 1—This option directs the
program to use one of two methods to calculate and print
maximum 10 concentration (or deposition) values. If this
option is .used, option ISW(12) must equal "0." The program
determines the maximum values and receptor locations from the
set of all receptors input. Method- 1: A value of "1" directs
the program to calculate and print the maximum 10 values and
respective receptors for each individual source and to
calculate and print the maximum 10 values and respective
receptors for the combined sources independently of each
other. The output for individual sources and combined sources
will in general show a different set of receptors. Method 2:
A value of "2" directs the program to first calculate and
print the maximum 10 values and respective receptors for the
combined sources (sum of sources) and then print the
contribution at each receptor of each individual source to the
combined sources maximum 10. This option can only be used if
one or more of the following conditions is met:
4-7
-------
Parameter
Name
Condition a - The run uses an output tape or data file
(user must specify NOFILE, if tape)
Condition b - The run uses an input tape or data file,
but has no input data card sources (all are
taken from tape; user must specify NOFILE,
if tape)
Condition c - The total number of input sources is less
than or equal to the minimum of I and J,
where
J = 300
and
I = (E - (Nx + Nv + 2NXV) - K-L-M) (4-1)
(Nse(NxNy + N»y»
E = the total amount of program data storage in BLANK
COMMON. The design size is 40,000.
Nx = Number of points in the input X-a:-:is of the
receptor grid system (NXPNTS)
Ny = Number of points in the input Y-axis of the
receptor grid system (NYPNTS)
Nxy= Number of discrete (arbitrarily placed) input
receptors (NXWYPT)
N,.= Number of seasons in the input meteorological data
(NSEASN)
K = N..(N,N,+N,y)
0 ; if ISW(4) = "0"
L = OR
N,Ny+Nxy; if ISW(4) = "1" or "-1"
OR
N,Ny+N,y; if ISW{25) = "1"
OR
2(N,Ny+NXy>; if ISW(4) and ISW(25) are
both non-zero.
M = 0; if ISW(4) = 1 or "-1" and ISW(ll) = 2 or
if ISW(7) = l,or NSEASN = l,or WGROUP = 0
N«Ny+Nxy; if ISW(4) = 0 or ISW(ll) ± 2 and if
ISW(7) = 1 and NGROUP = 0 and NSEASN = 1
4-8 12/87
-------
Parameter
Name
ISW(12)
ISW(13!
ISW(14)
ISW(15)
Maximum 10 Calculation Option 2—This option directs the
program to calculate concentration or deposition at a special
set of user supplied discrete (arbitrarily placed) receptor
points. If this option is used, option ISW(ll) must equal
"0." A value of "1" directs the program to expect to read
from 10 to 50 special receptors at which concentration or
deposition is to be calculated. If this option is selected
and 10 special receptors are input, both seasonal and annual
concentration or deposition values for individual sources and
combined sources are printed for the 10 user-specified
receptors. If more than 10 special receptors are input, the
program assumes the first 10 points are for season 1, the
second 10 points are for season 2, and the last 10 points are
for annual tables. This option requires the parameter NXWYPT
given below to be a multiple of 10. All input tape or data
file sources are recalculated with this option. Also, if an
input tape is being used, the receptor grid system, discrete
receptors and their elevations input from the tape are
discarded and the user inputs the new special set of receptor
points (with elevations if ISW(4) equals "1" or "-1" and
receptor heights above ground if ISW{25) = "1") via data card.
Print Output Unit Option—This option is provided to enable
the user to print the program output on a unit other than
print unit "6." If this value is not punched or a "0" is
punched, all print output goes to unit "6." Otherwise, print
output goes to the specified unit. Also, if this value is
punched non-zero positive, two end-of-file marks are written
at the end of the print file. If ISW(13) is a negative value,
the end-of-file marks are not written.
Optional File Input Unit Number—This option is provided to
enable the user to assign the unit number from which data are
read under ISW(5). If ISW(14) is not punched or is "0," the
program defaults to unit "2." If the input data are being
read from a mass-storage file, ISW(14) must be set to a
negative value. A positive value implies magnetic tape. Note
that ISW(14) is the internal file name used by the program to
reference the data file and must be equated with the external
file name used to assign the file (see Section 4.2.2).
Optional File Output Unit Number—This option is provided to
enable the user to assign the unit number to which tape or
output file data are written under ISW(5). If ISW(15) is not
punched or is "0", the program defaults to unit "3." If the
output data are being written to a mass-storage file, ISW(15)
must be set to a negative value. A positive value implies
magnetic tape. Note that ISW(15) is the internal file name
used by the program to reference the data file and must be
equated with the external file name used to assign the file
(see Section 4.2.2).
4-9
12/87
-------
Parameter
Name
ISW(16)
ISW(17)
ISW(18)
ISW(19)
ISW(20)
ISW(21)
ISW{22)
Print Output Paging Option—This option enables the user to
minimize the number of print output pages. A value of "1"
directs the program to minimize the output pages by not
starting a new page with each type of output table. If
this option is not punched or is "0". the program will
start each unrelated output table on a new page. The user
is cautioned not to exercise this option until familiar
with the output format because the condensed listing may be
confusing.
Lines Per Page Option—This option is provided to enable
the user to specify the number of print lines per page on
the output printer. The correct number of lines per page
is necessary for the program to maintain the output
format. If this value is not punched or is "0", the
program defaults to 57 print lines per page.
Optional Format for Joint Frequency of Occurrence—This
parameter is a switch used to inform the program whether it
is to use a default format to read the joint frequency of
occurrence of speed and direction (FREQ) or to input the
format via data card. If this option is not punched or is
"0", the program uses the default format given under FMT
below. If this option is set to a value of "1", the array
FMT below is read by the program.
Option to Calculate Plume Rise as a Function of Downwind
Distance—This option is applicable to all stack sources
and if set equal to "0" or not punched, the downwind
distance is not considered in calculating the plume rise.
If ISW(19) is set equal to "1", the plume rise calculation
is a function of downwind distance. ISW(19) is set to "0"
if the regulatory default option (ISW(22» is selected.
Option to Add the Briggs (1974) Stack-Tip Downwash
Correction to Stack Sources—This option is applicable to
all stack sources and if set equal to "0" or not punched,
no downwash correction is made. If ISW(20) is set equal to
"1", the Briggs (1974) downwash correction is applied to
the stack height for all stack sources. ISW(20) is set to
"1" if the regulatory default option (ISW(22)) is selected.
Buoyancy-Induced Dispersion Option—Allows the program to
modify the dispersion coefficients to account for
buoyancy-induced dispersion. A value of "0" directs the
program to modify the dispersion coefficients for
stack-type sources while a "1" directs the program to
bypass the modification. ISW(21) is set to "0" if the
regulatory default option (ISW(22)) is selected.
Regulatory Default Option—If chosen (this option is chosen
if ISW (22) = 0, otherwise ISW(22) should be set to 1), the
program will internally re-define some user defined input
4-10
-------
Parameter
Name
ISW(22)
Gont.
options to produce a simulation consistent with EPA regulatory
recommendations. The following features are ^.incorporated when
this option is selected:
1) Final plume rise is used at all downwind receptor locations.
2) Stack-tip downwash effects are included.
3) Buoyancy-induced dispersion effects are parameterized.
4) Default wind profile coefficients are assigned (.07, .07,
.10, .15, .35, .55, for the rural mode; and .15, .15, .20,
.25, .30, .30 for the urban modes).
5) Default vertical potential temperature gradients are assigned
(A:0.0, BrO.O, C:0.0, D:0.0, E:0.02, F:0.035 °K/m)
6) A decay half-life of 4 hours is assigned if SC>2 is modeled in
an urban mode; otherwise, no decay is assigned.
7) Revised building wake effects procedure is selected, which
uses either the method of Huber and Snyder, or that of
Schulman and Scire, depending on the stack height and
building dimensions (see Section 2.4.1.1.d).
Note that the model selects the appropriate urban or rural
mixing height, and that building downwash is calculated when
appropriate.
ISW(23) Pollutant Indicator Switch—If SC>2 is modeled the user should
set this option equal to "0". If a pollutant other than S02 is
modeled the user should set this option equal to "1". Note,
this switch is only used when ISW(22) = 0.
ISW(24) Input Debug Switch—If the user wants input data printed as soon
as it is entered set this option to "0", otherwise set this
option to "1". Note: any input data resulting from the
selection of ISW(6) will also be printed.
ISW(25) Above Ground ("flagpole") Receptor Option - Allows the user to
model receptor heights above local terrain. A value of "1"
directs the program to read user-provided receptor heights above
local terrain. The default value of "0" assumes no heights are
provided. This option is available regardless of the regulatory
default option setting.
NSOURC Number of Data Card Input Sources—This parameter specifies the
number of input card image sources. This includes card images
that specify a new source being entered and card images that
specify modifications or deletions to sources input from tape or
data file. If this value is not punched or is "0", the program
assumes all sources are input from tape or data file. Also, if
a negative value is punched for this parameter, the program will
continue to read source data card images until it encounters an
end-of-file or a negative source identification number in the
parameter NUMS below. There is no limit to the number of
sources the program can process when using tape output (see
(ISW(H)).
NGROUP Number of Source Combination Groups—This parameter is used to
select concentration (deposition) calculations for specific
4-11
12/87
-------
NGROUP sources or source combinations to be printed .under the
Cont. parameter ISW(8) above. A source combination consists of one
Parameter or more sources and is the sum of the concentrations
(deposition) calculated for those^sotlrces. If the user desires
only individual source output or only all sources combined or
both, the parameter NGROUP is not punched or is set equal to
"0" and ISW{8) is set according to which option the user
desires. Also, if NGROUP is not punched or is set equal to
"0", the parameters NOCOMB and IDSOR below are omitted from
the input data. However, if NGROUP is set greater than zero,
the program assumes the user desires to NGROUP restrict the
output of concentration tables to select individual sources or
select combinations of sources or both, depending on ISW(8).
The maximum value for NGROUP is 20. If more than 20 source
combinations are desired they must be produced in multiple
runs of ISCLT. This can be done by specifying an output tape
or data file on the first execution. The user would then use
this tape for input on subsequent runs to produce the remaining
desired source combinations. Also, only a few of the data
cards and values from the initial data deck are required on
subsequent runs. The parameter NGROUP cannot be used or
punched non-zero unless one or more of the following
conditions is met:
Condition a - The run uses an output tape or data file (user
must specify NOFILE, if tape)
Condition b - The run uses an input tape or data file, but has
no input data card sources (all are taken from
tape, NSOURC = "0") (user must specify NOFILE,
if tape)
Condition c - The total number of input sources (NSOURC +
input tape sources) is less than or equal to the
minimum of I and J, where
J = 300
and
I = [E - (N* + Ny + 2N» Ny) (4-2)
- K - L - M]/[N,.(N« Ny + Nxy)]
All of the variables in this equation except K are the same as those
defined under ISW(ll) above.
0 ; if ISW<8)=1 and
K = or
N,.(N,Ny+Niy); if ISW(8)*1 or ISW(11)=2
4-12 12/87
-------
Parameter
Nam'e
NXPNTS
X-Azis/Range Receptor Grid Size-This parameter specifies the
number of east-west receptor grid locations for the Cartesian
coordinate system X-axis, or the number of receptor grid
ranges (rings) in the polar coordinate system, depending on
which receptor grid system is chosen by the user under
parameter ISW(2). This is the number of X-axis points to be
input or the number of X-axis points to be automatically
generated by the program. A value of "0" (not punched directs
the program to assume there is no regular receptor grid being
used. The maximum value of this parameter is related to other
parameter values and is given by the equation
E > [N,+Ny+2N»y] + [(KN.+I)
(4-3)
where all variables in the above equation are the same as
those defined under ISW(ll) above except K and I, which are
defined as
1 ; if ISW(8)=1 and
K = or
2 ; if ISW(8)^1 or ISW(11)=2
0 ; if ISW(4)=0 (no terrain)
or
1=1 ; if ISW(4)=1 or "-1" or if ISW(25) = 1
or
2 ; if ISW(4) and ISW(25) are both non zero
This parameter is ignored by the program if tape or data file
input is being used.
NYPNTS Y-Axis/Azimuth Receptor Grid Size—This parameter specifies
the number of north-south receptor grid locations for the
Cartesian coordinate system Y-axis, or the number of receptor
azimuth bearings from the origin in the polar coordinate
system, depending on which receptor grid system is chosen by
the user under parameter ISW(2). If the parameter NXPNTS is
set non-zero, the parameter NYPNTS must also be non-zero. The
maximum value of this parameter is given by the equation under
NXPNTS above. The parameter NYPNTS is ignored by the program
if tape or data file input is being used.
NXWYPT Number of Discrete (Arbitrarily Placed) Receptors—This
parameter specifies the total number of discrete receptor
points to be input to the program. A value of "0" (not
punched) directs the program to assume no discrete receptors
are being used. This parameter must be set to a multiple of
10 if option ISW(12) is selected. Also, the maximum value of
this parameter is limited by the equation given under NXPNTS
above. This parameter is ignored by the program if input tape
or data file is being used, except in the case where the
ISW(12) option has been selected.
4-13
12/87
-------
Parameter
Name
NSEASN
NSTBLE
NSPEED
NSCTOR
NOFILE
Number- of Seasons—This parameter specifies the number of
seasons or months in the input meteorological data. A value
of "0" (not punched) defaults to "I". Also, if annual
meteorological data are being used, a value of "1" should be
specified. The maximum value of this parameter is "4". If
monthly STAR summaries and seasonal average mixing heights and
ambient air temperatures are used to calculate monthly
concentration or deposition values for each month of the year,
four separate program runs, each containing three "seasons"
(months), are required. This parameter is ignored by the
program if an input tape or data file is being used.
Number of • Paquill Stability Categories—This parameter
specifies the number of Pasquill stability categories in the
input joint frequency of occurrence of wind speed and
direction (FREQ). A value of "0" (not punched) causes the
program to default to "6" (maximum). This parameter is
ignored by the program if an input tape or data file is being
used.
Number of Wind Speed Categories—This parameter specifies the
number of wind speed categories in the input joint frequency
of occurrence of wind speed and direction (FREQ). A value or
"0" (not punched) causes the program to default to "6"
(maximum). This parameter is ignored by the program if an
input tape or data file is being used.
Number of Wind Direction Sector Categories—This parameter
specifies the number of wind direction sector categories in
the input joint frequency of occurrence of wind speed and
direction (FREQ). A value of "0" (not punched) causes the
program to assume the standard "16" (maximum) sectors are to
be used (see Section 2.2.1.2). This parameter is ignored by
the program if an input tape or data file is being used.
Tape Data Set File Number—This parameter specifies the output
tape file number or, if only an input tape is being used, the
input tape file number. This parameter is used by the ISCLT
program to position the tape at the correct file if multiple
passes through the data are required. This parameter must be
input if the user is using Condition a or Condition b under
ISW(ll) and/or under NGROUP. This parameter does not apply to
runs that use mass-storage (assumed one file) or runs that
satisfy Condition c under ISW(ll) and/or NGROUP. Also, the
user must position input and output tapes at the correct files
prior to executing the ISCLT program.
4-14
-------
Parameter
Name
NOCOMB Number of Sources Defining Combined Source Groups—This
parameter is not read by the program if the parameter NGROUP
above is zero or not punched. Otherwise, this parameter is an
array of NGROUP values where each value gives the number of
source identification numbers used to define a source
combination. The source identification numh_r is that number
assigned to each source by the user under the source input
parameter NUMS below. An example and a more detailed
discussion of the use of this parameter is given under IDSORC
below. A maximum of 20 values is provided for this array.
IDSORC Combined Source Group Defining Sources—This parameter is not
read by the program if the parameter NGROUP above is zero or
not punched. Otherwise, this parameter is an array of source
identification numbers that define each combined source group
to be output. The values punched into the array NOCOMB above
indicate how many source identification numbers are punched
into this array successively for each combined source output.
The source identification numbers can be punched in two ways.
The first is to punch a positive value directing the program
to include that specific source in the combined output. The
second is to punch a negative value. When a negative value is
punched, the program includes all sources with identification
numbers less than or equal to it in absolute value. Also, if
the negative value is preceded by a positive value in the same
defining group, that source is also included with those
defined by the negative number, but no sources with a lesser
source identification number are included. For example,
assume NGROUP above is set equal to 4 and the array NOCOMB
contains the values 3, 2, 1, 0. Also, assume the entire set
of input sources is defined by the source identification
numbers 5, 72, 123, 223, 901, 902, 1201, 1202, 1205, 1206, and
1207. To this point we have a total of 11 input sources and
we desire to see 4 combinations of sources taken from these
11. Also, the array NOCOMB indicates that the first 3 values
in the array IDSORC defines the first source combination, the
next 2 values (4th and 5th) in IDSORC define the second
combination, the 6th value in IDSORC defines the third
combination and the last combination has no defining (0)
sources so the program assumes all 11 sources are used.
Similarly, let the array IDSORC be set equal to the values 5,
72, -223, 1201, -1207, -902. The program will first produce
combined source output for source 5, and all so'urces from 72
through 223. The second combined source output will include
sources 1201 through 1207. The third will include source
numbers 1 through 902 and the last will- include all sources
input. Note that the source identification numbers in each
defining group are in ascending order of absolute value.
Also, if ISW{8> equals "2" (combined output only) and there
are groups with only one positive source number (individual
sources), the program logically treats these individual
sources as combined sources.
4-15
-------
Parameter
Name
FMT Optional Format for Joint Frequency of Occurrence—This
parameter is an array which is read by the program only if
ISW(18) is set to a value of "1". The array FMT is used to
specify the format of the }oint frequency of occurrences of
wind speed and direction data (FREQ, STAR summary,
fi.j.k.a in Table 2-4). The format punched, if used, must
include leading and ending parentheses. If ISW(18) is not
punched or is set to a value of "0", the parameter FMT is
omitted from the input deck and the program uses the default
format "(6F10.0)". This default format specifies that there
are 6 real values per card occupying 10 columns each,
including the decimal point (period), and the first value is
punched in columns one through ten. If the user has received
the STAR data from an outside source, the deck must also be
checked for the proper order as well as format.
b. Receptor Data These data consist of the (X,Y) or (range, azimuth)
locations of all receptor points as well as the elevations of the receptors
above mean sea level and heights of receptors above local terrain. The
minimum distance in meters between source and receptor for which calculations
are made is given by:
Stack Sources:
minimum distance =
1 ; no wake effects
or
MAX(1,3*HB) ; wake effects, squat building
or
MAX(1,3*HW) ; wake effects, tall building
Volume Sources:
minimum distance =
Area Sources:
minimum distance =
Where:
1 + 2.15*SIGYO
1 + 0.5*BW
HB = height of building (regular or direction specific)
HW = width of building (regular or direction specific)
SIGYO = standard deviation of the lateral source
dimension of building
BW = width of area source
4-16
12/87
-------
Parameter
Name
Receptor Grid System X-Axis or Range — This parameter is read
by the program only if the parameters NXPNTS and NYPNTS are
non-zero and only if an input tape or data file is not being
used. This parameter is an array of values in ascending order
that defines the X-axis or ranges (rings) (depending on
ISW(2)) of the receptor grid system in meters. If only the
first 2 values on the input card are punched and the parameter
NXPNTS is greater than 2, the program assumes the X-axis
(range) is to be generated automatically and assumes the first
value punched is the starting coordinate and the second value
punched is an increment used to generate the remaining NXPNTS
evenly-spaced points. If all receptor points are being input,
NXPNTS values must be punched.
Receptor Grid System Y-Axis or Azimuth — This parameter is read
by the program only if the parameters NXPNTS and NYPNTS are
non-zero and only if an input tape or data file is not being
used. This parameter is an array of values in ascending order
that defines the Y-axis or azimuth bearings (depending on
ISW(2)) of the receptor grid system in meters or degrees. If
only the first 2 values on the input card are punched (third
and fourth values are zero) and the parameter NYPNTS is
greater than 2, the program assumes the first value punched is
the starting coordinate and the second value punched is the
increment used to generate the remaining NYPNTS evenly-spaced
(rectangular or angular) points. If all receptor points are
being input, NYPNTS values must be punched. If polar
coordinates are being used, Y is measured clockwise from zero
degrees (north).
Elevation of Grid System Receptors — This parameter is not read
by the program if the parameter ISW(4) is zero or if an input
tape is being used or if NXPNTS or NYPNTS eguals zero. This
parameter is an array specifying the terrain elevation (feet
if ISW(4)=1, meters if ISW(4)=-1) above mean sea level at each
receptor of the Cartesian or polar grid system. There are
NXPNTS • NYPNTS values read into this array. The program
starts the input of values with the first Y coordinate
specified and reads the elevations for each X coordinate at
that Y in the same order as the X coordinates were input. A
new data card is started for each Y value and the NXPNTS
elevations for that Y are read. The program will expect
NYPNTS groups of data cards with NXPNTS elevation values
punched in each group. For example, assume we have a 5 by 5
Cartesian or polar receptor array. The values Zi through
Z$ are read from the first card group, the values Z«
through Zio from the second card group and Zzi through
from the last card group.
RHT Above Ground ("flagpole") Grid Receptor Heights - This
parameter is not read by the program if the parameter ISW(25)
is zero or if an input tape is being used or if NXPNTS or
NYPNTS equals zero. This parameter is an array specifying the
receptor heights (meters) above local terrain elevation at
each receptor of the Cartes ion or polar grid system. The
method of input is similar to that of Z described above.
4-17 12/87
-------
Parameter
Name
Rectangular
Z21 ,
Z6
Zl
Z22 .
Z7
Z2
Z23 ,
Z8
23
Z24 ,J
Z9
Z4
25
X
z
(Cont.)
- X5
- X4
- X3
- X2
- XI
4-18
-------
Parameter
Name
X Discrete (Arbitrarily Placed) Receptor X or Range—This
(Discrete) parameter is not read by the program if the parameter NXWYPT
is zero or if the program is using an input tape or data file
with the ISW(12) option set to zero. This parameter is an
array defining all of the discrete receptor X points. The
values are either east-west distances or radial distances in
meters, depending on the type of reference system specified by
ISW(3). NXWYPT points are read by the program.
Y Discrete (Arbitrarily Placed) Receptor Y or Azimuth—This
(Discrete) parameter is not read by the program if the parameter NXWYPT
is zero or if the program is using an input tape or data file
with the ISW(12) option set to zero-. This parameter is an
array defining all of the discrete receptor Y points in meters
and degrees. The values are either north-south distances or
azimuth bearings (angular distances) measured clockwise from
zero degrees (north depending on the type of reference system
specified by ISW(3). NXWYPT points are read by the program.
Z Elevation of the Discrete (Arbitrarily Placed) Receptors—This
(Discrete) parameter is not read by the program if the parameter ISW(4)
is zero or if the parameter NXWYPT equals zero of if an input
file is being used with the ISW(12) option equal to zero.
This parameter is an array specifying the terrain elevation
(feet if ISW(4)=1, meters if ISW(4)=-1) at each of the NXWYPT
discrete receptors.
RHT Above Ground ("flagpole") Discrete Receptor Heights - This
(Discrete) parameter is not read by the program if the parameter ISW(31)
is zero or if the parameter NXWYPT equals zero or if an input
file is being used with the ISW(12) option equal to zero.
This parameter is an array specifying the receptor heights
(meters) above local terrain elevation at each of the NXWYPT
discrete receptors.
c. Identification Labels and Model Constants. These data consist of
parameters pertaining to heading and identification labels and program
constants. These data, except for TITLE, are not read by the program if an
input tape or data file is being used.
Parameter
Name
TITLE Page Heading Label—This parameter is an array that allows up
to 80 characters of title information to be printed as the
first line of each output page.
UNITS Concentration/Deposition and Source Units Label—This
parameter is an array used for the optional input of two unit
labels. The first 40 characters of this array are provided
for an optional output units label for concentration or
deposition. This label is defaulted to "micrograms per cubic
meter" for concentration and "grams per square meter" for
deposition, if the parameter TK below is not punched or is
"0". The second 40 characters of this array are provided for
4-19 12/87
-------
Parameter
Name
UNITS an optional source input units label. This label is defaulted
Cont. to "grams per second" for concentration or "grams" for
deposition for stacks and volume sources and to "grams per
second per square meter" or "grams per square meter" for area
sources, if the parameter TK below is not punched or is "0".
ROTATE Wind Direction Correction Angle—This parameter is used to
correct for any difference between north as defined by the X,
Y reference grid system and north as defined by the weather
station at which the wind direction data were recorded. The
value of ROTATE (degrees) is subtracted from each
wind-direction sector angle (THETA). This parameter is
positive if the positive Y axis of the reference grid system
points to "the right of north as defined by the weather
station. Most weather stations record direction relative to
true north and the center of most grid systems are relative to
true north. However, some weather stations record direction
relative to magnetic north and the ends of some UTM (Universal
Transverse Mercator) zones are not oriented towards true
north. The user is cautioned to check the wind data as errors
in the wind direction distribution will lead to erroneous
program results. The default value of ROTATE is "0".
TK Model Units Conversion Factor—This parameter is provided to
give the user flexibility in the source input units used and
the concentration or deposition output units desired. This
parameter is a direct multiplier of the concentration or
deposition equation. If this parameter is not punched or is
set to a value of "0", the program defaults to "1 x 10s"
micrograms per gram for concentration and to "1" for
deposition. This default assumes the user desires
concentration in micrograms per cubic meter or deposition in
grams per square meter and the input source units are grams
per second or total grams for stack and volume sources and
grams per second per square meter or grams per square meter
for area sources, depending on whether the program is to
calculate concentration or deposition. Also, if the default
value for this parameter is selected, the program defaults the
unit labels in the array UNITS above. If the user chooses to
input this parameter for other units, he must also input the
units labels in UNITS above. This parameter corresponds to K
in Equations (2-51), (2-56), (2-57), and (2-58).
ZR Weather Station Recording Height—This parameter is the height
above ground level in meters at which the meteorological data
were recorded. If this parameter is not punched or has a
value of "0", the program defaults to "10" meters. This
parameter, corresponds to Z\ in Equation (2-1).
G Acceleration Due to Gravity—This parameter, which is used in
the plume rise calculations, is the acceleration due to
gravity. If this parameter is not punched or has a value of
4 20
-------
Parameter
Name
G "0", the program uses "9.8" meters per second squared as the
Cont. default value. This parameter corresponds to g in equation
(2-3).
DECAY Decay Coefficient—This parameter is the coefficient
(seconds'1) of time-dependent pollutant removal by physical
or chemical processes (Equations (2-20), (2-21)). If SOz is
modeled in an Urban Mode and the regulatory default option
(ISW(22)) is chosen, the program assigns a decay coefficient
corresponding to a half life of four hours. Otherwise,
pollutant decay is not considered.
d. Meteorological Data. These data are the meteorological input
parameters classified according to one or more of the categories of wind
speed, Pasquill stability, wind direction and season or annual. These
parameters are not read by the program if an input tape or data file is being
used.
FREQ Joint Frequency of Occurrence—This parameter array consists
of the seasonal or annual joint frequency of occurrence of
wind-speed and wind-direction categories classified according
to the Pasquill stability categories (STAR summary,
fi.j.k.a in Table 2-4). This parameter has no default and
must be input in "the correct order. The program begins by
reading the joint frequency table for season 1 (winter) and
stability category 1 (Pasquill A stability). The first data
card contains the joint frequencies of wind speed categories 1
through 6 (1 through NSPEED) for the first wind direction
category (north). The second data card contains the joint
frequencies of wind speed categories 1 through 6 for the
second wind direction category (north-northeast). The program
continues in this manner until the joint frequencies of the
last direction category (north-northwest) for stability
category 1, season 1 have been read. The program then repeats
this same read sequence for stability category 2 (Pasquill B
stability) and season 1. When all of the stability category
values for season 1 have been read, the program repeats the
read sequence for season 2, season 3, etc., until all of the
joint frequency values have been read. There are a total of
NSPEED«NSCTOR»NSTBLE»NSEASN data cards. If the total sum of
the joint frequency of occurrences for any season (or annual)
does not add up to 1, the program will automatically normalize
the joint frequency distribution by dividing each joint
frequency by the total sum. Also, the program assumes
stability categories 1 through 6 are Pasquill stabilities A
FREQfhrough F. Seasons 1 through 4 are normally winter,
spriny, summer and fall. See the parameter FMT above for the
format of these data.
-------
Parameter
Name
TA Average Ambient Air Temperature—This parameter array consists
of the average ambient air temperatures (Ta;i<,4 in
Table 2-4), classified according to season (or annual) and
stability category, in degrees Kelvin. One data card is read
for each season (1 to NSEASN) with the temperature values for
stability categories 1 through NSTBLE punched across the
card. When the program has completed reading these data
cards, it will scan all of the values in the order of input
and, if any value is not punched or is zero, the program will
default to the last non-zero value of TA it encountered.
HM Mixing Heights—This parameter array consists of the median
mixing layer height in meters (Hm/i,k,i in Table 2-4)
classified according to wind speed, stability and season (or
annual). The program begins reading the mixing layer heights
for season 1. The program reads the mixing layer height
•values for each wind speed category <1 to NSPEED) from each
card. There are NSTBLE (1 through NSTBLE) cards read for each
season. The program scans each value input in the order of
input and, for each season, if a zero or non-punched value is
found, the program defaults to the last non-zero value
encountered within the values for that season. The ISCLT
program automatically uses a mixing height value of 10000
meters for the E and F stability categories when the program
is run in the Rural Mode.
DPDZ Potential Temperature Gradient—This parameter array consists
of the vertical gradients of potential temperature (BQ/Bzl/k
in Table 2-4) classified according to wind speed and stability
category, in units of degrees Kelvin per meter. There are
NSTBLE (1 through NSTBLE) data cards read with the values for
wind speed categories 1 through NSPEED read from each card. A
value of 39/3z greater than zero indicates stable thermal
stratification and a value of 36/3z less than zero
indicates unstable thermal stratification. However, because a
blank input field is interpreted as zero, the program assumes
a zero input value means a default value is desired. Also,
because the same plume rise equation is used for adiabatic and
unstable conditions, a negative input value will direct the
program to use the plume rise equations for adiabatic or
unstable thermal stratification. If the first value on a data
card is not punched or is zero, a default value is used that
depends on the stability category. -If the stability category
is A, B, C or D, the value is left as a zero and the
adiabatic/unstable plume rise equation is used. However, if
the stability category is E or F, the value defaulted is
0.02 degrees Kelvin per meter for E and 0.035 degrees Kelvin
per meter for F stability. When any of the second through
sixth values of DPDZ on a data card are input as a zero or are
blank, the program will default to the previous value on the
data card. If the regulatory default option is selected
(ISW(22)=0) the default values will override any user input
values.
4-22
-------
Parameter
Name
UBAR
THETA
Wind Speed—This parameter array consists of the median wind
speeds in meters per second (ui in Table 2-4) for the wind
speed categories used in the calculation of the joint
frequency of occurrence of wind speed and direction (STAR
summary). There are NSPEED values read from this card. If
any value is not punched or is zero, the program defaults to
the following set of values: 1.5, 2.5, 4.3, 6.8, 9.5 and 12.5
meters per second.
Wind Direction—This parameter array consists of the median
wind direction angles in degrees for the wind-direction
categories used in the calculation of the joint frequency of
occurrence of wind speed and direction (STAR summary). There
are NSCTOR values read from 1 to 2 data cards and if the first
two values of this array are not punched or are zero, the
program defaults to the following standard set of values: 0,
22.5, 45, 67.5, 90, . . . , 337.5 degrees (N, NNE, NE
NNW). The wind direction is that angle from which the wind is
blowing, measured clockwise from zero degrees (north).
Wind Speed Power Law Exponent—This parameter array consists
of the wind speed power law exponent (p in Equation (2-1))
classified according to wind speed and stability categories 1
through NSTBLE. If the first value on any data card in this
set is not punched or is zero, the program defaults to the
value from the following set of values: Rural A = .07, B =
.07, C = .10, D = .15, E = .35, F = .55; Urban A = .15, B =
.15, C = .20, D = .25, E = .30, F = .30 depending on the
stability category A through F. Also, if any of the second
through last values on a card is not punched or is zero, the
value is defaulted to the previous value on the data card. If
a negative value is input, the result is a wind speed power
law exponent of zero. If the regulatory default option is
selected (ISW(22)=0) the default values will override any
user-input values.
e. Source Data. These data consists of all necessary information required
for each source. These data are divided into three groups: (1) parameters
that are required for all source types, (2) parameters that are required for
stack type sources, and (3) parameters that are required for volume sources
and area sources. The order of input of these parameters is given at the end
of this section.
4-23
-------
Parameter
Name
NUMS Source Identification Number—This parameter is the source
identification number and is a 1- to 5-digit integer. If this
number is negative, the program assumes NUMS is only a flag to
terminate the card source input data. Also, if NUMS is not
punched or is zero, the program will default NUMS to the
relative sequence number of the source input. This number
cannot be defaulted if source data are also being input from
tape or data file. Sources must be input in ascending order
of the source identification number.
DISP Source Disposition—This parameter is a flag that tells the
program what to do with the source. If this parameter is not
punched or has a value of "0", the program assumes this is a
new source for which concentration or deposition is to be
calculated. Also, if the program is using an input tape or
data file, this new source will be merged into the old sources
from file or will replace a file source with the same source
identification number. If the parameter DISP has a value of
"1", the program assumes that the file input source having the
same source identification number is to be deleted from the
source inventory. The program removes the source as well as
the concentration or deposition arrays for the source. If the
parameter DISP has a value of "2", the program assumes the
source strengths to be read from data card for this source are
to be used to rescale the concentration or deposition values
of the tape input source with the same source identification
number. The new source strengths input from card replace the
old values taken from the input tape and the concentration or
deposition arrays taken from tape are multiplied by the ratio
of the new and old source strengths. The DISP option equal to
"2" can only be used if QFLG equals zero and the tape input
source has QFLG equal to zero.
TYPE Source TYPE—This parameter is a flag that tells the program
what type of source is being input. If this parameter is not
punched or is "0", the program assumes a stack source. If
this parameter has a value of "1", the program assumes a
volume source. Similarly, if this parameter has a value of
"2", an area source is assumed.
QFLG Source Emission Option—This parameter is a flag that tells
the program how the input source emissions are varied. If
this value is not punched or is "0", the program assumes the
source emissions vary by season (or annual) and only NSEASN
values are read by the program. If this parameter has a value
of "1", the program assumes the source emissions vary by
stability category and season. If this parameter has a value
of "2", the program assumes the source emissions vary by wind
speed category and season. If this parameter has a value of
"3", the program assumes the source emissions vary hy wind
speed category, stability category and season. The order of
input of the source strengths under each of these options is
discussed under the parameter Q below.
4-24
-------
Parameter
Name
DX Source X Coordinate—This parameter gives the Cartesian X
(east-west) coordinate in meters of the source center for
stack and volume sources and the southwest corner for area
sources (X in Table 2-6) relative to the origin of the
reference grid system being used.
DY Source Y Coordinate—This parameter gives the Cartesian Y
(north-south) coordinate in meters of the source center for
stack and volume sources and the southwest corner for area
sources (Y in Table 2-6) relative to the origin of the
reference grid system being used.
H Height of Emission—This parameter gives the height above
ground in meters of the pollutant emission. For volume
sources, this is the height to the center of the source.
ZS Source Elevation—This parameter gives the terrain elevation
in meters above mean sea level at the source location and is
not used by the program unless receptor terrain elevations are
being used.
Q Source Emission—This parameter array gives the emission rate
of the source for each category specified by QFLG above. If
QFLG above is "0", NSEASN values are read from one data card.
IF QFLG is "I", NSEASN data cards are read with the source
emission values for stability categories 1 through NSTBLE read
from each card. If QFLG is "2", NSEASN data cards are read
with the source emission values for wind speed categories 1
through NSPEED read from each card. If QFLG is "3", NSPEED (1
through NSPEED) source emission values are read from each data
card and there are NSTBLE (1 through NSTBLE) data cards read
for each season. There are no default values provided for the
parameter Q and the program assumes "0" is a valid source
emission, the input units of source emission are:
PARAMETER Q
Source Type
Stack or
Volume
Area
*Default units
Concentration
Deposition
mass per unit time
(g/sec)*
mass per unit time
per unit area
(g/sec*m2))*
total mass
(g)*
total mass per unit
area
(g/m2)*
4-25
-------
Parameter
Name
NVS
VS
FRQ
GAMMA
Number of Particulate Size Categories—This parameter gives
the number of particulate size categories in the particulate
distribution used in calculating ground-level deposition or
concentration with deposition occurring. If ground-level
deposition (ISW{1) = "2") is being calculated, this parameter
must be punched and has a maximum value of 20. Also, if the
program is calculating concentration and this value is punched
greater than zero, concentration with deposition occurring is
calculated. If the parameter NVS is greater than zero, the
program reads NVS values for each of the parameter variables
VS, FRQ and GAMMA below.
Settling Velocity—This parameter array is read only if NVS
above is greater than zero. This parameter is the settling
velocity in meters per second for each particulate size
category (1 through NVS). No default values are provided for
this parameter.
Mass Fraction of Particles—This parameter is read only if NVS
above is greater than zero. This parameter is the mass
fraction of particulates contained in each particulate size
category (1 through NVS). No default values are provided for
this parameter.
Surface Reflection Coefficient—This parameter array is read
only if NVS above is greater than zero. This parameter is the
surface reflection coefficient for each particulate size
category (1 through NVS). A value of "0" indicates no surface
reflection (total retention). A value of "1" indicates
complete reflection from the surface. The reflection
coefficient range is from 0 to 1 and no default values are
provided.
Stack Source
Parameters
TS
VEL
Stack Gas Exit Temperature—This parameter gives the stack gas
exit temperature (T$ in Table 2-6) in degrees Kelvin. If
this parameter is zero, the exit temperature is set equal to
the ambient air temperature. If this parameter is negative,
its absolute value is added to the ambient air temperature to
form the stack gas exit temperature. For example, if the
stack gas exit temperature is 15 degrees Celsius above the
ambient temperature, enter TS as -15 {the minus sign is used
by the program only as a flag).
Stack Gas Exit Velocity—This parameter gives the stack gas
exit velocity in meters per second.
Stack Diameter—This parameter gives the inner stack diameter
in meters and no default is provided.
4-26
-------
Stack-Source
Parameters
HB Building Height—This parameter, gives the height above ground
level in meters of the building adjacent to the stack. This'
parameter and BW below control the wake effects option. If HB and
BW are punched non-zero, wake effects for the respective source
are considered. A negative value of HB (or the selection of the
regulatory default option) instructs the program to use the
revised building wake effects procedures, which uses either the
methods of Huber and Snyder or those of Schulman and Scire,
depending on the stack height to building height ratio (see
Section 2.4.1.1.d).
BW Building Width—This parameter gives the width in meters of the
building adjacent to the stack. If the building is npt square,
input the dimension of a square building of equal horizontal
area. If HB is not punched or is zero, this value should not be
punched. The effective width used by the program is the diameter
of a circle of equal area to the square of the side length BW.
Regulatory applications generally require the use of the "maximum
projected width". This can be accomplished by setting BW = 0.886
MPW where MPW is the maximum projected width.
WAKE Supersquat Building Wake Effects Equation Option—This option is
used to control the equations used in the calculation of the
lateral virtual distance (Equations (2-37) and (2-38)) when the
effective building width to height ratio (BW/HB) is greater than
5. If the parameter is not punched or has a value of "0" and the
width to height ratio is greater than 5, the program will use
Equation (2-37) to calculate the lateral virtual distance
producing the upper bound of the concentration or deposition for
the source. If this parameter has a value of "1", the program
uses Equation (2-38) producing the lower bound of the
concentration or deposition for the source. The appropriate value
for this parameter depends on building shape and stack placement
with respect to the building (see Section 2.4.1.1.d).
DSBH Direction Specific Building Height—This parameter array is read
only when the Schulman-Scire wake effects method is used (see
Section 2.4.1.1.d). The values are building heights for aaqh wind
sector flow vector starting at north and proceeding clockwise to
north-northwest. A total of NSCTOR values are read. Negative
values of DSBH provide the same option as WAKE on Card Group 17.
DSBW Direction Specific Building Width—This parameter array is read
only when the Schulman-Scire wake effects method is used (see
Section 2.4.1.1.d). The values are building widths for each wind
sector flow vector starting at north and proceeding clockwise to
north-northwest. A total of NSCTOR values are read.
Volume Source
Parameters
SIGYO Standard Deviation of the Crosswind Distribution This
parameter gives the standard deviation of the crosswind
distribution of the volume source (a
yo
in Table 2-6) in
4-27
12/87
-------
SIGYO meters. See Section 2.4.2.3 to determine the correct value
(Cont.) for this parameter. No default value is provided.
SIGZO Standard Deviation of the Vertical Distribution—This
parameter gives the standard deviation of the vertical
distribution of the volume source (ozo in Table 2-6) in
meters. See Section 2.4.2.3 to determine the correct value
for this parameter. No default value is provided for this
parameter.
Area Source
Parameters
XO Width of Area Source—This parameter gives the width of the
area source (x0 in Table 2-6} in meters. This parameter
should be the length of one side of the approximately square
area source. No default is provided for this parameter.
f. Source Data Input Order. There are from one to four data input card
groups of one or more cards each required to input the source data. The data
cards and parameters required depend on the source type (TYPE) and on the para-
meters DISP, QFLG, NVS and the concentration/deposition option parameter
ISW(l). Card Group 17 is always included in the input deck for each source
input (1 to NSOURC). Card group 17a through 17c are included only if NVS on
Card Group 17 is non-zero. Card Group 17ca and 17cb are included only if HB
on Card Group 17 is negative or if the regulatory default mode (ISW(22)) has
been selected. Card Group 17d is included only if DISP on Card Group 17 equals
"0" or "2". The order of input of these source cards is Card Group 17 followed
by those used from 17a through 17d for each successive source input. DO NOT
group all of 17 together, all of 17a together, etc. or the program will
terminate in error.
Source Input
Card Group 17
Required Source Parameters for Card Group 17—The parameters read from
the first data card for each source and their order are:
Stack Sources — NUMS, DISP, TYPE, QFLG,, DX, DY, H,
ZS, TS, VEL, D, HB, BW, WAKE, NVS
Volume Sources - NUMS, DISP, TYPE, QFLG, DX, DY, H.
ZS, SIGYO, SIGZO, NVS
Area Sources NUMS, DISP, TYPE, QFLG, DX, DY, H.
ZS, XO, NVS
4-28 12/87
-------
If the parameter DISP on this card is set to value of "0",
all parameters on this card are expected to have the correct
value- and the program may read Card Groups 17a, 17b and 17c
(depending on NVS), 17ca and 17cb (depending on either
ISW(22) = "0".or HB < 0 ), and will read Card Group 17d. If
DISP is set to a value of "1", only the parameters MUMS and
DISP are referenced (required) on this card, the program
assumes it is to delete an incoming tape or data file source
and only this data card is read for this source. If DISP is
set up to a value of "2", only the parameters MUMS, DISP and
QFLG are referenced (required) on this card because the
program assumes it is to read the source strengths from Card
Group 17d and to rescale the concentration or deposition of
an incoming tape or data file source. Parameters not
referenced on this first data card are set from tape or date
file source data by the program.
Source Input
Card Groups
17a, 17b,
and 17c
Source Particulate Distribution Data—This card group
consists of three sets of one or more data cards each and is
read by the program only if DISP is set to "0" and the
parameter NVS is set to a value greater than zero for
concentration calculations with deposition occurring or for
deposition calculations. The first data card(s) contains the
values of the parameter array VS, the second contains the
values of the parameter array FRQ and the third contains the
values of the parameter array GAMMA. A total of NVS values
are read from each set of cards.
Source Input
Card Groups
17ca and 17cb
Source Input
Card Group 17d
Direction Specific Building Dimensions - This card group
consists of two sets of cards.each and is read by the program
if HB is negative on the source card or if the regulatory
default mode (ISW(22) = "0") has been selected. The first
set of cards contains the values of the parameter array DSBH
and the second card set contains the values of the parameter
array DSBW. A total of NSCTOR values are read from each set
of cards.
Source Emissions—the last input card group for a source
contains the source emission values for the source. This
card group consists of one or more data cards and is read
4-29
12/8:
-------
only if the parameter DISP is not equal to "1". . The number
of cards required and the order of values input depends on
the parameters QFLG and is given under the source strength
parameter 0_ above.
4.1.3. Output Information
The ISCLT program generates five categories of program output. Each
category is optional to the user. That is, the user controls what output
other than warning and error messages the program generates for a given run.
In the following paragraphs, each category of output is related to the
specific input parameter that controls the output category. All program
output are printed except for magnetic tape or data file output.
a. Input Parameters Output. The ISCLT program will print all of the input
data except for source data if the parameter ISW(6) is set equal to a value of
"1" or '3". An example of this output is shown in Appendix D.
b. Source Parameters Output. The ISCLT program will print the input card
and tape source data if the parameter ISW(6) is set to a value of "2" or "3".
An example of the printed source data is shown in Appendix D.
4-29a
-------
c. Seasonal/Annual Concentration or Deposition. The parameter ISW(l)
specifies whether the program is to calculate concentration or deposition and
the parameter NSEASN specifies if seasonal or annual input meteorological data
is being used. The option ISW(7) is used to specify whether seasonal output
or annual output or both is to be generated. If the input meteorological data
are seasonal (winter, spring, summer, fall), the program can be directed to
produce tables of seasonal as well as annual concentration or deposition by
setting the parameter ISW(7) equal to "0" or "3". Also, only seasonal tables
are produced if ISW{7) equals "1". If the parameter NSEASN is set equal to a
value of "1" and only annual output is selected (ISW(7)="2"), the program
labels the output concentration or deposition as annual calculations.
However, if seasonal output is selected with NSEASN equal to "1", the output
tables are labeled seasonal. Also, all seasonal output is labeled season 1,
season 2, etc., requiring the user to keep track of the actual meteorological
season. Example Annual output tables are shown in Appendix D.
d. Concentration or Deposition Printed for the Maximum 10 and/or All
Receptor Points. The ISCLT program is capable of printing the concentration
or deposition calculations for each receptor point input to the program or
printing only the maximum 10 of those receptors or both. The parameter
ISW(IO) is used to determine which calculations are to be printed. Examples
of output tables giving the calculations at all points and the maximum 10 are
given in Appendix D.
e. Magnetic Tape or Data File Output. The ISCLT program will write all
input data and all concentration (deposition) calculations to magnetic tape or
data file. These data are written to the logical unit number specified by the
parameter ISW(15). This tape or data file must be assigned to the run prior
to the execution of the ISCLT program, positioned to the correct file and must
be equated to the logical unit number given in ISW(15). ISW(15) must be a
4-30
-------
positive value for magnetic tape or a negative value for mass storage. If
seasonal meteorological input data are used, the program saves only seasonal
concentration (deposition) on the output file and if input is annual, only
annual calculations are saved. This output file can be read back into the
ISCLT program to print tables not output in the original run and/or to modify
the source inventory for corrections or updates in the source emissions.
4.2 User's Instructions for the ISCLT Program
4.2.1 Program Description
The ISC long-term (ISCLT) program is designed to calculate average
concentration or total deposition values produced by emissions from multiple
stack, volume and area sources. The concentration or total deposition values
can be calculated on a seasonal (monthly) or annual basis or both for an
unlimited number of sources. The program is capable of producing the seasonal
and/or annual results for each individual source input as well as for the
combined (summed) seasonal and/or annual results from multiple groups of
user-selected sources. The program calculations of concentration or
deposition are performed for an input set of receptor coordinates defining a
fixed receptor grid system and/or for discrete (arbitrarily placed) receptor
points. The receptor grid system may be a right-handed Cartesian coordinate
system or a polar coordinate system. In either case, zero degrees (north) is
defined as the-positive Y axis and ninety degrees (east) is defined as the
positive X axis and all • points are relative to a user-defined hypothetical
origin (normally X=0, Y=0), although the Universal Transverse Mercator (UTM)
coordinates may be used as the Cartesian coordinate system).
4-31 12/87
-------
The ISCLT computer program is written in ANSI FORTRAN-77 and is designed
to execute on most medium to large scale computers with minimal or no
modifications. The program requires approximately 80,000 words (UNIVAC 111C)
of executable core for instruction and data storage. The program design
assumes a minimum of 32 bits per variable word and a minimum of four character
bytes per computer word. The program also requires from two to four
input/output devices, depending on whether the tape input/output options are
used. Input card image data is referenced as logical unit 5 and print output,
which requires 132-character print columns, is referenced as logical unit 6.
The optional tape or data file input is referenced as logical unit 2 and the
output is referenced as logical unit 3. The user has the option of either
using the default logical unit numbers given here or specifying alternate
logical unit numbers. The computer program consists of a main program (ISCLT)
and 22 subroutines as shown in Appendix. F. The FORTRAN source code for the
entire model is given in Appendix B.
4.2.2 Data Deck Setup
The card image input data required by the ISCLT program depends on the
program options desired by the user. The data may be partitioned into five
major groups as shown in Figure 4-1. The five groups are:
1. Title Record (1 data card)
2. Program Option and Control Records (2 to 5 Records)
3. Receptor Data Records (the number of records included in this
group depends on the parameters ISW(4), ISW(5), ISW(12),
ISW(25), NXPNTS, NYPNTS and NXWYPT)
4. Meteorological Data (only if ISW(5) is less than or equal to 1)
5. Source Data Cards (this record group is included only if NSOURC
is greater than zero)
4-32 12/87
-------
(5)
NUHS, DISP, etc. (this deck consists
of all source data cards (Card
Group 17) and is included in the
data deck only if NSOURC > 0).
(3)
(4)
FMT (this deck consists of parameter
card groups FMT (group 9} through
parameter card group P (group 16)
and is included in the data deck
only if ISW(5) <_ 1)
|
XDIS,YDIS,ZDIS,RHT (discrete receptors;
|
RHT (grid system receptor height deck)
f
Z (grid system elevations deck)
Y (grid system Y-axis deck)
X (grid system X-axis deck)
| UNITS (read only if ISW(5) <_ 1)
| IDSORC (read only if NGROUP > 0)
| NOCOMB (read only if NGROUP > 0)
(NSOURC, NGROUP, NXPNTS, etc.
\
ISW
(I)
TITLE
FIGURE 4-1. Input data deck setup for the ISCLT program.
4-33
-------
4.2.3 Input Data Description
Section 4.1.2 provides a summary description of all input data parameter
requirements for the ISCLT program. This section provides the user with the
FORTRAN format and order in which the program requires the input data
parameters. The input parameter names used in this section are the same as
those introduced in Section 4.1.2. Two forms of data may be input to the
program. One form is card image input data (80 characters per record) in
which all required data may be entered. The other form is magnetic tape or
mass storage. Both forms of input are discussed below.
a. Card Input Requirements. The ISCLT program reads all card image input
data in a fixed-field format with the use of a FORTRAN "A", "I" or "F" editing
code (format). Each parameter value must be punched in a fixed-field on the
data card defined by the start and end card columns specified for the
variable. Table 4-4 identifies each variable by name and respective card
group. Also, Table 4-4 specifies the card columns (fixed-field) for the
parameter value and the editing code ("A", "I" or "F" for alpha-numeric,
integer and real variables, respectively) used to interpret the parameter
value.
Card Group 1 in Table 4-4 gives the print output page heading and is
always included in the input data deck. Any information to identify the
output listing or data case may be punched into this card. If the card is
left blank, the heading will consist of only the output page number or the
heading will be taken from the input tape or data file, if used.
Card Group 2 gives the values of the program option array ISW. This card
is always included in the input data deck. However, the values of ISW(l)
through ISW(4) are automatically set by the program if you are using an input
(source/concentration or deposition inventory) tape. The options on this card
that determine whether or not some card groups are included in the input data
4-34
-------
s
CO
p/
W
U
2 2
Q; O
< M
a. H
Oi
£«
O* CJ
sa
a
§ Q
< 2
CJ IK.
S H
8i
OS fri
cu
EH
J
U
I-H
C
o
JJ
cu
Ll
U
en
cu
Q
CO
<. O JJ
2 o
CU
JJ
13 01
•H >,
Li 01
0*
'O
Ll -H
O Li
JJ 01
Cu
(V Li
U O
CO JJ
Li CU
CU
CU 0
JJ CU
s u
•H CU
T3 JJ
O CJ
o •*
u -o
U
C 0
to O
•H U
01
CU Ll
JJ 10
tO O
O QJ
II II
i— i r^
Li
0
0
X
c
(0
1— 1
JQ
rH
l-l
^<
£
«••*
(M
2
00
M
•0
cu
U
1— 1
Ql
^
,— 4
•H
a)
Li
JJ
A
'Tj
rtj
cu
JJ
cu
Li
U
01
• H
•o
to
•H
OI
cu
JJ
m
U
II
1— 1
t_l
0
o
X
c
-H -H
•"H *Tj (TJ
a, > >
cu cu
C rH rH
•H CU 0)
(0
Li C C
Ll "H -H
cu to to
1 1 (j JJ
t-l t-l
O cu cu
C JJ JJ
H M ii
O i— 1 i— 1
1
s^
o
X
C
(0
rH
jQ
rH
M
00
4C
JIB*.
^t
«•*>
3
CO
M
cu
Cu
fO
JJ
JJ
Qj >1
JJ rH
3 C
O 0
Li 0)
O Cu
<0
^ 1 1 1
3
CU JJ
C 3
JJ
O 3
C O
II II
O i-H
Li
O
X
C
(0
rH
rH
M
0
rH
^*
m
^
CO
M
,
01
CU
0J
(0
JJ
JJ
3
jJ
§
>i
"c "c
o to
0) JJ
Oj 2
m cu
JJ C
•rl
Qj JJ
G O
•H jQ
II II
(M ro
tO
JJ
€ >,
rH
•a jj c
cu 3 O
JJ Qj
c c to
• H -H JJ
S-i fT3
CU O> T3 (0
U JJ
jj sj jj to
O 3 3 -O
C O Cu
01 C JJ
Q) -H 3
Li JJ Qj
(0 3 0) C
43 CJ -H
ro t-i
JJ i-H 3 rH
tO r- 1 O r- 1
T3 fO 01 rO
3 C C C
Pl -H -H -H
C L, S-l S-I
•H Qj CU Qj
II II II II
O rH CM fO
S-i
o
X
c
ro
rH
rH
M
fM
rH
.••*
vD
3
CO
M
•H
*-•*
n
u
O
rM
II
in
'*.-•
[^
CO
i— i
"^
4J
Oj
£
• H
CU
Qj
ra
jj
u-i
•H
T3
CU
C
(0
^^
U
0)
_Q
1 i
o
c
m
u
T3
C
m
'U
£
S-i
01
O
Cu
CU
.c
JJ
£*1
XJ
r-H
. — 1
U
•H
4J
rfl
O
JJ
3
(0
JJ
cu
01
cu
Li
rO
01
Li
CU
JJ
cu
£ .
tO T3
L, a)
Qj 3
CU 01
W C
CU -(
jC cu
H .a
*
4-35