United States
Environmental Protection
Agency
Region 10
1200 Sixth Avenue
Seattle WA
Alaska
Idaho
Oregon
Washington
Environmental Services Division July 1986
EPA-910/9-86-144
c/EPA User's Guide for POSTZ
A Post-Processor for the SHORTZ Air
Quality Model
-------
EPA-910/9-86-144
July 1986
USER'S GUIDE FOR, POSTZ
A POST-PROCESSOR FOR THE
SHORTZ AIR QUALITY MODEL
. Prepared by
Kirk D. Winges
EPA Contract No. 68-02-3886
Project Officer
Robert B. Wilson
U. S. Environmental Protection Agency, Region 10
1200 Sixth Ave.
Seattle, Washington 98101
TRC Environmental Consultants, Inc.
15924 22nd Ave. SE
Mill Creek, Washington 98012
(206) 485-2992
-------
DISCLAIMER
This report has been reviewed by Region 10, U. S.
Environmental Protection Agency, and approved for publication.
Approval does not signify that the contents necessarily
reflect the views and policies of the U. S. Environmental
Protection Agency, nor does mention of trade names or
commercial products consitute endorsement or recommendation
for use.
-------
USER'S GUIDE FOR POSTZ
A POST-PROCESSOR FOR THE SHORTZ AIR QUALITY MODEL
TABLE OF CONTENTS
1 .0 INRODUCTION 1
1 .1 General 1
1.2 About SHORTZ 2
1.3 Running SHORTZ to Create an Output
Tape for POSTZ 3
2 . 0 POSTZ USER' S INSTRUCTIONS 5
2.1 . Program Options.' 5
2.1.1 Output Tables 5
2.1.2 Averaging Times 6
2.1.3 Calms Policy 6
2.1.4 Running or Block Averages 7
2.1.5 Restricting Time Limits. . 7
2.1.6 Restrict Receptors 7
2.1.7 Restrict Sources 8
2.1.8 Specify Background 8
2.1.9 Scale Concentrations 8
2.2 Model Parameters and Definitions 8
2.3 Format for the Input File 14
2.4 Executive Control language (ECL) for
the POSTZ Program 14
3 . 0 SAMPLE INPUT AND OUTPUT FILES 27
APPENDIX A COMPLETE FORTRAN LISTING FOR THE
POSTZ PROGRAM (UNIVAC Version)
-------
1.0 INTRODUCTION
1.1 General
The POSTZ computer program is a post-processor for the SHORTZ
computer program. The SHORTZ air quality model is designed
with the capability for producing a sequential tape output of
hourly concentrations at each receptor for each source. The
SHORTZ program can use this tape to "re-do" a particular run,
with a few changes. Essentially, the program saves the time
of repeating all the analyses by reading the results of a
previous analysis, and only making changes where requested by
the user, such as adding a new source. However, , it is also
possible to use this sequential tape to do other analyses,
such as:
o printing summary tables of concentrations for
different averaging times,
o restricting the analysis to certain receptors or
to a shorter portion of the meteorological data
base than the original SHORTZ run,
o restricting the analysis to certain sources, or
scaling the contribution from particular
sources,
o specifying separate background concentrations
for individual receptors,
o changing the manner of computation for averages
from block averages to running averages, and
o implementing the EPA calms policy.
Until now, no computer program was available to read the
SHORTZ output tape and conduct these analyses. The POSTZ
program is specifically designed to provide this capability.
This user's guide has been organized into 3 major sections,
including this introduction. In the introduction, information
is provided on the SHORTZ Model and the necessary information
which must be provided to the SHORTZ Model to prepare the
output tape for the POSTZ program. In Chapter 2.0, specific
user instructions for the POSTZ program are provided. In
Chapter 3.0 sample input and output streams are provided.
Appendix A has been included which presents the entire FORTRAN
listing for the POSTZ program.
-------
1.2 About SHORTZ
The SHORTZ air quality model was developed by the H. E. Cramer
Company. The program and it capabilities are discussed in
detail in the "User's Instructions for the SHORTZ and LONGZ
Computer Programs, Volumes I and II", EPA Document Numbers
EPA-903/9-82-004a and b. The SHORTZ program was written
specifically to run .on the Sperry/UNIVAC computer. The
program in the form on the UNIVAC, is not transportable to
other computer systems, .since it requires two specific UNIVAC
assembly language routines to operate. The routines allow the
program to define dynamic mass storage on the fixed disk
system. The information written to mass storage is not saved
by the system, but rather used only in the operation of the
computer program and deleted after a successful run.
In addition to the mass storage, the user can' elect to write
output to a magnetic tape. The magnetic tape output is very
comprehensive, and includes essentially "all the detailed
information of the sources, receptors and model parameters,
including specification of all the input parameters for the
SHORTZ run. The program also writes to tape a sequential
listing of all the concentrations computed for each hour,
including the concentrations by each source for each receptor
for each hour.
The POSTZ Model was originally developed on a micro computer
with a Microsoft Disk Operating .System (MS-DOS) . One of the
early steps in the POSTZ program development was the creation
of an MS-DOS version of the SHORTZ computer program. The
program generally uses the same FORTRAN code as the UNIVAC
version with several distinct differences. The mass storage
subroutines used by SHORTZ on the UNIVAC were replaced by
dynamic file opening control used by the MS-DOS FORTRAN. The
mass storage is performed on the default disk drive for the
system. Another major change to the program from the UNIVAC
version concerns the tape input and output routine. The MS-
DOS version of the program writes an output file very similar
in format to the output tape written by the UNIVAC version,
but the output media is the default disk drive, rather than a
tape unit. This necessitated major changes to the output
handling routine of the SHORTZ program. Also the capability
to read input from tapes was not germane to the POSTZ program
development and was not included in the capabilites of the MS-
DOS version of SHORTZ. There were a number of other minor and
cosmetic changes to the SHORTZ program from the UNIVAC version
to the MS-DOS version.
The MS-DOS version uses roughly 234 kilobytes of memory on an
IBM-PC, IBM-AT or compatible computer system. Each word
consists of 16 bits on the IBM-PC and compatible computers as
-------
compared to the UNIVAC where each word consists of 36 bits.
1.3 Running SHORTZ to Create an Output Tape for POSTZ
The instructions for running the SHORTZ Model are well
described in the User's Instructions. To write an output
tape, a value of 2 or 3 must be given for the parameter KSW
(column 1, Card 1). For the MS-DOS version, only a value of 2
can be used, since there is no capability for tape input with
the MS-DOS version of the program.
For the UNIVAC version, there are additions required in the
Executive Control Language (ECL) statements which appear at
the start of the program. Unless directed otherwise, the
program writes the tape output to unit 3. The following
statements must appear in the ECL:
@ASG,options 'output-tape-file.,type,reel-number
@USE 3,output-tape-file
The options are described in the User's Instructions for the
SHORTZ program. They concern the nature of the tape, whether
it is a temporary tape,r a high d-ensity tape, an unlabeled
tape, etc. Since writing will occur, it is always necessary
to include a /W in the option list. .Perhaps the most common.
use is an internal tape with a reel number assigned by the
NCC. The options should then be F/W. The F option allows the
tape to be written even if a label has been previously written
to the tape, but the user does not know the label, or cannot
remember the label. It requires only that the user remember
the reel-number, and of course that the user's identification
number be authorized to write to the tape. The output tape
file is a name given by the user for the output file to be
written. An example would be SHZOUT.RUN1. The type specifies
the density of the tape. Most tapes are 1600 or 6250 bits per
inch and the specification here should be T. The reel-number
identifies the tape within the NCC system and is assigned by
the tape librarian when the user calls and requests a tape.
The value of 3 in the USE statement above indicates that the
output is to be written to unit 3, the default. It should be
noted that the user should change this if internal
specifications in the input stream for the model have changed
the output unit number.
The output files for the SHORTZ Model, depending on the length
of the meteorological data base, can be quite long. It is
possible to write a tape file that is longer than will fit on
a typical tape. The SHORTZ Model has the capability to write
more than one tape, but this capability has not been put into
the POSTZ program. The user should ensure that the SHORTZ run
is restricted enough so that more than one tape is not
-------
written. The SHORTZ User's Instructions provide discussions
concerning the output length for tapes. This will typically
not be a problem for most runs of either program. Given the
possibility for long output, a tape should be dedicated to the
SHORTZ output. However, it is possible to write the output to
a single file within a tape, using the MOVE statement. The
format for this statement is:
@MOVE output-tape-file.,n
In this above statement, the value of n is the number of end
of file marks to skip on the tape to which the output is to be
written. A similar statement is needed in the ECL for the
POSTZ model which will read the tape (discussed later). For
more information on running SHORTZ on the UNIVAC computer,
consult the User's Instructions.
The MS-DOS version of the SHORTZ Model is much simpler to run.
An executable version of the program is generally obtained on
a floppy diskette, and called SHORTZ.EXE. The program can be
executed by typing d: SHORTZ where d is the disk drive where
the SHORTZ model is currently located. The user should be
logged into the drive where the input files are located. The
program will prompt for the name of the card input file name,
which is, an ASCII -disk file in exactly the same format as the
card-image input file for the UNIVAC version. It is suggested
that the user create the input file with a screen or line
editor program. The SHORTZ program will also prompt for an
output file name where the printed output is to be written.
Finally, the SHORTZ program will prompt for the tape output
file name. After completion of the SHORTZ run, the output can
be copied to the line printer by the command COPY
outputfilename LPT1:. Do not attempt this with the tape
output file since the tape output is not in ASCII, but rather
in unformatted form. No job control statements are required
in the input stream for the MS-DOS version of the SHORTZ
program.
Whether running on the UNIVAC, or on a micro-computer, once
the SHORTZ output tape (or file) is created, the user is ready
to run POSTZ. The following chapter discusses the preparation
of inputs for the POSTZ program.
-------
2.0 POSTZ USER INSTRUCTIONS
2.1 Program Options
The POSTZ computer program is designed to analyze an output
tape from the SHORTZ air quality model. The SHORTZ program
should be run first. Some guidance is given in Chapter 1.0 on
the running of SHORTZ, but the user is encouraged to refer to
the SHORTZ User's Instructions ("User's Instructions for the
SHORTZ and LONGZ Computer Programs, Volumes I and II", EPA
Document Numbers EPA-903/9-82-004a and b).
2.1.1 Output Tables
The POSTZ program's major function is to produce a series of
output tables. There are three separate output tables that
can be produced by the POSTZ program. Each of these tables is
described as follows:
o High-5 Table. The POSTZ program produces a
.table of the highest five concentration
estimates at each receptor. For example, if the
model has been instructed to compute 3-hour
average concentrations (see Section 2.1.2 for a
discussion of averaging time options) for a full
year of hourly concentration data produced by
SHORTZ, the POSTZ program will use the hourly
concentrations produced by SHORTZ to compute 3-
hour averages for every 3-hour period in the
year (see Section 2.1.3 concerning option of
running versus block averages). The POSTZ
program will also sort all these 3-hour averages
to determine the five highest values at each
receptor location. When printing the results,
the model also prints the ending hour and day
for each concentration printed.
o Top-50 Table. The POSTZ program sorts all the
computed average concentrations (for a
particular averaging time) to determine the
highest 50 values regardless of receptor. The
program prints the ending hour, day and the
receptor coordinates for each concentration
printed.
o Exceedance Table. The POSTZ program sorts and
records all values above a user specified
criteria value. For example, the user may
-------
desire to print out all 24-hour average
concentrations above 37 micrograms per cubic
meter (24-hour PSD Increment for TSP). A
maximum of 100 values may be printed in this
table. The values are sorted in order of
concentration, very similar to the Top-50 table.
The POSTZ program also prints the ending hour,
day and receptor coordinates for each
concentration printed.
2.1.2 Averaging Times
The user can request the POSTZ program to print any one, any
two or all three of the above tables for each of 6 separate
averaging times. The averaging times are:
"o one-hour average concentrations,
o _ three-hour average concentrations,
o eight-hour average concentrations,
o . '24-hour average concentrations,
o annual-average concentrations (in reality, this
option is for any long-term average, since the
number of days per year is an input specified by
the user, thus if 30 day averages are required,
the user can specify that there are 30 days per
year), and
o an additional averaging time to be specified by
the user. The length of this averaging time may
be any value from one hour to 24 hours.
The switches at the start of the program enable or disable any
of the above tables.
2.1.3 Calms Policy
The POSTZ program can be requested to implement the EPA calms
policy when computing averages. If the calms policy is
requested, when values of wind speed are less than 1.0 meters
per second, the program considers the concentrations as
invalid and averages on the basis of the other concentrations
within the time period. A full discussion of the calms policy
as it is implemented in the POSTZ program is contained in the
Guideline on Air Quality Models, (EPA-450/2-78-027R) . It
should be noted that the use of the calms policy is an option
of the POSTZ program, and is entirely at the discretion of the
user.
-------
2.1.4 Running or Block Averages
The user can specify whether running or block averages are to
be computed. Running averages mean that the average is
computed every hour for the previous n hours where n is the
averaging time selected. For example, if there are 24 hours
in the data base and running 3-hour average concentrations are
to be computed, the first period will cover hours 1 through 3,
while the second period will cover hours 2 through 4. Twenty-
two different but overlapping 3-hour averages will be
calculated from the 24-hour data base.
Block averages mean that the individual averaging periods do
not overlap. For example, in the above case if block averages
were to be computed rather than running averages, the first
period would include hours 1 through 3 while the second period
would include hours 4 through 6. Eight non-overlapping 3-hour
averages will be calculated from the 24-hour data base.
The POSTZ program requires that all averaging times are
controlled by the same switch for running or block averages,
so that it is not possible in the same run to compute running
averages for 3-hour concentrations and block averages for 24-
hour concentrations. The one exception to the above rule is
for annual averages. Annual averages (or any averaging period
chosen to be processed as an annual average) must be computed
with block averages. Even if running averages are used for
the remainder of the averaging times, block averages will be
automatically used for the annual averages.
2.1.5 Restricting Time Limits
The POSTZ program allows the user to select a specific portion
of the sequential data tape for analysis. Thus, if a SHORTZ
run was made for a full year and a sequential concentration
tape subjected to POSTZ analysis, the user may restrict the
POSTZ analysis and all the tables produced in the analysis to
a smaller period of time than the full year. For example, a
POSTZ analysis might be run for just the summer months. The
user restricts the period of time by providing to the program
a start date and time for the analysis to begin, and ending
date and time for the program to stop.
2.1.6 Restrict Receptors
Similar to the above analysis, the user may restrict the POSTZ
analysis to consider only a portion of the receptor network
for the SHORTZ run. Both grid and discrete receptors may be
selected for analysis. The default condition for the program
is to assume all receptors are to be included. An initial
switch tells the program that receptor restriction will be
used, and the program then assumes all receptors are to be
-------
excluded. One by one, then, the user instructs the program
which receptors are to be included.
2.1.7 Restrict Sources
Similar to the receptor restriction, it is possible for the
program to exclude the concentrations from certain sources.
The program initially begins by assuming all sources are to be
included. Once the user tells the program (through a switch)
that source restriction is to be performed, the program then
assumes that all sources are to be excluded, until .the user
selects, one by one, the sources to be included in the
analysis.
2.1.8 Specify Background
The user can add a background concentration to each receptor.
The program allows a "blanket" value of the background to be
assigned to all receptors, ,as well as a specific format for
separate background concentrations for each receptor.
2.1.9 Scale Concentrations
The concentrations resulting from specific sources can be
scaled. The individual source is identified ,and a scaling
factor is used for the contribution to the concentrations
computed by SHORTZ for that source.
2.2 Model Parameters and Definitions
In this section the major model parameters are discussed.
Specifically, the model parameters are defined, and the format
for the input file is presented in the next section.
Following are descriptions of the input parameters:
PTITLE This is a run title for the POSTZ run. It is
not the same as the run title for the SHORTZ
run. POSTZ prints both the run title for the
SHORTZ run and PTITLE on the first page of the
POSTZ output. The title may be used for record
keeping or to label the output.
ISW A total of 25 separate switches are included in
the array ISW(25). They assume values of 0, 1
or 2 depending on the switch and desires of the
user. In general a value of 0 is used to
decline a particular option (leaving the space
blank also declines the option). Placing a 1 or
2 in the switch selects the option and often
requires subsequent cards containing more
information which would be left out if the
switch were not selected. The first 5 switches
8
-------
control the various program options. ISW(l)
specifies whether t,he time period for the
analysis is to be restricted to particular hours
or days. If this option is selected, the user
is required to provide information later in the
runstream which tells the program how to
restrict the time period. ISW(2) determines if
the analysis is to be restricted to certain
receptors. If the switch is selected, the user
is required to provide information on which
receptors to restrict the analysis to on a later
card. ISW(3) determines if the analysis is to
be restricted to certain sources. If selected
the user must provide a list of the sources to
include later in the runstream. ISW(4)
specifies that background concentrations are to
be added to fe'ceptors. Depending on whether the
value for the switch is 1 or 2 the user may add
a uniform background, to all receptors, or he may
specify a specific background for each receptor.
ISW(5) determines if concentrations from certain
sources are to be scaled. If selected the user
must provide information later in the runstream.
Switches 6 through 23 select the various output
tables to be prepared. In general, the switches
are grouped in threes for each averaging time.
Thus switches 6-8 refer to 1-hour averages,
switches 9-11 refer to 3-hour averages, switches
12-14 refer to 8-hour averages, switches 15-17
refer to 24-hour averages, switches 18-20 refer
to annual averages, and switches 21-23 refer to
the user specified averaging time. Each group
of three switches determines whether a high-5,
top-50 or exceedance table are to be prepared.
If the any of the switches pertaining to the
user specified averaging time are selected
(switches 21-23) the user is required to provide
information on the averaging time desired later
in the runstream. ISW(.24) determines whether
running or block averages are to be computed for
all averaging times except annual (annual
averages are always block averages). If the
switch is 0, block averages are used, if the
value is 1 running averages are used. Finally,
ISW(25) determines whether the EPA calms policy
is to be implemented or not. If the switch is
0, the calms policy is not implemented, if the
value is 1 then the calms policy is implemented.
JULST One item of information not passed to the tape
from the SHORTZ program is the Julian day number
of the start of the meteorological record.
-------
IYR
IDYST
IYRST
IDYEND
IYREND
NXINCL
Therefore, the user is required on all runs to
specify the Julian day number for the start of
the sequential analysis. JULST
number (from 1 to 366) for the
meteorological data that
run to create the tape
current run of POSTZ.
is the day
start of the
was used in the SHORTZ
being analyzed by the
Like'JULST, IYR is the year for the start of the
meteorological data. The year numbers can be
provided using values four digits long, such as
1986 or it can be provided using values with two
digits, such as 86, as long as the user is
consistent in all date information, using the
same format for all dates provided to the POSTZ
program.
If the user elects to restrict the analysis to
certain days, he specifies the day the analysis
iS' to start on with IDYST. If IDYST is less
than JULST, IDYST is set equal to JULST (unless
IYRST is greater than IYR).
If the user elects to restrict the analysis to
certain days, he specifies the year the analysis
is to start in with IYRST. If IYRST is less
than IYR, IYRST is set equal to IYR.
If the user elects to restrict the analysis to
certain days, IDYEND is used to specify the last
day for the analysis. If IDYEND is selected so
that the ending day is beyond the end of the
sequential tape, IDYEND is reset to reflect the
end of the tape.
If the user elects to restrict the analysis to
certain days, IYREND is used to specify the year
in which the last day of the analysis occurs.
If IYREND is selected so that the ending day is
beyond the end of the sequential tape, IYREND is
reset to reflect the end of the tape.
If the
If the user elects to restrict the analysis to
certain receptors, the parameter NXINCL is used
to indicate the number of X grid points which
will be included in the analysis. The SHORTZ
program allows the specification of two types of
receptors: grid receptors and discrete
receptors. This parameter concerns only the
grid receptors and specifies the number of X
grid points which will be included in the
analysis. If there are only discrete receptors
10
-------
NYINCL
XX
YY
NXYINC
XDISCR
YDISCR
NSELIM
K
in this run, NXINCL should be set to 0.
If the user elects to restrict the analysis to
certain receptors, the parameter NYINCL is used
to indicate the number of Y grid points which
will be included in the analysis (see note under
NXINCL concerning type of receptors). If the
value of NXINCL is set to 0, there are no grid
points, and NYINCL is not read.
The X-grid values to include in the analysis if
the analysis is to be restricted to certain
receptors. It is important to note that these
X-grid values must have been in the original
SHORTZ run. NXINCL values of XX are read.
The Y-grid values to include in the analysis if
the analysis is to be restricted to certain
receptors. It is important to note that these
Y-grid values must have been in the original
SHORTZ run. NYINCL values of YY are read.
If the user elects to restrict the analysis to
certain receptors, the parameter NXYINC
indicates the number of discrete receptors to
include in the analysis (see note under NXINCL
concerning the type of receptors). If there are
no discrete receptors, this value should be set
to 0.
If the user elects to restrict the analysis to
certain receptors, the individual X coordinate
values for the discrete receptors to include in
the run are stored in XDISCR. It is important
to note that the discrete receptors specified
here must have been in the original SHORTZ run.
NXYINC values of xdiscr are read.
If the user elects to restrict the analysis to
certain receptors, the individual Y coordinate
values for the discrete receptors to include in
the run are stored in YDISCR. It is important
to note that the discrete receptors specified
here must have been in the original SHORTZ run.
NXYINC are read.
Contrary to the implied name of this parameter,
it is used to specify the number of sources to
include in the run if the user elects to
restrict the analysis to certain sources.
If the user elects to restrict the analysis to
11
-------
certain sources, the individual source
identification numbers as specified in the
original SHORTZ run are stored in K. NSELIM
values of K are read.
BACK
NUMDIS
XI
If the user elects to assign a uniform
background concentration to all the 'receptors,
the value is stored in BACK.
If the user elects to specify a different
background concentration for specific receptors,
the number of receptors to be specified in this
manner is entered in NUMDIS.
If the user elects to specify a different
background concentration for specific receptors,
the values ,are input into the program by first
identifying the receptor by coordinates, and
then specifying the value to use for background
at that receptor. The receptor coordinates must
have been in the original SHORTZ run. XI is
used to specify the X-coordinate value for
background identification. NUMDIS values of XI
are read.
Yl
VAL
NSCALE
ISCALE
If the user elects to specify a different
background for specific 'receptors, (see note
under XI above) Yl is used to specify the Y-
coordinate value for background identification.
NUMDIS values of Yl are read.
If the user elects to specify a different
background concentration for specific receptors,
the value to be used at a specific receptor
whose coordinates have been identified by XI and
Yl is stored in VAL.
If the user elects to scale the concentrations
for a specific source, the number of sources to
scale in this manner is specified by NSCALE.
If the user elects to scale the concentrations
for a particular source, the individual source
identification number as stated in the original
SHORTZ run is specified in ISCALE. NSCALE
values of ISCALE are read.
SCALE
If the user elects to scale the concentrations
for a particular source, the scaling value is
specified by SCALE. To double the
concentrations from a particular source, the
value of SCALE would be set to 2. NSCALE values
12
-------
of SCALE are read.
NDYYR
IOTHER
ONEEX
THREX
EIGEX
TWYEX
ANNEX
OTHEX
If annual concentrations are to be computed (as
determined by the setting to a value of 1 for
any of the switches ISW(18)
ISW(20)), the number of days to
one year are specified by NDYYR.
noted that if the user desires
average of less than one year, but
ISW(19) or
be considered
It should be
to compute an
greater than
one day (e.g. 30-day average), the NDYYR value
may be set to any number desired. There is no
requirement that the "annual" concentrations
actually be one year averages. There is a
requirement, however, that any averages computed
using the annual options must be block averages.
If, the user desires to specify an averaging time
other than the ones established, the value of
IOTHER allow the specification of the number of
hours for a different average. The restrictions
are that the value must be at least one hour,
and that the value cannot be greater than 24
hours.
I.f the user elects to print out exceedance
tables, the criteria for the exceedance of one-
hour average concentrations is specified in
ONEEX.
If the user elects to print out exceedance
tables, the criteria for the exceedance of
three-hour average concentrations is specified
in THREX.
If the user elects to print out exceedance
tables, the criteria for the exceedance of
eight-hour average concentrations is specified
in EIGEX.
If the user elects to print out exceedance
tables, the criteria for the exceedance of 24-
hour average concentrations is specified in
TWYEX.
If the
tables
annual
ANNEX.
user elects to print out exceedance
the criteria for the exceedance of
average concentrations is specified in
If the user elects to print out exceedance
tables, the criteria for the exceedance for the
user specified averaging time period is
13
-------
specified in OTHEX.
2.3 Format for the Input File
The POSTZ program requires two types of input. The first is
the SHORTZ output tape, which has been discussed in Section
1.0. The second type of input 'is the specification of all
the parameters listed in Section 2.2. The format for input of
these parameters is an ASCII "card-image" file.. In actuality,
very little processing is currently done with cards. However,
the format for the input to all computer models is usually
based on the concept of 80-column records, simulating card
input. The records are typically stored in a disk file. For
the UNIVAC version of the model, the card-image input file may
be stored in a different file than the Executive Control
Language (ECL) command file which controls the POSTZ run, or
they may be placed between the ECL statements at the proper
location. In Section 2.4, the format for the ECL to use with
the POSTZ Model is discussed. Here the format for the card-
image input file is presented in Table 2-1 . The MS-DOS
version of the computer program is much simpler to use. The
card-image input file is located in a disk file with no
control statements. When the program is initiated by the
command POSTZ, the program prompts the user for the name of
the input file. Other prompts are for the name of the SHORTZ
output file, and for the name of a file to store the output
from the POSTZ run.
2.4 Executive Control Language (ECL) for the POSTZ Program
To run the POSTZ Model on the UNIVAC computer it is necessary
to create a batch input command file, instructing the UNIVAC
as to the location of the input files, the location of the
program to be executed (POSTZ) and the destination for any
output. The assumption has been made here that the user is
familiar with the UNIVAC operating system, and can create the
input streams with an editor program, or alternatively to
create them off-line on a micro-computer and "upload" the
input stream to the UNIVAC. The command file must contain
certain Executive Control language (ECL) statements for the
POSTZ program to work. The statements and the order they
should appear are discussed as follows:
@RUN,priority jobid/95,acctnum/userid, projectid/uid,time
where priority is the run priority, jobid is a six character
name for the run, acctnum is the user's account number, userid
is the user identification code, projectid is the project
identification code, uid is the ADP identifier and time is the
maximum runtime in minutes.
14
-------
Table 2-1
User's Instructions for Creating the Card Image
Input File
Card Card
Group No. Parameter Format Description
1 1 PTITLE A80 Run title for the POSTZ run,
different from the run title
for the SHORTZ•run used to
create the output tape
2 1 ISW(l) 2512 Switch to determine if time
(col 2) period of analysis is to be
restricted. If =0, use
entire time period on the
output tape from the SHORTZ
run. If =1, include Card
Group 4 to specify time
period restriction
ISW(2) (col 4) Switch to determine if
analysis is to be restricted
to .certain receptors. If
=0, all receptors will be
included in the POSTZ run.
If =1, restrict analysis
only to those receptors
specified in Card Group 5.
ISW(3) (col 6) Switch to determine if the
analysis is to be restricted
to only certain sources. If
=0, all sources will be
included. If =1, restrict
analysis to only those
sources specified in Card
Group 6.
ISW(4) (col 8) Switch to determine if
background concentrations
are to be added to the
SHORTZ Model results. If
= 0, no background
concentrations are included;
if =1, uniform background is
to be added to all receptors
as specified in Card Group
7; if -2, separate
15
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card Card
Group No. Parameter Format Description
background is to,be.added to
each receptor as specified
in Card Group 7.
ISW(,5) (co.l 10) Switch to determine if
contributions from
particular sources are to be
scaled. If =0, no source
scaling is done; if =1,
scaling will be performed
for sources and values
specified in .Card Group 8.
ISW(6) (col 12) Switch to determine if a.
• high-5 output table is to be
prepared for 1-hour average
concentrations. If =0, no
, ' high-5 table .for 1-hour
concentrations.. If '= 1 ,
prepare a high-5 table by
receptor for 1 hour average
concentrations.
ISW(7) (col 14) Switch to determine if a
Top-50 table is to be
prepared for 1-hour average
concentrations. If =0, no
top-50 table for 1-hoilr
concentrations. If =1,
prepare a top-50 table for
1-hour average
concentrations.
ISW(8) (col 16) Switch to determine if an
exceedance table is to be
prepared for 1-hour average
concentrations. If =0, no
exceedance table for 1-hour
concentrations. If =1,
prepare an exceedance table
for 1-hour average
concentrations. Exceedance
16
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card
No. Parameter Fo_r_mat DesjcrijDt_ic)n
value is specified in Card
Group 11 (ONEEX).
ISW(9) (col 18) Switch to determine if a
high-5 output table is to be
prepared for 3-hour average
concentrations. If =0, no
high-5 table for 3-hour
concentrations. If =1,
prepare a high-5 table by
receptor for 3 hour 'average
concentrations.
ISW(IO) (col 20) Switch to determine if a
Top-50 table is to be
prepared for 3-hour average
concentrations. If =0, no
. ' top-50 table for 3-hour
concentrations. If =1,
prepare a top-50 table for
3-hour average
concentrations.
ISW(ll) (col 22) Switch to determine if an
exceedance table is to be
prepared for 3-hour average
concentrations. If =0, no
exceedance table'for 3-hour
concentrations. If =1,
prepare an exceedance table
for 3-hour average
concentrations. Exceedance
value is specified in Card
Group 11 (THREX).
ISW(12) (col 24) Switch to determine if a
high-5 output table is to be
prepared for 8-hour average
concentrations. If =0, no
high-5 table for 8-hour
concentrations. If = 1 ,
prepare a high-5 table by
17
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card Card
Group No. Parameter Fo_r_majt Description
receptor for 8 hour average
concentrations.
ISW(13) (col 26) Switch to determine if a
Top-50 table is to be
prepared for 8-hour average
concentrations. If =0, no
top-50 table for 8-hour
concentrations. If =1,
prepare a top-50 table for
8-hour average
concentrations.
ISW(14) (col 28) Switch to determine if an
exceedance table is to be
prepared for 8-hour average
concentrations. If =0, no
exceedance table for 8-hour
concentrations. If =1,
prepare an exceedance table
for 8-hour average
concentrations. Exceedance
value specified in Card
Group 11 (EIGEX).
ISW(15) (col 30) Switch to determine if a
high-5 output table is to be
prepared for 24-hour average
concentrations. If =0, no
high-5 table for 24-hour
concentrations. If =1,
prepare a high-5 table by
receptor for 24 hour average
concentrations.
ISW(16) (col 32) Switch to determine if a
Top-50 table is to be
prepared for 24-hour average
concentrations. If =0, no
top-50 table for 24-hour
concentrations. If =1,
prepare a top-50 table for
18
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card Card
No. Parameter Fp_rma_t Description
24-hour average
concentrations.
ISW(17) (col 34) Switch to determine if an
exceedance table is to be
prepared for 24-hour average
concentrations. If =0, no
exceedance table for 24-hour
concentrations. If =1,
prepare an exceedance table
for 24-hour average
concentrations. Exceedance
value specified in Card
Group 11 (TWYEX).
ISW(18) (col 36) Switch to determine if a
high-5 output table is to be
prepared for annual average
concentrations. If =0, no
high-5 table for annual
concentrations. If =1,
prepare a high-5 table by
receptor for annual average
concentrations. Must input
the number of days per year
in Card Group 9.
ISW(19) (col 38) Switch to determine if a
Top-50 table is to be
prepared for annual average
concentrations. If =0, no
top-50 table for annual
concentrations. If =1,
prepare a top-50 table for
annual average
concentrations.
ISW(20) (col 40) Switch to determine if an
exceedance table is to be
prepared for annual average
concentrations. If =0, no
exceedance table for annual
19
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card
No. Parameter Format Description
concentrations.. If =1,
prepare an exceedance table
for annual average
concentrations. Exceadance
value specified in Card
Group 11 (ANNEX).
ISW(21) (col 42) Switch to determine if a
high-5 output table is to be
prepared for user specified
average concentrations. If
=0; no high-5 table for user
specified concentrations.
If =1, prepare a 'high-5
table by receptor for user
specified average
concentrations.
ISW(22) ' (col 44) Switch to determine if a'
Top-50 table is to be
prepared for user specified
average concentrations. If
.=0, no top-50 table for user
specified concentrations.
If =1, prepare a top-50
table for user specified
average concentrations.
ISW(23) (col 46) Switch to determine if an
exceedance table is to be
prepared for user specified
average concentrations. If
=0, no exceedance table for
user specified
concentrations. If =1,
prepare an exceedance table
for user specified average
concentrations. Exceedance
value specified in Card
Group 11 (OTHEX).
ISW(24) (col 48) Switch to determine whether
20
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card Card
Group No.. Parameter Format
ISW(25) (col 50)
JULST
15
Description
block averages or running
averages are computed in all
analyses. If =0, then block
averages are- computed. If
=1, then running averages
are computed.
Switch to determine if EPA
calms policy is to be
implemented. If =0, the
calms policy' is not
implemented, if =1, the
calms policy is implemented.
Julian day for the start of
the data on the tape from
the SHORTZ run.
IYR
15
Year for the start of the
data on the tape from the
SHORTZ run. It can either
be of the form 1986 or of
the form 86 -- the program
will use either.
[This card is only read if ISW(1)=1]
IDYST
IYRST
IDYEND
IYREND
15
15
15
15
The Julian day for the start
of the analysis.
The year for the start of
the analysis.
The Julian Day for the end
of the analysis.
The year for the end of the
analysis.
[This card group is only read if ISW(2)=1]
1 NXINCL 15 The number of X grid points
to included in the analysis.
21
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card
No. Parameter Format Description
2 [This card is not read if NXINCL=0]
[This card can be repeated as often as necessary
to include all the X grid points]
XX 8F10.0 The X grid points to include
in the analysis. Note, each
X value to be used must be'
specified here. It is not
possible to specify the
starting x value and the
increment for other values
as it is in SHORTZ.
3 [This card is. not read if NXINCL=0]
NYINCL 8F10.0 The number of Y grid points
to include in the analysis.
4 [This card is not read if NXINCL=0]
[This card can be repeated as often as necessary
to include all the Y grid points]
YY 8F10.0 The Y grid points to be
included in the analysis.
(See note on format for X
grid points above).
5 NXYINC 15 The number of discrete
receptors to include in the
analysis.
6 [This card is not read if NXYINC=0]
[This card can be repeated as often as necessary
to include all the discrete receptors]
XDISCR 8F10.0 The x values for the
discrete receptors.
7 [This card is not read if NXYINC=0]
[This card can be repeated as often as necessary
to include all the discrete receptors]
YDISCR 8F10.0 The y values for the
discrete receptors.
6 [This card group is only read if ISW(3)=1]
1 NSELIM 15 The number of sources to
22
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card
No. Parameter Format Description
include in the run
2 [This card is repeated NSELIM times]
K 15 The source identification
number of the source to be
included in the run.
[This card group is not read if ISW(4)=0]
If ..ISW(_4)=1 the following format applies
1 BACK F10.0 The uniform background
concentration to apply to
all receptors.
IfISWJ4)=2 the following format applies
1 NUMDIS 15 The number of receptors for
which a background value is
to be specified
2 [This card is repeated NUMDIS times]
XI F10.0 The x-coordinate of the
receptor for background
specification
Yl F10.0 The y-coordinate of the
receptor for background
specification
VAL F10.0 The background concentration
to use for the identified
receptor.
[This card group is not read if ISW(5)=0]
1 NSCALE 15 The number of sources to
scale
2 [This card is repeated NSCALE times]
ISCALE 15 The source identification
number of the source to be
scaled.
SCALE F10.0 The scaling value •
23
-------
Table 2-1 (continued)
User's Instructions for Creating the Card Image
Input File
Card
No. Parameter Format Description
[This card
ISW(20) =1]
group is only read if ISW(18), ISW(19). or
1 NDYYR • 15 The number of days per year,.
(Note this value can be any
desired number less than
366. For example, if the
user desired a 30 day
average value, the annual
can be set to '30 days with
this card.)
10 [This card group is only read if ISW(21), ISW(22) or
ISW(23) =1]
1 IOTHER 15 The averaging time to be
specified by the user.
(Number of hours)
11
1 [This card is not read if ISW(8)=0]
ONEEX F10.0 The exceedance value for 1-
hour concentrations
[This
THREX
card is not
F10.0
read if ISW(11)=0]
The exceedance value for
hour concentrations
3-
[This card is not read if ISW(14)=0]
EIGEX F10.0 The exceedance value for 8-
hour concentrations
[This card is not read if ISW(17)=0]
TWYEX F10.0 The exceedance value for 24-
hour concentrations
[This card is not
ANNEX F10.0
[This
OTHEX
card is not
F10.0
read if ISW(20)=0]
The exceedance value for
annual concentrations
read if ISW(23)=0]
The exceedance value for
user specified average
concentrations
24
-------
@CAT,P print-file(+l) .
where print-file is a name specified for the output
destination. This is an optional card which allows the output
from the POSTZ run to be stored in a file on the disk for
later access with an editor program. This structure is ideal
for remote batch jobs.
@BRKPT PRINTS,print-file
This card is only included if the @CAT card is included. It
redirects the output for the run to the file named in the foiCAT
command.
@ASG/options output-tape-file.,type,reel-number
@USE 2,output-tape-file
These two cards specify the SHORTZ output tape. .The parameter
output-tape-file is a 'name used for the' purposes of this run
to identify the file on the tape used fpr input to the POSTZ
run. The values of type and reel-number are the same as for
the SHORTZ run which created the tape (,see Chapter 1.0). As
with the SHORTZ run, it is possible to use a specific file on
a tape, rather than the entire tape as the above command would
imply. If a specific file on the tape is to be read, the
additional command shown below is needed:
@MOVE output-tape-file.,n
where n is the number of end-of-file marks to skip on the tape
of the output file from the SHORTZ run.
@ASG,A program-file.
This card locates the program file and assigns it to the run.
The actual name for the program file will be determined by the
user and the system he is working on. • To inquire regarding
the means.to obtain a copy of POSTZ or the proper program file
name if using the National Computer Center, the user should
contact the nearest regional office of EPA.
@ASG,A card-input-file.
This card is optional, depending on whether the card-image
input file is included in the command runstream, or in a
separate file. If in a separate file, the above command must
be included, and card-input-file is the specification of the
filename for the card-image input file.
@XQT program-file.ABS
This card actually begins execution of the program. The AES
25
-------
element for the file containing the program is assumed to be
the relocatable (executable) program produced by the FORTRAN
compiler. The user may have this file named differently than
is assumed above. If the input • data is located in a file,
then the following statement is required:
@ADD card-input-file.
If the cards are to be included in the command runstream, the
previous statement must not be included' and the cards
themselves occur here in the runstream in the order and format
indicated in Table 2-1.
®BRKPT PRINTS
@FREE print-file.
These cards are only included' if the output has been
redirected to a print file defined by the previous @CAT
statement.
@FIN
This card terminates the run.
2.5 M-S DOS Operation of POST2
The MS-DOS version of the program is generally obtained on a
floppy diskette. The program uses over 330 kilobytes of
memory and is called POSTZ.EXE. Execution of the MS-D'OS
version of the program is accomplished by typing:
d:POSTZ
where d is the drive designator for the location of the
program. The program will prompt the user for the name of the
card-image input file (Table 2-1) and the SHORTZ output taps
file. It will also ask for a new disk file name to be given to
the output from the POSTZ program. The POSTZ program informs
the user on the progress of the run by advising on the day and
hour of data currently being processed. At the completion of
the POSTZ run, the user can transfer the output to the line
printer by typing COPY output-filename LPT1:
26
-------
3.0 SAMPLE INPUT AND OUTPUT FILES
Sample input and output files have been taken from the User's
Instructions for the SHORTZ Model. Essentially the same input
file used in the test case for the SHORTZ Model has been used
here, with two exceptions: the output tables from the test
SHORTZ run are quite voluminous and not germane to the POSTZ
development, so. the control options have been shut off for the
output, and the output tape selection has been made by setting
KSW to 2 in Column 1 of Card 1 . The input file from the
SHORTZ run is included as Figure 3-1. The output file from
the SHORTZ run has been included as Figure 3-2.
The POSTZ test run uses the output tape written by the SHORTZ
program for the test case above as an input. The card-image
input stream is depicted in Figure 3-3. The output file from
the POSTZ run is depicted in .Figure 3-4. These Figures 'Were
made from runs of the POSTZ Model on a COMPAQ DESKPRO
computer, using MS-DOS Version 2.1. The Program was also
tested on the NCC UNIVAC with this same test case and produced
identical output values. The output results also match the
values given in the User's Instructions for the SHORTZ Model.
27
-------
Figure 3-1 Sample SHORTZ Input File
-------
2EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ
17 15 1
0 24
602000
607000
612000
605167
4487000
4491500
4489107
358
387
375
230
366
320
250
274
381
241
274
366
288
250
366
338
280
360
354
320
296
335
338
329
326
338
335
311
280
363
305
300
396
3 24
0
603000
607500
4488000
4492000
347
305
230
317.
259
250
308
250
296
241
308
274
280
290
320
323
317
369
344
399
244
274
5
311
604000 604500
608000 608500
4488500
4493000
317
280
226
296
230
241
253
280
290
329
232
265
226
323
226
320
229
363
226
349
229
290
4489000
4494000
341
381
229
247
226
305
226
363
226
293
229
311
232
315
235
379
238
360
244
305
250
271
605000
609000
4489500
4495000
375
300
274
280
262
286
238
329
232
347
229
355
229
354
258
372
268
354
286
283
305
290
605500
609500
4490000
4496000
396
256
349
366
282
326
268
293
229
354
290
354
282
311
244
372
235
290
296
283
274
286
606000
610000
4490500
4497000
332
286
360
354
320
372
314
347
238
360
323
347
296
314
293
277
286
320
229
293
268
317
606500
611000
4491000
351
305
305
373
370
335
262
314
244
372
280
335
335
290
326
344
305
369
299
360
296
320
-------
286
323
326
317
, 232
323
349
256
344
329
280
344
274.
.683
1 1
116
116
117
118
0100
0200
0300
0400
0500
0600
0700
0800
0900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
0000
1
117
170
190
21,0
220
245
255
255
250
250
250
250
250
250
250
255
260
260
265
260
250
250
240
260
270
318
318
1260
5.
6.
10.
9.
8.
9:
9.
10.
9.
8.
8.
7.
9.
7.
6.
6.
7.
6.
7.
6.
5.
6.
6.
5.
349 232
314 323
363 364
286 235
317 323
232 232
347 320
302 280
366.7
1 1
118
.15
.15
.00
4017
6878
0316
7744
2311
2600
7744
2889
0028
4883
2311
2022
2600
7167
6878
1733
7167
6878
7167
6878
9161
1733
1733
9161
-117 -118
6043804488740
6043804488740
6023304487800
953
1068
1184
1299
1415
1530
1645
1598,
1551
1504
1457
1410
1363
1316
1269
1221
1174
1127
1080
1033
986
939
892
845
226
274
265
293
308
229
235
229
67.1
62.5
229.0
283
284
285
285
283
282
280
280
280
279
279
279
279
278
278
277
276
276
275
275
275
275
274
274
229
229
229
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
.17
285
290
238
296
229
268
258
226
472
444
411
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
293
288
232
31L
241
323
311
254
160.98
162.14
881.46
.1051
.1051
.1051
.1051
.1051
, .1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
299 329
296 308
226 226
274 290
280 267
308 338
286 280
314 335
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
.0735
2.60
1.85
3.20
.1051
. .1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
.1051
-------
Figure 3-2 Sample SHORTZ Ouput File
-------
SHORT! (VERSION 82326)
AN AIR QUALITY DISPERSION MODEL IN
SECTION 2. NON-GUIDELINE MODELS,
IN UNAHAP (VERSION 5) DEC 82.
SOURCE: FILE 23 ON UNAMAP MAGNETIC TAPE FROM NTIS.
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORT!
DATE
, CASE 1, PA6E
TABLE
- GENERAL INPUT DATA -
NUMBER OF INPUT SOURCES
NUMBER OF X GRID COORDINATES
NUMBER OF Y GRID COORDINATES
TOTAL NUMBER OF HOURS IN EACH DAY
NUMBER OF DAYS OR CASES
NUMBER OF CONCENTRATION REPORTS (SOURCE COMBINATIONS)
NUMBER OF DISCRETE CALCULATION POINTS
MET DATA INPUT CARD RATE (1=HOURLY,2=2 HOURLY,ETC)
IS CONCENTRATION CALCULATED AT BASE RATE PRINTED
NO. OF HOURS IN FIRST AVERAGE CONCENTRATION PRINTED
NO. OF HOURS IN SECOND AVERAGE CONCENTRATION PRINTED
NO. OF HOURS IN THIRD AVERAGE CONCENTRATION PRINTED
ARE TERRAIN ELEVATION HEIGHTS USED
IS WIND SPEED TERRAIN FOLLOWING
ARE CONCENTRATIONS AVERAGED OVER DAYS OR CASES
IS THE FORMAT FOR SOURCE DATA READ
IS COORDINATE SYSTEM CARTESIAN (=0) OR POLAR (=1)
ARE DISCRETE RECEPTORS CARTESIAN (=0) OR POLAR (=1)
ARE SOURCE COORDINATES CARTESIAN (=0) OR POLAR (=1)
SIGEPU SIGAPU FOR ALL SOURCES OPTION
RURAL/URBAN MODE OPTION (RURAL=0),(UR8AN=1)
MODEL UNITS CONVERSION FACTOR
ACCELERATION OF GRAVITY
HEIGHT'OF MEASUREMENT OF 'WIND SPEED, ETC
ENTRAPMENT COEFFICIENT FOR UNSTABLE ATMOSPHERE
ENTRAPMENT COEFFICIENT FOR STABLE ATMOSPHERE
DISTANCE OVER WHICH RECTILINEAR PLUME EXPANSION OCCURS
DECAY COEFFICIENT FOR PHYSICAL OR CHEMICAL DEPLETION
ANGULAR DISPL OF GRID SYSTEM FROM TRUE NORTH
ELEVATION OF BASE OF WEATHER STATION
X ORIGIN OF POLAR COORDINATES
v ORIGIN OF POLAR COORDINATES
(NSOL'RC) =
(NXPNTS) ~-
(NYPNT3) =
(NHOURS) =
( NOAYS) =
(NGROUP) =
(NXWYPT) =
(!SW(1)) =
(ISW(2)) =
(ISW(3)) =
(ISW(4)) --
(ISW(5)) =
USW(7')) --
(!SW(9)) =
(ISW(11))=
(ISW(12))=
(ISW(13))=
(ISW(14))=
(ISW(15))=
(ISW(16))=
(!SW(17))=
( TK ) --
( 6 ) =
( ZR ) =
(GAMMA!) =
(GAMMA2) =
i (, XRY ) =
(DECAY ) =,
(ROTATE) =
(HA) --
(UTMX) =
(UTMY) =
3
17
IK
24
1
c
1
1
0
24
0
n
1
0
0
0
0
0
0
n
0
.10000000E+07
9.8000
5.0960
.600
.560
50.0000
.OOOOOOOOE+00
.883
366.70
.00
.00
(HETERS/SEC**2)
(METERS)
(METERS)
(METERS)
(METERS)
(METERS)
*-* COORDINATE SYSTEM X AXIS (METERS) *-*
.502000QOE+OS,
.59700000EiQ6,
.5'200000E+06,
.44S70000E+07,
J4915000E+07,
.5C300000E+06, .50400000E+06, .5045000CE+Q6, .50500000E-05.
.60750000E+06, .50800000E+06, .60850000E+06, .6Q900000E4-06,
.60950000EA05,
*-* COORDINATE SYSTEM Y AXIS (METERS )*-*
.44880000E+07, .44885000E+07, .44890000E+07,
.44920000E+07, .44930000E+07, .44940000E+07,
J48950QQE+07, .44900000E+07, J4905QOCE^-07,
.44950000E+07, .14960000E+07, .
-------
. EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ 9ATE . CASE 1, °AGE
TABLE 1 (CONT)
- GENERAL INPUT DATA - '
NUMBER OF SOURCES IN EACH CONCENTRATION REPORT GROUP 1 TO 5 - 1, 1. 1. 1, 1,
SOURCE NUMBERS FOR EACH OF THE ABOVE GROUPS = 116, 117, 118, -117, -113,
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ
TA8LE 2
- GENERAL 'INPUT DATA -
*-* GRID SYSTEM TCRRA!N HEIGHTS (METERS'* *-*
- GRID SYSTEM X AXIS (METERS) -
DATE
CASE 1 PAGE 3
"
AXIS (METERS
4497000.000
4496000.000
4495000.000
4494000.000
4493000.000
4492900.900
4491500.000
4491000.900
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487000.000
329.0000000
349.0000000
317.0000000
286.0000000
305.0000000
311.0000000
326.0000000
335.0000000
354.0000000
338.9000000
288.0000000
241.0000000
250.0000000
230.0000000
358.0000000
347.0000000
317.0000000
353.0000000
349.0000000
244.0000000
344.0000000
317.0000000
320.0000000
280.0000000
308.0000000
296.0000000
308.0000000
259.0000000
230.0000000
347.0000000
320.0000000
323.0000000
364.0000000
232.0000000
229.0000000
226.0000000
229.9000000
226.0000900
225.0000000
232.0000000
290.0000000
253.0000000
230.0000000
226.0000000
317.0000000
- HEIGHT
235.0000000
308.0000000
255.0000000
225.0000000
250.0000000
244.9000000
238.9000000
235.0000000
232.0000000
229.0000000
226.0000000
226.0000000
226.0000000
229.0000000
341.0000000
-
258,
,229,
'238,
285,
305.
286.
268
258,
229,
229.
232.
238.
262,
274,
375.
.0000000
.0000000
.0000000
,0000000
.0000000
.9000000
.0000000
.0000000
,0000000
.0000000
,0000000
,0000000
.0000000
.0000000
.0000000
311.0000000
241.0000000
232.0000000
293.0000000
274.0000000
296.0000000
235.0900000
244.0000000
282.0000000
290.0000000
229.0000000
258.0000000
282.0000000
349.0000000
396.0000000
285.0000000
280.0000000
226.0000000
299.0000009
258.9000900
229.9009000
285.0000090
293.0000000
296.0000000
323.0000000
238.0000900
314.0000000
320.0000000
360.0000000
332.0000000
280.
267.
226.
329.
296.
299
205.
325.
325.
230.
244.
262.
370.
3^5
351.
0000900
9000000
9009000
0090000
0000000
0090999
nnnnnnn
9000900
0000000
0000000
0000000
0900900
0000000
9000000
0000900
280.0000090
255.0900001
222.0000001
*3 o *5 onnono1"'
200. 0990091
TOO nnnnnol
'•-•••-•••-"-" I
no nnnnnn"
•w w' '. . J w *J „ J - .
,.,:<,. .".. Jj'.'w-
320.00CCCC']
280 0090991*
250.0000000
274.0000001
274.000000|
365.0090000
387.00QC001
- GRID SYSTEM X AXIS (METERS) -
AXIS (METERS
4497000.000,
4496000.000
flflQcnnn non
4494000.000
i493009.000
/M G o n n o nnn
- 4 _ i. _• u v . U 'J v.1
4491590.000
4491000 'JOO
•H9Q500.QOO
4490000.000
4489500.000
4489000.000
4488500.000
4438000.000
4487000.000
507500.000
)
3C2. 0000000
232.0000000
''96 0000000
314.0000000
274.9000000
399.9000900
359.0000000
323.00Q990Q
290.0090000
274.9000000
241.0000000
250.0000000
250.0000000
317.0000000
305.0000000
508000.000
280.0000000
232.0000000
235.9000000
323.0000000
290.0000000
349.0000000
353.0000000
320.0000009
323.0000000
265.0000000
329.0000000
280.0000000
241.0000000
296.0000000
280.0000000
508500.000
229.0000000
229.0000000
293.9000000
274.0000000
271.0000000
205.0000900
360.9090000
379.9000000
215.QOOOOOO
311.0000000
293.0000000
363.0000000
305.0000000
247.0000000
381.0000000
509000.900
- HEIGHT
226.0000000
258.0000000
296.0000000
290.0000000
299.0000000
283.0000000
351.9000000
372.0000000
35*1.0000000
355.0000000
347.0000000
329.0000000
286.0000000
280.0000000
300.0000000
509500.000
-
254,
323,
311.
288.
286
282
290
372
31!
354
354
293
326.
366.
256.
,0000000
.0000000
.0000000
.0000000
.9900900
.0000000
.0000000
.0000000
.0009000
.0000000
.0000000
.0000000
.0000000
.0000000
.0000000
610000.000
314.9000009
308.0000000
274.0COOOOQ
296.0000000
317.9900000
293.0000009
320.0000000
277.0000009
.314.0000000
347.0000000
350.0000000
347.0000000
372.0000000
354.0000000
286.0000000
611000.000
90s o on noon
338.0000000
290.0090000
308.0000000
320.9900090
250.0000009
369.9000000
344.0000900
290.Q9090CC
335.0000000
372.0000000
314.0000000
335.0000000
373.0000000
305.0000000
sunnn nnn
i. i t. J 0 •„ . -j j v
1/1/1
•^ 1 T .
344 .
•211
326.
395
263.
335.
32°
295.
350.
355.
356.
381.
320.
375.
9090000
nnfloonn
0900090
0000090
Q900090
n n g n n n A
9099000
.0000000
9COOOOO
0090900
OOOOOQC
0000009
0000000
0000000
0000000
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ DATE CASE 1, PAGE
TABLE 2 (CCNT)
- GENERAL INPUT DATA -
*-* DISCRETE POINT TERRAIN HEIGHTS (METERS) *-*
X Y HEIGHT X Y HEIGHT X Y HEIGHT
(METERS) (METERS ) (METERS) (METERS ) (METERS) (METERS )
505167.0 4489107.0 274.0000000
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORT!
DATE
CASE 1, 5AGE 5
TABLE
- SOURCE INPUT DATA -
- SOURCE INVENTORY -
C T SOURCE T SOURCE X Y HEIGHT IF TYPE=0 IF TYPE-9 ANGLE'
A A NUMBER Y STRENGTH COORDINATE COORDINATE ABOVE TEMP (OEG K) VOL. EMISS. T0
R ? P(GRAMS/SEC) (METERS) (METERS ) GROUND IF TYPE=10R2 RT. M**3/SEC LONG
D E
E
(METERS) LENGTH SHORT IF TYPE=10R2 SIDE
SIDE (MTRS) LENGTH LONG (DEG)
SIDE (MTRS)
STACK ELEVATION
INNER AT
RADIUS STACK
(METERS) BASE
(METERS)
A T c n T c T 9 r PI i f
SETTL!NG
VCI ATI TV
OF
(METERS/SEC) OCCURRENCE
(FRACnCN)
X
X
X
115 0
1170
118 0
318.150
318.150
1260.000
604380.00
504380.00
602330.00
4488740.00
4488740.00
4487800.00
67.10
52.50
229.00
472.000
444,
411
.000
.000
160.980
162.140
881.460
.0
.0
.0
2.600
1.850
3.200
229.00
229.00
229.00
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORT!
DATE
CASE
PAGE
TABLE
- METEOROLOGICAL INPUT DATA -
HOUR WIND
DIRECTION
(DEGREES)
THETA
100 170
200 190
300 210
400 220
500 245
600 255
700 255
800 250.
arm 2^n
1900 250,
1100 250.
1200 250
1300 250
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
OQQO
.0000
.0000
.0000
.0000
1400 250.0000 •
1500 255
1600 260
1700 250
1800 265
1900 260
2000 250,
2100 250,
2200 240,
2300 260,
0 270.
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
,0000
.0000
WIND LAYER
cpccn DFPTH
(MTR/SEC) (METERS)
UBAR HM
5.4017
6.6878
10.0316
9.7744
8.2311
9.2600
9.7744
10.2889
9.Q028
3.4883
8.2311
7.2022
9.2600
7.7167
6.6878
6.1733
7.7157
. 6.6878
7.7167
6.5878
5.9151
6.1733
5.1733
5 Q151
953.000
1068.000
1184.000
1299.000
1415.000
1530.000
1545.000
1598.000
1551 .000
1504.000
1457.000
1410.000
1363.000
1316.000
1269.000
1221.000
1174.000
1127.000
.1080.000
1033.000
985.000
939.000
892.000
845.000
AMBIENT
TEMP
(OEG K)
TA
283.000
284.000
285.000
285.000
283.000
282.000
280.000
280.000
280.000
279.000
279.000
279.000
279.000
278.000
278.000
277.000
275.000
276.000
275.000
275.000
275.000
275.000
274.000
274.000
VERT GRAD STAB WND SPD
OF POT TMP I1 ITY POWER LAW
(DEG K/M) CAT. EXPONENT
DPDZ ISTBLE P
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.0000
.1700
.1700
.1700
.1700
.1790
.1700
.1700
.1700
.1700
.1700
.1790
.1700
.1700
.1700
.1709
.1700
.1700
.1700
.1700
.1700 '
.1700
.1700
.1790
.1700
STD DEV EL
ANGLE, SOR
TYPE 0
SIGEPU(RAD)
.0735000
.0735000
.0735000
.0735000
.9735000
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000
' .0735000
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000'
.0735000
STD DEV AZ
ANGLE, SOR
TYPF 0
SIGAPU(RAD)
.1051000
.1051000
.1051000
.1051000
.1051000
.1051000
J051000
.1051000
.1051000
.1051000
.1051000
.1051DOO
.1051000
.1051000
.1051000
.1051000
.1051000 •
.1051000
.1051000 •
.1051000
.1051000
.1051000
.1051000
.1051000
CTO ncv Fl
ANGLE, SOR
TYPE 10R2
SIGEPL(RAO)
.0735000
.0735000
.9735000
.0735000
.0735000
.0735090
.0735000
.0735000
.0735000
.0735090
.0735000
.0735000
.0735000
.0735000
.0735000
.0735000
.'0735000
:0735009
.0735000
.073500C
.9735000
.9735909
.0735090
.0735000
crn new 47
ANGLE, SOR
TYPE ioR2
SIGAPL(RAO)
.1051000
.1051000
.1051000
.1051000
.1051000
.1051000
.1051000
.1051000
.1051000
. 1051009
.1051000
.1051900
.1051900
.1051009'
.1051000
.1051009
.1051000
.1051000
.!9S1000
.'95199C
.1951900
.195:000
.1951090
.1051000
1 ,|JTC&4>
OIFFL'S!CN
COEFFICIENT
ALPHA
.9090
.9000
.9000
.9000
.9000
.9000
.9009
Q o n n
Oiv-.n
on [if.
.9000
. .9009
.9000
.9009
.9900
.9009
.9009
on no
.9090
qniif
.9909
.90CO
.9900
onpj
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ
DAT
, CASE 1, PAGE
TABLE
5
24 HOUR GROUND LEVEL CONCENTRATION (M1CROGRAMS/CU8IC METER) FROM SOURCES 115
- HOUR(S) 0 TO 2300 -
502000.000
Y AXIS (METERS
- GRID SYSTEM X AXIS (METERS) -
(THE MAXIMUM CONCENTRATION IS 703.9803000 AT X= 605000.0, Y=4489000.0)
503000.000 504000.000 504500.000 505000.000 505500.000 505000.000
- CONCENTRATION -
606500.000 507009.009
4497000.000
4496000.000
4495000.000
4494000.000
4493000.000
4492000.000
4491500.000
4491000.990
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487000.000
AXIS (METERS
4497999.000
4496000.000
4495009.009
4494090.900
4493990.900
•U92999.900
4491599.900
4*191900.900
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487999.099
.5925970
.3250352
.1003984
.0107155
.9001398
.0000000
.0000000
.0000000
.0000000
.0000000
.0000000
.9000000
.0000000
.0000000
.0000000
507500 900
1.
i
2.
o
5044350
9542720
3557710
1938170
1.0225120
(THE
1291423
0104869
0001106
0000000
0000000
0000000
0000000
0000000
0000000
0000000
MAXIMUM
598000.000
2
4
6
9
10
4
.2215371
.3437609
.5860771
.9932383
.0061160
.5009410
.9260990
.5898230
.3085100
.0858920 '
.0795788
.0000000
.0000000
.0000000
.0000000
- GRID
CONCENTRATION
.2029481
.2930555
.4359421
.5869537
1.2503640
2.5044530
4.0643040
5.8793920
12.7106900
24.4354800
26.5601800
2.0146020
.0000000
.0000000
.0000000
.8092649
1.2350180
, 1.9995920
3.4535380
6.0568220
9.5372980
9.9648290
7.4840510
3.7723870
23.5652100
54.1037000
703.9803000
.0004179
.0000000
.0000000
1.6497090
2.2151770
2.8550010
3.3057270
2.6153250
1.1373510
2.2496440
9.0951679
18.4789900
23.2355400
31.3449190
490.3597000
1.1893320
.0900000
.0000000
1.5834970
1.5885710
1.2547629
.5900617
.6121413
4.5135740
8.0388390
19J509100
10.1667199
8.7312210
175.5237009
222.5909000
5.6338900
.0000022
.0000009
2.
s
P
5
2
.7528550 .19576^
.4777890 .1430093
0/1 1 T3/10 t 1 < 001 1
. L *» > J J S vj . H I 1 lj tl ^M
.5326683 1.24087«
.7725079 3.292!:^
o c -i 1 p c ft i o c (\ *j .1 o 0
03001(10 '; -50011! 1
.5114570 1 . 1*25'!
.4862359 11.5944300
35.0789500 71.55649(1
192.
95
7.
.4516990 123.25499'!
.9382000 57.1597500
.7552729 5.659445"
.9905453 .9129!5|
.9099999 099999*
SYSTEM X AXIS (METERS) -
IS 793.9393009 AT X=
508500.900 609000. 000'
699500.000
605000.0, Y=4489000.0)
619990.900
611099.900
512099.990
) - CONCENTRATION -
.1205282
.4177227
1.3120880
2.2535730
2.9139369
2.2827479
.7158939
4.4914030
22.4247200
37.2355190
77.3923500
34.9214000
5.9596820
.0545519
.0000009
1
1
2.
i
9.
36.
73.
55.
'24.
c
3690108
9779585
5332830
3939370
5960250
5253664
3991579
5342490
0857500
9470800
7816800
8581000
3525480
1355394
0090009
i
i
i
1
i
14
45
55
39
18
5
.7620249
.2349750
.3753450
.9455149
.2115910
.9919550
.7386799
.8854800
.3380200
.4491100
.3820900
.4492000
.1360030
.2301350
.0000092
.9542982
1.0731000
1.3519499
1.5292850
.3447785
2.4378300
7.4591309
21.3618900
44.9431500
40.5044500
29.2389300
13.7789500
4.4392390
.3511621
.0000041
.8712346
.9768595
1.3455440
.7422271
.2388307
4.1883640
19.3277990
25.7578700
33.2155300
30.7198900
21.9905100
10.5927400
3.9732790
.4552012
.0000330
.7558341
1.9377049
.9755375
.2529563
.7725705
5.7801799
14.00RQ°nn
28.3180990
30.4879400
24.1052800
16.8919300
8.5643200
3.5169310
.5462135
.0001559
.3099805
5559574
1 0 1 J 0 C p
TQQCTOQ
2.3539259
0 0 0 Q 0 q c q
10 iio()onn
23.1287190
19.2423300
15.8505300
10.5752900
5.7936580
2.7542090
.5486255
.0914541
1
2
". 'i
n
1 q
13
19
7
4
2
,» 7 Q c o f> c
.'539801
1 •} ,1 1 •} -3 o
.9955530
T 0 T c c .j (•;
'-' •) •) p ,1 n n
co^ncni
.9452599
.2881190
.9232909
.1162459
.2427380
.245C21C
.5844387
.0958393
-------
EXAMPLE SHORT TERM PROBLEM. - PROGRAM SHORT! DATE , CASE 1, PAGE 3
TABLE 5 (CONT)
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/CUBIC METER) FROM SOURCES 116
- HO'JR(S) 0 TO 230G -
- DISCRETE POINT RECEPTORS -
(THE MAXIMUM CONCENTRATION IS, 719.1874000 AT X= 605167.0, Y=44891Q7.0)
X Y CONCENTRATION X Y CONCENTRATION X Y CONCENTRATION
(METERS) (METERS ) (METERS) (METERS ) (METERS) (METERS )
605167.0 4489107.0 719.1874000
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORT!
DATE
TA8LS
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/CUBIC METER) FROM SOURCES
- HOUR(S)
0 TO 2300 -
602000.000
Y AXIS (METERS
- GRID SYSTEM X AXIS (METERS) -
(THE MAXIMUM CONCENTRATION IS 216.5484000 AT X= 505500.0, Y=4489000.0)
503000.000 504000.000 504500.000 505000.000 505500.000 506000.000
- CONCENTRATION '-
506500.000 507000.000
4497000.000
4496000.000
4495000.000
4494000.000
4493000.000
4492000.000
4491500.000
4491000.000
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487000.000
AXIS (METERS
4497000.000
4496000.000
4495000.000
4494000.000
4493000.000
'-492000.000
4491500.000
4491000.000
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4437000.000
.5840322
.3203222
.0988631
.0105380
.0001374
.0000000
.0000000
.0000000
.0000000
.0000000
.0000000
.9000000
.0000000
.0000000
.0000000
507500.000
/
. 1150303
.4096686
1.2995060
2.2527070
2.8433530
2.2813640
.5442634
3.2392370
18.4208900
74.0384400
59.5277000
25.0575900
4.3904240
.0421785
.0000000
1.4828680
1.9250540
2.3213540
2.1602910
1.0023560
.1268465
.0102633
.0001079
.0000000
.0000000
.0000000
.0000000
.0000000
.0000000
'.0000000
(THE MAXIMUM
508000.000
1.
4.
6,
q
9,
3
.2146315
.3346512
.5734432
.9705455
,9604090
.3821580
,7157480
,3189390
.7613190
,7404190
,0689247
.0000000
,0000000
.0000000
.0000000
- GRID
CONCENTRATION
.1615749
.2346455
.3430317
.5254368
.9395726
1.8483790
2.7222210
4.1925190
6.4423120
7.9690680
1.8308670 '
.0000014
.0000000
.0000000
.0000000
SYSTEM X AXIS
.6288304
.9492511
1.5264980
2.6471320
4.5258500
5.9935610
6.9393590
4.8390230
2.3080330
13.3092200
16.7915200
41.7861600
.0000538
.0000000
.0000000
(METERS)
IS 216.5484000 AT X=
508500.000 609000.000 509500.000
1.
1
2.
2.
i _
i
7.
16.
19.
12.
215.
_
605500
2922060
7073390
1777650
5434830
9679530
9255843
9615770
8298220
3799900
8600300
9612100
5484000
5913.732
0000000
0000000
1.2360030
1.2343310
.9569063
.5392794
.5648772
4.1450890
7.5753530
9.9173630
9.5169190
5.1210470
109.0088000
155.3864000
3.8970960
.0000015
.0000000
.5874423
.3717956
.1953373
.5245304
"i T> g o t; ') n
5J48836Q
6.7775340
5.4485790
1.9119500
27.0237300
138.98380'!"
63.3541500
6.0713320
.0004690
.0000000
.1554S5|
.1279575
.453S50|
1.334732!
'5 0 S 0 Q 0 1 r
A - Q i •! n n p
" ' "* " '' 1m
3 3970171
"'9579441
8.7907050
en ^ooo^nB
OC iopoc.nl
19.5558900
4.336996|
.009472!
.cooooor
.0, Y=4489000.C)
610000.000
611000.000
51200Q.OOO
- CONCENTRATION -
.3562494
.9619373
1.5006310
1.3383940
i cCQlIRn
.4903683
1.4357500
7 5257510
32.0906000
52.9497000
46.1479800
18.7355900
4.0098480
.1047905
.0000000
1,
1
i
i
•3
12.
40
48.
31.
14.
3,
.7521905
,2152960
.3666539
,9192530
. 1378810
.5732105
.5918760
,9623100
,5405300
,5968200
,5282300
,5463000
,9984690
,1766849
,0090002
.9521827
1.0632390
1.3435700
1.5154350
.3303909
1 .8319390
6.1586300
19.2873900
40.7618100
35.5159700
23.8167400
10.8235500
3.4552500
.2744120
.0000032
.3637392
.9753856
1 342059"
.7354209
.2248232
3.2802450
3.3994750
24.3718700
34.1396900
25.4697500
17.8261500
8.2445820
3.1392980
.3728145
.0000258
i
1"
12.
25.
27.
20.
13.
• 5.
2.
7544968
0347220
9668295
2452025
5794514
iq?oi 1 n
5591900
2739100
0921400
4153800
5433500
7851320
8085160
4362422
OOC131 1
.8096063
.5547487
1 Q p 1 0 '1 1
21 pOOTC
1 p c o c n i o
0 QVV'P'in
17.4358900
20 3775100
16 .5419800
13.0948200
8.4997210
4.5713420
2.1907600
5200550
.0011653
. 479(1340
•, c o ? 7 1 q
i n o •"> Q c £
•} 1 11) ,1 n Q
•3 -i /I '7 '; p ' '':
1 1 '.i o '! o PJ o n
IS. 9378500
1/1 i •? ceo on
, r . [ J •*• ;. J \}
11.2396900
3.9330140
5.5953920
3.3774980
1.7971560
.5455423
.0046919
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHCRTZ DATE , CASE 1, PAGE 'C
TABLE 6 (CONF)
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/CU8IC METER) FROM SOURCES '17
- HOUR(S) 0 TO 2300 -
- DISCRETE POINT RECEPTORS -
(THE MAXIMUM CONCENTRATION IS, 227.1075000 AT X= 505167.0, Y=«89107.0)
X Y CONCENTRATION X Y CONCENTRATION X Y CONCENTRATION
(METERS) (METERS ) (METERS) (METERS ) (METERS) (METERS ) '
605157.0 4489107.0 227.1076000
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHOSTZ DATE , CASE !, PAGE ':
TASLE 7
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/CUBIC METER) FROM SOURCES 118
- HOUR(S) 0 TO 2300 -
- GRID SYSTEM X AXIS (METERS) -
(THE MAXIMUM CONCENTRATION IS 65.3336300 AT X= 606000.0, Y=4489000.0)
502000.000 503000.000 504000.000 504500.000 505000.000 605500.000 505000.000 505500.000 5C7000.00C
Y AXIS (METERS ) - CONCENTRATION -
4497000.000
4496000.000
4495000.000
4494000.000
4493000.000
4492000.000
4491500.000
4491000.000
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487090.000
.3135576
.4302749
.5746663
.7817338
1.2206500
1.3315820
2.3091350
2.6643630
2.9568150
1.4772880
.0589084
.0000054
.0000000
looooooo
•.0000000
1.1191820
1.5187040
2.3235350
3.3899610
3.7246600
5.4932910
5.8106370
4.8580150
1.3521610
.5725184
.5232643
.5161535
.0001355
.0000000
.0000000
2.4490670
2.5337910
2'.3928300
1.2991400
.5328915
.7938274
2.0385550
3.7374350
3.9000510
3.3751590
2.1860000
.2956836
2.7162780
.2453971
.0000000
1.3756790
1.0875440
".5726251
.2389899
.9622607
3.5516540
4.4751840
4.2439250
4.3490660
1.7218000 •
.4017198
4.4699620
12.3209900
1.8873240
.0000166
.4973253
.2505524
.2453905
.9875341
3.2909870
4 4052410
4.5382950
4.7512300
1.8716770
.4834537
3.7743600
23.8142900
28.5942000
7.3827760
.0016518
.1553836
.2450550
.8452041
2.4715770
2.4518530
4.3713300
, 3.7920220
1.3459250
.5002166
3.7874290
'12.9385800
50.7727500
'35.0744900
17.9980300
.0189031
.2596258
.7950175
1.8052000
2.3515600
3.0045340
3.1392580
1'. 7306620
.5377139
2.3015290
11.1454100
29:6031100
55.3336300
39.9559200
18 (12384UO
.0523170
.7094654
•1.5225510
2. 1242720
2.4758090
7 IKQIflSn
i c i o TJ c n
.1855125
1. 4548860
6 5225160
1 8. 0556000
<3 ipannn
44.5825300
41.1220700
12.8005000
.1456552
1.26976ll
i 7 o 7 7 7 7P
1.82215CC
2.530550B
1 '>01RC'lB
. J -- •
1 1 ">. 1 ^ 0 •">
.383226fa
^H
: . 97 121 V
1 P 0 ' 0 1 r 0*
10 is T i c n f i
;17 cojjrugH
•59 95l2?cP
21.505990C
U.423570B
OHkC QsH
607500.000
Y AXIS (METERS )
- GRID SYSTEM X AXIS (METERS) -
(THE MAXIMUM CONCENTRATION IS 55.3336300 AT X= SOSQOO.'O, Y=4439000.0)
608000.000 508500.000 509000.000 609500.900 510000.000 511000.000
- CONCENTRATION -
612COO.QQO
4497000.000
1496000 ^00
4495000.000
'.491000.000
4193000.000
4492000.000
4191500.000
4491000.000
4190500.900
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4187000.000
1
1
i
2
i
i
5
1?
37
« 1
T 1
30
19
11
.5469320
.5470840
.3253250
.1053330
.0285170
.5079088
.5406560
.6531310
.9165200
.3819400
.7538000
.9935600
.5042800
.0185500
.3544337
1.3850730
1. 3814250
1.3545220
i K, it 1 o /i 1 n
.3159035
'! .5147780
4.5690850
9.5975520
24.2241200
39.5715300
42.7579900
28.7169900
16.4967900
9.3390420
.4713520
i
i
i
i
6
14.
19
40.
32
. 28.
16
7.
.1510450
.4539680
.7132410
7510139
,2472275
.7550530
.5648190
.7412200
.3172900
,3504600
.9927500
.5403300
.3192200
.5187000
.7421083
1
1.
i
4.
g.
19,
34,
n
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ , DATE , CASE 1, PAGE :2
TABLE 7 (CONT)
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/CU8IC METER) FROM SOURCES 1.^2
- HOUR(S) 0 TO 2300 -
- DISCRETE POINT RECEPTORS -
(THE MAXIMUM CONCENTRATION IS, 32.5627500 AT X= 605157.0, Y=d48910?.fl}
X Y CONCENTRATION X Y CONCENTRATION X Y CONCENTRATION
(METERS) (METERS ) (METERS) (METERS ) . (METERS) (METERS )
605167.0 4489107.0 32.5627500 .
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ CATE , CASE 1, °AGE 13
TASLE 8
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/CUSIC METER) FROM SOURCES 'M6 -117
- HOUR(S) 0 TO 2300 -
- GRID SYSTEM X AXIS (METERS) -
(THE MAXIMUM CONCENTRATION IS 745.7665000 AT X= 605000.0, Y=4489000.C)
602000.000 503000.000 604000.000 604500.000 605000.000 605500.000 606000.000 505500.000 507001.000
Y AXIS (METERS ) - CONCENTRATION -
4497000.000
4496000.000
4495000.000
4494000.000
4493000.000
4492000.000
4491500.000
4491000.000
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487000.000
AXIS (METERS
449700C.QQ9
•1495000.000
4495000.000
4:194000.000
/i/io-jnQQ inn
.ua9nnn'n(jn
4491500.000
4491000.000
4490500.000
4490000.000
4489500.000
4489000.000
4483500.000
4488000.000
4437000.000
1.1766290
.5453575
.1992516
.0212535
.0002772
.0000000
.9000000
.0000000
.0000000
.0000000
.0000000
.0000000
.0000000
.0000000
.doooooo
507500.000
•\
i
.2355535
.8273913
2.5115930
a. 5162800
5.7527990
4.5541110
1.3601570
7.5906390
40.8456100
151.2740000
137.0201000
59.9789900
10.3501100
.0968305
.0000000
2.9873030
3.8793260
4.6771250
4.3541080
2.0248780
.2559888
.0207502
.0002185
.0000000
.0000000
• .0000000
.0000000
.0000000
.0000000
.0000000
(THE MAXIMUM
508000.000
.7352602
1.9399060
3.2339130
o 7 g o o o n n
5. 15535 1C
1 .0157350
3.3349050
17.2099900
68.1773400
136.8968000
101.9296000
43.5936900
9.3623940
.2403299
.0000000
1
1
3
8
13
19
20
.4361686.
.6784121
.1595200
.9637840
.9665250
.3831010
.6418500
.0087500
.0698300
.8253110
.1485035
.0000000
.0000000
.0000000
.0000000
- GRID
CONCENTRATION
.3645229
.5277020
.7789738
1.2123910
2.1900370
4.4528310
6.7865250
11.0719100
19.1530000
32.4045500
.28.3910400
2.0146030
.0000000
.0000000
.0000000
1.4380950
2.1842790
3.5260900
5.1006700
10.6826700
16.5308500
16.9041900
12.3230700
6.0804200
35.8754300
70.8952200
745.7665000
.0004718
.9000000
' .9000000
2.9419150
3.9225160
5.0327660
5.8492090
4.5832890
2.0629360
4.2112210
15.9259900
34.8589300
43.0955700
44.3061200
706.9080000
1.7807050
.0000000
.0000000
2.
2_
2.
1.
1
0
15.
20.
19.
14.
285.
379.
o_
8195000
8230010
2116590
2293410
1770191
5596630
7151900
3682800
6836300
8522700
5375000
4773000
5309860
0000038
0000000
1
i
c
1(1
n
in
4
C1
^ -J
331
150
13
.3403080 .3512211
.3495857 .2709561
rt 0 7 P 7 1 7 Q 0 0 ^ 7 1 £
.0572999 ' 3.6756521
c ni o c o ft p. c. ,1 •) « 1 Q |
.i 1 1 rt ') ft n o c ft 7 .i rt Q n
.7105000 :.''852':7'|
9500400 2.'nn5"l
* J " 1
. 1027500 130 ^737001
.1 1 <; i P n n 910 n c •; n n o J
. ^ -j -J ^ u u '• _ i r . \t •. j -• u . •
2Q9/jnoo Q7 "•'>1::^^''> "
59SRPO'1 1 ' .i^7/l,(ilp
. 0 1, w J \j V. ' * i . t '. . '- *t v, \j
im 1 1 jg n9 MOO]
.0000009 .9000901
SYSTEM X AXIS (METERS) -
IS 745.7565000 AT X^
6085nn.nno snonoo.nnn
1
2
2
3
i)
1
3
27
85
104
70
33
9
.5142150
.4502710
.7430080
.8553560
1 0 0 (! "7 1 0
.5742550
.4305550
.8477900
.9785500
.0459000
.9103200
.0955000
.1344720
.4068209
.0000004
snusoo.fion
- CONCENTRATION -
1.9164810
2.1353890
2.5956190
3.0447210
.6751694
4.2597700
13.6138100
40.5492800
85.7049800
76.2204100
53.0556700
24.6025200
7.8954890
.5255742
.0000073
1.9522450
2.6885030
1.4777480
C1 •JI;KOQ
7.4586090
C 1 1 1QT £ 00
72.3553200
57.1895500
39.3166500
18. 8374300
7.1125770
.8390155
.0000538
505000.0, Y=4439000
610000.000
1 cip?-)in
2.0724260
1 .9423570
.4982588
1 o <; 9 n 1 9 n
10.5773300
25.5561900
53.5920100
57.5800700
44.5206700
30.5352700
15.3494500
6.325447C
.9824557
0002971
m
. \. /
611000 000
1
i.
f
18.
36.
14.
35.
28.
19.
10.
4.
1
6!QSfJ7"
•5 1 Q 0 n P p
383752?
5093515
0 1 P. 0 9 7 0
o, 1 1 7 1 n o
5649790
0953200
3843599
9453500
(1750200
3559000
9449700
1585829
9025299
s 1 o n 0 n n o n
•i
•j
.>,-
r.
•o
24
10
10
. (_
7
4
1
'; n e •• c 7 o
o -j 7 •:• T o c
Qoioort n
0 ,1 " •: ••> c n
0 T i 1 n 0
<: •} ,1 o •) n i%
•jliT'nn
^ o 7 o n o n
Q P 1 o Q no
Q 1 1 PA on
p,9fW}c n
Q ( 1 1 « 1 J
.229981?
ft 1 Q C '1 1 0
I
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORT! DATE , CASE 1, PAGE
TASLE 8 (CONT)
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/C'JBIC METER) FROM SOURCES 115 -117
_ uoiiocc'i n TO Tjnn _
i > V v I \ ^ •—• y \* IV (.-JWV
- DISCRETE POINT RECEPTORS -
(THE MAXIMUM CONCENTRATION IS, 946.2949000 AT X= 505167.0, Y=4489107.0)
X Y CONCENTRATION X Y CONCENTRATION X Y CONCENTRATION
(METERS) (METERS ) (METERS) .(METERS ) (METERS) (METERS )'
605167.0 4489107.0 945.2949000
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORTZ DATE , CASE 1, PAGE 15
TABLE 8
24 HOUR GROUND LEVEL CONCENTRATION (MICROGRAMS/CU8IC METER) FROM SOURCES 115 -118
- HOUR(S) 0 TO 2300 -
- GRID SYSTEM X AXIS (METERS) -
(THE MAXIMUM CONCENTRATION IS 769.5807000 AT X= 6Q5000.0, Y=4489000.C)
602000.000 503000.000 504000.000 504500.000 505000.000 505500.000 505000.000 506500.000 S07COO.OCO
Y AXIS (METERS ) - CONCENTRATION -
4497000.000
4495000.000
4495000.000
4494000.900
4493000.000
4492000.000
4491500.000
4491000.000
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487000.000
1
1
1
1
2
2
2
i
.4901870
.0755320
.7739279
.8029872
.2209370
.8315820
.,3091350
.5643630
.9568150
.4772880
.0689084
.0000054
.0000000
.0000000
.0000000
507500.000
AXIS (METERS
4497000.000
4495000.000
4495000.000
•1494000.000
4492000.000
4192000.000
4491500.000
4491000.000
4490500.000
4490000.000
4489500.000
4489000.000
4488500.000
4488000.000
4487000.000
}
1
2
d
K
6
<;
3
'4
56
198
178
90
29
1 1
4.1064860
5.3980290
7.0006590
7.7440680
5.7495380
5.7492790
5.8313880
4.8582350'
1.8521610
.6725184
.5232643
.5161535
.0001355
.0000000
.0000000
(THE MAXIMUM
508000.000
2.8852360
3.2122040
3.5523500
3.2529240
4.4994160
9.5759280
15.5804000
22.7462000
23.9598800
11.2024700
2.3345040
.2956836
2.7162780
.2453971
.0000000
- GRID
CONCENTRATION
1.7412020
1.5152460
1.3515990
1.5013800
3.1522970
8.1144850
11.2517100
15.3158400
23.5020700
34.1263500
28.7927700
5.4845650
12.3209900
1.8873240
.0000166
1.9354210
2.4348310
3.7714800
7.0883040
13.9735600
20.9361000
21.4424800
17.0743000
7.9520980
37.3588800
74.6695800
769.5807000
28.6946800
7.8827750
.0015618
3.
4.
5.
8.
. 8.
6.
8.
18.
35.
46.
57.
757.
36.
17.
0982990
1675710
8779710
3207860
0351470
4372550
0032420
7719100
4591900
8830000
2448000
5808000
8551900
9980300
0189031
3.0791260
3.5190190
4.0158690
4.0809010
4.1816530
1 1 .7990200
17.4958600
20.9059900
21.9852600
25.9975800
315.1409000
444.8110000
49.4869000
18.0238400
.0523170
2
2
2
••>
9
i 1
14
11
11
31
374
204
.54
12
.0497730 1. 5209821
.3722400 2. 0587441
.5520440 2. 75252^
c o » 1 n o n < •'.' n c ? 1 ••
. -' ^' H , « u -. J . ~ » * ; • *M
77iip,jn c Oldftngjj
o?Q?5^ ^ Q*; i oopn
lOR^HIft n £ 7 /Wl A •
. , ^_ ... „ ,..,•,,.,•
i1 1 d Q 0 0 0, fi H7107oB
.TUT.4..JW W.wlt'vl-J
.0207100 31.3042300
'!5°3600 160 ''4*5004
.7249000 255.149000|
.9750000 137.08590cl
9486800 35 0934300
flniRinn i/i /u c icnl
. 9VIS..VW .-, .->-- . Jv|
.1456552 .3055951
SYSTEM X AXIS (METERS) -
IS 769.5807000 AT X=
608500.000 509000.000
509500.000
505000
.0, YM489000.0)
510000.000
611000.000
512000.000
- CONCENTRATION -
.7824900
.3744750
.4379190
.9225130
.7914150
.1720200
.9008140
.3437700
.7521400
.5560000
.7738000
.9726600
.9543800
"53900
.3644837
2.1203330
3.3213310
5.0984340
5.4242720
5.5012530
2 .5305130
7.9039910
26.9075400
92.4014600
175.5683000
144.5877000
72.3106800
25.8591800
9.5793710
.4713520
2.6582600
3.9142390
4.4612490
4.5177800
2.5457000
4.3403220
14.9953700
42.5890100
115.8258000
144.3964000
103.9031000
61.5358300
25.4537000
7.9255210
.7421088
3.0758480
3.6631620
3.8334720
3.3325130
1.2546840
8.3285240
22.3257300
50.0235100
120.0477000
114. .1141000
82.9730500
47.5352400
21.3157200
7.9752100
.7505735
3.0042030
3.2553250
3.2715640
1.5400540
1.7653810
12.3707100
29.9277600
74.7041500
104.5034000
89.0978300
65.3772400
37.4559900
19.3997500
8.4880530
.7901954
2.
2.
2.
7550700
9227390
1825450
7884794
3.5234380
17
40.
77.
87.
71.
52.
32.
17.
7.
5753400
7653100
4123300
4419900
1388700
6987800
5167500
5595500
7687810
92525!1
1. 5423740
.5453243
1 .5539020
8.0529610
30.5239800
56.5495700
59.0581300
58.3218300
48.3125400
35.2530100
22.7806800
13.2539600
5.9364420
1.0.974830
•1
i
1 0
41
52
Slj
41
35
26
17
10
5
i
- ,1 C 77 0 ;1
-1569'I5
.3927575
.3100500
••; s P fi c n n
FflOOK00
POO 'snn
.557790C
.5529600
.3335100
.3083400
.5959300
.9959300
Q 1 1 ft oor)
1945000
-------
EXAMPLE SHORT TERM PROBLEM - PROGRAM SHORT! DATE , CASE 1, PAGE 15
TABLE 9 (CONT)
24 HOUR C-ROUNO LEVEL CONCENTRATION (M1CR06RAMS/CUSIC METER) FROM SOURCES 115 -118
- HOUR(S) 0 TO 2300 -
- DISCRETE POINT RECEPTORS -
(THE MAXIMUM CONCENTRATION IS, 978.8577000 AT X= 605157.0, Y=4489107.C)
X Y CONCENTRATION X Y CONCENTRATION X Y CONCENTRATION
(METERS) (METERS ) (METERS) (METERS ) (METERS) (METERS )
1
505157.0 4439107.0 978.8577000
-------
Figure 3-3 Sample POSTZ Input File
-------
THIS IS A
01121
1 1986
3
602000.
4
4494000.
1
605167.
4489107.
2
116
117
5
602000.
605167.
602000.
603000.
604000 .
2
116 2.
117 2. ,
TEST OF THE POSTZ PROGRAM
0000
603000.
4495000 .
4494000.
4489107.
4495000.
4496000 .
4497000 .
00000100000000
604000.
4496000. 4497000.
1000.
1000.
1000.
1000.
500.
1 0
-------
Figure 3-4 Sample POSTZ Output File
-------
POSTZ - VERSION1 1.0
A POST PROCESSOR FOR THE
SHORT! AIR QUALITY MODEL
JULY, 1986
-------
POSTZ - A POST PROCESSOR FOR THE SHORTZ MODEL
POSTZ RUN TITLE: THIS IS A TEST OF THE POSTZ PROGRAM
SHORTZ RUN TITLE: EXAMPLE SHORT TERM PROBLEM - PROGRAM SHQRTZ
POSTZ - °AGE NO.
IS«(1) RESTRICT TIME LIMITS (1=YES,0=NO)
ISW(2) LIMIT RECEPTORS TO ANALYZE (1=YES,0=NO)
ISW(3) LIMIT SOURCES TO ANALYZE (1=YES,0=NO)
ISW(4) SPECIFY BACKGROUND CONCENTRATIONS (0=NO,1=UNIFORM,2^BY RECEPTOR)
ISW(5) SCALE CONCENTRATIONS FOR SPECIFIC SOURCES ,(1=YES,0=NO)
1-HOUR AVERAGE ANALYSIS:
!SW(5) HIGH-5 TABLE PREPARED (1=YES,0=NO)
ISW(7) TOP 50 TABLE PREPARED1 (1=YES,0=NO)
!SW(S) EXCEEDANCE TABLE PREPARED (1=YES,0=NG)
3-HOUR AVERAGE ANALYSIS:
!SW(9) HIGH-5 TABLE PREPARED (f=YES;0=NO)
ISW(IO) TOP 50 TABLE PREPARED (1=YES,0=NO) •
ISW(II) EXCEEDANCE TABLE PREPARED (1=YES,9=NO)
3-HOUR AVERAGE ANALYSIS:
!SW(12) HI6H-5 TABLE PREPARED (1=YES,0=NO)
!SW(13) TOP-50 TABLE PREPARED (1=YES,0=NO)
ISW(H) EXCEEDANCE TABLE PREPARED (1=YES,0=NO)
2 i-HOUR AVERAGE ANALYSIS: '
ISW(15) HIGH-5 TABLE PREPARED (1=YES,Q=NO)
ISW(16) TOP-50 TABLE PREPARED (1=YES,0=NO)
ISW(17) EXCEEDANCE TABLE PREPARED (1=YES.O=NO)
ANNUAL AVERAGE ANALYSIS:
ISW(18) HIGH-5 TABLE PREPARED (1=YES,Q=NO)
!SW(19) TOP-50 TABLE PREPARED (1=YES,0=NO)
ISW(20) EXCEEDANCE TABLE PREPARED (l^YES.O^NO)
/>YCRAG!NG TI^C ,4NALVSIS
pn uTCu.c TAOI c oocoioen M-VCC n-Mfi'i
^c ' f i.iWi. -j inoLfc. ' Au_i rt-\uJ \ • ~ . Lw, w~i>*j/
no-sn T4Pi c opcoApcn CI-YC': n-\m^
ui •j w irtwt.— LUi_i.i(Si-L/ \ i ~ 1 taW t -j >\^j /
ISW(23) EXCEEDANCE TABLE PREPARED (1=YES,0=NO)
!SW(2d) BLOCK OR RUNNING AVERAGES (DEBLOCK, 1-RUNNING)
PAI M<; oni irv n=vcc n=Mn^
i i«
JULIAN DAY FOR START OF ANALYSIS
YEAR FOR START OF ANALYSIS
JULIAN DAY FOR END OF ANALYSIS
YEAR FOR END OF ANALYSIS
1986
1986
-------
RESTRICTING ANALYSIS TO CERTAIN RECEPTORS
THE FOLLOWING RECEPTORS WERE INCLUDED
POST! - pi.GE MO
(
*.
(
1
\
!
\
502000.,
503000
504000.,
505167.,
4494000.)
4195000.)
4496000.)
4489107.)
i
\
i
\
(•
i
503000. ,
504000. ,
502000.,
4424000.)
4495000.)
4497000.)
' 504000
f 502000.,
( 503000.,
4494000.)
4495000. N
4497000.)
( 5Q200n
( 503000.,
( 604000.,
f t a c n p n \
«/. ORnno >
r ^ «' u w w w > /
4497000.)
-------
POSTZ - PAGE NO. 3
THE ANALYSIS HAS SEEN CONFINED TO CERTAIN SOURCES
THE FOLLOWING SOURCES (SY IDENTIFICATION NUMBER) WERE INCLUDED
•MR
-------
I
POSTZ - PAGE NO.
SEPARATE BACKGROUND CONCENTRATIONS HAVE SEEN
SPECIFIED FOR EACH RECEPTOR
(X-COORD,Y-COORD) VALUE
( 602QOO., 4494090.) 1000. ( 605167., 4489107.) 1000. ( 502000., 44Q5000.) 1000. ( 603000., 4495000.)
504000., -U97000.) 500.
-------
POSTZ - PAGE NO. 5
THE CONCENTRATIONS FROM CERTAIN SOURCES WERE
SCALED SY THE FOLLOWING AMOUNTS
(IDENTIFICATION NUMBER,SCALING VALUE)
( 115, 2.00) ( 117, 2.00) ( 118, 1.00') (
-------
HIGH FIVE TA2LE FOR 24 HOUR AVERAGES
°OSTZ - PAGC NO
METERS
X
502000.
503000.
504000.
502000.
503000.
604000.
602000, '
603000.
504000
502000
503000.
504000.
605157.
V
4494000.
4494000.
4494000.
4495000.
4495000.
4495000.
4496000.
4496000.
4496000.
4497000 .
4497000.
4497000.
4489107.
CLEV.
286.
349.
232.
317. .
353.
364.
349."
317..
323.
32°
347.
320,.
274.
HIGHEST
1000. 04( 1
8.71( 1
3.93( 1
1000. 40( 1
9.35( 1
2.32( 1
1.29( 1
1007. 76( 1
135'' 1
2.35( 1
5.970 1
500. 87( 1
2892. 59( 1
VALUE(OAY,
2ND HIGH
,24)
, 2)
, 2)
,24) •
, 2)
, 2)
r ' i
,24)
2)
1 N
, 2)
,24)
,24)*
958. 38(
8.71(
3.93(
958. 73(
9.35(
2.32(
1.29(
965. 09(
''.36(
2.35(
5.97(
430. 04(
2850. 92(
1,23)
1, 3)
1, 3)
1,23)
1, 3)
1, 3)
1, 2)
1,23)
!, 3)
1 1N
1 , '-I
1, 3) '
1 23)
1,23)*
ENDING ^O'JR'1
3RD HIGH
915. 71( 1
8.71( 1
3.93( 1
917. 07( 1
9.35( 1
2.32( 1
1.29( 1
924. 43( 1
,1.35( 1
2.35( 1
5.97( 1
459. 21( 1
2807. 47( 1
,22)
, 0
, 4)
,22)
, 4)
, *)
, 3)
,22)
n
, 3)
, 4)
22^
^22)
4TH HIGH
875. 04( 1,21)
8.71( 1/5)
3.93( 1, 5)
875. 40( 1,21)
o 3<5( 1 <^
. . J J ^ , , V / ,
0 it! 1 -\
(..*•.{ , , « i
1.29( 1, 4)
882. 76( 1,21)
;.35( 1,5)
235' 1 4 "*
5.97( 1, ,5)
438. 37( 1,21)
2594 41 ' 1 '1 ^
qTU LI 1C
00-3 oof
u w w' . 0 3 k
3.?1(
i QO/
-'•--' \
833. 73(
9.35(
2.32(
!.29(
0/11 pfl/
1.25(
•1 -5 C ,'
C Q "t I
fit S :\ 1
2586. 46(
u
1,20)
1, 5)
1, 5)
' r - J /
^ 5)
1. 1 5 )
s \
> t •* /
1 .20)
, ^ \
- i
; C \
• • J
•i on ^
-. ••*n\
I
- DENOTES PEAK VALUE
-------
APPENDIX A
COMPLETE FORTRAN LISTING FOR THE POSTZ PROGRAM
(UNIVAC Version)
-------
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
PROGRAM POSTZ
POSTZ - A POST PROCESSOR FOR THE SHORTZ MODEL. POSTZ WAS WRITTEN
BY KIRK D. WINCES, TRC ENVIRONMENTAL CONSULTANTS, INC.,
15924 22ND AVE. SE, MILL CREEK, WA, 98012 (206) 485-2992
THE SHORTZ MODEL MUST BE RUN FIRST AND THE OPTION SELECTED TO
WRITE AN OUTPUT TAPE (CARD 1, COLUMN 1 SHOULD BE A "2").
THE SHORTZ MODEL HAS THE CAPABILITY TO USE A LARGE NUMBER
OF RECEPTORS. IF POSTZ IS TO BE USED TO ANALYZE THE OUTPUT
TAPE FROM A SHORTZ RUN, THE NUMBER OF RECEPTORS MUST BE LIMITED
TO A TOTAL OF NO MORE THAN 300 (GRID AND DISCRETE'RECEPTORS).
THE POSTZ PROGRAM HAS TWO INPUT STREAMS: THE OUTPUT TAPE FROM THE
SHORTZ MODEL AS NOTED ABOVE, AND A "CARD IMAGE" FILE SPECIFYING
HOW THE POST PROCESSING IS TO BE ACCOMPLISHED. THERE ARE MANY
OPTIONS FOR,THE POSTZ ANALYSIS OF AN OUTPUT TAPE. IN GENERAL,
THE OPTIONS TO BE USED ARE SPECIFIED WITH A'SERIES OF SWITCHES
ON CARD 2. THE REMAINING CARDS TO BE READ ARE DETERMINED BY
WHICH SWITCHES HAVE BEEN SET. THE FOLLOWING DESCRIPTION IS
INTENDED TO PROVIDE AN OUTLINE FOR THE PROGRAM .OPTIONS AND
REQUIRED INPUT VARIABLES IN THE CARD IMAGE FILE.
USER'S INSTRUCTIONS FOR CREATING THE CARD,IMAGE
INPUT FILE
CARD CARD
GROUP NO.
' 1 1
PARAMETER FORMAT
PTITLE ' A80
ISW(l)
2412
ISW{2)
ISW(3)
DESCRIPTION
RUN TITLE FOR THE POSTZ RUN,
DIFFERENT FROM THE RUN TITLE
FOR THE SHORTZ RUN USED TO
CREATE THE OUTPUT TAPE
SWITCH TO DETERMINE IF TIME
PERIOD OF ANALYSIS IS TO BE
RESTRICTED. IF =0, USE
ENTIRE TIME PERIOD ON THE
OUTPUT TAPE FROM THE SHORTZ
RUN. IF =1, INCLUDE CARD
GROUP 4 TO SPECIFY TIME
PERIOD RESTRICTION
SWITCH TO DETERMINE IF
ANALYSIS IS TO BE RESTRICTED
TO CERTAIN RECEPTORS. IF
=0, ALL RECEPTORS WILL BE
INCLUDED IN THE POSTZ RUN.
IF =1, RESTRICT ANALYSIS
ONLY TO THOSE RECEPTORS
SPECIFIED IN CARD GROUP 5.
SWITCH TO DETERMINE IF THE
ANALYSIS IS TO BE RESTRICTED
TO ONLY CERTAIN SOURCES. IF
-------
C =0, ALL SOURCES WILL BE
C INCLUDED. IF =1, RESTRICT
C ANALYSIS TO ONLY THOSE
C SOURCES SPECIFIED IN CARD
C GROUP 6.
C
C ISW(4) SWITCH TO DETERMINE IF
C BACKGROUND CONCENTRATIONS
C ARE TO BE ADDED TO THE
C SHORTZ MODEL RESULTS. IF
C =0, NO BACKGROUND
C CONCENTRATIONS ARE INCLUDED;
C IF =1, UNIFORM BACKGROUND IS
C TO BE ADDED TO ALL RECEPTORS
C AS SPECIFIED IN CARD GROUP
C 7; IF =2, SEPARATE
C BACKGROUND IS TO BE ADDED TO
C EACH RECEPTOR AS SPECIFIED
C IN CARD GROUP 7.
C
C ISW(5) SWITCH TO DETERMINE IF
C CONTRIBUTIONS FROM
C PARTICULAR SOURCES ARE TO BE
C SCALED. IF =0, NO SOURCE
C SCALING IS DONE; IF =1,
C SCALING WILL BE PERFORMED
C FOR SOURCES AND VALUES
C SPECIFIED IN CARD GROUP 8.
C
C ISW(6) SWITCH TO DETERMINE IF A
C HIGH-5 OUTPUT TABLE IS TO BE
C PREPARED FOR 1-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C HIGH-5 TABLE FOR 1-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A HIGH-5 TABLE BY
C RECEPTOR FOR 1 HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW(7) SWITCH TO DETERMINE IF A
C TOP-50 TABLE IS TO BE
C PREPARED FOR 1-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C TOP-50 TABLE FOR 1-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A TOP-50 TABLE FOR
C 1-HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW(8) SWITCH TO DETERMINE IF AN
C EXCEEDANCE TABLE IS TO 3E
C PREPARED FOR 1-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C EXCEEDANCE TABLE FOR 1-HOUR
C CONCENTRATIONS. IF =1,
-------
C PREPARE AN EXCEEDANCE TABLE
C FOR 1-HOUR AVERAGE
C CONCENTRATIONS. EXCEEDANCE
C VALUE SPECIFIED IN CARD GROUP
C 11 (ONEEX).
C
C ISW(9.) SWITCH TO DETERMINE IF A
C HIGH-5 OUTPUT TABLE IS TO BE
C PREPARED FOR 3-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C HIGH-5 TABLE FOR 3-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A HIGH-5 TABLE BY
C RECEPTOR FOR 3 HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW(IO) SWITCH TO DETERMINE IF A
C TOP-50 TABLE IS TO BE
C PREPARED FOR 3-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C TOP-50 TABLE FOR 3-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A TOP-50 TABLE FOR
C 3-HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW(ll) SWITCH TO DETERMINE IF AN
C EXCEEDANCE TABLE IS TO BE
C PREPARED FOR 3-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C EXCEEDANCE TABLE FOR 3-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE AN EXCEEDANCE TABLE
C FOR 3-HOUR AVERAGE
C CONCENTRATIONS. EXCEEDANCE
C VALUE SPECIFIED IN CARD GROUP
C 11 (THREX).
r;
C ISW(12) SWITCH TO DETERMINE I? A
C HIGH-5 OUTPUT TABLE IS TO BE
C PREPARED FOR 8-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C HIGH-5 TABLE FOR 8-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A HIGH-5 TABLE BY
C RECEPTOR FOR 8 HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW(13) SWITCH TO DETERMINE IF A
C TOP-50 TABLE IS TO BE
C PREPARED FOR 8-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C TOP-50 TABLE FOR 8-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A TOP-50 TABLE FOR
-------
C 8-HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW(14) SWITCH TO DETERMINE IF AN
C EXCEEDANCE TABLE IS TO BE
C PREPARED FOR 8-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C EXCEEDANCE TABLE FOR 8-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE AN EXCEEDANCE TABLE
C FOR 8-HOUR AVERAGE
C CONCENTRATIONS. EXCEEDANCE
C VALUE SPECIFIED IN CARD GROUP
C 11 (EIGEX).
C
C ISW(15) SWITCH TO DETERMINE IF A
C HIGH-5 OUTPUT TABLE IS TO BE
C PREPARED FOR 24-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C HIGH-5 TABLE FOR 24-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A Hl'GH-5 TABLE BY
C RECEPTOR FOR 24 HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW(16) SWITCH TO DETERMINE IF A
C TOP-50 TABLE IS TO BE
C PREPARED FOR 24-HOUR AVERAGE
C CONCENTRATIONS. IF =0, NO
C TOP-50 TABLE FOR 24-HOUR
C CONCENTRATIONS. IF =1,
C PREPARE A TOP-50 TABLE FOR
C 24-HOUR AVERAGE
C CONCENTRATIONS.
C
C ISW{17) SWITCH TO DETERMINE I? AN
C EXCEEDANCE TABLE IS TO BE
C PREPARED FOR 24-HOUR AVERAGE
C CONCENTRATIONS. IF =0 NO
C EXCEEDANCE TABLE FOR 24-HOUR
C CONCENTRATIONS. IF =1
C PREPARE AN EXCEEDANCE TABLE
C FOR 24-HOUR AVERAGE
C CONCENTRATIONS. EXCEEDANCE
C VALUE SPECIFIED IN CARD GROUP
C 11 (TWYEX).
C
C ISW(18) SWITCH TO DETERMINE IF A
C HIGH-5 OUTPUT TABLE IS TO BE
C PREPARED FOR ANNUAL AVERAGE
C CONCENTRATIONS. IF =0, NO
C HIGH-5 TABLE FOR ANNUAL
C CONCENTRATIONS. IF =1,
C PREPARE A HIGH-5 TABLE BY
C RECEPTOR FOR ANNUAL AVERAGE
-------
C CONCENTRATIONS. MUST INPUT
C NUMBER OF DAYS PER YEAR IN
C CARD GROUP 9.
C
C ISW(19) SWITCH TO DETERMINE IF A
C TOP-50 TABLE IS TO BE
C PREPARED FOR ANNUAL AVERAGE
C CONCENTRATIONS. IF =0, NO
C TOP-50 TABLE FOR ANNUAL
C CONCENTRATIONS. IF =1,
C PREPARE A TOP-50 TABLE FOR
C ANNUAL AVERAGE
C CONCENTRATIONS.
C
C ISW(20) SWITCH TO DETERMINE IF AN
C EXCEEDANCE TABLE IS TO BE
C PREPARED FOR ANNUAL AVERAGE
C CONCENTRATIONS. IF =0, NO
C EXCEEDANCE TABLE FOR ANNUAL
C CONCENTRATIONS. IF =1,
C PREPARE AN EXCEEDANCE TABLE
C FOR ANNUAL AVERAGE
C CONCENTRATIONS. EXCEEDANCE
C VALUE SPECIFIED IN' CARD GROUP
C 11 (ANNEX).
C.
C ISW(21) SWITCH TO DETERMINE IF A
C HIGH-5 OUTPUT TABLE IS TO BE
C 'PREPARED 'FOR USER SPECIFIED
C AVERAGE CONCENTRATIONS. IF
C =0, NO HIGH-5 TABLE FOR USER
C SPECIFIED CONCENTRATIONS.
C IF =1, PREPARE A HIGH-5
C TABLE BY RECEPTOR FOR USER
C SPECIFIED AVERAGE
C CONCENTRATIONS.
C
C ISW(22) SWITCH TO DETERMINE I? A •
C TOP-50 TABLE IS TO BE
C . ' PREPARED FOR USER SPECIFIED
C AVERAGE CONCENTRATIONS. IF
C =0, NO TOP-50 TABLE FOR USER
C SPECIFIED CONCENTRATIONS.
C IF =1, PREPARE A TOP-50
C TABLE FOR USER SPECIFIED
C AVERAGE CONCENTRATIONS.
C
C ISW(23) SWITCH TO DETERMINE IF AN
C EXCEEDANCE TABLE IS TO BE
C PREPARED FOR USER SPECIFIED
C AVERAGE CONCENTRATIONS. IF
C =0, NO EXCEEDANCE TABLE FOR
C USER SPECIFIED
C CONCENTRATIONS. IF =1,
C PREPARE AN EXCEEDANCE TABLE
-------
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c.
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
FOR USER SPECIFIED AVERAGE
CONCENTRATIONS. EXCEEDANCE
VALUE SPECIFIED IN CARD GROUP
11 (OTHEX).
ISW(24) SWITCH TO DETERMINE WHETHER
BLOCK AVERAGES OR RUNNING
AVERAGES ARE COMPUTED IN ALL
ANALYSES. IF =0, THEN BLOCK
AVERAGES ARE COMPUTED. IF
=1, THEN RUNNING AVERAGES
ARE COMPUTED.
ISW(25) SWITCH TO DETERMINE IF EPA
CALMS POLICY IS TO BE
IMPLEMENTED. IF =0, CALMS
POLICY IS NOT IMPLEMENTED, IF
=1, CALMS POLICY IS
IMPLEMENTED.
JULST 15 JULIAN DAY FOR THE START OF
THE DATA ON THE TAPE FROM
THE SHORTZ RUN.
IYRI5 YEAR FOR THE START OF THE
DATA ON THE TAPE FROM THE
SHORTZ RUN. IT CAN EITHER
BE OF THE FORM 1986 OR OF
•THE FORM 86 — THE PROGRAM
WILL USE EITHER.
[THIS CARD IS ONLY READ IF ISW(1)=1]
IDYST 15 THE JULIAN DAY FOR THE START
OF THE ANALYSIS.
IYRST
IDYEND
IYREND
15
15
15
THE YEAR FOR THE START OF
THE ANALYSIS.
THE JULIAN DAY FOR THE END
OF THE ANALYSIS.
THE YEAR FOR THE END OF THE
ANALYSIS.
[THIS CARD GROUP IS ONLY READ IF ISW(2)=1]
1 NXINCL 15 THE NUMBER OF X GRID POINTS
TO INCLUDED IN THE ANALYSIS.
2 [THIS CARD IS NOT READ IF NXINCL=0]
[THIS CARD CAN BE REPEATED AS OFTEN AS NECESSARY
TO INCLUDE ALL THE X GRID POINTS]
XX 8F10.0 THE X GRID POINTS TO INCLUDE
IN THE ANALYSIS. NOTE, EACH
X VALUE TO BE USED MUST BE
SPECIFIED HERE. IT IS NOT
-------
C POSSIBLE TO SPECIFY THE
C STARTING X VALUE AND THE
C INCREMENT FOR OTHER VALUES
C AS IT IS IN SHORTZ.
C
C 3 [THIS CARD IS NOT READ IF NXINCL=0]
C NYINCL 8F10.0 THE NUMBER OF Y GRID POINTS
C TO. INCLUDE IN THE ANALYSIS.
C
C 4 [THIS CARD IS NOT READ IF NXINCL=0]
C [THIS CARD CAN BE REPEATED AS OFTEN AS NECESSARY
C TO INCLUDE ALL THE Y GRID POINTS]
C YY 8F10.0 THE Y GRID POINTS TO BE
C INCLUDED IN THE ANALYSIS.
C (SEE NOTE ON FORMAT FOR X
C GRID POINTS ABOVE).
C
C 5 NXYINC 15 THE NUMBER OF DISCRETE
C RECEPTORS TO INCLUDE IN THE
C ANALYSIS.
C
C 6 [THIS CARD IS NOT READ IF NXYINC=0]
C [THIS CARD CAN BE REPEATED AS OFTEN AS NECESSARY
C TO INCLUDE.ALL THE DISCRETE RECEPTORS]
C XDISCR 8F10.0 THE X VALUES FOR THE
C DISCRETE RECEPTORS.
C
C 7 [THIS CARD IS' NOT READ IF NXYINC=0]
C [THIS CARD CAN BE REPEATED AS OFTEN AS NECESSARY
C TO INCLUDE ALL THE DISCRETE RECEPTORS]
C YDISCR 8F10.0 THE Y VALUES FOR THE
C DISCRETE RECEPTORS.
C
C 6 [THIS CARD GROUP IS ONLY READ IF ISW(3)=1]
C 1 NSELIM 15 THE NUMBER OF SOURCES TO
C INCLUDE IN THE RUN
C
C 2 [THIS CARD IS REPEATED NSELIM TIMES]
C K 15 THE SOURCE IDENTIFICATION
C NUMBER OF THE SOURCE TO 3E
C INCLUDED IN THE RUN.
C
C 7 [THIS CARD GROUP IS NOT READ IF ISW(4)=0]
C
C IF ISW(4)=1 THE FOLLOWING FORMAT APPLIES
C 1 BACK F10.0 THE UNIFORM BACKGROUND
C CONCENTRATION TO APPLY TO
C ALL RECEPTORS.
C
C IF ISW(4)=2 THE FOLLOWING FORMAT APPLIES
C 1 NUMDIS 15 THE NUMBER OF RECEPTORS FOR
C WHICH A BACKGROUND VALUE IS
C TO BE SPECIFIED
C
C 2 [THIS CARD IS REPEATED NUMDIS TIMES]
-------
c
c
c
c
c
c
c
c
c
G
C
c
c
c
c
c
c
,c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
r-
w
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
XI
Yl
VAL
F10.0 THE X-COORDINATE OF THE
RECEPTOR FOR BACKGROUND
SPECIFICATION
F10.0 THE Y-COORDINATE OF THE
RECEPTOR FOR BACKGROUND
SPECIFICATION
F10.0 THE BACKGROUND CONCENTRATION
TO USE ,FOR THE IDENTIFIED
RECEPTOR.
[THIS CARD GROUP IS NOT READ IF ISW(5)=0]
1 NSCALE 15 THE NUMBER OF SOURCES TO
SCALE
2 [THIS CARD IS REPEATED NSCALE TIMES]
ISCALE 15 THE 'SOURCE IDENTIFICATION
NUMBER OF THE SOURCE TO BE
SCALED.
SCALE
F10.0
THE SCALING VALUE
10
[THIS CARD GROUP IS ONLY READ IF ISW(18), ISW(19) OR
ISW(20) =1]
1 NDYYR 15 THE NUMBER OF DAYS PER YEAR,
(NOTE THIS VALUE CAN BE ANY
DESIRED NUMBER LESS THAN
366. FOR EXAMPLE, IF THE
USER DESIRED A 30 DAY
AVERAGE VALUE, THE ANNUAL
CAN BE SET TO 30 DAYS WITH
THIS CARD.)
[THIS CARD GROUP IS ONLY READ IF ISW(21), ISW(22) OH
ISW(23) =1]
1 IOTHER 15 THE AVERAGING TIME TO BE
SPECIFIED BY THE USER
(NUMBER. OF HOURS)
11
[THIS CARD IS NOT READ IF ISW(8)=0]
ONEEX F10.0 THE EXCEEDANCE VALUE FOR
HOUR CONCENTRATIONS
1-
[THIS CARD IS NOT READ IF ISW(11)=0]
THREX F10.0 THE EXCEEDANCE VALUE FOR 3-
HOUR CONCENTRATIONS
[THIS CARD IS NOT READ .IF ISW(14)=0]
EIGEX F10.0 THE EXCEEDANCE VALUE FOR 8-
HOUR CONCENTRATIONS
[THIS CARD IS NOT READ IF ISW(17)=0]
TWYEX F10.0 THE EXCEEDANCE VALUE FOR 24-
HOUR CONCENTRATIONS
-------
c
C
C
C
C
C
C
C
C
C
C
[THIS CARD IS NOT READ IF ISW(20)=0]
ANNEX F10.0 THE EXCEEDANCE VALUE FOR
ANNUAL CONCENTRATIONS
[THIS CARD IS NOT READ IF ISW(23)=0]
OTHEX F10.0 THE EXCEEDANCE VALUE FOR
USER SPECIFIED AVERAGE '
CONCENTRATIONS
YDISC(400) ,ZZ(400)
,DMYISW) , (NUMSQ( 1) , DMYNSQ)
INTEGER TITLE, TYPE, HOUR, DATS
COMMON /DIM/ NSOURC,NGROUP,NXPNTS,NYPNTS;NXWYPT,NHOURS,NDAYS
1 , IOVRSN
COMMON /IO/ NOGO,IT1REC, IT2REC , NINTP , NOTTP , NINFL ( 3 ) , NQTFL ( 3 )
COMMON /FP1/ ROTATE, TK,G,ZR, HA, GAMMA1 , GAMMA2 , XRY , DECAY
1 , UTMX , UTMY
COMMON /MET/ THETA , UBAR , HM , TA , DPDZ , ISTBLE , P , SIGEPU , SIGAPU , SIGEPL ,
1SIGAPL, ALPHA, HOUR , NUMSQB ( 300 ), QB ( 300 ) ,TSB(300) ,VOLB(300)
COMMON /SOR/ NUMSQ(SOO) ,TYPE(300) ,Q(300) ,DX(300) ,DY(300) ,H(300) ,
1HS(300) ,TS(300) , VOL (300) , DTK (300) ,RDS(300) ,NS(300) , VS ( io , 300 ) ,
2FREQ(10,300) , JFLG(300) ,H10(300) ,H10SQ(300) ,VSQ(300) ,W(300)
3 , IFTAP ( 300 ) , CHGQ ('300 )
COMMON /RSI/ ISW(20) , TITLE (20) , KUNR ( 6 ) , KFNR ( 3 ) , DD1 (3) , 1 3, DATS (2). ,
1IH, IB, IPRT( 100) ,NLINES,KSW,NGROPS, IHOUR( 1464) , MUNT ( 4 ) , ITABLE , IPAGE
2,ND,LINE, IPC,IBLNK ;
COMMON /FILES/ IUNT,JUNT
' DIMENSION X(400) ,Y(400) ,XDISC(400)
EQUIVALENCE ( NSOURC , DMYNSO ) , (ISW( 1
IUNT=5
JUNT=6
C WRITE THE OUTPUT HEADER
WRITE( JUNT,9005)
9005 FORMAT ( ' I1 ,40X, ' POSTZ - VERSION 1.0' ,/,
1 40X,'A POST PROCESSOR FOR THE'
2 40X,'SHORTZ AIR QUALITY MODEL1
3 40X, ' JULY, 1986 ' )
C THE FOLLOWING STATEMENTS READ IN THE MAJOR "FIXED" DATA FROM
G' THE SHORTZ OUTPUT TAPE. THE VARIABLES, SUCH AS HOURLY CONCEN-
C TRATIONS ARE READ IN SUBROUTINE SEQANL .
CALL INPOUP (-1, DMYNSO, 8)
CALL INPOUP (-1 ,DMYISW,49)
C NXPNTS IS THE NUMBER OF X GRID POINTS, NYPNTS IS THE NUMBER OF
C Y GRID POINTS, NXWYPT IS THE NUMBER OF DISCRETE RECEPTORS.
C NSOURC IS THE NUMBER OF SOURCES, NHOURS IS THE NUMBER OF HOURS
C PER DAY OF METEOROLOGICAL DATA, NDAYS IS THE NUMBER OF DAYS OF
C METEOROLOGICAL DATA IN THE SHORTZ RUN.
J 1 =NXPNTS+NYPNTS+ 2 * NXWYPT
CALL INPOUP (1,ZZ,J1)
IF (NXPNTS. EQ.O. OR. NYPNTS. EQ.O) GOTO 10
C THE FOLLOWING ROUTINES STORE THE RECPTOR COORDINATES IN THE X ANY Y
C ARRAYS. DISCRETE COORDINATES ARE STORED IN XDISC AND YDISC ARRAYS.
DO 11 1=1, NXPNTS
X(I)=ZZ(I)
/
/
-------
11 CONTINUE
DO 12 I=1,NYPNTS
Y(I)=ZZ(NXPNTS+I)
12 CONTINUE
10 IF(NXWYPT.EQ.O) GOTO 13
DO 14 I=1,NXWYPT
XDISC(I)=ZZ(NXPNTS+NYPNTS+I)
YDISC(I)=ZZ(NXPNTS+NYPNTS+NXWYPT+I)
14 CONTINUE
C NXXYY IS THE TOTAL NUMBER OF RECEPTOR POINTS
13 NXXYY=NXPNTS*NYPNTS+NXWYPT
CALL INPOUP(1,ZZ,NXXYY)
C NIJ IS THE NUMBER OF GRID RECEPTOR POINTS
NIJ=NXPNTS*NYPNTS
CALL INPOUP(1,ROTATE,11)
CALL INPOUP(1,DMYNSQ,33*300)
CALL SEQANL(NSOURC,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,XDISC,YDISC,
1ZZ,TITLE,NUMSQ,NHOURS,NDAYS,IUNT,JUNT,NIJ)
WRITE(JUNT,9015)
9015 FORMAT('I1)
END
-------
SUBROUTINE SEQANL(NSOURC,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,XDISC,
1YDISC,ZZ,TITLE,NUMSQ,NHOURS,NDAYS,IUNT,JUNT,NIJ)
C
C SUBROUTINE SEQANL -- PART OF POSTZ
C THIS SUBROUTINE IS THE MAJOR PORTION OF THE POSTZ PROGRAM. ALL
C THE INPUT AND OUTPUT IS MANAGED FROM THIS ROUTINE. SUBROUTINE
C SORT IS THE MAJOR ANALYTICAL TOOL, WHILE SEVERAL, ROUTINES ARE
C CALLED TO WRITE OUTPUT TABLES
C
INTEGER HOUR,PGCT,PTITLE(20),TITLE(20)
COMMON ONE(300,5),
*IONEDY(300,5),IEIGSK(8),ITWYSK(24),
2IONETM(300,5) , ONE50 ( 50 ) -, ION50D( 50 ) ,ION50T(50) ,ION50R(50) ,
30NEE(100),IONED(100),IONET(100),IONER{100),THR(300,5),
4ITHRDY(300,5),ITHRTM(300,5),THR50(50),ITH50D(50),ITH50T(50), ,
, THRE(IOO),ITHED(IOO),ITHET(IOO),ITHER(100),
IEIGDY(300,5),IEIGTM(300,5),EIG50(50),IEI50D(50),
,IEI50R{50) , EIGE(IOO') , IEIED(100) ,IEIET(100) ,IEIER(100)
ITWYDY(300,5) ,ITWYTM(300,5) ,TWY50(50) , ITW50D'(50) ,
, ITW50R(50),TWYE(100),ITWED(100),ITWET(IOO),ITWER(IOO)
5),IANNTM(3QO,5),ANN50(50),IAN50D(50),
,ANNE(100),lANED(lOO),IANET(100),IANER(100)
,.5) ,IOTHTM(300,5) ,OTH50(.50) ,IOT50D(50) ',
, OTHE(IOO),IQTED(IOO),IDTET(IOO), IpTER(lOO)
5ITH50R(50)
6EIG(300,5)
7IEI50T(50)
8TWY(300,5)
9ITW50T(50)
AANN(300/5)
BIAN50T(50)
COTH(300,5)
DIOT50T(50)
IANNDY(300
IAN50R(50)
IOTHDY(300,
IOT50R(50)
ETHRC(300),ANNC(300),
FEIGC(300),TWYC(300),OTHC(300),ISW(25),IELIM(300),XDISCR(300),
GYDISCR(300),ISELIM(300),BCKG(300),SCLE(300),CON(300),CONT(300),
HONEC(SOO),BACKOT(300) . , .
COMMON /MET/ THETA , UBAR', HM , TA , DPDZ , ISTBLE , P ,SIGEPU , SIGAPU , SIGEPL ,
1SIGAPL,ALPHA,HOUR,NUMSQB(300),QB(300),TSB(300),VOLB(300) • .
DIMENSION X(300),Y(300),XDISC(300),YDISC(300),ZZ(300),NUMSQ(300).
1XX(300),YY(300),ISHOLD(8),HOLD(300,24)
/
USED
DATA HOLD/7200*0
C CLEAR THE SCLE ARRAY
C SOURCE
DO 7 1=1,NSOURC
SCLE(I)=1.
7 CONTINUE
'PGCT=1
LNCT=4
C WRITE TITLE FOR OUTPUT
WRITE(JUNT,9020) PGCT .
9020 FORMAT(IX,'1',105X,'POSTZ -
WRITE(JUNT,9021)
9021 FORMAT(43X,'POSTZ - A POST PROCESSOR
C READ THE TITLE CARD FOR THE POSTZ RUN
READ(IUNT,9006)
9006 FORMAT(20A4)
C READ THE SWITCHES
READ(IUNT,9007)
9007 FORMAT(25I2)
C READ THE JULIAN DAY
C SHORTZ RUN METEOROLOGICAL DATA
READ(IUNT,9008) JULST,IYR
9008 FORMAT(I5,15)
TO SCALE CONCENTRATIONS FOR A PARTICULAR
PAGE NO.
(PTITLE(J),J=l,20)
FOR THE SHORTZ MODEL1,/)
(ISW(I),1=1,25)
(JULST) AND THE YEAR FOR THE START OF THE
-------
C WRITE THE SWITCH INFORMATION IN THE OUTPUT FILE
WRITE(JUNT,9022) (PTITLE(I),1=1,20),(TITLE(I),1=1,20),
1(ISW(J),J=l,11)
9022 FORMAT(20X,'POSTZ RUN TITLE: ',20A4,/,20X,'SHORTZ RUN TITLE: ',
120A4,//,20X,'ISW(l) RESTRICT TIME LIMITS (1=YES,0=NO)',40X,
2I1,/,20X,'ISW(2) LIMIT RECEPTORS TO ANALYZE (1=YES,0=NO)',
334X,I1,/,20X, 'ISW(3) LIMIT SOURCES ,TO ANALYZE (1=YES,0=NO) ' ,
436X,I1,/,20X,'ISW(4) SPECIFY BACKGROUND CONCENTRATIONS (0=NO,1=UNI
5FORM,2=BY RECEPTOR) ',, 9X, II ,/, 20X, ' ISW( 5 ) SCALE CONCENTRATIONS',
6' FOR SPECIFIC SOURCES (1=YES,0=NO)',19X,II,//,
715X,'1-HOUR AVERAGE ANALYSIS:',/,20X,'ISW(6) HIGH-5 TABLE1,
8' PREPARED (1=YES,0=NO)',39X,II,/,20X,'ISW(7) TOP 50 TABLE1,
9' PREPARED (1=YES,0=NO)',39X,II,/,20X,'ISW(8) EXCEEDANCE TABLE1,
A1 PREPARED ( 1=YES , 0=NO )', 35X , II ,// ,.15X ,' 3-HOUR AVERAGE ANALYSIS:1,
B/, 20X, ' ISW(9) HIGH-5 TABLE PREPARED ( 1=YES , 0=NO ) ' , 39X , 11 , / ,.20X ,
C'ISW(IO) TOP 50 TABLE PREPARED (1=YES,0=NO)',38X,II,/,20X,
D'ISW(ll) EXCEEDANCE TABLE PREPARED (1=YES,0=NO)',34X,I1,/)
; WRITE(JUNT,9023) (ISW(J),J=l2,25)
9023 FORMAT(15X,'8-HOUR AVERAGE ANALYSIS:',/,20X,
,1'ISW(12) HIGH-5 TABLE PREPARED (1=YES,0=NO)',38X,II,/,20X,_
2'ISW(13) TOP-50 TABLE PREPARED (1=YES,0=NO)',38X,I1,/,20X,
3'ISW(14) EXCEEDANCE TABLE PREPARED (1=YES,0=NO)',34X,II,//,15X,
4'24-HOUR AVERAGE ANALYSIS:',/,20X,
5'ISW(15) HIGH-5 TABLE PREPARED ( 1=YES , 0=NO ) ' , 38X ,,11 , / , 20X , '
6'ISW(16) TOP-50 TABLE PREPARED (1=YES,0=NO)',38X,I1,/,20X,
7'ISW(17) EXCEEDANCE TABLE PREPARED ( 1=YES , 0=NO ) '', 34X , 11 , // , 1 5X ,
8'ANNUAL 'AVERAGE ANALYSIS:',/,2OX,
9'ISW(18) HIGH-5 TABLE PREPARED ( 1=YES , 0=NO ) ,' , 38X , 11 , / , 20X ,
A'ISW(19) TOP-50 TABLE PREPARED (1=YES,0=NO)',38X,I1,/,20X,
B'ISW{20) EXCEEDANCE TABLE PREPARED (1=YES,0=NO)',34X,I1,//,15X,
C'USER-SPECIFIED AVERAGING TIME ANALYSIS',/,20X,
D'ISW(21) HIGH-5 TABLE PREPARED (1=YES,0=NO)',38X,I1,/,20X,
E'ISW(22) TOP-50 TABLE PREPARED (1=YES,0=NO)',38X,I1,/,20X,
F'ISW(23) EXCEEDANCE TABLE PREPARED (1=YES,0=NO)',34X,II,//,20X,
G'ISW(24) BLOCK OR RUNNING AVERAGES (0=BLOCK,1-RUNNING)',27X,I1,
H/,20X,'ISW(25) CALMS POLICY (1=YES,0=NO)',47X,I1)
IYRST=IYR.
IDYST=JULST
IYREND=IYR
IDYEND=JULST+NDAYS-1
401 ITEST=IYREND/4
ITEST2=ITEST*4
IMAX=365
IF(ITEST.EQ.ITEST2) IMAX=366
IF(IDYEND.LE.IMAX) GOTO 400
IYREND=IYREND+1
IDYEND=IDYEND-IMAX
GOTO 401
400 IF(ISW(1).NE.1) GOTO 10
C THIS ROUTINE ALLOWS THE USER TO SELECT A PORTION OF THE TIME
C PERIOD IN THE SEQUETIAL SHORTZ TAPE FOR ANALYSIS. IT BASICALLY
C RESTRICTS THE ANALYSIS TO A STARTING JULIAN DAY AND YEAR AND AN
C ENDING JULIAN DAY AND YEAR.
READ(IUNT,9009) IDYST,IYRST,IDYEND,IYREND
9009 FORMAT(4I5)
-------
IF(IYR-IYRST) 1,2,3
3 IYRST=IYR
2 IF(JULST.GT.JDYST) IDYST=JULST
1 IF(IYR-IYREND) 10,5,6
6 IYREND=IYR
5 IF(JULST.GT.IDYEND) IDYEND=JULST
10 WRITE(JUNT,9024) IDYST,IYRST,IDYEND,IYREND
9024 FORMAT(//,20X,'JULIAN DAY FOR START OF.ANALYSIS',46X,13,/,
120X,'YEAR FOR START OF ANALYSIS',SOX,15,/,20X,'JULIAN DAY ',
2'FOR END OF ANALYSIS',48X,13,/,20X,'YEAR FOR END OF ANALYSIS',
352X,I5)
IF(ISW(2) .EQ.'O) GOTO 20
C THIS ROUTINE ALLOWS THE USER TO RESTICT THE ANALYSIS TO ONLY CERTAIN
C SELECTED RECEPTORS. AN ARRAY CALLED IELIM IS INITALIZED AT 1 TO
C INDICATE ALL RECEPTORS ARE TO BE ELIMINATED. FOR THOSE TO BE
C THE ARRAY ELEMENT IS LATER CHANGES TO 0
DO 11 I=1,NXXYY
IELIM(I)=1
11 CONTINUE
READ(IUNT,9011) NXINCL
C IF THE NUMBER OF'X GRID ELEMENTS TO INCLUDE IS 0, THERE ARE NO
C GRID RECEPTORS TO INCLUDE
IF(NXINCL.EQ.O) GOTO 12
C READ THE X AND Y RECEPTORS TO INCLUDE
READ(IUNT,9012) (XX(I),1=1,NXINCL)
READ(IUNT,9011) NYINCL ' . •
READ(IUNT,9012) (YY(I),1=1,NYINCL)
9011 FORMAT(15)
C THIS ROUTINE CHECKS TO ENSURE THE RECEPTORS THE USER HAS SPECIFIED
C WERE INDEED IN THE ORIGINAL SHORTZ RUN
DO 13 11=1,NXINCL
DO 13 JJ=1,NYINCL
IGO=0
DO 14 I=1,NXPNTS
DO 14 J=1,NYPNTS
IF(X(I).NE.XX(II).OR.Y(J).NE.YY(JJ)) GOTO 14
IGO=1
K=(J-1)*NXPNTS+I
IELIM(K)=0
GOTO 13
14 CONTINUE
IF(IGO.EQ.O) CALL EROUT(2)
13 CONTINUE
12 READ(IUNT,9011) NXYINC
C IF NEITHER GRID NOR DISCRETE RECEPTORS ARE SELECTED, AN ERROR HAS
C DETECTED IN THE INPUT STREAM
IF(NXINCL.EQ.O.AND.NXYINC.EQ.O) CALL EROUT(1)
IF(NXYINC.EQ.O) GOTO 25
C THE DISCRETE RECEPTORS ARE READ
READ(IUNT,9012) (XDISCR(I),1=1,NXYINC)
9012 FORMAT(8F10.0)
READ(IUNT,9012) (YDISCR(I),1=1,NXYINC)
C THIS ROUTINE CHECKS TO ENSURE THAT SELECTED DISCRETE RECEPTORS WERE
C INDEED IN THE ORIGINAL SHORTZ RUN
DO 17 1=1,NXYINC
-------
IGO=0
DO 18 J=1,NXWYPT
IF(XDISC(J).NE.XDISCR(I),OR.YDISC(J).NE.YDISCR(I)) GOTO 18
IGO = 1
IELIM(NIJ+J)=0
18 CONTINUE
IF(IGO.NE.l) CALL EROUT(3)
17 CONTINUE
C PGCT IS A PAGE COUNT INDEX USED TO NUMBER PAGES
25 PGCT=PGCT+1
WRITE(JUNT,9020) PGCT
WRITE(JUNT,9025)
9025 FORMAT(45X,'RESTRICTING ANALYSIS TO CERTAIN RECEPTORS',/,
148X,'THE FOLLOWING RECEPTORS WERE INCLUDED ',//, IX, ''
3 ----------------- _ --------------------- . )
NOPNTS=0
IF(NXPNTS.EQ.O.AND.NYPNTS.EQ.O) GOTO 26
DO 27 I=1,NYPNTS
DO 27 J=1,NXPNTS '
K=(t-l) *NXPNTS+J
IF(IELIM(K) .EQ.l) GOTO 27
NOPNTS=NOPNTS+1
XX(NOPNTS)=X( J)
YY(NOPNTS)=Y(I)
27 CONTINUE
26 DO 28 I=1,NXWYPT
K=NIJ+I
IP(IELIM(K) .EQ.l) GOTO 28
NOPNTS=NOPNTS+1
XX(NOPNTS)=XDISC(I)
YY(NOPNTS)=YDISC( I)
28 CONTINUE
WRITE (JUNT, 9026) (XX(K) ,YY(K) ,K=1,NOPNTS)
9026 FORMAT ( 10X,I(I/F10.0,'/1,F10.0,1)',5X,1(I;F10.0,I,I/F10.0
20 IF(ISW(3) .EQ.O) GOTO 30
C THIS ROUTINE ALLOWS THE USER TO ELIMATE CERTAIN SOURCES FORM
C INCLUSION IN THE POSTZ RUN
C THE PARAMETER NSELIM IS THE NUMBER OF SOURCES TO INCLUDED IN
C THE POSTZ RUN
READ (IUNT, 9011) NSELIM
DO 21 I=1,NSOURC
ISELIM(I)=1
21 CONTINUE
DO 22 1=1, NSELIM
READ ( IUNT, 9011) K
IGO=0
C THIS ROUTINE CHECKS TO ENSURE THE SOURCE TO BE INCLUDED WAS INDEED
C IN THE ORIGICAL SHORTZ RUN
DO 23 J=1,NSOURC
IF(NUMSQ( J) .NE.K) GOTO 23
IGO=1
ISELIM( J)=0
23 CONTINUE
-------
IF(IGO.EQ.O) CALL EROUT ( 4 }
22 CONTINUE
PGCT=PGCT+1
WRITE(JUNT,9020) PGCT
WRITE( JUNT,9027)
9027 FORMAT (4 IX, 'THE ANALYSIS HAS BEEN CONFINED TO CERTAIN SOURCES'
135X, 'THE FOLLOWING SOURCES (BY IDENTIFICATION NUMBER) WERE ',
2'INCLUDED1 ,//, ' ----------------------------------- •
NOSO=0
DO 24 I=1,NSOURC
IF(ISELIM(I) .EQ.l) GOTO 24
NOSO=NOSO-I-1
ISHOLD ( NOSO ) =NUMSQ ( I )
IF(NOSO.LT.S) GOTO 24
WRITE (JUNT, 9028) ( ISHOLD ( J ), J=l , NOSO )
9028 FORMAT( 10X ,,8 ( 15 , 6X) )
NOSO=0
24 CONTINUE
WRITE (JUNT ,9028) ( ISHOLD ( J ), J=l, NOSO )
30 IF(ISW(4)-1) 40,31,32
C THIS ROUTINE ALLOW THE SPECIFICATION OF BACKGROUND VALUES TO BE
C , ADDED TO THE CONCENTRATIONS ON THE TAPE. BACKGROUND CAN BE
C SPECIFIED IN TWO WAYS:. A UNIFORM BACKGROUND IS USED FOR ALL RECEPTORS
C IN ALL ANALYSES, AND SEPARATE BACKGROUND CONCENTRATIONS ARE USED FOR
C INDIVIDUAL RECEPTOR
31 READ (IUNT, 9013) BACK
9013 FORMAT(FIO.O) . ' .
DO 33 I=1,NXXYY
BCKG(I)=BACK • ,
33 CONTINUE
PGCT=PGCT+1
WRITE( JUNT, 9020) PGCT
WRITE( JUNT, 9029) BACK
9029 FORMAT(42X, ' A UNIFORM BACKGROUND CONCENTRATION OF ',F10.2,/,
137X, 'MICROGRAMS PER CUBIC METER. HAS BEEN ADDED TO ALL RECEPTORS')
GOTO 40
32 PGCT=PGCT+1
C THIS ROUTINE ALLOW THE DIFFERENT BACKGROUND CONCENTRATIONS TO 3E USED
C FOR EACH RECEPTOR
WRITE (JUNT, 9020) PGCT
WRITE( JUNT, 9030)
9030 FORMAT ( 4 3X, 'SEPARATE BACKGROUND CONCENTRATIONS HAVE BEEN ',/,
153X, 'SPECIFIED FOR EACH RECEPTOR ',/, 55X ,'( X-COORD , Y-COORD ) VALUE1,
21 \ ___________ _ _________________________________________________
/ , ---
o __________________________________________________________________
4' )
34 READ(IUNT,9011) NUMDIS
C THE NUMBER OF RECEPTORS TO HAVE A BACKGROUND SPECIFICATION IS
C READ AS NUMDIS
IF (NUMDIS) 49,49,45
45 NBACK=0
DO 46 1=1, NUMDIS
C THE COORDINATES FOR EACH RECEPTORS TO HAVE A BACKGROUND SPECIFIED
C ARE READ ALONG WITH THE BACKGROUND VALUE
-------
READ(IUNT,9012) X1,Y1,VAL
IGO=0
DO 35 II=1,NXPNTS
DO 35 J=1,NYPNTS
IF(X1.NE.X(II).OR.Yl.NE.Y(J)) GOTO 35
IGO=1
K=(J-1)*NXPNTS+II
NBACK=NBACK+1
C XX, YY AND BACKOT ARE USED STRICTLY FOR OUTPUT
XX(NBACK)=X(II)
YY(NBACK)=Y(J)
BAGKOT(NBACK)=VAL
BCKG(K)=VAL
GOTO 46
35 CONTINUE
DO 47 J=1,NXWYPT
IF(X1.NE.XDISC(J).OR.Yl.NE.YDISC(J)) GOTO 47
IGO=1
K=NIJ+J
NBACK=NBACK+1
XX(NBACK)=XDISC(J)
YY(NBACK)=YDISC(J)
BACKOT(NBACK)=VAL
BCKG(K)=VAL
GOTO 46
47 CONTINUE
IF(IGO.EQ.O) CALL EROUT(12)
46 CONTINUE
49 NHOLD=0
IF(NBACK.GT.200) NHOLD=NBACK-200
NBACK=NBACK-NHOLD
37 WRITE(JUNT,9031) (XX(I),YY(I),BACKOT(I),1=1,NBACK)
9031 FORMAT(2X,'(',F9.0,',',F9.0,')',F9.0,2X,'(',F9.0,',',F9.0,')',
1F9.0.2X, ' ( ' ,F9.0, ' , ' ,F9.0, ' ) ' ,F9.0,2X, ' ( ' ,F9.0, ' , ' ,F9.0, ' ) ' ,
2F9.0)
IF(NHOLD.EQ.O) GOTO 40
PGCT=PGCT+1
WRITE(JUNT,9020) PGCT
NBACK=NBACK-l-NHOLD
WRITE(JUNT,9031) (XX(I),YY(I),BACKOT(I),1=201,NBACK)
40 IF(ISW(5).EQ.O) GOTO 50
C THIS ROUTINE ALLOWS THE CONTRIBUTION FROM ANY INDIVIDUAL SOURCE TO BE
C SCALED BY A PARTICULAR VALUE, SCALE
READ(IUNT,9011) NSCALE
DO 41 1=1,NSCALE
READ(IUNT,9014) ISCALE,SCALE
9014 FORMAT(15,F10.0)
IGO=0
DO 42 J=1,NSOURC
IF(NUMSQ(J).NE.ISCALE) GOTO 42
IGO=1
SCLE(J)=SCALE
42 CONTINUE
IF(IGO.EQ.O) CALL EROUT(5)
41 CONTINUE
-------
PGCT=PGCT+1
WRITE(JUNT,9020) PGCT
WRITE(JUNT,9032)
9032 FORMAT(44X,'THE CONCENTRATIONS FROM CERTAIN SOURCES WERE',/,SOX,
I1SCALED BY THE FOLLOWING AMOUNTS',/,47X,'(IDENTIFICATION NUMBER/SC
2ALING VALUE) ',//,'
O I \
j — — — — — — — — — — — — — — — — -.____ — — — — — — ^
WRITE(JUNT,9033) (NUMSQ(I),SCLE(I),1=1,NSOURC)
9033 FORMAT(5X, ' ( ' ,15, ' , ' ,F7.2, ' ) ' ,5X, ' ( ' ,15, ' , ' ,F7.2, ' ) ' ,
15X,'(',15,',',F7.2,')',5X,'(',15,',',F7.2,')',
, 25X,'(',15,',',F7.2,')',5X,'(',15,',',F7.2,')')
C
C THE FOLLOWING ROUTINES CHECK TO ENSURE THERE IS ENOUGH MET DATA
C IN THE ORIGINAL SHORTZ RUN TO COMPLETE THE ANLAYSES REQUESTED BY
C THE USER
50 IF(ISW(9).EQ.O.AND.ISW(10).EQ.O.AND.ISW(11).EQ.0) GOTO 60
IF(NHOURS.LT.3) CALL EROUT(6)
60 IF(ISW(12).EQ.O.AND.ISW(13).EQ.O.AND.ISW(14).EQ.O) GOTO 70 ' '
IF(NHOURS.LT.S) CALL EROUT(7)
70 IF(ISW(15).EQ.O.AND.ISW(16).EQ.O.AND.ISW(17).EQ.O) GOTO 80
IF(NHOURS.LT.24) CALL EROUT(8)
80 IF(ISW(18).EQ.O.AND.ISW(19).EQ.O.AND.ISW(20).EQ.O) GOTO 90
C IF AN ANNUAL CONCENTRATION IS TO BE COMPUTED, THIS ROUTINE ASKS
C THE NUMBER OF DAYS PER YEAR IN THE ANNUAL RUN
READ(IUNT,9011) NDYYR
NHRYR=NDYYR*24
IF(NDAYS.LT.NDYYR) CALL EROUT(9)
90 IF(ISW(21).EQ.O.AND.ISW(22).EQ.O.AND.ISW(23).EQ.O) GOTO 100
C OTHER IS USED WHEN THE USER WANTS TO SPEICFIY AN AVERAGING OTHER
C THAN THE STANDARD VALUES (E.G. 1-HOUR, 3-HOUR, 8-HOUR, 24-HOUR, AND
C ANNUAL.
READ(IUNT,9011) IOTHER
IF(NHOURS.LT.IOTHER) CALL EROUT(IO)
100 JDAY=0
JULD=JULST-1
C IF AN EXCEEDANCE TABLE IS TO BE PREPARED, THE CRITERIA FOR THE
C EXCEEDANCE IS SPECIFIED IN THE NEXT LINES
IF(ISW(8).EQ.1) READ(IUNT,9012) ONEEX
IF(ISW(11).EQ.1) READ(IUNT,9012) THREX
IF(ISW(14).EQ.1) READ(IUNT,9012) EIGEX
IF(ISW(17).EQ.1) READ(IUNT,9012) TWYEX
IP(ISW{20).EQ.l) READ(IUNT,9012) ANNEX
IP(ISW(23).EQ.l) READ(IUNT,9012) OTHEX
C
C MAIN LOOP - OVER THE TOTAL NUMBER OF DAYS
C
DO 1000 IDAY=1,NDAYS
JDAY=JDAY+1
JULD=JULD+1
C TEAT TO DETERMINE IF LEAP YEAR
ITEST=IYR/4
IF(ITEST-IYR/4) 110,111,110
110 IF(JULD.GT.365) GOTO 112
GOTO 120
111 IF(JULD.GT.366) GOTO 112
-------
GOTO 120
112 IYR=IYR+1
JULD=1
C TEST TO DETERMINE IF BEYOND END YEAR OR DAY, OR BEFORE START
C YEAR OR DAY.
120 IF(IYR-IYREND) 130,121,2000
121 IF(JULD.GT.IDYEND) GOTO 2000
C
C LOOP OVER NUMBER OF HOURS PER DAY
C
130 DO 900 IHR=1,NHOURS
C READ THE INPUT TAPE (FROM THE SHORTZ RUN) TO OBTAIN THE MET DATA
CALL INPOUP(1,THETA,1213)
C ISKIP IS AN PART OF THE CALMS POLICY IMPLEMENTATION
ISKIP=0
C THE FOLLOWING ARE COUNTERS FOR THE NUMBER OF "VALID" (AS DEFINED
C BY THE CALMS POLICY) AND TOTAL HOURS PROCESSED
I3CT=I3CT+1
IF(I3CT.LE.3) GOTO 132
I3CT=1
132 I8CT=I8CT+1
IF(I8CT.GT.8)
IF(UBAR.GE.1.)
GOTO 134
133 'I8CT=1
I8VAL=0
IF(UBAR.GE.l.)
134 I24CT=I24CT+1
IF(I24CT.GT.24
IF(UBAR.GE.1.)
GOTO 136
135 I24CT=1
I24VAL=0
IF(UBAR.GE.1
136 IANNCT=IANNCT+1
IF(IANNCT.GT.NHRYR) GOTO 138
IF(UBAR.GE.1.) IANVAL=IANVAL+1
GOTO 139
138 IANNCT=1
IANVAL=0
IF(UBAR.GE.1.) IANVAL=1
139 IOTHCT=IOTHCT+1
IF(IOTHCT.GT.IOTHER) GOTO 137
IFJUBAR.GE.1.) IOTVAL=IOTVAL+1
GOTO 140
137 IOTHCT=1
IOTVAL=0
IF(UBAR.GE
140 IF(ISW(25)
DO 151 1=1
CONT(I)=0.
151 CONTINUE
GOTO 133
I8VAL=I8VAL+1
I8VAL=1
GOTO 135
I24VAL=I24VAL+1
) I24VAL=1
1.) IOTVAL=1
EQ.1.AND.UBAR.LT.l.) ISKIP=1
NXXYY
C LOOP ON SOURCES
-------
c
DO 180 ISRC=1,NSOURC
C READ.THE CONCENTRATION DATA FROM THE TAPE
CALL INPOUP(1,CON,NXXYY)
IF(IYR-IYRST) 180,150,160
150 IF(JULD.LT.IDYST) GOTO 180
160 IF(ISELIM(ISRC).EQ.1) GOTO 180
IF(ISKIP.EQ.l) GOTO 180
DO 170 I=1,NXXYY
IF(ISELIM(ISRC).EQ.l) GOTO 170
IF(IELIM(I).EQ.1) GOTO 170
CONT(I)=CONT(I)+CON(I)* SCLE(ISRC)
170 CONTINUE
180 CONTINUE
IF(IYR-IYRST) 900,181,182
181 IF(JULD.LT.IDYST) GOTO 900
182 IF(ISKIP.EQ.l) GOTO 191
DO 190 I=1,NXXYY
C CONT IS THE TOTAL CONCENTRATION FROM ALL SOURCES INCLUDING BACKGROUND
CONT(I)=CONT(I)+BCKG(I)
190 CONTINUE
C
C THE MAIN LOOP ON RECEPTORS
•C '
191 DO 300 IR=1,NXXYY
' IF(IELIM(IR).EQ.1) GOTO 300
DO 290 1=1,23
J=25-I
C HOLD IS USED FOR RUNNING AVERAGES. THE PREVIOUS 24-HOURS OF
C TOTAL CONCENTRATIONS BY RECEPTOR ARE STORED IN HOLD
HOLD(IR„J)=HOLD(IR,J-1)
290 CONTINUE
HOLD(IR,1)=CONT(IR)
C THE SORT SUBROUTINE IS THE MAJOR ANALYSIS TOOL FOR THE POSTZ
C PROGRAM. FOR EACH AVERAGIN TIME, SORT KEEPS TRACK OF THE
C.HIGH-5, TOP-50 AND EXCEEDANCE TABLE INFORMATION. IN THIS FIRST
C CALL, THE THE HIGH-5 FOR THE 1-HOUR ANALYSIS ARE STORED IN THE
C ARRAY ONE, WHILE THE TOP 50 FOR ALL RECEPTORS ARE STORED IN
C ONE50 AND THE EXCEEDANCE TABEL DATA IS STORED IN ONEE.
C SUBSEQUENT CALLS USE SIMILAR NOMENCLATURE
. CALL SORT(1,ISW(6),ISW(7),ISW{8),CONT,IR,ONEC,HOLD,1,1,JULD,
1IHR,ONE,IONEDY,IONETM,ONE50,ION50D,ION50T,ION50R,ONEE,IONED,
2IONET,IONER,ONEEX,ISW(24))
CALL SORT(3,ISW(<9) , ISW( 10) , ISW( 11) , CONT , IR , THRC , HOLD , I3CT,3, JULD,
1IHR.THR,ITHRDY,ITHRTM,THR50,ITH50D,ITH50T,ITH50R,THRE,ITHED,
2ITHET,ITHER,THREX,ISW(24))
IT=0
C THIS ROUTINE IS PART OF THE CALMS POLICY IMPLEMENTATION
DO 200 1=1,7
IEIGSK(I)=IEIGSK(I+1)
IT=IT+IEIGSK(I)
200 CONTINUE
IEIGSK(8)=ISKIP
IT=8-IT-ISKIP
IF(IT.LT.6) IT=6
-------
CALL SORT(8,ISW(12),ISW(13),ISW(14),CONT,IR,EIGC,HOLD,I8CT,IT,
1JULD,IHR,EIG,IEIGDY, IEIGTM , EIG50 , IEI50D , IEI50T , IEI50R,EIGE,
2IEIED,IEIET,IEIER,EIGEX,ISW(24))
C THIS ROUTINE IS PART OF THE CALMS POLICY IMPLEMENTATION
IT=0
DO 210 1=1,23
ITWYSK(I)=ITWYSK(1+1)
IT=IT+ITWYSK(I)
210 CONTINUE
ITWYSK(24)=ISKIP
IT=24-IT-ISKIP
IF(IT.LT.ia) IT=18
CALL SORT(24,ISW(15),ISW(16),ISW(17),CONT,IR,TWYC,HOLD,I24CT,IT,
1JULD,IHR,TWY,ITWYDY,ITWYTM,TWY50,ITW50D,ITW50T,ITW50R,TWYE,ITWED,
2ITWET,ITWER,TWYEX,ISW(24))
IF(ISW{25).NE.l) IANVAL=IANNCT
CALL SORT(NHRYR,ISW(18),ISW(19),ISW(20),CONT,IR,ANNC,HOLD,IANNCT,
1IANVAL,JULD,IHR,ANN,IANNDY,IANNTM,ANN50,IAN50D,IAN50T,IAN50R,ANNE,
2IANED,IANET,IANER,ANNEX,0)
CALL SORT(IOTHER,ISW(21),ISW(22),ISW(23),CONT,IR,OTHC,HOLD,
1IOTHCT,IOTHER,JULD,IHR,OTH,IOTRDY,IOTRTM,OTH50,IOT50D,IOT50T,
2IOT50R,OTHE,IOTED,IOTET,IOTER,OTHEX,ISW(24))
300 CONTINUE
900 CONTINUE
•1000 CONTINUE
C
c THE'MAJOR LOOPS HAVE BEEN CLOSED, THE NEXT CALLS WRITE THE
C OUTPUT TABLES
c
2000 IF(ISW(6).EQ.l) CALL H5OUT(1,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,
,1XDISC,YDISC,ZZ,PGCT,LNCT,ONE,IONEDY,IONETM,IELIM)
IF(ISW(7).EQ.l) CALL T500UT(1,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,ONE50,ION50D,ION50T,ION50R)
IF(ISW(8).EQ.1) CALL EXOUT(1,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,ONEE,IONED,IONET,IONER,ONEEX)
IF(ISW(9).EQ.1) CALL H50UT{3,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,THR,ITHRDY,ITHRTM,IELIM)
IF(ISW(10).EQ.1) CALL T500UT(3,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,THR50,ITH50D,ITH50T,ITH50R)
IF(ISW(11).EQ.1) CALL EXOUT(3,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,THRE,ITHED,ITHET,ITHER,THREX)
IF(ISW(12).EQ.l) CALL H50UT(8,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,EIG,IEIGDY,IEIGTM,IELIM)
IF(ISW(13).EQ.l) CALL T500UT(8,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,EIG50,IEI50D,IEI50T,IEI50R)
IF(ISW(14).EQ.l) CALL EXOUT(8,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,EIGE,IEIED,IEIET,IEIER,EIGEX)
IF(ISW(15).EQ.l) CALL H50UT(24,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,TWY,ITWYDY,ITWYTM,IELIM)
IF(ISW(16).EQ.1) CALL T500UT(24,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,TWY50,ITW50D,ITW50T,ITW50R)
IF(ISW(17).EQ.l) CALL EXOUT(24,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,TWYE,ITWED,ITWET,ITWER,TWYEX)
IF(ISW(18).EQ.1) CALL H50UT(NHRYR,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,ANN,IANNDY,IANNTM,IELIM)
-------
IF(ISW(19).EQ.l) CALL T500UT(NHRYR,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,ANN50,IAN50D,IAN50T,IAN50R)
IF(ISW(20).EQ.l) CALL EXOUT(NHRYR,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,ANNE,IANED,IANET,IANER,ANNEX)
IF(ISW(21).EQ.1) CALL H50UT(IOTHER,NXPNTS,NYPNTS,NXWYPT,NXXYY,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,OTH,IOTHDY,IOTHTM,IELIM)
IF(ISW(22).EQ.l) CALL T500UT(IOTHER,NXPNTS,NYPNTS,X,
1Y,XDISC,YDISC,ZZ,PGCT,LNCT,OTH50,IOT50D,IOT50T,IOT50R)
IP(ISW(23).EQ.l) CALL EXOUT(IOTHER,NXPNTS,NYPNTS,X,Y,
1XDISC,YDISC,ZZ,PGCT,LNCT,OTHE,, IOTED , IOTET , IQTER, OTHEX)
RETURN
END
-------
SUBROUTINE SORT(N,II,12,13,CONT,IR,CAL,HOLD,ICT,NVAL,JULD,IHR,H5,
1IH5DY, IH5TM,T50, IT50DY, IT50TM, IR50,E, IEDY, IETM, IRE, EX, IBR)
C
C SUBROUTINE SORT — PART OF POSTZ
C THIS SUBROUTINE KEEPS 'TRACK OF TOP 5, TOP 50 AND EXCEEDANCE
C TABLE DATA FOR EACH AVERAGING TIME. SUBROUTINE SEQANL CALLS SORT
C FOR EACH OF THE AVERAGING TIMES SELECTED BY THE USER
C
DIMENSION CONT(300),CAL(300),HOLD(300,24),H5(300,5),IH5DY(300,5),
1IH5TM(300,5) , T50 ( 50 )., IT50DY( 50 ) ,IT50TM(50) ,E(100) , IEDY( 100) ,
2IETM(100),IRE(100),IR50{50)
C IF THIS AVERAGING TIME HAS NOT ,BEEN SELECTED, RETURN TO SEQANL
IF(I1.EQ.O.AND.I2.EQ.O.AND.I'3.EQ.O) RETURN
C TEST TO SEE WHETER RUNNING OR BLOCK AVERAGES
IF(IBR.EQ.l) GOTO 10
C ROUTINE TO COMPUTE BLOCK AVERAGES
CAL(IR)=CAL(IR)+CONT(IR)
C FOR BLOCK AVERAGES,' RETURN IF THE COUNTER DOES NOT EQUAL THE
C AVERAGING TIME
IF(ICT;LT.N) RETURN -
GOTO 30
C ROUTINE TO COMPUTE RUNNING AVERAGES
10 Cl=0.
DO 20 1=1,N
C1=C1+HOLD(IR,I)
20 CONTINUE •••'•'
C1=C1/NVAL
GOTO 36 ' ' '
C Cl IS THE AVERAGE CONCENTRAION TO BE EXAMINED TO DETERMINE
C IF IT NEEDS TO BE ADDED TO THE TOP 5. TABLE FOR THE RECEPTOR, OR
C TO THE TOP-50 TABLE FOR ALL RECEPTORS
30 C1=CAL(IR)/NVAL
36 IF(Il.NE.l) GOTO 40
C TEST FOR THE HIGH-5 TABLE
DO 31 1=1,5
IF(C1.GT.H5(IR,I)) GOTO 32
31 CONTINUE
GOTO 40
32 IF(I.EQ.S) GOTO 33
DO 34 J=I,4
K=5+I-J
C MAKE ROOM IN THE ARRAYS FOR VALUE, DAY AND TIME
H5(IR,K)=H5(IR,K-1)
IH5DY(IR,K)=IH5DY(IR,K-1)
IH5TM(IR,K)=IH5TM(IR,K-1)
34 CONTINUE
33 H5(IR,I)=C1
IH5DY(IR,I)=JULD
IH5TM(IR,I)=IHR
40 IF(I2.NE.1) GOTO 50
DO 41 1=1,50
IF(C1.GT.T50(I)) GOTO 42
41 CONTINUE
GOTO 50
C THE ROUTINE FOR THE TOP-50 TABLE
-------
42 IF(I.EQ.SO) GOTO 43
DO 44 J=I,49
K=50+I-J
T50(K)=T50(K-1)
IT50DY(K)=IT50DY(K-1)
IT50TM(K)=IT50TM(K-1)
IR50(K)=IR50(K-1)
44 CONTINUE
43 T50(I)=C1
IT50DY(I)=JULD
IT50TM(I)=IHR
IR50(I)=IR
C THE ROUTINE FOR THE EXCEEDANCE TABLE
50 IF(I3.NE.1) GOTO 100
IF(Cl.LE.EX) GOTO 100
DO 51 1=1,100
IF(C1.GE.E(I)) GOTO 52
51 CONTINUE
GOTO 100
52 IF(I.EQ.IOO) GOTO 53
DO 54 J=I,99
K=100+I-J
E(K)=E(K-1)
IEDY(K)=IEDY(K-1)
IETM(K)=IETM(K-1)
IRE(K)=IRE(K-1)
54 CONTINUE
53 E(I)=C1
IEDY(I)=JULD
IETM(I)=IHR
IRE(I)=IR
100 IF(IBR.EQ.l) RETURN
CAL(IR)=0.
RETURN
END
-------
SUBROUTINE H50UT(N,NX,NY,NXY,NTOT,X,Y,XD,YD,Z,PGCT,LNCT ,
1H5,IDAY,ITIME,IELIM)
C
C SUBROUTINE H50UT — PART OF POSTZ
C THIS SUBROUTINE WRITES THE HIGH-5 OUTPUT TABLE
C
INTEGER PGCT
COMMON /FILES/ IUNT,JUNT
CHARACTER*! HSTAR(300,2)
DIMENSION X(300),Y(300),XD(300),YD(300),Z(300),H5(300,5),
1IDAY(300,5),ITIME(300,5),IELIM(300)
PK1=0.
PK2=0.
IPK1=0
IPK2=0
DO 50 1=1,NTOT
HSTAR(I,1)=' '
IF(H5(I,1).LT.PK1) GOTO 60
PK1=H5(I,1)
IPK1=I
60 HSTAR(I,2)=' '
IF(H5(I,2).LT.PK2) GOTO 50
PK2=H5(I,2)
IPK2=I
50 CONTINUE
HSTAR(IPK1,!)='*'
HSTAR(IPK2,2)='*'
PGCT=PGCT+1
WRITE(JUNT,9000) PGCT
9000 FORMAT(IX,'1',105X,'POSTZ - PAGE NO. ',I3,///)
WRITE(JUNT,9001) N
9001 FORMAT(46X,'HIGH FIVE TABLE FOR ',15,' HOUR AVERAGES1,//,
12X,'RECEPTOR COORDINATES',42X,'HIGHEST FIVE CONCENTRATIONS',/
28X,'METERS',5IX,'VALUE(DAY,ENDING HOUR)',//,
35X,'X',11X,'Y',7X,'ELEV.',9X,'HIGHEST',12X,'2ND HIGH1,12X,
4'3RD HIGH',12X,'4TH HIGH',12X,'5TH HIGH1,/,'
6 --------------------------------------------- . )
LNCT=11
IF(NX.EQ.O.OR.NY.EQ.O) GOTO 20
DO 10 1=1, NY
DO 10 J=1,NX
K=(I-1) *NX+J
IF(IELIM(K) .EQ.l) GOTO 10
WRITE (JUNT, 9002) X(J),Y(I),Z(K),H5(K,1) ,IDAY(K, 1) , ITIME (K, 1)
1HSTAR(K, 1) ,H5(K,2) , IDAY(K, 2) , ITIME(K, 2) ,HSTAR(K, 2) ,
2(H5(K,L) , IDAY(K,L) , ITIME (K,L) ,L=3,5)
9002 FORMAT ( 1X,F10.0,2X,F10.0,2X,F6.0,2X,F10.2, ' ( ' ,13, ' , ' ,12, ' ) ' ,
LNCT=LNCT+1
IF(LNCT.LT.60) GOTO 10
PGCT=PGCT+1
WRITE (JUNT, 9000) PGCT
WRITE( JUNT, 9003) N
-------
9003 FORMAT(46X,'HIGH FIVE TABLE FOR ',15,' HOUR AVERAGES (CONT.)',//,
12X,'RECEPTOR COORDINATES',42X,'HIGHEST FIVE CONCENTRATIONS',/,
28X,'METERS',51X,'VALUE(DAY,ENDING HOUR)',//,
35X,'X',11X,'Y',7X,'ELEV.'>9X,'HIGHEST',12X,'2ND HIGH
4'3RD HIGH',12X,'4TH HIGH',12X,'5TH HIGH1,/,1
,12X,
LNCT=11
10 CONTINUE
20 IF(NXY.EQ.O) GOTO 30
DO 40 1=1,NXY
',K=NX*NY+I
IF'(IELIM(K) .EQ.l) GOTO 40
WRITE(JUNT,9002) XD(I),YD(I),Z(K),H5(K,1),IDAY(K,1),ITIME(K,1),
1HSTAR(K, 1) , H5(K,2).,IDAY(K,2) , ITIME(K,2) ,HSTAR(K,2)
2(H5(K,L),IDAY(K.L),ITIME(K,L),L=3,5)
LNCT=LNCT+1
IF(LNCT.LT.60) GOTO 40
PGCT=PGCT+1
WRITE(JUNT,9000) PGCT
WRITE(JUNT,9003) N
40 CONTINUE
30' WRITE(JUNT,9004)
9004 FORMAT(/,IX,'* -'DENOTES PEAK VALUE1)
RETURN ' . '
• END '
-------
SUBROUTINE T500UT(N,NX,NY,X,Y,XD,YD,Z,PGCT,LNCT,
1T50,IDAY,ITIME,IR)
C
C SUBROUTINE T500UT --PART OF POSTZ
C THIS SUBROUTINE WRITES THE TOP-50 OUTPUT TABLE
C
INTEGER PGCT
COMMON /FILES/ IUNT,JUNT
DIMENSION X(300),Y(300),XD(300),YD(300),Z(300),T50(50),IDAY(50),
,1ITIME(50),IR(50) '
PGCT=PGCT+1
WRITE(JUNT,9000) PGCT
9000 FORMAT(IX,'1',105X,'POSTZ - PAGE NO. ',I3,///)
WRITE(JUNT,9001) N
9001 FORMAT(4IX,'TOP 50 TABLE FOR ',15,' HOUR AVERAGE CONCENTRATIONS',
1//,33X,'RECEPTOR COORDINATES',32X,'ENDING1,/,7X,'RANK',8X,'VALUE',
21IX, 'X1 , 14X, 'Y1 ,11X, 'ELEVATION' , 6X, 'DAY1 , 6X, 'TIME' ,•/, '
3 L
4 .)
LNCT=9
NT=NX*NY
DO 10 K=l,50
IF(IDAY(K).EQ.O) RETURN
IF(IR(K).GT.NT) GOTO 20
' J=(IR(K)-1)/NX.+ 1 '
I=IR(K)-(J-l)*NX
X1=X(I)
Y1'=Y(J) ' ' • '
GOTO 30
20 I=IR(K)-NT
X1=XD(I)
Y1=YD(I)
30 WRITE(JUNT,9002) K,T50(K),XI,Yl,Z(IR(K)),IDAY(K),ITIME(K)
9002 FORMAT(6X,I5,5X,F10.2,5X,F10.0,5X,F10.0,5X,F10.0,5X,I5,5X,I5)
LNCT=LNCT+1
IF(LNCT.LE.60) GOTO 10
PGCT=PGCT+1
WRITE(JUNT,9000) PGCT
WRITE(JUNT,9004) N
9004 FORMAT(4IX,'TOP 50 TABLE FOR ',15,' HOUR AVERAGE CONCENTRATIONS',
1' (CONT.)',
2//,33X,'RECEPTOR COORDINATES',32X,'ENDING',/,7X,'RANK',8X,'VALUE',
31IX,'X',14X,'Y',11X,'ELEVATION1,6X,'DAY',6X,'TIME',/,'
4
• 5 .)
LNCT=9
10 CONTINUE
RETURN
END
-------
SUBROUTINE EXOUT(N,NX,NY,X,Y,XD,YD,Z,PGCT,LNCT,
,1E,IDAY,ITIME,IR,EX)
C
C SUBROUTINE EXOUT — PART OF POSTZ
C THIS SUBROUTINE WRITES THE EXCEEDANCE OUTPUT TABLE
C
INTEGER PGCT
COMMON /FILES/ IUNT,JUNT
DIMENSION X(300),Y(300),XD(300),YD(300),Z(300),E(100),IDAY(IOO),
IITIME(IOO),IR(IOO)
PGCT=PGCT+1
WRITE(JUNT,9000) PGCT
9000 FORMAT(IX,'1',105X,'POSTZ - PAGE NO. ',I3,///)
WRITE(JUNT,9001) N,EX
9001 FORMAT(39X,'EXCEEDANCE TABLE FOR ',15,' HOUR AVERAGE CONCENTRATION
IS',/,48X,'VALUES ABOVE ',F10.2,' ARE PRINTED1,
1//,33X,'RECEPTOR COORDINATES',32X,'ENDING',/,7X,'RANK',8X,'VALUE',
21IX,'X1,14X,'Y',11X,'ELEVATION1,6X,'DAY',6X,'TIME',/,'
2 ,
4 . ,
LNCT=10
NT=NX*NY
DO 10 K=l,100
IF(IDAY(K).EQ.O) RETURN
IF(IR(K).GT.NT) GOTO 20
J=(IR(K)-1)/NX+1
I=IR(K)-(J-l)*NX
X1=X(I)
Y1=Y(J)
GOTO 30
20 I=IR(K)-NT
X1=XD(I)
Y1=YD(I)
30 WRITE(JUNT,9002) K,E(K),XI,Yl,Z(IR(K)),IDAY{K),ITIME(K)
9002 FORMAT(6X,I5,5X,F10.2,5X,F10.0,5X,F10.0,5X,F10.0,5X,I5,5X,I5)
LNCT=LNCT+1
IF(LNCT.LE.SO) GOTO 10
PGCT=PGCT+1
WRITE(JUNT,9000) PGCT
WRITE(JUNT,9004) N
9004 FORMAT(39X,'EXCEEDANCE TABLE FOR ',15,' HOUR AVERAGE CONCENTRATION
IS (CONT.)',/,48X,'VALUES ABOVE ',F10.2,' ARE PRINTED',
1//,33X,'RECEPTOR COORDINATES',32X,'ENDING',/,7X,'RANK',8X,'VALUE',
21IX,'X1,14X,'Y',11X,'ELEVATION',6X,'DAY',6X,'TIME',/,'
4
LNCT=10
10 CONTINUE
RETURN
END
-------
SUBROUTINE EROUT(N)
C
C SUBROUTINE EROUT -- PART OF POSTZ
C THIS SUBROUTINE WRITES THE ERROR MESSAGES
C
COMMON /FILES/ IUNT,JUNT
GOTO(1,2,3,4,5,6,7,8,9,10,11,12) ,N
1 WRITE(JUNT,9000)
9000 FORMAT ( IX, '*'*ERROR** NO RECEPTORS HAVE BEEN INCLUDED')
STOP
2 WRITE(JUNT,9001)
9001 FORMAT(IX,'**ERROR** GRID RECEPTOR TO BE INCLUDED IS NOT1,/,
1 IX,1 IN ORIGINAL SHORTZ RUN1)
STOP
3 WRITE(JUNT,9002)
9002 FORMAT ( IX, ' **ERROR** DISCRETE RECEPTOR TO BE INCLUDED IS NOT',/,
1 IX,1 IN ORIGINAL SHORTZ RUN')
STOP
4 WRITE(JUNT,9003)
9003 FORMAT(IX,'**ERROR** SOURCE NUMBER TO BE INCLUDED IS NOT IN',/,
1 IX,1 ORIGNIAL SHORTZ RUN1)
STOP
5 WRITE(JUNT,9004)
' 9004 FORMAT(IX,'**ERROR** SOURCE TO BE SCALED IN NOT IN ORIGINAL1,/,
1 IX,' SHORTZ RUN')
STOP
6 WRITE(JUNT,9005)
9005 FORMAT(IX,'**ERROR** THERE ARE TOO FEW HOURS TO COMPUTE A',/,
1 IX,1 3-HOUR AVERAGE')
STOP
7 WRITE(JUNT,9006)
9006 FORMAT(IX,'**ERROR** THERE ARE TOO FEW HOURS TO COMPUTE AN',/,
1 IX,1 8-HOUR AVERAGE1)
STOP
8 WRITE(JUNT,9007)
9007 FORMAT(IX,'**ERROR** THERE ARE TOO FEW HOURS TO COMPUTE A1,/,
1 IX,' 24-HOUR AVERAGE')
STOP
9 WRITE(JUNT,9008)
9008 FORMAT(IX,'**ERROR** THERE ARE FEWER DAYS IN THE TAPE FILE',/,
1 IX,' THAN HAVE BEEN SPECIFIED AS THE NUMBER1,/,
2 IX,' OF DAYS PER YEAR1)
STOP
10 WRITE(JUNT,9009)
9009 FORMAT(IX,'**ERROR** THERE ARE TOO FEW HOURS TO COMPUTE THE1,/,
1 IX,1 USER-SPECIFIED AVERAGING TIME1)
STOP
11 WRITE(JUNT,9010)
9010 FORMAT(1X,'**ERROR** THE GRID RECEPTOR FOR WHICH BACKGROUND',/,
1 IX,1 HAS BEEN SPECIFIED IS NOT IN THE ORIGINAL',/,
2 IX,' SHORTZ OUTPUT FILE')
STOP
12 WRITE(JUNT,9011)
9011 FORMAT(IX,'**ERROR** THE DISCRETE RECEPTOR FOR WHICH BACKGROUND',/
1, IX,' HAS BEEN SPEICIFED IS NOT IN THE ORIGINAL',/,
-------
2 IX,1 SHORTZ OUTPUT FILE')
STOP
RETURN
END
-------
SUBROUTINE INPOUP(I,ADR1,ADR2)
C
C SUBROUTINE INPOUP -- PART OF POSTZ,
C THIS SUBROUTINE IS TAKEN FROM THE SHORTZ MODEL AND IS USED TO
C READ THE INPUT TAPE. IT CAN ALSO BE USED TO WRITE THE OUTPUT
C TAPE, AND THAT IS ONE OF ITS PURPOSES IN THE SHORTZ MODEL
C
INTEGER TITLE,DATS,ADR2
DIMENSION XEND(14),ADR1(1),BUFIN(2000),BUFOT(2000,2),IXEND(14)
COMMON /DIM/ NSOURC,NGROUP,NXPNTS,NYPNTS,NXWYPT,NHOURS,NDAYS
1,IOVRSN '
COMMON /IO/ NOGO,IT1REC,IT2REC,NINTP,NOTTP,NINFL(3),NOTFL(3) •
COMMON /RST/ ISW(20),TITLE(20),KUNR(6),KFNR(3),DD1(3),I3,DATS(2),
1IH , IB ,,IPRT ( 100 ) , NLINES , KSW, NGROPS , IHOUR ( 1464 ) , MUNT ( 4 ) , ITABLE , IPAGE
2,ND,LINE,IPC,IBLNK
COMMON /FILES/ IUNT,JUNT
EQUIVALENCE (XEND,IXEND),(JOVRSN,XOVRSN),(IVSN,XVSN)
DATA JOVRSN/9903/
DATA IXEND/0541600000000,13*0/,NI,NO,ISREC/3*!/,IT1REC,IT2REC,NOB,
1IEND,LSW,IWRD/6*0/,NIB,NWRD/2*2000/,XFT/0.0/
NOGO = 0
NINTP=1
NINFL(1)=2
IF (I .EQ. 2) GO TO 200
IF (I .EQ. '3) GO TO 350 .
IF (LSW .EQ. 0) GO TO 20
IT1REC = IT1REC+1
CALL NTRAN$'(NINFL(NI) ,22,2,ADR2,ADR1,IEOF)
10 IF (IEOF .GT. 0) GO TO 370
IF (IEOF .EQ. -1) GO TO 10
IF (IEOF .EQ. -2) GO TO 180
GO TO 190
20 N = 0
30 N = N+l
IF (N .GT. ADR2) GO TO 370
40 NIB = NIB+1
IF (NIB .GT. IWRD) GO TO 80
ADR1(N) = BUFIN(NIB)
IF (I .GT. 0) GO TO 30
IF (ADR2 .NE. 8) GO TO 50
IF (N .LT. 8) GO TO 30
XVSN = ADR1(N)
IF (ADR1(N) .EQ. XOVRSN) GO TO 30
IVSN = 9902
NIB = NIB-1
IOVRSN = JOVRSN
GO TO 30
50 IF (IVSN .EQ. JOVRSN) GO TO 30
IF (ADR2 .NE. 49.0R.N .LT. 49) GO TO 30
NIB = NIB-13
K = 2
CALL FFDASC(K,TITLE(15),KFNR(1))
K = 4
CALL FFDASC(K,TITLE(11),KUNR(1))
DO 60 K=l,10
-------
60 TITLE(21-K) = TITLE(ll-K)
K = 10
CALL FFDASC(K,TITLE(11).TITLE(1
DO 70 K=16,20
70 TITLE(K) = IBLNK
GO TO 30
80 JSW = 1
IEND = 0
IT1REC = ITi:
90
100
NIB
CALL
IF (
IF (
IF (
IF (
IF (
= 0
NTRAN$(:
IEOF
IEOF
IEOF
JSW
IEND
.GT
.EQ
.EQ
.EQ.
.EQ
110
120
130
140
150
160
170
180
NINFL(NI) ,22,2/NWRD,BUFIN/IEOF)
0) GO TO 130
-1) GO TO 100
-2) GO TO 110
2) GO TO 150
1) GO TO 160
WRITE (JUNT,9001) NINFL ( NI )
LINE = 200
GO TO 190
IF (JSW .EQ. 2) GO TO 150
IEND = IEND+1
IF (IEND . GT. 1) GO TO 120
IT1REC = IT1REC-1
LE.
IT1REC
IF (IT1REC
GO TO 90
WRITE (JUNT
LINE = 200
GO TO 180
IWRD = IEOF
IF (IT1REC
IF (IEOF .NE
IF (JSW .EQ.
IF (IEND .EQ
IF (BUFIN(l)
IF (IEOF .NE
JSW = 2
GO TO 90'
CALL NTRAN$
NI = NI+1
IF (NI .GT.
GO TO 80
WRITE (JUNT
LINE = LINE+2
CALL NTRAN$(NINFL(NI
GO TO 180
WRITE (JUNT,9004)
LINE = LINE+2
NI = 1
GO TO 80
NOGO = 1
IT1REC = 0
LSW = 0
IWRD = 0
NIB = 2000
GO TO 370
0) GO TO 80
9002) NINFL ( NI ), IT1REC
GT. 1) GO TO 140
NWRD) LSW = 1
2) GO TO 160
0) GO TO 40
.NE. XEND(1)) GO TO 160
14) GO TO 160
NINFL(NI),22,11,22,22,22)
NINTP) GO TO 170
9003) NINFL(NI),IT1REC
) ,22,7,-1,22,22)
-------
190 NOGO = 2
IF (IT2REC .LE. 0) GO TO 370
GO TO 360
200 JSW = 1
N = 0
210 N = N+l
IF (N .GT. ADR2) GO TO 370
220 NOB = NOB+1
IF (NOB .GT. NWRD) GO TO 230
BUFOT(NOB,2) = ADR1(N)
GO TO 210
230 IT2REC = IT2REC+1
240 IF (NOB .GT. NWRD) NOB = NWRD
IF (NOB .LT. 4) NOB = 4
CALL NTRAN$(NOTFL(NO) ,22,1,NOB,BUFOT(1,2),IEOF)
250 IF (IEOF .GT. 0) GO TO 330
IF (IEOF .EQ. -1) GO TO 250
IF (IEOF .EQ. -2) GO TO 270
260 WRITE (JUNT,9005) NOTFL(NO),IT2REC
LINE = 2,00
GO TO 190
270 IF (NO+1 .LE. NOTTP) GO TO 280
WRITE '(JUNT,9006)
LINE = LINE+2
NOGO = 2
GO TO 370
280 IF (IT2REC .GT. 2) GO TO 310
WRITE (JUNT,9009) NOTFL(NO)
LINE = LINE+2
NO = NO+1
IF (NO .GT. NOTTP) GO TO 270
IF (IT2REC .LE. 1) GO TO 300
CALL NTRAN$(NOTFL(NO),22,1,NWRD,BUFOT(1,1),IEOF)
290 IF (IEOF .GT. 0) GO TO 300
IF (IEOF .EQ. -1) GO TO 290
IT2REC = 1
GO TO 260
300 IT2REC = 0
GO TO 230
310 CALL NTRAN$(NOTFL(NO),22,7,-2,22,22)
CALL NTRAN$(NOTFL(NO),22,9,22,22,22)
CALL NTRAN$(NOTFL(NO),22,1,14,XEND,IEOF)
CALL NTRAN$(NOTFL(NO),22,9,9,11,22)
IT2REC = IT2REC-2
WRITE (JUNT,9007) NOTFL(NO),ISREC,IT2REC
LINE = LINE+2
NO = NO+1
IF (NO .GT. NOTTP) GO TO 270
XFT = FLOAT(NWRD)*0.005625+0.75
IT2REC = IT2REC+1
ISREC = IT2REC
CALL NTRAN$(NOTFL(NO),22,1,NWRD,BUFOT(1,1),IEOF)
320 IF (IEOF .GT. 0) GO TO 230
IF (IEOF .EQ. -1) GO TO 320
GO TO 260
-------
330 XFT = XFT+FLOAT(NWRD)*0.005625+0.75
NOB = 0
IF (JSW .NE. 1) GO TO 360.
DO 340 J=1,NWRD
340 BUFOT(J,1) = BUFOT(J,2)
GO TO 220
350 JSW = 2
IF (NOB .GT. 0) GO TO 230
360 CALL NTRAN$(NOTFL(NO),22,9,9,8,-1)
XFT = (XFT+12.0')/12.0
WRITE (JUNT,9008) NOTFL(NO),ISREC,IT2REC,XFT
LINE = LINE+1
370 CONTINUE
RETURN
9001 FORMAT (24H1*** READ ERROR ON UNIT ,I2,11H AT RECORD ,I6/)
9002 FORMAT (25H1*** END OF DATA ON UNIT ,I2,2H, ,I6,13H RECORDS READ/)
9003 FORMAT (25HO*** END OF FILE ON UNIT ,I2,2H, ,I6,13H RECORDS READ/)
9004 FORMAT (86HO*** WARNING - MORE INPUT REELS THAN UNITS ASSIGNED PRO
1G. GOING TO FIRST UNIT ASSIGNED/)
9005 FORMAT ( 25H1 *** WRITE ERROR ON UNIT , 12', 11H AT RECORD ,I6/)
9006 FORMAT (71HO*** ERROR - MORE OUTPUT REELS REQUIRED THAN UNITS AS'SI
1GNED, PROG STOPS)
9007 FORMAT (32HO*** END OF OUTPUT' REEL ON UNIT,,12,9H RECORDS ,I6,9.H T
1HROUGH ,I6,8H WRITTEN/)
9008'FORMAT (32H * * * END OF OUTPUT DATA ON'UNIT ,I2,9H RECORDS ,I6,9H T
1HROUGH ,I6,10H WRITTEN, ,F8.3,18H FEET OF TAPE USED/)
9009 FORMAT (47HO*** WARNING - NOT ENOUGH ROOM ON REEL ON UNIT •, I2,44H
1PROG. STARTS FIRST OUTPUT REC. ON NEXT REEL/)
END
-------
80779*101
REPORT DOCUMENTATION
PAGE
1. REPORT NO.
EPA-910/9-86-144
3. Recipient1 • AccMslon No.
4, T»l* and Subtltl*
USER'S GUIDE FOR POSTZ
A POST-PROCESSOR FOR THE
SHORTZ AIR QUALITY MODEL
S. Rtpert Oat*
July 1986
7. Author**)
Kirk D. Winges
B. Performing Organization R*pt. No:
3481-R81
9. Performing Organization Nam* and Addr***
TRC Environmental Consultants,
15924 22nd Avenue SE .
Mill Creek, Washington 98012
10. Project/Task/Work Unit No.
Inc.
11. Contract(C) or Grant(G) No.
(0 68-02-3886
(G)
12. Sponsoring Organization Nam* and Address
U. S. Environmental
Region 10
1200 Sixth Avenue
Seattle. Washington
Protection Agency
98101
13. Typ* of Report & Period Covered
Final
14,
IS. Supplementary NotM
IS. Abstract (Limit: 200 word*)
This report provides instructions for the use of the POSTZ computer
program, a post-processor for the SHORTZ air quality model. The
POSTZ program is designed to read a sequential output tape'of
concentrations produced by the SHORTZ model. The POSTZ program
provides summary tables of the highest concentrations for different
averaging times. The POSTZ program allows the user to restrict the
analysis to specific receptors, sources or periods of meteorological
data. In addition, the user can scale the concentrations from a
particular source if desired. The user can select between running
or block averages, and can implement the EPA calms policy if
desired. The output tables can be in the form of the highest five
concentrations for each receptor, the top 50 concentrations for all
receptors, or all concentrations computed to be above a user
specified criteria value.
17. Document Analyses a. Descriptor*
Air Pollution, Meteorology, Turbulent Diffusion
b. Identifiers/Open-ended Terms
Dispersion Modeling, Complex Terrain, SHORTZ, Post-processor
c. COSATl Field/Group
18. Availability Statement
Release unlimited
19. Security Class (This Report)
Unclassified
20. Security Class (This Page)
Undas siiie_d_
21. No. of Pages
96
22. Price
(SeeANSI-239.18)
See Instructions on Reverse
OPTIONAL FORM 272 (4-77)
(Formerly NTIS-35)
Department of Commerce
------- |