EPA-650/3-73-001
September 1972
Ecological Research Series
-------
EPA-650/3-73-001
USER'S MANUAL
FOR THE APRAC-1A
URBAN DIFFUSION MODEL
COMPUTER PROGRAM
by
R.L. Mancuso and F.L. Ludwig
Stanford Research Institute
Menlo Park, California 94025
Contract No. CPA 22-69-0064
Project Element No. 1A1009
EPA Project Officer: Dr. Warren B. Johnson, Jr.
Meteorology Laboratory
National Environmental Research Center
Research Triangle Park, N.C. 27711
Prepared for
COORDINATING RESEARCH COUNCIL
30 ROCKEFELLER PLAZA
NEW YORK, NEW YORK 10020
CONTRACT CAPA-3-68(l-69)
and
OFFICE OF RESEARCH AND DEVELOPMENT
U.S. ENVIRONMENTAL PROTECTION AGENCY
WASHINGTON, D.C. 20460
September 1972
-------
This report has been reviewed by the Environmental Protection Agency and
approved for publication. Approval does not signify that the contents
necessarily reflect the views and policies of the Agency, nor does
mention of trade names or commercial products constitute endorsement
or recommendation for use.
11
-------
CONTENTS
LIST OF ILLUSTRATIONS IV
LIST OF TABLES V
LIST OF SYMBOLS yi
I INTRODUCTION 1
II SUMMARY OF MODEL DEVELOPMENT AND
EVALUATION PROGRAM 3
III BASIC MODEL 7
A. Traffic Data and Emission Rate 7
B. Intraurban Diffusion 13
C. Extraurban Diffusion 17
D. Local Street Diffusion 17
E. Transport Wind, Mixing Depth, and Stability
Index 20
IV COMPUTER PROGRAM • 21
A. Program Usage 21
B. Data Input 26
C. Data Output and Program Capabilities 33
REFERENCES 39
APPENDICES
A. SRI CDC 6400 VERSION OF APRAC-1A COMPUTER PROGRAM 41
B. EPA IBM 360/50 VERSION OF APRAC-1A COMPUTER PROGRAM 69
C. EXAMPLE OF INPUT DATA 107
D. EXAMPLE OF OUTPUT DATA 117
iii
-------
ILLUSTRATIONS
Figure 1 Computer Display of Traffic Links
for Chicago 8
Figure 2 Hourly Distribution of Traffic for Two
Facility Types in St. Louis 10
Figure 3 Portion of a Typical Traffic Map 11
Figure 4 Diagram of Segments Used for Spatial
Partitioning of Emissions 14
figure 5 Vertical Diffusion According to Gaussian
Formulation 15
Figure 6 Schematic of Cross-Street Circulation
Between Buildings 18
Figure 7 Specifications for Leeward and Windward Cases
on the Basis of Receptor Location, Street
Orientation, and Wind Direction 19
Figure 8 APRAC-1A Flow Chart 22
Figure 9 Meteorological Inputs to the Model with Observed
and Calculated CO Concentrations at the St. Louis
CAMP Station 15
Figure 10 Calculated St. Louis CAMP Station CO Concentration
Frequency Distribution for 1965 Traffic Conditions:
0800, 1200, and 1800 Hours 36
Figure 11 Calculated St. Louis Concentration Patterns for Two
Grid Sizes 37
IV
-------
TABLES
Table 1
Table 2
Table 3
Table 4
Table 5
Table 6
Facility Codes and Car Speeds for
St. Louis ,
Values of a for Cars Produced After 1970
9
10
Basic Input Information: Cards A
through M 27
Traffic Input Data: Cards N and 0 31
Meteorological Input Data: Cards P,
Q,, and R 32
Central Processor Times on the SRI CDC 6400
Computer 38
-------
LIST OF SYMBOLS
Text
Symbols
a, b
c
c
Program
Symbols
A, B
PF1, P01
CCAL
CCAL
EXTRAQ
E
F
H
h
i
J
K> Lo
P
t
—
FUEL
Z3
HT
I
IJ or JI
CK, XL0
PT12, PT34,
PT6, PTSAT,
PTSUN
u
W
X
WS
WWST
R or YR
Definition
Parameters to define o
Vehicle emission constants
CO concentration at receptor
CO concentration computed by
extraurban and intraurban models
CO concentration computed by extra-
urban model
CO concentration computed by street
model
Vehicle emission rate
Annual consumption of fuel within a
22.5-degrees angular sector from
32 to 1000 km upwind of receptor
Average building height
Mixing depth
Upwind area segment index
Stability class index
Empirical constants
Hourly traffic factor
Emission rate for area segment
Average vehicle speed on
traffic link
Vertical dispersion coefficient
Airport wind speed
Transport or rooftop wind speed
Width of street
Upwind distance from receptor
VI
-------
I INTRODUCTION
*
The APRAC-1A diffusion model has been developed as a versatile and
practical model for computing the concentrations of pollutants at any
point within a city. It is the result of the studies of Ludwig et al.
t
(1970, 1972) and Johnson et al. (1971) and includes the most recent
modeling features of these studies. The model calculates pollutant
contributions from diffusion on various scales, including:
• Extraurban diffusion, mainly from sources in upwind cities.
• Intraurban diffusion from freeway, arterial, and feeder
street sources.
• Local diffusion of emissions within a street canyon.
Currently, the model treats only carbon monoxide (CO), a relatively inert
gas in the atmosphere but an important pollutant in terms of health.
Motor vehicles are the major source of this gas.
A brief summary of the APRAC-1A model development and evaluation
program is given in Section II. The basic model formulation is briefly
described in Section III. More detailed descriptions of both the theory
and the numerical techniques are given in reports by Ludwig et al. (1970,
1972) and Johnson et al. (1971). These reports should be used to sup-
plement the information in this manual.
*
The acronym APRAC is derived from the initial letters of the Air
Pollution Research Advisory Committee, under whose auspices this
research has been conducted. The members of this committee are
drawn from the Coordinating Research Council (CRC) and the Environ-
mental Protection Agency (EPA). The designation 1A refers to the
present version of the model; with future improvements the designa-
tion will be changed to APRAC-1B, APRAC-1C, and so on.
t
References are listed at the end of the body of the report.
-------
A detailed description of the APRAC-1A computer program and its
use is given in Section IV. The computer program can be used to make
calculations of the following types:
• Synoptic model: hourly concentrations as a function of time,
for comparison and verification with observed concentrations
and for operational applications.
• Climatological model: the frequency distribution of concen-
trations, for statistical prediction of the frequency of
occurrence of specified high concentrations in connection
with planning activities.
• Grid-point model: concentrations at various locations in
a geographical grid, providing detailed horizontal con-
centration patterns for operational or planning purposes.
A complete listing of the program and listings of input and output
data of an actual run are given in Appendices A, B, C, and D.
-------
II SUMMARY OF MODEL DEVELOPMENT AND EVALUATION PROGRAM
The APRAC-1A urban diffusion model was developed to simulate CO
concentrations from readily available meteorological and traffic data.
The model is based on existing experimental data and previous research
results. The CO concentrations calculated with the model were initially
compared with measured data from Continuous Air Monitoring Program (CAMP)
stations ; the calculated and observed values often differed signif-
icantly in magnitude, although they had similar trends.
An extensive measurement program was undertaken in San Jose,
California to determine the causes of the discrepancies between calcu-
lated and observed concentrations. The measurements showed that roof-
level winds blowing across a street canyon cause a helical circulation
in the canyon. The resulting street-level CO concentrations differed
by as much as a factor of three from one side of the street to the
other. One of the principal accomplishments of the research in San
Jose was the development of a new submodel to describe these street-
canyon effects. The submodel substantially improved the agreement
between observations and calculations. The San Jose program also un-
covered and corrected other shortcomings of the original model. The
resulting changes resulted in more realistic specification of atmospheric
stability and turbulent diffusion in urban areas.
San Jose is a moderate-sized city and the question arose whether the
San Jose results would be generally applicable to larger cities with
taller buildings. To answer this question another extensive measurement
program was undertaken in St. Louis. One of the primary concerns of the
St. Louis research was evaluation of the performance of the street-canyon
-------
submodel in street canyons deeper than those studied in San Jose. To
test the generality of the San Jose observations, two street canyons, with
height-to-width ratios of about 1.5 and 2, were instrumented in St. Louis,
so that CO concentrations could be measured on both sides of each street
at five heights, from 4 m to roof level. Concentrations were also measured
in midstreet at 7 m and 35 m (roof level). Winds in the street canyon
were measured on either side, at 4.5 m and at roof level. Larger-scale
airflow and CO concentrations in the area were monitored with instru-
mentation up to a height of 130 m on a television transmitting tower on
top of a building at the intersection of the two streets. These data
*
are available on magnetic tape from the National Climatic Center.
The data collected show that a single-helix circulation is found in
the deep street canyons of St. Louis and that the simple model developed
from the San Jose data is fundamentally correct for the deeper canyons.
Some slight modifications were required to account for the entrainment of
recirculating polluted air in the downward-flowing part of the helical
circulation. There had been evidence of this entrainment in the San
Jose data also. The data indicate that the helical circulation develops
when the roof-level winds are at an angle of more than 30 degrees to the
street direction. When the winds are more nearly parallel to the street,
cross-street gradients were found to be small. For winds parallel to the
street, the street-canyon submodel describes the vertical gradients as
an average of the two expressions used when the winds are blowing across
the street. The small changes that were made in the street-canyon sub-
model have improved the model's ability to predict the CO gradients in
street canyons. Observations made at street level with mobile equipment
*
National Climatic Center
National Oceanic and Atmospheric Administration
Federal Building
Asheville, N.C., 28801.
-------
indicate that the street-canyon submodel is applicable through most
of the block, at least to within about 10 m of the intersection.
The submodel used to calculate atmospheric stability was revised
to give results that are more consistent with the fluctuations of wind
direction observed on the television tower. It was found that for a
given stability type, there is appreciably greater fluctuation in wind
direction when the air has an urban fetch than when it has a nonurban
fetch. This fact lends support to the revisions made during an earlier
phase of the program that effectively increased diffusion rates in the
urban situation as compared to the rural.
Model calculations of mixing depth were compared with lidar (laser
radar) observations of the aerosol layer, and with radiosonde measure-
ments of the temperature profile near the downtown center. From these
comparisons it was concluded that the mixing-depth submodel does as
well as is now possible using routinely available data.
Helicopter and van measurement of CO concentrations around the
downtown area were combined with wind speed measurements so that a mass
budget analysis could be performed to estimate the rate of CO emissions
by traffic in the study area. The emission submodel was applied to the
same area, and the results were compared. Uncertainties in the wind
field, possible changes in CO emission rates during the measurement
periods, and uncertainties in traffic amounts all contribute to the
difficulties in making reliable comparisons, but the results were
adequate to show any serious deficiencies in the model. On the average,
the three cases analyzed in St. Louis and the four earlier cases from
San Jose show agreement within a factor of two. With the data currently
available there does not seem to be sufficient justification for changing
the emission submodel.
-------
The revisions and additions that have been made in the model since
it was originally formulated have substantially improved its performance.
When the revised model was applied in St. Louis and compared with observed
hour-average CO concentrations, the root-mean-square difference between
the calculated and observed values ranged from 2.6 to 3.9 ppm, depending
on the particular observation site. This is about half the uncertainty
of the original model when it was applied to this same city. If the
calculated and observed data are fit by linear regression, the correspon-
ding differences are reduced to values between 1.6 and 3.3 ppm. Also, the
correlation coefficients between calculated and observed CO have been
improved substantially. They are now in the range 0.4 to 0.7, as opposed
to the 0.2-to-Q.4 range found before revision. The ability of the model
to specify frequency distributions of concentration is good. Median and
90-percentile concentrations are specified within about 2 and 3 ppm of the
observed values, respectively. Use of regression relationships derived
from the observed concentrations and those calculated with the model
reduces the error in specifying median and 90-percentile concentrations
to about 1.3 ppm.
The APRAC-1A model is now sufficiently accurate that it can be used
for planning purposes. Some improvements and extensions are still de-
sirable, including better specification of emissions and a new submodel
to describe the effects that take place in the immediate vicinity of
a freeway.
-------
Ill BASIC MODEL
A. Primary Traffic Data and CO Emission Rate
The APRAC-1A diffusion model is based primarily on the CO emissions
from a network of traffic road segments or links. An example of the
network of traffic links for Chicago is shown in Figure 1. Each of these
links is assigned an average daily traffic volume (expressed in vehicles
per day) based on historical, current, or forecast data obtained from
appropriate traffic agencies. Each link is identified by the geographical
locations of its end points and is designated as a particular road type,
such as freeway arterial, or local street (see Table 1). To calculate
the traffic volume for a given hour, an hourly factor (P.), which is
illustrated in Figure 2, is first applied to the daily traffic volumes.
Values for this factor can also be obtained from traffic agencies. Then
the CO emission rate, E (g-CO vehicle-mi ), is estimated from the mean
vehicle speed. S (mi h ) by using an empirical equation of the form
E = aS~0 (1)
0' and (3 are constants that depend on the characteristics of the emission
control devices installed and the mixture of old and new model cars on
the road. For current CO calculations ot is taken to be 700 and (3 to be
0.75. These values are appropriate to a mixture of about half pre-1968
(1966 for California) and half newer cars (Johnson et al., 1971). For
cars produced since 1968, the value of p has been 0.48. Existing and
potential legislation requires Q< to decrease with time, as shown in
Table 2. For future years, the effective values of & and p for use in
Eq. (1) have to be determined on the basis of the fraction of the total
-------
TA-7874-78
FIGURE 1 COMPUTER DISPLAY OF TRAFFIC LINKS FOR CHICAGO
cars represented by each model year. The model values of S are shown in
Table 1. These values depend on the type of road and the time of day,
that is, whether peak or off-peak traffic hours. The total CO emission
from a given traffic link is found by multiplying E by the length of the
link (mi) and by the hourly traffic volume.
-------
Table 1
FACILITY CODES AND CAR SPEEDS FOR ST. LOUIS
Traffic
Link Type
Downtown
freeway
arterial
Suburban
freeway
arterial
Local
Street model
Code
1
2
3
4
5
6
7
Average Car Speed
in Off-Peak Hours*
(mi h-1)
43
9
53
20
12
9
5
Peak-hour car speeds are set at 0.85
of off-peak-hour values.
Although it is possible to use emission data averaged over 1- or
*
2-mi grid squares, this procedure is believed to be undesirable because
of the loss of spatial resolution. Carbon monoxide concentrations are
highly variable over short distances in the vicinity of roadways (Ott,
1972), and much of this variability is reflected in the calculated con-
centrations when the input data for individual links are used. However,
sample calculations using smoothed emission fields indicate that by this
means the variability is reduced to about half that obtained by using
link input data (Coventry and Ruff, 1972).
This modification of the model has in fact been mpde by D. Coventry
and R. Ruff of the Environmental Protection Agency, Research Triangle
Park, North Carolina.
-------
12
10
c
I-
>-
_l
<
I- 4
o
a:
RADIAL EXPRESSWAY
CIRCUMFERENTIAL ARTERIAL
I
oo
M
08 12
HOUR OF DAY
16
LST
20 24
TA-7874-15
FIGURE 2 HOURLY DISTRIBUTION OF TRAFFIC FOR
TWO FACILITY TYPES IN ST. LOUIS
Table 2
VALUES OF Q' FOR CARS PRODUCED AFTER 1970
Model Years
1972-1974
1975-1979
After 1980
3
160
16
8
A typical form in which traffic data are made available by traffic
agencies is shown in Figure 3, a segment of the traffic map of the St.
Louis area. As is often the case, the downtown traffic must be obtained
from a separate, more detailed map of the same general form. Figure 3
10
-------
1970
TRAFFIC MAP
OF
ST. LOUIS
ST. LOUIS COUNTY
NOTE- Figures represent estimated average annual two-way weekday traffic
volumes.
SOURCE: Prepared by the Missouri State Highway Department Division of
Planning in cooperation with the U.S. Department of Transporta-
tion Federal Highway Administration.
FIGURE 3 PORTION OF A TYPICAL TRAFFIC MAP
TA-8563-128
11
-------
also illustrates another common problem, that of different jurisdictions
within the same metropolitan area. In this case, East St. Louis (Illinois)
traffic must be obtained from another map, and that map may be of a
different scale, requiring some coordinate transformation.
When traffic maps have been obtained from the appropriate agencies,
the x and y coordinates must be measured for the end points of each of
the links. These coordinates, the daily traffic volume, and the street
type must be punched on cards. The street type should probably be assigned
by a traffic engineer or someone of similar background. The rest of the
operation can be accomplished by a data aide. Typically, it takes one
ti two weeks for a data aide to convert the data for an area like St.
Louis from a form like that shown in Figure 3 to punchcards suitable for
input to the model . The more maps of different scales that are required,
the longer will be the time.
Not all the city's traffic is represented by the primary links
shown in Figure 3. The number of vehicle-miles traveled on streets not
represented by the primary network is computed from an estimate of the
total vehicle-miles traveled in the area (based on such a measure as
total fuel consumption) less the total vehicle-miles traveled on the
links of the primary network. The number thus computed is distributed
over the study area by estimating the relative density of local streets
as opposed to parks, open spaces, or streets already coded, for each 4-
2
mi area of a 2-mi-by-2-mi grid covering the area. The emissions from
the local street travel in a given square are assumed to emanate uniformly
from that square. Although the emission per vehicle mile is high on the
local streets because of low speeds, the overall contribution is small
because of the small number of vehicles.
Most urban areas in the United States have completed an area-wide
transportation study to determine traffic demands and transportation
12
-------
facility needs for future years, for example, 1980 or 1990. Such a study
is required for participation in a Federal Highway Administration program,
An important result of such a study is a design of,a future traffic net-
work for an area and a forecast of the traffic volumes on the links of
this network. The procedure for conducting these studies has been highly
developed and partially standardized. The emission inventory components
of the APRAC-1A diffusion model are designed so that the network descrip-
tion (including link length and facility types), link volume, and link
speed data of the widely used traffic-planning computer programs can
serve as input for the diffusion model. In most cases, the only manual
step required to use forecast traffic conditions will be the measurement
and coding of node coordinates for the network.
Traffic forecasts include travel on both primary and secondary net-
works. The primary network links are usually represented in the traffic
forecast analysis much as they appear in Figure 3. However, local or
feeder streets are represented in the analysis as connectors between the
assumed center of population of a traffic zone, where all traffic in that
zone is assumed to originate, and points on the primary network. The
vehicle-miles on connector links therefore approximate those expected on
local streets, but the traffic is concentrated on a few fictitious links,
rather than spread over a broad area. This can be compensated for in the
model by averaging the emissions from these fictitious secondary links
over the 2-mi by 2-mi background grid, so that the connector links never
explicitly appear in the calculations as do the individual primary net-
work links.
B. Intraurban Diffusion
The intraurban diffusion calculation uses a number of area segments
spaced at logarithmic upwind-range intervals from a receptor point, as
shown in Figure 4. These area segments are oriented normal to the
13
-------
DC
O
H
Q.
LU
O
HI
CC
CO
Z
O
^
5
ai
O
Z
O
.
O Z
5 ID
5 u.
UJ go
> HI c
< I "J
CL Z (_ U
X Z ^ lil
uj < 2
-------
direction of the transport wind, and they overlay the traffic network.
The traffic links and portions of links falling within each area segment
are identified, and the emissions from the individual links are accumu-
lated. Emissions are assumed to be released uniformly over the entire
area of the segment. To save computer time, the emissions within the
four segments farthest from the receptor are calculated by using a grid-
point technique described by Ludwig et al. (1970). The CO contributions
from each of the nine area sources are computed individually and then
added to find the total intraurban concentration of CO at the receptor
poi nt .
A "Gaussian-plume" diffusion formulation is used for the calcula-
tions. The vertical concentration profile from a crosswind line source
(such as a road) is assumed to be Gaussian in shape, as shown schemati-
cally in Figure 5. The spread of this vertical concentration distribution
02 DEPENDS UPON
• TRAVEL DISTANCE
• ATMOSPHERE STABILITY
HEIGHT
GAUSSIAN
VERTICAL
CONCENTRATION
PROFILE
LINE
SOURCE
DISTANCE
TA-8563-49
FIGURE 5 VERTICAL DIFFUSION ACCORDING TO GAUSSIAN FORMULATION
is characterized by its standard deviation, a . On the basis of experi-
z
mental data, a is taken to have the form
' z
15
-------
= ax , (2)
where x is the downwind distance, and a and b are parameters that depend
on the atmospheric stability (see Johnson et al., 1971, Section V-D).
The contributions from one of the upwind area sources to the CO concen-
tration (C) at the receptor are then given by
•I -J\
I '- • ± \ 1 " . i ~ j I » V ^ /
-2 -1
where Q is the average area emission rate (g m s ) from the segment,
and u is the transport wind speed. The subscript i denotes the different
segments; x. is the end point and x. the beginning point of Segment i.
The subscript j denotes different stability classes. The model is so
formulated that the values of a and b need not be the same for all seg-
ments. However, the values currently used in the computer program are
functions only of stability class, j, as suggested by urban tracer ex-
periments (Johnson et al., 1971).
A simple "box" model,
x - v
M a i+1
V
is applied for distant segments when there is a limiting mixing depth
(h) determined by the vertical temperature stratification. Under these
conditions, pollutants tend to be distributed uniformly in the vertical
after sufficient travel has taken place. A change from the Gaussian
model to the box model is made at the distance where the two would give
equal values of concentration if applied to a line source.
16
-------
C. Extraurban Diffusion
The extraurban diffusion calculation is simpler than the intraurban.
On the basis of the box-model concept, the CO contributed by extraurban
sources (C ) is estimated from the equation
e
5.15 x 10~UF ...
C = ; , (5)
e uh
-1
where F is the annual consumption of fuel (gal yr ) within a 22.5-degree
angular sector extending from 32 to 1000 km upwind of the receptor loca-
tion. The input for a given city includes a table of 16 values of F that
are associated with 16 wind direction categories. These segment fuel
consumptions can be calculated from information available in Federal
Highway Administration publications. Only one value of C is computed
for each day, by using a speed and a direction that are representative
of the strongest wind observed during the day.
D. Local Street Diffusion
Evidence oi a helical air circulation in street canyons, as illus-
trated in Figure 6. lias been observed (see Johnson et al., 1971). Re-
ceptors on the leeward side of a building (to the right side as shown in
Figures 6 and 7) are exposed to substantially higher concentrations than
are those on the windward (left) side because of the reverse flow component
across the street, near the surface. Thus, the concentration (C) at a
receptor may be considered as having two superimposed components. One
component is the concentration (C ) of the air entering the street canyon
from above. (It is assumed that the concentration computed by the extra-
urban and intraurban diffusion models represents C, ) . The other component
b
(AC) arises from the locally generated CO emissions within the street.
Hence, we have
C = Cb + AC (6)
17
-------
MEAN
WIND
(u)
BACKGROUND
CO CONCENTRATION
PRIMARY RECEPTOR
VORTEX
TRAFFIC
LANE
-W-
TA-8563-67R
FIGURE 6 SCHEMATIC OF CROSS-STREET CIRCULATION
BETWEEN BUILDINGS
Equations for calculating the £C components on both the leeward side
C/\C ) and the windward side (Ac,.,) were derived by Johnson et al . (1971)
Li 'V
and modified by Ludwig et al. (1972). The leeward component is calculated
AC =
I 2 2\1/2 1
(u + 0.5) MX + z ) + L0J
(7)
In this equation, K is an empirical nondimensional constant (K =• 7) ; L
is a dimension representing the vehicle size (L — 2 m); and x and z are
o
the horizontal and vertical distance of the receptor relative to the
center of the traffic lane (see Figure 6). Also, u is the rooftop wind
speed (ms ) generally estimated from the airport wind speed U by a
a
regression relationship, and Q is the average rate of emission of CO
18
-------
LEEWARD
(3 + 210° SINGLE 0 + 150'
STREET
TA-7874-22R
FIGURE 7 SPECIFICATION FOR LEEWARD AND WINDWARD
CASES ON THE BASIS OF RECEPTOR LOCATION,
STREET ORIENTATION, AND WIND DIRECTION
(g m s ) in the street. Values of Q are computed by multiplying the
street's daily traffic volume by the hourly traffic factor (P ) and the
vehicle emission rate (E), and by•converting the units. The windward-side
component (&C ) is calculated by
AC... =
K Q (H - z)
W W(u + 0.5)H
(8)
where W is the street width, and H is the average building height. When
the wind direction is such that neither a leeward nor a windward case is
appropriate, an intermediate concentration (AC,) is calculated by combining
the above two equations,
=1/2
(9)
19
-------
E, Transport Wind, Mixing Depth, and Stability Index
The meteorological input variables required for the model are the
transport wind direction and speed, the mixing depth, and the atmospheric
stability type. The model is designed to be generally applicable to any
city where conventional airport weather observations might be the only
observations available. Thus, the airport surface wind speed and direc-
tion are used to estimate the transport wind. Special methods were de-
veloped to calculate the mixing depth and stability index from the available
meteorological observations. The method used for the mixing-depth cal-
culation is based on the U.S. Weather Service's nearest morning (1200
G" ""/ , upper air temperature sounding. This sounding, together with the
maximum afternoon temperature at the surface, permits the afternoon or
maximum mixing depth to be calculated. The morning or minimum mixing
depth is calculated by using a simple urban model and an empirical re-
lationship involving city size and urban and rural nighttime temperatures.
Hourly mixing depths are then interpolated on the basis of the observed
hourly surface temperatures for the daylight and premldnight hours;
mixing depth is assumed to be constant between the hours of midnight and
dawn. The method used to determine the stability index depends on pre-
vailing insolation and wind speed during daylight hours and on cloud
amount (opaque) and wind speed during nighttime hours.
20
-------
IV COMPUTER PROGRAM
A flow diagram of the computer program of the APRAC-1A diffusion
model is shown in Figure 8. The program was originally written for a
CDC 6400 computer, and a CDC FORTRAN (Version 2.1) listing compatible
with that computer is given in Appendix A. The program was revised by
D. H. Coventry and R. E. Ruff of the EPA so that it could be run on that
Agency's IBM 360/50 computer. The revised IBM-compatible program listing
is given in Appendix B. Both these listings contain some programming
statements that have been used for special purpose calculations but are
by-passed in the current configuration.
The program has been organized in the form of several subprograms.
Some of these subprograms are designed to convert conventional meteoro-
logical observations to the requirements of the model, that is, to
stability index or mixing depth. By using subprograms for these purposes,
changes can be made with minimum disruption of the rest of the program.
Thus, if future research provides better methods for determining stability
index, the user can incorporate these new methods relatively easily.
Similarly, the local street-canyon effects are treated in a separate
subprogram that can be changed or replaced as knowledge of the phenomenon
increases . The program has been designed to be as practical as possible
in regard to the data that it requires and the ease with which it can
be modified to incorporate the results of continuing research.
A . Program Usage
Three versions of APRAC-1A are included in this program: the synop-
tic, climatological, and grid-point models.
21
-------
SUBROUTINE INDAT
Reads in basic input data.
Run information (run type and output form).
City information (population, latitude, center location, and gasoline consumption rates);
Traffic information (hourly frequencies, peak hours, and car speed code);
Receptor information (no. of points and locations);
Starting and ending day and holiday information.
SUBROUTINE BASIC
Computes basic quantities used throughout program.
SUBROUTINE LINKS
Reads in and stores historical primary traffic link emission and secondary grid emission data.
ITYP = 2
CLIMATOLOGICAL
MODEL
SUBROUTINE STORE
Computes and stores the C/Q values for nine upwind segments.
five stability classes, and seven mixing depth categories.
Computes and stores the Q values for nine upwind segments,
36 wind directions, and NPT receptor locations
{The computations of C/Q and Q are made by calling subroutines CALXOQ and CALQUE.)
SUBROUTINE RAOBHMM
* Reads in the morning raob and mm-max surface temperatures.
Determines the mm-max mixing depth.
TA-8563-1290
FIGURE 8 APRAC-1A FLOW CHART
22
-------
SUBROUTINE SFCOBS1
Reads m first hour of surface observations (0000 LST)
SUBROUTINE SFCOBS2
Reads in the remaining surface observations for the day
SUBROUTINE EXTURB
Determines the extraurban concentration.
HOUR LOOP 1
DO 200, I - 1 TO 24;
SUBROUTINE MINWIN
Assigns wind speed a minimum value of 1 ms~ and uses the last observed wind direction
for calm conditions
SUBROUTINE STABIN
Determines the stability class.
SUBROUTINE DEPTH
Determines the mixing depth. Also reads in at sunset the next day's raob. min-max
temperature, and first hour of surface data (by calling subroutine RAOBHMM and SFCOBS1).
TA-aS63-129b
FIGURE 8 APRAC-1A FLOW CHART (Continued)
23
-------
Yes
ITYP = 2
CLIMATOLOGICAL
MODEL
No
SUBROUTINE LOCXOQ
Locates in storage the precompiled
values of C/Q for each of the
9 segments.
SUBROUTINE CALXOQ
Computes the C/Q values for each
of the 9 segments
POINT LOOP
DO 190 N = 1 TO NUMBER
OF POINTS TO BE TREATED (NPTI
Yes
ITYP « 2
CLIMATOLOGICAL
MODEL
7
No
SUBROUTINE LOCQUE
Locates in storage the precompiled CO
emission values (Q) for each of the
9 segments.
SUBROUTINE CALQUE
Computes the CO emission values
(Q) for each of the 9 segments
SUBROUTINE CALCON
Computes the background CO concentrations at one of the selected city locations, based on
the equation
9
C - £ (C/Q.) 'Q
No / ISM
STREET MODEL
TA-8563-129C
FIGURE 8 APRAC-1A FLOW CHART (Continued)
24
-------
To Point Loop
To Hour Loop
SUBROUTINE STREET
Computes the street profile of CO concentrations.
SUBROUTINE PPDATA
Writes and punches out the desired results.
ARE THERE MORE
DAYS OF DATA
TA-8563-129d
FIGURE 8 APRAC-1A FLOW CHART (Concluded)
25
-------
The synoptic and climatological models use the hour-by-hour values
of meteorological and traffic input data to make hourly CO computations
for a limited number (one to ten) of receptor locations within the city.
In the climatological model the major computations are made once and are
then stored in data arrays. The climatological model is more efficient
than the synoptic model when the number of CO computations per receptor
exceeds about 50 (or about two days of hourly computations). The synoptic
model is useful in an operational or evaluational sense, whereas the
climatological model is designed primarily as a tool for use in street
and freeway planning operations.
The grid-point model permits CO computations for as many as 625
points spaced selectively throughout the city. If this model is used,
the computations are restricted to one specified hour of the day, and
the street profile computations are automatically bypassed. Thus, the
grid-point model has limited use, but it is valuable for portraying a
horizontal pattern of CO concentration over an entire city, which would
be closely representative of background or rooftop values.
B. Data Input
Various basic information is specified on cards denoted as A through
M, as given in Table 3 and also in Appendix C. Card A is used to specify
values for ITYP, ISM, NPT, IOUT(1), IOUT(2) and IOUT(3). The run type
is determined by giving ITYP a value of 1 (synoptic), 2 (climatological),
or 3 (grid point). The street model is used if ISM is set equal to 1
and not used if it is set to 0. Also, if line-printer, punch-card, or
magnetic tape output is desired, the corresponding variable, IOUT(1),
IOUT(2), or IOUT(3), is set equal to 1; if not, it is set to 0. The
number of receptor points (NPT) for which concentrations are to be cal-
culated is also specified on Card A. However, the receptor locations
and information pertaining to the streets on which the receptors are
26
-------
CB <
B Q
§ 5
g J
H a
O
'*.
Z
B
a.
z
c
0
•H
•P
C.
CJ
(U
Q
cu
CO
w
CO
O
D.
^
Cfl
•P
E
0>
3
n
^
Ul
P
C
rH
0
.Q
CO
(H
O
fe
E
3
O
O
•o
CO
U
TJ
QJ
CO
OJ
f*
(U
-°
0
fl) Q) (0 "1
•0 T3 T3 P
o •- o o c
E -H IH E E H
(U O O
• • Q CO O CO QJ 11 +•* 3
rHTJCfgUCJCUrH 3 •• D.
O 0 H -H in ft O O. -P P
PEbDP'OPP-P -P 3 p
« o c c in in a 3 a o
O O rH -H -,-( V O P
•H -H O O QJ 4) U 3 Q)
"DPPD.rH-O'OO (H O Q.
C G. CO | Q) 3 £ h Qj a
•P. O E TI TJ rH rH -P J3 P
c .H r( o o u «H c tno v> i tn
a en £ be * * 5 ° a. ° ** |§!£>£C>>
P II H II QJ II II 0) t It 11 1 II II P>»O 0 0
"• H H
^ O r-H_^4_O O O
1 1 1 III
1 1 1 III
~ * ^
* * i I I
~ w jz; www
O O O O O O
•-* «-l rH] r-J rH rH
»-• rH r-( w IH rH
O O O O O
o CN n Tp m (o
^ 11 iii
1 (-1 i-H
I-H
in
O
O
V
9
i
o
o
L4
-H
g
O_
O
^
p
x
"1
o
ft,
W
O
r-f
£
W
(0
c
o
-H
a)
0
rH
>,
X
O
a
0)
u
(£
^^
c
M
rH
O
>,
ea
•H
p
P
i
rH~
II
rH
a
o
^H
^
G
0
4J
HI
U
U
•U
cn
o
IH
0
C
0
•H
-Direct
o
r~
rH
o
f*)
0
o
£
o
c
Cfl
a>
V
u
be
E-
co
•*
7?
O
tu
O
n
t
«
•e
P
a
g
n
fH
O
recepi
(D
X]
-P
C
0
^
T3
h
: of vehicles
Numbei
o
o
8
o
A!
(d
'O
rH
a>
Ct,
01
0)
rH
O
•H
0)
§
u
•a"
O
rH
O
1
rH
CO
tor street
recepi
*
a
•a
:or street wi
a.
cu
o
0)
a:
-T
«
Q
A
IA
0}
s
^
S
^
«
o
r-t
O
1
?
V
V
£
•p
n
!•
o
+*
o.
0
*•
o
1
c
11
T!
O
U
O
Traffi
*
00
0
*••
IH
13
^,
g
U
t— i
O
rH
tH
O
1
,— J
H>J
0)
(1)
•p
U)
o
4-t
O.
0)
0
rH
3
O.
C
U
.H
4H
01
•o
01
43
o
h O
P. B
SH -H
0 -P
at
0 ro
•H ^
tfl CJ
IH 0)
>
CQ
•P -P
27
-------
•o
HI
3
c
0
H
p
a
•H
0
0)
p
CD
CO
rH
CO
O
•H
P.
EH
cn
-p
•H
E
, •]
0
3
ca
»
+j
•r(
C
^
ri
CO
B
0
a
e
3
rH
O
U
•o
h
ca
o
^
o
•H U
"H -H
«H > -U
r| r-. 0
d r< ca
•O ca +J
C B a
0 -H C
Or. O
0) a u
vj
+j » >> E>1 r*l CO 3 3
•H -H tH -H CT ft W
O O O U CD > >
in
o
o
CM
m"
• o m
CO CO rH • , r-
• CO CO O O
00 • Si • O O
CO N rH in ["- 1
m o
§000
rH r-H rH
0 O 0 O
+i +j *j v O O
o o o o A V
o
•O CD
3 rH
£ a
•rl 0
ca a
rH
Oi 0 B C
CD H CD
t* rH rH CJ
60 rH O ,,
CD -H • 0 II
Q E O 0, ii
r
5 ^ frf 3 ^ g
W & X U (i &<
O •
O O rH O O O
-• rH fa rH rH rH
Em fe CM fc fe fc
o o o o o
M||l ||
| rH rH rH rH rH
rH rH CM ^ IT) CO
U
c
O 01
£ o
O. •*-» CQ
e o 4»
3 >
C -P
O rl
o ca T3
rH CO
(Dp O
C & rl
•H C
r-l CO -P
O «Q
in co be
M rH -H
be 0)
a> o o
£ -*-1
o a
bC H 3
C -P
rl U JH
> C 0
H 3 «M
be (H
cn
cn ca -o
•a 4»
r-. in O
ca m a.
o w
a>
£ ca A
H ** <->
C7 rH
o it on
l~l '"•' "^ ^" '
0 0
A A
tn
c
o
f— t
rH
CO
bi
O
CO JH
c ca
O CD rH
rH £
H CD -H
SO. E
x-. -4 Ot>
rH M W rH
0 0
00 00
o o
00 QO
rH rH
O3
n
•D
C
n)
rH
P W
rl U
3 -H
O «H
ca
o p
ca
ca
L« a
a» «M
.c o
be o
c
-P TO i<
CO CD 3
o a o
H >C
•D CO rH
-rl "- O CO
H .C 0-
cn a
•o >> ^* i
ij-t
at .c n n
-
U O H
frj r" Tl
(H u ca
ca o
S C
•H ca
bfl O r.
01 «M >>
•a
>, i*
0) rH
J! '« O
to oo
CN
r- n
o
• r-t
O ^~•
rH
o
p
o
c
o
-H
U
ca
rl
[>H
CM
^ rH
rH II
O
§8
i
i
rH
N
O CO
« U
rH^ *O
^ C
O c»
in
a
L,
ca
rt
o
-H
0)
o
X
CD
CD
>
*
CO
If)
in
c
d) ca
a *,
W ^
00
to ii
o
O •*. '
o
o
c
o
•rl
•p
o
rH
^s
CO II
&4 rH
o
rH
fe
S
1
tN
ffi CO
^ ffi
^ -a
v^ C
X CO
OJ
•o
o
6
^j
0)
0)
tn
-P
Cfl
(^
o
a.
cv
u
X
(U
(D
>
s
ca
(A
ca
Or-
el
W
ID 00
r- n
o
* rH
O w
o
o
c
o
•p
u
ca
h
s
rH rH
(O 1)
£ rH
O
rH
S8
i
rH
CM
rH C*)
•. I-H
r7 *D
^r« C
M CO
C
CO
to
ca
•D
rl
•p
ca
CQ
f->
o
«H
p
a
CD
u cn
X CD
CD O.
•p"
CD
0 ca
J2 0
CO rl
W rH
CO rH
ca
01
ca o
W *M
X
to n
0
o -^
o
•p
o
c
0
•H
•p
u
ca
rl
rH
CO II
O
rH
§8
8
rH
CN
rj C*3
^ Hj
rH "O
^^ Q
HS ca
•a
c
ca
cn
ca co
1 1
h T3
O CO
!H
o
* ,
ca ca
T3
ca o
CO rC
00
X
rH ||
o ^^
rH
O
0
c
0
+•»
u
ca
rl
,—v M
rH
CO II
B M,
O
rH
Sg
a
i
rH
N
X CO
-S
rH T3
^ c
« a
28
-------
c
c
•H
a
H
c
CO
0
3
CO
CO
CJ
•H
a
H
tn
H
£
•H
CD
CO
>
cn
H
C
0
.Q
-H
CO
s
o
c
s
o
Q
•o
0
•
S-*. *"-N
CD CN
ca
•a CD
be co
c -o
•H
C bo
c c
•>-( -H
bo *D
CD C
CO W
t£ to
00
O O
0 t>
W
i
1 M
1 Al
O CJ
cj «
a a
o
W u
f P
a a
O 0
I-H r-<
O
O N
<-t 1
1 rf
l-l l-l
J
CD
CO
•o
be
c
•rl
C
•H
bC
(D
P
^
O
«H
0)
CD
0)
+j
o
CO
Q
^
CO
T3
tn
3
£_,
•W
[^
0
r^
01
1
o
ra
D
g
)_(
O
o
CO
1
CN
^j
G
•g
Q.
1
^
^
be
c
•H
•o
CD
OS
CD
•M
CD
CO
o
•p II
OJ p
•D H
*M ^-^
O
^ 0
3 -0
O O
K E
X
oa
o
rH
3
S
1
t-t
o
o
1
r-l
CO
U]
ca
•o
r-t
0
o
4)
XI
E
3
55
r-4
O
Al
CD
3
S5
X
O
O
1
TP
in
ca
TJ
r-i
0
o
x:
J
h- -^/
0
01
H P
Q Q
CD
-H
CO
Q
HH
S
o
o
i
«
j~
Q ii
s -
29
-------
located are specified on B cards (one card for each location). Informa-
tion pertaining to the city, such as population, latitude, city-center
location, secondary emission amount, and emission constants (a and p),
is given on Card C. Gasoline consumption rates for 16 different angular
sectors are given on Cards D(l) and D(2) . These data are used to compute
the extraurban contribution to the CO concentration. Vehicle-speed,
peak-hour, and hourly traffic distribution information is specified on
Cards E through K, as given in Table 3. The beginning and ending days
of the period to be analyzed and the number of holidays during the period
are specified on Card L. The dates of the holidays are given on the M
-ards, one card for each holiday.
The primary and secondary traffic data are specified on the N and
0 cards, as shown in Table 4. The primary traffic data are given by
individual links, with a separate card for each link. The read-in of
the traffic link data cards is terminated when a card with a 9 in Column
1 is detected; however, the number of these cards should not exceed
1200. The secondary traffic data are given for grid points spaced 2 mi
apart. There is a separate 0 card for each grid point; each 0 card con-
tains the x and y location of the grid point and the percentage of the
2
city's total secondary emission that falls within the 4-mi area surround-
ing the point. The total amount of secondary traffic is set equal to
the percent of the total primary traffic (CLE) that is specified on Card
C. The read-in of this secondary emission data is also terminated when
a card with a 9 in Column 1 is detected.
Meteorological data are specified on Cards P, Q, and R, as shown
in Table 5. Card P, a header card, contains the city's name, the date,
and the daily maximum and minimum temperatures. This card also contains
an indicator of whether it is a day when there is a time change. The
program initially assumes local standard time (LST). Thus, if the first
day is during a daylight saving period, Card P should always include
30
-------
V
o
n
H
w c
§ s
H c
g V
z -o
rt Q
c
o
•H
a.
u
o
Q
3
CO
>
i-l
CO
CJ
a.
£
-t-j
•H
E
•H
^
CD
3
ca
in
.-I
c
i-i
I
C.Q
CO
E
O
fM
c
E
3
O
O
•o
ca
U
0)
43
be
c
-H
•H
be
tn
•o
ca
o
o
in
&
3
C
0)
+J
,_<
c
H
0)
•a
c
c
<
CN
03
^
in
01
CO
m
o
o
-H
o
H
E
O
O
£
^
rH
X
0
in
N
Q
I
rH
1-4
x
x
c
•H
I-H
o
+J
o
c
o
0
n
o
0
c
0
CO
o
0
^
^
It-
CD
•o
c
I-H
^
1-4
II
M
O
T3
C
O
O
O
43
0
c
o
•H
TO
CJ
O
c
>> i-H
^ at
x 43
o
en
^
m
CO
O
0
o
•H
o
o
J
CN
s
N
^
O
tn
PM
O
CO
1
r-l
CJ
C
r-l
C
0
CO
^
a
0) XI
rH 6
CJ 3
H C ^
43 C
0 O -H
> *o -«
o
'K O
tr* [JL( J
i> w _i
o
«j in m
ft, I-H (— <
^ ,— t (O
CO 'f ^
1 1 1
CO CO *f
0
c
•H
— (
(A
•D
0
0
h
a>
XI
3
0)
-*V
•H
C
H
HH
Or-
"D
C
H
C
<
o
o
in
^
o
o
CO
"
m
o
o
-*->
o
-H
E
O
o
^
t~t
X
I-H
o
r-l
CN
O
CO
1
i-H
I-l
*
i
•a
0)
o
CO
a
to
(Q
c
o
a
•H
IM
b£
«M
O
C
0
H
CO
u
o
^
^
h
O
•o
c
r- 1
^
rH
II
Ul
CO
IH
0)
C
H
0)
i-H
g
1
-p
ca
JH
CO
t3
C
o
(A
ca
•H
o
en
>i
H
U
0)
-*J
o
c
o
o
(H
0
CN
O
o
0
o
c
0)
^
£
N
M
O
rH
t-t
O
1
rH
CO
(I)
-*->
C
•H
43
-P
C
'C
^
y
o
(0
•H
+->
ca
43
•H
'H
ca
v
c
g.
•H
a
u>
c
•H
•o
c
3
O
rH
3
in
CQ
2
CO
N
•H
|
*
31
-------
c
o
T-t
+J
o.
H
escr
Q
O
P
eo
ca
o
•H
H
tn
P
i-i
(9
>
to
•H
C
1-4
s,
CO
-H
CO
e
£
B
B
3
i-i
O
U
•o
CO
u
•a a>
» § 5
3 -H 4->
O +>
•H cn ^
> -a MX
OJ rl B r-1
a co TH
B a in -H
2+> -rl
0) *J C
CO CO «H .C "H
T3 13 O M
OJ *> .H TJ
O O B U 'S. B
^ 01
CTi C rH CO O CO
r-l 3 3 OJ 0) -^
•M CO CO -H 4-» 01 t; E
p v a> ex B a -r>
M a a
3 B B o -i £ a -o
< 9 Ci) X CJ co ^
E CD O T3 O CO -O
C8CNBBWB C
a ^ 3 3 co II u co
E E 0 II v
XOJ>HH •-! CO 01
co t- « m r-i
-
o
1 1 1 1 ->
1 1 1 1 1
0) b
•P 3
I co pu h 0
1 Q o o I
< H E- H 0
H < X fe CO
co Q a i-i o
r-l r-l S * r-l
0 0
r-i r-i m in m
o in o in
O CM M CO CO
r-l 1 1 1 1
1 r-l -1 ID -1
r-l r^ W N CO
0.
rl CJ
O £
0)
4J • CO
bo co en TJ
B TJ r-l
•H eg a
> bo > E
•H B V
M .ri rJ O
•o o
01 c -o in
TJ 3 rl
IH O M O
03 Cn TJ J3
U C -H
"M h *• E
O 3 01 01
„ 8*f
tp OJ C
^ rl CO CO TJ
BO. CJ O
B TJ O CO h
q co a
01 3 O -H B
+> .rt B OJ
.ri gj IH h tp
B rl ^1 CJ £
•H 3 e *>
IH v bo m
01 CO ^ C to
"D rl 01 YH J3
B 01 |
•rt a £ Tl TD
B -r> 0 h
B CJ O OJ CO
< *> O h O
r-l
O
in • n
«• 8«
i m ^
§
i
1 Al
to
rl
at
O
^
U -H
o 55
r-t M
E- ft
r-t r-l
[».
M 5
1 1
co o
a>
13
B
<-<
ft
M II
5 M
C
o
•H
ca
e
o
n
a CM
o
a
H§
n o
1°
0 B
ca ^
CO O CD
CJ r* CO
CO 0
*•* v +> o
0 0 ! 0 Al
01
CJ
rl W
DDrl
•S e
-------
an indication of a +l-hour change from standard to daylight time. Cards
Q and R contain the 1200 GMT radiosonde (raob) data and the hourly sur-
face data. The raob data include both significant and standard levels
from Card Deck Format 505 of the National Climatic Center (see footnote
of Section II) . There is a separate Q card for each level, and the last
card must always be for the 500-mb level, since this terminates the
reading of data. The surface data are specified for each hour of the
day, beginning at 0000 LST and ending at 2300 LST; there is a separate
R card for each hour. The program reads in a separate set of the P, Q,
and R cards for each day until the CO computations have been completed
through the specified ending date. The P and Q cards and the first R
card (surface data for hour zero) for the day following the ending date
must be included at the end of the data deck.
C. Data Output and Program Capabilities
The calculated CO concentrations can be output on the line printer,
on punch cards, or on magnetic tape, depending on the values assigned to
IOUT(1), IOUT(2), and IOUT(3). An example of the synoptic model's CDC
6400 line-printer output is shown.in Appendix Df (This output was pro-
duced by the data input shown in Appendix C). The line-printer output
for the climatological model and that for the grid-point model have basi-
cally the same format; however, one would normally desire only punch-card
or magnetic-tape output when using these two models. The computed CO
concentrations (ppm) always include the background value (C ). When ISM
is set equal to 1, the low-level (3.65-m) and high-level (22.8-m) con-
centrations are calculated for both the right and left sides of the street
(facing in the street direction specified on Card B). The outputs also
include the hourly values of cloud amount (tenths), surface temperature
(°K), wi '." direction (degrees), wind speed (ms ), stability index, and
mixing depth (m).
33
-------
Punch-card or magnetic-tape output from this program can be used as
the input data into the computer to produce graphic displays. Graphic
displays of the various outputs are very useful. For example, in evaluat-
ing the performance of the synoptic model, extensive comparisons were
made of the calculated hourly concentrations (C ) with the values observed
b
at CAMP stations in Chicago, St. Louis, Denver, Cincinnati, and Washington,
B.C. One of the graphic comparisons for St. Louis is shown in Figure 9.
Also shown in Figure 9 are the graphic displays of the mixing depth,
stability index, cloud amount, and wind direction and speed. Frequency
distributions of CO concentrations for various different times of day at
the St. Louis CAMP station are shown by graphic display in Figure 10.
l.iese results were obtained from the output of the climatological model.
An example of computer contouring based on 625 grid-point values is
presented in Figure ll(a) . The road or link network that was used is
shown as an underlay. Figure 1Kb) illustrates the telescoping grid or
"zoom" capability of the grid-point model; that is, the grid spacing was
reduced by a factor of ten to depict the detailed concentration pattern
of downtown St. Louis.
As noted earlier, the APRAC-1A program was originally written in
CDC FORTRAN Version 2.1 for use on the SRI CDC 6400 computer. It has
been modified for use on an IBM 360/50, but further modification of the
program may be necessary if it is to be used with other types of compu-
ters. The central processing times on the CDC 6400 computer for the
synoptic, climatological, and grid-point models are given in Table 6.
The storage requirement of the program is 45,000 words of memory, although
this would need to be increased for cities larger than St. Louis.
34
-------
o
UJ
Z
o
o
o
o
Q.
UJ
O
18
16
14
12
10
8
6
4
2
0
400C
2000
i I
OBSERVED
CALCULATED
O
P 400
o
a
UJ
UJ
a.
§
•o .c
.= I/I
o o
co o
zoo
0
20
10
0
10
5
0
\!
^ <-> 0 20 40 60 80 100 120 140 160 HOURS
[— MON .|. TUES 4— WED -4- THURS-4— FRI -|- SAT -|- SUN —j
ST. LOUIS, MO. DATA (19-25 OCTOBER. 1964)
TA-7874-49R
FIGURE 9 METEOROLOGICAL INPUTS TO THE MODEL WITH OBSERVED
AND CALCULATED CO CONCENTRATIONS AT THE ST. LOUIS
CAMP STATION. The calculated concentrations shown here are
generally underestimates because they do not include the contributions
from the local street.
35
-------
«>u
_, 50
>
IT
£ 40
z
30
in
_j
o
\ 20
H
Z
UJ
o 10
cr
Q.
0
60
NTERVAL
4> 01
o 0
in 30
05
?20
z
o 10
UJ
a.
0
60
INTERVAL
* 2
0 0
in
3 30
o
PERCENT /
_ ru
ooo
1 ,.,.,.. 1 ..,..,., , • .
-
-
—
! rr-r-
(a) 0800 HOURS
•
—
i 'rrri.i i
i ,.,...., , . . , , , . .
-
-
-
•
i . •
i .iii
(b) 1200 HOURS
-
-
•
—
1 t t . . . i i
1 . . | • •••! | . • | . . t-| | i .
-
-
-
.1 0.2 0.5 1
CO CONCENTRA
-------
12
10
8
6
4
2
0
-2
-4
-6
-8
-10
-12
1500-1600 COT
15 OCTOBER 1964
WIND 310"/1.6 m I-1
MIXING DEPTH 1670m
UNSTABLE
-12 -10-8-6-4-2 0 2 4 6 8 10 12
DISTANCE EAST OF CAMP STATION — miles
TA-7874-26
(a) 1-MILE (1.6 km) GRID SPACING
-1.2
-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1.0 \.2
DISTANCE EAST OF CAMP STATION — miles
TA-7874-24
(b) 0.1-MILE (0.16 km) GRID SPACING
FIGURE 11 CALCULATED ST. LOUIS CONCENTRATION
PATTERNS FOR TWO GRID SIZES
37
-------
Table 6
APRAC-1A CENTRAL PROCESSOR TIME (SECONDS) ON THE SRI CDC 6400 COMPUTER*
(For computing CO concentrations in St. Louis)
Process
Compilation
computation
Basic
computation
(per recep-
tor per
hour)
Model
Synoptic
One
Receptor
(NPT = 1)
10.5
6.1
-1
1.14 X 10
More Than
One
Receptor
(NPT > 1)
10.5
5.9
-1
2.25 X 10
Climatological
10.5
5.9 + 7.5 X NPT
-3
3.7 x 10
Grid-Point
10.5
5.9
-1
2.85 x 10
This computer has an add time of 1.1 microsecond and a cycle time of
1.0 microsecond.
38
-------
REFERENCES
Coventry, D. H., and R. E. Ruff, 1972: personal communication.
Johnson, W. B., W. F. Dabberdt, F. L. Ludwig, and R. J. Allen, 1971:
"Field Study for Initial Evaluation of an Urban Diffusion Model
for Carbon Monoxide," Comprehensive Report, Contract CAPA-3-68(1-69),
Stanford Research Institute, Menlo Park, California, 240 pp.
Ludwig, F. L., W. B. Johnson, A. E. Moon, and R. L. Mancuso, 1970: "A
Practical Multipurpose Diffusion Model for Carbon Monoxide," Final
Report, Contracts CAPA-3-68 and CPA 22-69-64, Stanford Research
Institute, Menlo Park, California, 184 pp.
Ludwig, F. L., and W. F. Dabberdt, 1972: "Evaluation of the APRAC-1A
Urban Diffusion Model for Carbon Monoxide," Final Report, Contract
CAPA-3-68(l-69), Stanford Research Institute, Menlo Park,
California, 167 pp.
Ott, Wayne, 1972: "An urban survey technique for measuring the spatial
variation of carbon monoxide concentrations in cities." APCA
Paper No. 72-17, 65th Annual Meeting of the Air Pollution Control
Association, Miami Beach, Florida.
39
-------
Appendix A
SRI CDC 6400 VERSION OF APRAC-1A COMPUTER PROGRAM
PROGRAM APRAC1A(INPUT,OUTPUT . PUNCH ,TAI>E41)
C
C AIR POLLUTION DIFFUSION MUDEL
C
C THIS PROGRAM COMPUTES THE CU CONCENTRATIONS AT VARIOUS CITV
C LOCATIONS bASED UN TRAFFIC AND METEOROLOGICAL INPUT DATA.
C
C SET ITYP EQUAL TO
C 1 FUR SYNUPTIC MUDtL
C 2 FOR CLIMATOLOGICAL MODEL
C 3 FOR GRID PT MODEL.
C
COMMON/CITY/ ITYP , 1SM,IOUT{IOJ
COMMON/CPTS/ NPT,XPn&25),YPT(625)
COMMON/CSTM/ VCARl10),AST<10J,WWST(101,ICST(10),NOCD(10J
COMMON/CHANG/ HANG,ACV
COMMON /CORD/ XXC.YYC
COMMON /DAY/ IOW.NH,IDAHO,IHW,IDWT
COMMON/CSTA6/ SP130.CP1UO
COMMON/CFUEL/ FUELt16),NHUL(10)
COMMON/CPFC/ PF1,PF2,P01.P02,S(9),KT(24)
COMMON /MIX/ ISTA.1YR.1MO,I DA,DAY,IDS,I HO,1ST,1UATt,IDATE I,IUATE9
COMMON /HMM/POP4,P(25),T(25),MAXT,MI NT,HMAX.HMIN,NSL,SAP
COMMON /SFC/ IHRI2«»,1CHTI2*),IHDI2*),WSI2*J . ICLDI 24) , I TEMPI 2<»)
COMMON /CAR/ PT12(24),PT3<»(24),PT6124),PTSATC24),PTSUN(24)
COMMON /UUE/ ULUO) ,I3M(10) ,0.12(10) ,U3*(10) ,U6UO) ,XOU( 10)
COMMON /PNT/ IC,RT2,IJS<24),HS(24),HR(168),CCAL(1200)
COMMON/CSAV/ J5,L7 ,1 9,K36,HM1,SXOQI320J,SQ12(32^J),SQ3*(32^0),
1 SQLI3240)
COMMON /CLOQ/ VCL (2500),VCM(2500),Cb< 3),iB13),NRM,NCL,MRM,MCL,CLE,
IGSP,RGS,£OX
COMMON/CLDAT/II-K1 1200) , XII 1200) .X2U200) ,YH1200i , Y2« 1200 ) , E (1200 )
COMMON /ALLC/ YSt351),IYNl351)
1,XR(10),YR(101,ARtlO),EFAC{9),ACR,LR*.LR5.LR9,LK10,AG,YC1,YC2,XKC,
2YRC.XC2.8Y2,BY1,612,811,CRI,IT,IJ,XG,YG,SA,CA,NN,C«,LRL
I FORMAT 1IH1,55X*SYNOPTIC AIR POLLUTION MODEL*/)
2 FORMAT!1H1,52X,*CLI MATOLOG1CAL AIR POLLUTION MODtL*/)
3 FORMAT (1H1,53X,*URIO POINT AIR POLLUTION MODEL*/)
C
C IN^UT THE BAMC PARAMETERS.
CALL INDAT
IF (ITYP.EQ.l) PKINT 1
IF ( IIYP.fcO.2) PRINT 2
IF (ITYP.fcU.3) PKINT J
C
C COMPUTE THE EMJiilON FACTOR AND OTHER INITIAL UUANTlTItS.
CALL UAilC
C
C INPUT THE MAJOR AND MINOR TRAFFIC EMISSION DATA.
CALL LINKS
C
C PRECOMPUTE AND STORE ALL THE X/U AND U VALUES (IF 1TYPE = 2).
IF (1TYP.E0..2) CALL STORE
C
C INPUT THE HEADER CARD ANU RAOB OBSERVATIONS AND
C DETERMINE THt MIN-MAX MIXING DEPTH FUR THE FIRST DAY.
41
-------
CALL RAOBHMM
C
C INPUT THE FIRST SFC DBS CARD.
CALL SFCOBS1
C
C INPUT THE REMAINING SFC OBS FOR THE DAY.
100 CALL SFC08S2
C
C DETERMINE THE EX1RAURBAN CO CONTRIBUTION.
CALL EXTURB (EXTRAO)
C
C BEGIN HUUR LOOP
DO 200 1*1,24
C
C IN GRID POINT MODEL THE COMPUTATION IS DONE FOR HUUR IHD ONLY.
IF (ITYP.EQ.31 1MHO
C
C DETERMINE THE MINIMUM MIND SPEED AND DIRECTION.
CALL MINWIN IDS IWI'IMDd)
C
C DETERMINE THE STABILITY INDEX.
CALL STABLEII.SAS.U)
C
C DETERMINE THE MIXING DEPTH AND AT SUNSET INPUT THE HEADER
C CARD. THE RAOB DATA, AND THE FIRST FOUR HOURS OF SFC DATA
C FOR THE NEXT DAY.
CALL DEPTH tI,SAS,HT)
C
C DETERMINE THE X/U VALUES USING THE GAUSSIAN AND BOX MODELS.
IF (ITYP.EQ.2) 140,150
C
C SELECTS APPROPRIATE X/U VALUES FROM STORAGE.
140 CALL LOCXOOIHT,IJJi GU TO 160
C
C COMPUTE APPROPRIATE X/Q VALUES.
150 CALL CALXOQIHT.UI
C
C BEGIN POINT LOOP
160 DO 190 N=l,NPT
C
C DETERMINE THE ti VALUES FROM TRAFFIC EMISSION DATA.
IF UTYP.EQ.2J 170,160
C
C SELECTS APPROPRIATE ti VALUES FROM STORAGE.
170 CALL LOCJUE(1HI,N)t GO TO 189
C
C COMPUTE APPROPRIATE Q VALUES.
180 CALL CALQUEI1WI,N)
C
C DETERMINE THE TOTAL BACKGROUND CO CONCENTRATION.
139 CALL CALCON (I,N,HT,EXTRAti)
C
C DETERMINE THE STREET PROFILE CO CONCENTRATIONS IF ISM IS
C NOT EQUAL TO ZERO.
IF (ISM.Nt.O) CALL STREET(I.N)
C
42
-------
190 CONTINUE
IF (ITYP.EQ.3) GO TO 210
200 CONTINUE
210 CONTINUE
C
C PRINT AND PUNCH THE OtSIREO PARAMETERS.
CALL PPOATA
C
C DAY OF THE WEEK CHECK.
IF (NH.EU.O) GO TO 300
IF (IDAHO.Eg.0) GO TO
IOAHO=0$ NH=NH-1
295 IF uoATEi.EQ.NHOLiNH>)
300 IDW=IDW*1
IF ( IOW.LT.S) GU TO 400
IUW=1HW-1» IC=0
C
C CHECK FOR IHt END OF THE TIME PEKIUD.
<»00 IF UOATE1.LE.IOATE9* GU TO 100
500 STOP * END
43
-------
SUBROUTINE INUAT
THIS SUBROUTINE READS IN BASIC DATA
COMMON/CITY/ ITYP,ISM,IUUTi10)
COMMON/CPTS/ NPT,XPT{625J.YPT1625)
COMMON/CSTM/ VCAR(10).ASTC10),MWST(10).ICSTt10),NUCD(10)
COMMON /CORD/ XXC.YYC
COMMON/CSTAB/ SP180,CP180
COMMON /DAY/ 1 DM,NH,IDAHO,IHU,IDWT
CONMON/CPFC/ PFl,PF2,P01,P02tSm,KTI2*)
COMMON/CFUEL/ FUEL<16),NHOLl10)
COMMON /MIX/ lSTA,irR,lMO,IDA,UAY, 1DS,IHO,IST,IDATE,IDATE1,IDATE9
COMMON /GRAF/ 11(7) , I A(2*),1B(16,8).1V<5,8),ICS
COMMON /HMM/POP*,PI25),T(25),MAXT,MINT.HMAX,HMIN,NJ>L,SAP
COMMON /CAR/ PT12(2*),PT3*(2*),PI6<2*),PTSAM24).PTSUNI2*)
COMMON /PHI/ IC.RT2,IJS(2*),HS(2*),HR(168),CCAL(1200)
COMMON /CLOQ/ VCL12500),VCM(2500),Cb(J),SB(3),NftM,NCL.MRM,MCL,CLE,
1GSP,RGS,ZQX
COMMON/CN2*/ N2*
1 FORMATI8HO)
3 FORMATI8F10.*)
* FORMAT(!jFlO.*,2I10)
5 FORMATI2F10.*)
READ 1,ITYP,ISM,NPT,I I OUT(I),I = 1,5)
IF IITYP.GT.2) GO TO 15
READ *,(XPT(I),YPTII).AST(I),VCARtI).MMSTtI),ICSTiI),NOCO(I),
2 1=1,NPT)
GO TO 17
15 ISM=0
READ 5,UPT(ll ,YPTII),I = l,NPT)
17 READ 3,SLAT,POPS,XXT,YYT,CLE,PF1,POI
3,(FUEL( 1,1=1,16)
3,(S( I),
1, UT( i)
3,JPT12(
3,(PT3*(
3,(PT6(1 ,1=1,2*)
3,(PTSATU),1=1,2*)
3,(PTSUNI1),1=1,2*)
1,IDATE,IDATt9,IDW.lHD.NH
INH.GT.O) READ 1,(NHOL(I),1=1,NH)
(ITYP.LT.2.AND.NPT.GT.10)
19
20
30
1 = 1, 2*,)
>.]
,1 =
READ
READ
READ
READ
READ
READ
READ
READ
READ
READ
IF
IF
N2*=2*$ IF (ITYP.LT.J) GO TO 19$ N2*=l* 1SM=0
DO 20 1=1,9
IF (S( D.LT.1.0) Sll) = 99.0
XXC=-XXTt YYC=-YYT
DO 30 I=1,NPT$ XPTl1)=XPT(I)»XXC
YPT(I) = YPT( I )*YYC
POPS=POPS*1.0E6$ POP*=SQRTtSOHTtPOPS>)
IGS=IHM$ IC=IHW-1» 1DAHO=0* SAP=-1.0
ACR=0.017*533
PL180=ACR*SLAT * SP180=S1N(PL180) $ CP180=COS(PL180)
PF1=PF1/3600.0
RETURN* END
44
-------
SUBROUTINE BASIC
C THIS SUBROUTINE COMPUTES BASIC QUANTITIES UStU IMMUUUHOUl PMOliMAM
C
COMMON/CHANG/ HANG,ACV
COMMON/CPFC/ PFl,PF2,POi,P02,S(9),Kn24)
COMMON /CAR/ PT12<24),PT34(24),PT6(24),Pr»An24I.PrSt,yR(10),AR(10),EFACI9),ACK,LR^,L«5,LRV.L«lD,Ali.irCl.tfL2.XrtC.
2YRC,XC2,BY2,BYI,B12.BI1,CRI,IT,IJ,XG,YG,SA.CA.NN,CK.LIU
DATA ACR,CR»HANG,LRL/0.0174533.62.137,0.196349625,$/
C
C COMPUTES BASIC FACTORS FOR EMISSION MODEL
C KTtI)=PEAK OR NON PEAK TRAFFIC AT ITH HOUR
C PT(I)=FRACTION OF DAILY TRAFFIC DURING ITH HOUR
C EFACII1=EMISS10N FACTOR FOR ITH FACILITY TYPE
C S(I)=SPEE» FOR ITH FACILITY
C P01 AND PF1=INPUT CONSTANTS.
TN1125=TANlHANli) * RR=7.7671 S CKl = IO.O/CM * ACV-CRI*CRI«1.0£*
LR4=* I LR5=L«4+1 S LK9-9 $ LR10-LR9»1 S XP-0.«5»*POl
00 50 1=1,2*
IF (KT(I).Eg.U 45,50
<»5 Pfl2(I )=PT12( I )*XP » PT3*( I I=PTJ
-------
SUBROUTINE LINKS
IF
C THIS SUBROUTINE READS IN DATA FOR MAIN AND LOCAL LINKS.
C
COMMON/CRTS/ NPT.XPM625), YPT( 625)
COMMON/CSTM/ VCAR(IO),AST!ID)tWWSI<10),ICSI{10)»NOCUI10)
COMMON /CORD/ XXQ.YYC
COMMON/CHANG/ HANG.ACV
CQMMON/CLOAT/ IFKU200) , XI ( 1200) , X2 11200 J ,Yl ( 1200),Y2I 1200) ,6(1200)
COMMON/A'LLC/ YS(351),IYNI351)
1 ,XR(10) vYRUO) ,AR(10),EFAC(9) ,ACR(LR4,LR5,LR9.LR10. Aoi YCl, YC2, XRC,
2YKCrXC2,BY2,BYl,BI2tBIl.CRIiIT,IJ.XGiYGiSAfCAfNNiCRtLKL
COMMON /CLOU/ VCL(2500),VCM(2500),l.B<3l,SBm,NRW,NCL,MRW,NCLfC(.E,
IGSP.RGS.ZOX
DATA GSP,NCL.NRW/100.0,49,<»9/
1 FORMATIIl,9X,4F5.0tF6.0,215)
2 FORMAT(I1,9X,3I10)
3 FORMAT (1H ,59X,*NO. OF LINKS =*tI5,/)
MRW=NHW*0.5 $ MCL = NCL*0.5 * LT=NRW*NCL t ARI»1.0/«,SP*GSP*ACV)
RGS=100.0/GSPt DGS=,05*RGS* ZJX=0.01*RuS> CRAD2=YC1*YCI/O.96*1000
WL=0.0$ NCO'=0» L = l * A75=0.833*HANG» NCD1 = NUCO(1)
XG=XPT(1)$ YG=YPT(1)
CB(1)=1.0 » CB(2)=CU(3)=COS(A75)
S8(l)=0.0 $ SB(2J=SINJ-A75) » SB«3)=-SB<2)
DO 75 1=1,LT
75 VCMII)=VCLII)=0.0
C
C READS IN DATA FOR MAIN LINKS
C (LINK DATA STORED IN ARRAYS XI,Yl,X2»Y2,AND E).
IF (LRL.GE.LR9) 100,120
100 READ 1,M1,X1L,YIL,X2L,Y2L,VV,IF.LL
IF (M1.LT.9) 110,170
110 VL=VV*LL*O.Ot » SVL=SVL»VL t Et=fcFAC(IF)«VL » IFK(L)=IF
XD=X1L-X2L i YD=Y1L-Y2L
XI U)*X1L+XXC * YI(L)=Y1L*YYC
X2(L)=X2L+XXC t Y2(L) =Y2L»YYC
EJL)*EE*EE/(XD*XD*YO*YD) * L="L»1 $ GO TO 100
120 READ 1,M1,X1L,Y1L,X2L,Y2L,VV,IF,LL
IF IMI.GT.8) GO TO 170$ NCO=NCD*l
130 VL=VV»LL*0.01 $ SVL*SVL*VL * Efc»EFAC(IF)*VL
XD=X2L-X1L $ YD=Y2L-YIL
DD=XD*XO+YD*YO $ X1L=X1L»XXC * YIL=YIL*YYC
C
C THE MAIN LINK DATA ALSO USED TO COMPUTE EMISSION VALUES AT GRID POINTS
C (VALUES STORED IN VCM, NO. OF GRID POINTS = NRWXNCL).
IL2=USQRT(DU)*DGS» ZII = 1.0/ILZ * ZIL~EE*Z1I
YM=YD*ZQX*ZlI S XM=XD*ZOX«ZII
YJ=Y1L*ZQX-YM*0.5*MRW S XJ=X1L*ZUX-XM*0.5*MCL
DO 135 1 = 1, ILZ
VJ=YJ*YM» XJ=XJ+XM S JY=INTIYJ) * JX=1NTIXJ)
IF (JY.LT.l.OR.JY.Gt.NRW) GO TO 135
IF (JX.LT.l.OR.JX.GE.NCL) GO TO 135
OY=YJ-JY$ DX=XJ-JX* 1XY=(JX-1)*NKM*JY
DYl = l.O-DY $ 0X1 = 1.0-DX » IXYl»IXY»l ( IN=UY*NRW i INlMN + l
VCMIIXY)=ZIL*DY1*DX1»VCM(IXY» * VCM
-------
2 $ ILM1=1L-1 S XDD=XD/ILMl S ¥DD=YD/ILM1
IF (NPT.EQ.l) 142iUl
1<*0 IF (NCD1.EQ.NCD) NOCO(U=L
m EU) = EE*EE/DO t XKL) = X1L $ YltLI = YU $ JFMU = IF
X2(L)=X2L*XXC $ Y2(L»=Y2L+YYC S L*L*l * GO TO 120
1*2 XTL=X1L-XG* YTL=Y1L-YG
00 160 1=1 tlL
IF ( (XTL*XTL+YTL*YTL).LT.CRAD2) GO TO 140
150 XTL=XTL+XDD
160 YTL = YU+YDD $ GO TO 120
170 NN=L-l * rZS=CON=0.0
PRINT 3,NN
C
C KbAOS IN LOCAL EMISSION OATA AND TRANSFORMED IT ONTO A OR10
C (GRID VALUES STORED IN ARRAY VCD.
250 READ 2, Ml, IX, IY, 12
IF (Ml.GT.8l GO TO 275
CX=( I X*XXC 1*0.01* CY=( IY»YYC)*O.Ol
YJ=CY*RGS*MRW $ XJ=CX*RGS+MCL S JY=INTIYJt i JX=INT(XJJ
IF I JY.LT. l.OR.JY.GE.NRW) GO TO 270
IF ( JX.LT.l.OR. JX.GE.NCLJ GO TO 270
IXY=(JX-U*NRH + JY * IXYl = IXYH S IXYN=IXY*NRW » UYDX=UY*OX
IXYN1=IXYN»1 S VCL(IXY) = U *VCHIXY) » VCL ( I XY 1 J = U »VCLUXY1)
VCL1IXYN)=U +VCHIXYN) S VCL ( IXYNlls I i »VCLIIXYN1)
270 CONTINUE
GO TO 250
275 I2S=11S*4 t IF (US.EQ.OJ GU TO 280
CON = 0.01*CLE*SVL*ARI*EFAC(5I/US
280 DO 290 L=1,LT
VCMJL)=ARI#VCM(LI
290 VCL«L)=CON*VCL(L1
RHTU«N $ END
47
-------
SUBROUTINE STORE
THIS SUBROUTINE COMPUTES AND STORES X/Q AND Q VALUES (IF ITYP EO 2),
COMMON/CPTS/ NPT.XPTI625J, VPU625I
COMMON/CSAV/ J5,L7,I9,K36,HMI,SXOQ1320),SU12(3240).
1 SQL(32^0)
COMMON /UUE/ OH10J,QM(10»,m2<10) ,Q3^I10J,g6(10J ,XOg(iO)
J5=5$ L7=7* 19-9$ K36=36$ HMI=37.5
DO 100 J=1,J5$ jaS*(J-l)*L7*I9S HM=HMI
DO 100 L=lfL7$ LBS=(L-11*19*JBS» HM=HM»HM
CALL CALXOU(HM.J)
DO 100 1 = 1.19* IBS = ULBS
100 SXOQ(IBS)^XOOIl)
DO 200 N=l,NPTt NBS=(N-1)*K36*1,9
UU 200 K=l,K36t KBS=(K-l)*!9tNBS» KDR=K»10
CALL CALQUE(KDR,N)
DO 200 1=1.J9» IBS'1»KBS
SQ12llBS)=Q12tI)S SQ3*iIBS)=Q3V(II
200 SQL(IBS)-OLd)
RETURN* END
48
-------
SUBROUTINE RAOBHMM
C
C THIS SUBROUTINE READS IN A HtADfcR CARD ANU A KAUB SOUNDING
C AND THE MAX AND WIN VALUtS OF THE MIXING OEPTtH ARE CALCULATED.
C
DIMENSION THETA125)
COMMON /MIX/ ISTA,IYR,IMO.IDA,DAY,IDS,IHDt1ST,IDATt, lUATt i , I UA Tt9
COMMON /HMM/PUP4.P125) ,T125) ,MAXT,Ml NT,MMAX,HMIN,NSL,SAP
I FORMAT (/,1H ,*HM1N DENOM IS ZERO*)
5 FORMATlIA1U.1110.3J5)
6 FORMAT (22X, IF-,. 1.13X, IF5.1 J
C
C KEAO IN HEADER CAKU AND RAOB DATA.
READ 5,1STA,I DAT,MAXT,MINT.IDSC
H)S=IUS+1DSC* 1DATEIMDAT
IYR = IDAT*O.OOOl» JMU= IOAT-IYR*10000» IMO*JMU*0.01
IDA=JMO-lMO*lOO* UAY=30.5* »27 J. /
00 2i 1=1,100
READ 6,nUtC(L)» I I L ) =T < L ) «-27 J.2» IF { P ( L ) , tU. bOO.O J t.0 IU il
25 CONTINUE
27 NSL=L
C
C JETERM1NE JHt PRESSURE LEVtL AT WHICH THt POTENTIAL TEMPERATURE
C (THtlA) IS oRtATER THAN UK LUUAL TU THE SFC MAX POTLNTIAL TEMPERATURE.
60 IF (FLQATIMAXTI.LE.T(I)) 70,80
70 HMAX=0.0 $ GO TO 120
80 TtTMX=(UOOO.O/PI1))*«0.286)*MAXT
DO 105 I-ltNSL.
THETAtI ) = TI I)*l1000.0/P(I•)**0.206
IF (THETAII).LT.TETMX) 105,100
100 IS=I S GO TO 110
105 CONTINUE
1,0 Tl) 115
C
C CALCULATION OF MAX MIXING DEPTH.
110 lSl=li-I
PM=Pmi)MPl 1SJ-P1 ISin*(TbTMX-THETAlI<»m/lTMkTA( JSJ-THETAUSl) )
HMAX= L^t.7«lMAXT+Tl IS ) ) *ALUG( P ( IJ/PMI
IF (HMAX.GT.400U.O) 115,120
115 HMAX =
-------
SUBROUTINE SFCObSl
C
C THIS SUBROUTINE READS IN THE FIRST SURFACE UflS CARD WHICH
C CONSISTS OF THE 2<»00, 0100, 0200, AND 0300 HOURS DBS.
C
COMMON /SFC/ IHK(2<»), ICHT(2<») , IHDU*),W$12<»J ,
t> FORMAT {6X, 14,215, 5X, 15, F5.0)
READ 5,1HR(24) ,ICLD(24I ,'ITEMP(24) ,IWDl£4l ,WS(24)
L=2A$ IHR(L)=L
ITEMPIL) =
RETURN i END
50
-------
SUBROUTINE SFCOBS2
C
C THIS SUBROUTINE READS IN THE KEST OF IHt SURFACE OBSERVATIONS
C FOR THE DAY.
C
COMMON/CPTS/ NPT,XPT(625).rPT«625)
COMMON /SFC/ lHRl2*l,lCMTU*ltlWD(2<>).WJ>(2*)ilCLU(24),ITEMPl2'»)
COMMON /PNT/ ICiRT2tIJS(24)tHSl24).Him68»fCCAlU2pO»
5 FORMAT (6XtI4,2l5t5X,I5.F5.0J
READ 3, (lHR(L),ICLD(L»,ITEMP(L»,IWO(C),WSa»,L=l,2JJ
DO 150 L-li23
IWD(L)=10.0*1WD(L)$ HSIL)«0.5154*HS(L)
150 ITEMP
-------
SUBROUTINE tXTURB (EXTRAQ)
C THIS SUBROUTINE DETERMINES THE MAX WIND SPEED FOK THE JAY AND
C OBTAINS THE WIND DIRECTION VECTORl ALLY. IT THtN CALCULATES THE
C EXTRAURBAN CO CONTRIBUTION ON THE BASIS OF GASOLINE CONSUMPTION IN
C THE StCTOR.
C
C THE FUEL-ARRAY CONTAINS THE RATE OF GAS CONSUMPTION IN THE
C SEGMENTS AS A FUNCTION OF WINP DIRECTION.
C
COMMON/CFUEL/ FUEL ( 16 ) ,NHOU 10)
COMMON /HMM/PUP4,P(25) , T (25) , MAXT, Ml NT, HMAX, HMIN.NSL, SAP
COMMON /SFC/ lHRt24),ICHT(24),IHDI24),WSI2<»)tICLDt2<»),lTEMP124)
1 FORMAT 1/ilH ,*MAX WIND IS CALM.*)
C
C DETERMINE THb 24-HOUR MAX HIND SPCtO.
WSMAX=0.0
DO 100 1 = 1, 24
IF (WSl I ).GT.WSMAX) riSMAX=WS(I)
100 CONTINUE
C
C JETERMINE THE COMPONENT* OF THE AVfcRAtit WIND DIRECTION FOR
C THE MAX HIND SPEED.
NWD=XWS=YWS=0.0
DO 105 1=1.24
IF (WS( D.NE.HSMAX) GO TO 105
WDMAX=IWOt I 1*0.0174533 $ NWD=NWD*l
S
C
C DETERMINE THE AVERAGE WIND DIRECTION ACCORDING TO 16 D1KECT1UNAL
C POINTS IN RADIANS.
IF (NWU.EQ.U GU TO 135
IF IYWS.NE.0.0) GU TO 130
110 IF 1XWS) 120,115,125
115 PRINT 1 t STOP
120 WDMAX^.71238dSI t GO TO 135
125 WDMAX=1. 5707963 » GO TO 135
130 WUMAX=ATANl XWS/ YKS )
IF (YWS.LT.0.0) WDMAX=«DMAX»3. 1415927
IF (WUMAX.LT.0.0) HDMAX=rtUMAX*6. 28 31 853
135 I01R = 2.5't64791*WDMAX*0.5
IF ( I JIR.kO. J) 1DIR=16
C
C CALCULATE THE EXTRAURBAN CONCENTRATION.
EXTRAU=(FUEL( ID1 R ) *0. 02900 / I WSMAX*HMAX) »*2.4t-4
RETURN t END
52
-------
SUBROUTINE M1NWIN UTI
C THIS SUBROUTINE CHECKS THE HOURLY MIND SPEED AND SETS EACH MIND
C SPEED LESS THAN 1.0 M/S TO 1.0 ( FOR WIND SPEEDS EQUAL 10 0 THE
C MIND DIRECTIONS ARE SET EQUAL TO THAT OF THE PREVIOUS HOUR).
C
COMMON XSFC/ IHKU4l,ICHTU«J,lNO(2*l«MS(2*ltlCLOU«)t ITEMP(2<»)
MI*MS(IT)
IF (MI.CE.1.0) 00 TO 105
IF (MI.LE.0.0) IWOtlTJMDP* MSUTI-1.0
105 IDP-IWOUT)
RETURN * END
-------
SUBROUTINE STABLE I I .SAL , U )
C THIS SUBROUTINE DETERMINES A STABILITY INDEX THRUSH A SERIES
C OF CHITERIAS CONCERNING CLOUD COVERt WIND SPEED, AND SOLAK
C ELEVATION (INSTABILITY INOEXI.
C
DIMENSION IX(25J«HCOi(,25)
COMMON/CITY/ I TYP , I SM, IOUT ( 10 )
COMMON/CSTAB/ SP180.CP180
COMMON /MIX/ ISTA.IYR, IHO, I DA, DAY, IDS, IHOt I ST , IDATt , IOAIE1 • I OATE9
COMMON /SFC/ 1HR(24J , lCHT(2
-------
c
C CALCULATION OF NIGHTTIME STABILITY
305 IF IWSP.GT.6.0) GO TO 310
IF (CC.GE.0.5.AND.WSP.GT.3.01 GO TO 310
IJ*5
310 IJSUJ = IJ
RETURNS END
55
-------
SUBROUTINE DEPTH IIT.SAS.HT)
C THIS SUBROUTINE CALCULATES THE ATMOSPHERIC MIXING JEPTH FOR EACH
C HOUR. AT SUNSETi SUBROUTINES RAOBHHM AND SFCOBS1 ARE CALLED INITIATING
C CALCULATIONS OF NIGHT-TIMt MIXING DtPTHS.
C
DIMENSION IWKI7)
COMMON/CITY/ ITYP,1SM,IOUT(10)
COMMON /DAY/ IDW.NH,IDAHO,IHH,IDHT
COMMON /MIX/ ISrA,lYR,IMO,lDA,DAY,IDS,IHU,lST,lDATE,lDATEl,IDATE9
COMMON /HMM/POP<»,P(25),T(25),MAXT,MINT,HMAX,HMIN,NSL,SAP
COMMON /SFC/ IHRI24) ,ICHT(24),IWDI24),WS(24),ICLD(24t,1TEMPI24)
COMMON /PNT/ IC,RT2,IJS(2<»),HSt24),HR(168)fCCALU200)
DATA HSUN,HRSUN/1000.0,18.0/
DATA IIWKIL),L=1,7> /5H(MON),5HITUt) , 5HIWkU) ,5HITHU ) ,5H(FR1),
15H(SAT),SHI SUN)/fIHOL/5HIHOL)/
1 FORMAT I/.1H ,20X, *DATE = *, IX,16,2X,A5,*X,*NO. OF
1 RAOB LEVELS =*13,*X*SFC PKtSS =*F7.1,4X*SFC MAX TEMP =*I*,^X*SFC
2MIN TEMP =*14,/1
v FORMATUH ,63X,Ald,/)
ITT=ITEMP(IT)t IF (ITT.LT.MINTJ ITT=M1NT
C
C IS IT DAYTIME tSAS UKEATtR THAN 01.
IF (SAS.GT.0.0) GU TO 100
C
C IS IT SUNbtT
IF (SAP.GT.0.0) 105,120
C
C LINEAR TEMPERATURE INTbRPOLAT 1L)N IS USED TO CALCULATE THE
C DAYTIME MIXING UbPTH.
100 HT=HM1N»(ITT-MINT)*IHMAX-HMIN)/IMAXT-MINT) $ GU TO 132
105 HT = HMIN»(ITT-MINT)* IHMAX-HMIN)/(MAXT-M1NT) * HSUN=HT » HRSUN=IT
C
C AT SUNSET, PK1NT THE BASIC STATION DATA. READ THE NEXT SUUNUINU
C AND CALCULATE THE MAX AND MIN MIXING DEPTHS, AND HEAD THE NEXT
C DAYS FIRST FOUR HOURS OF SFC DATA.
108 IF (lOUT(l).LT.l) GO TO 115
IF (ITIM.NE.U PRINT 9.ISTA
IF (IDAHO.EO.O) GO TO 110
PRINT 1, 1DATE1,IHOL,NSL,PID.MAXT.MINT S GU TO 115
110 PRINT 1, 1DATE1,IHK
-------
SUBROUTINE CALXOQ (HT.JI)
C THIS SUBROUTINE CALCULATES THE X/U VALUES FOR tACH OF THt
C NINE SECTOR SEGMENTS USING THE GAUSSIAN AND SOX MODELS.
C
DIMENSION RUO),A19,6),B19,6),AA(54|,BB<5*)
COMMON /QUE/ QL 110) fQM( 10) ,Q12( 101 ,0341 10) ,06(10) , XOOl 10)
COMMON /ALLC/ YS(351),IYN(351)
ItXRUOJ.YR(lO) ,ARtlO),EFACI9),ACR.LR4,LR5,LR9,LRlO.AG,YCl,VC2,XRC,
2YRC,XC2,Br2,BYl,Bl2,8Il,CRI,IT,IJ,XGirG,SA,CA,NN,CR,LRL
EQUIVALENCE
-------
SUBROUTINE CAIQUE UWD.M)
C THIS SUBROUTINE COMPUTES THE CO EMISSION WITHIN EACH SECTOR SEGMENT.
C
DIMENSION Y(20),X(20),IY(20),JY<20)
COMMON/CITY/ ITYP , ISM,IOUT(10»
COMMUN/CPTS/ NPT,XPT(625).YPT(625)
COMMON/CSTM/ VCAR(10),AST(10),WHST110),ICST(10),NOCO(10)
COMMON /QUE/ UL( 10),UMUO) ,Q12(10) ,034(10) ,46(10) .XOQ(IO)
COMMON/CLDAT/IFK(1200),X1(1200).X2(1200),Y1(1200),Y211200),E(1200)
COMMON /ALLC/ YS(351),IYNt351)
1,XR(10),YR(10),AR(10),EFAC(9),ACR,LR4,LR5,(.R9.LRIO.AI,,YC1,YC2,XRC,
2YRC,XC2,BY2,BYl,BI2,flll,CRl,IT,U,X(;,YG,SA,CA,NN,C«,LRL
XG=XPT(M)$ YG=YPT(M)
MCU=NOCD(M)$ IF { ISM.EQ.O) MCO=0
ANG=ACR*IWD $ CA=CUS(ANG) S SA=SIN{ANG)
Y12=0.0625*YRC
DO 210 1=1,9
Q12m = Q34(I)=>O.U
210 CONTINUE
C
C BEGINNING OF LOOP FOR COMPUTING EM1SS1UN WITHIN tACH SECTOR SEGMENT.
DO 300 NC=1,NN
IF (MOD.EQ.NCI GO TO 300
C
C TRANFURMATION OF COORDINATES.
221 XD1 = XHNC)-XG »YD1=Y1INC)-YG »XD2=X2(NC)-XG »YD2=Y2(NC1-YG »1YT=0
YSI=XD1*SA*YIU*CA$ YS2=XD2*SA+Y02»CA$ IF (YS1.1T.YS2) GU TO 216
YSS=YSlt YS1=YS2* YS2=YSS * IYf=l
216 IF IYS2.LT.O.UR.YS1.GT.YC1) GO TO 300
JY(1)=JY(2)=0
XSI=XD1*CA-YD1*SA» XS2=XD2*CA-YD2*SA
IF IIYT.EQ.O) GU TO 222
XSS=XS1$ XS1=XS2$ XS2=XSS
C
C CHECK TO SEE IF LINK LIES WtTHIN SECTOR.
222 Bl=B2=BY2t IF (YSl.GT.YC2) GO TO 22*
B1=BY1J IF (YS2.LT.YC2) 62=BYl
22<» IF (ABS(XSl).LT.Bl»YSl) GU TO 230
IF (ABS(XS2).LT.B2*YS2) GO TO 232
IF (XS1*XS2.LT.O) GO TO 2J*
IF (YS1.GT.YC2.UR.YS2.LT.YC2) GU TO 300
IF (ABSIXS1).LT.XC2.0R.ABS«XS2).LT.XC2) 234,300
2JO If (ABS(XS2).LT.B2*YS2) JY(2)-1* JY(1)=1$ GO TO 234
2J2 JY(2)=1
234 X(1)=XS1$ X(2)=XS2» Y(l)=YSli Y«2)=YS2
YD=YS1-YS2*0.0001 S XD=XS1-XS2
BL=XD/YDt AL=XS1-BL*YS1
C
C LOCATES THE SEGMENT WITHIN WHICH THt END POINTS OF THE LINK LIE.
IF (YSl.GT.O) GU TO 236
Y(l)=YSl=O.Ot X(1)=XS1=AL
IF (XS1.EQ.O) JY1U = 1
236 IYl=YSl*CRI*l * IY2=YS2*CRI*1
IF (IY2.GT.3iO) IY2=350
IY11)=IYN( IY1IS IF (YSl.GE.YS(IYl)) 1Y(1) = IY(1)* 1
IY(2)=IYN(IY2)» IF (YS2.GT.YS(IY2)) IY(2) = IY(2)* 1
58
-------
c
C LOCATES THE INTERSECTIONS OF LINK WITH SECTOR OIVIUING LINES.
N*2t IF ( lYm.EQ.W2)) GO TO 250$ L=IY(1)
238 YP=YR(LIS XP=AL+BL*VP* If ( ABS ( XPJ .GE.XRI LI ) GO TO .240
X(N)=XP» Y(N1=YP$ JY(N)«I Y(N)=L
IF (L.E=034(NS1)*QT
290 JY(N2)=0.0 * GO TO 294
292 CONTINUE
294 JY(N1)=0
295 CONTINUE
C
C END OF MAIN LOOP
300 CONTINUE
C
C COMPUTES CO EMISSION WITHIN SEGMENTS FROM MINOR LINK DATA.
CALL CALLOCIMJ
KETURN S END
59
-------
SUBROUTINE CALLOC IN)
C
C THIS SUBROUTINE USES THE GRID POINT VALUES (VCL AND VCN) TO
C COMPUTE THE AVERAGE CO EMISSIONS WITHIN EACH SECTOR SEGMENT. THE
C RESULTS ARE STORED IN QL AND QM.
C
COMMON/CPTS/ NPT,XPT(625»,YPTI625I
COMMON /QUE/ QL UO) . QM( 10) ,Q121 10) . Q3M10) ,Q6{ 10) , XOUl 10)
COMMON /CLOQ/ VCL(2500),VCM(2500),CB(31,SBi3)tNRH,NCL,MRW,MCLtCLfc,
IGSPfRGSfZQX
CQMHUN/ALLC/ YS «J51) , I YN(3!>1)
l,XR<10),YR(10),AR(10),£FAC(9),ACRiLR*,LR5tt»<9tLRlOtAGiYClfYC2,XftCt
2YRC,XC2,BY2iBYl,BI2i8ll,CRl,IT.U,XGtYG,SA,CA,NNtCRiLRL
XG=XPT1N)$ YG=YPTJN)
00 30 L*lt9
YRL-YR(L) t SEVL=SEVM*XDS=0.0 $ KTN*LO=1
IF U.CT.5) GO TO 5 ( KTN=3 $ LD=L-3
5 WCT=0.25/LD
00 20 JzltLD
FXDS=1.0-WCT*(X.O*2.0*«LO-J) ) t YRR=YRL*FXOS > SSfcVM«SSEVL=«0.0
00 10 K=1,KTN
mvl-YRR/CB(KI
YM=ID1VI*(CA*CB(K)-SA*S8(K))+YG)*ZQX+MKW
XM*(01VI*(SA*CBIK)*CA*SB1K.))»XG)*ZQX«-MCL
IYM=INTIYMJ* IXM=1NT(XM)
IF (IYM.LT.I.OK.IYM.GE.NRW) GO TO 10
IF (IXM.LT.l.OR.IXM.GE.NCU GO TO 10
DYM=YM-IYM* L)XM=XM-IXMi IXY= ( IXH-1 )*NRW»IYM
DYM1=1.0-OYM$ OXM1=1.0-DXM$ IXYl«IXY*l
TTEVL= (OYMl*VCHIXY)+DYM*VCL(IXY1))*OXH1
I *(DYM1*VCL(IXY* NRW»*DYM*VCL«lXYl*NRtO)*DXM
SSEVL=SSEVL+TTEVL
IF (L.LE.LRL) GO TO 10
TTEVM= (DYM1*VCM{IXY)*UYM*VCMJIXY1))*DXMI
1 +10YM1*VCM
-------
SUBROUTINE LOCXOU(HT.J)
C THIS SUBKOUTINt LOCATtS THE PROPER X/Q VALUES FROM STUKAGE.
C
COMMON/CSAV/ J5 tL7,19,K36, HMl , SXOW 13201 ,SQ12< 32*0 J ,
1 SULU240)
COMMON /QU£/ QLl 10) ,QMUO) ,012(10) .Q3M10I ,06110) ,XOU( 10J
HB zHMI+HMI+HMl
DO 10 t*ltL7» IF tHT.LT.HB) CO TO l
-------
SUBROUTINE LOCQUE(KD.N)
C
C THIS SUBROUTINE LOCATES THE PROPER Q VALUES FROM STORAOE.
C
COMMON/CSAV/ Jt>,L7,l9,K36,HMI,SXOq<320) , SO 12(3240) . SU3<»(
I !>UL(.i240l
COMMON /JUE/ ULllO),QM(lOI.Jl^ilOJ,034(10),U6(10l,XOgiU)
K=0.l*KO*0.5$ tF (K.LT.l) K = K36$ KBS= ( K-l) *l 9* «N-l) *K 3t>* 19
DU 100 1=1,19$ IBS=I*KBS
012(1 t = SQ12(IBS)S U3'»III=SQJ
-------
SUBROUTINE CALCON ( I ,N,HT, EXTRAU)
C THIS SUBROUTINE COMPUTES THE CO CONCENTRATION AT THE
C RECEPTOR POINT BASED ON THE VARIOUS MOUfcLS. THE OiURNAL
C PATTERNS FOR WEEKDAYS, SATURDAYS, SUNDAYS, AND HOLIDAYS ARE INCLUDED
C IN THE CO CALCULATIONS.
C
COMMON/CITY/ I T YP . I SM, IOUT » 10)
COMMON /CPTSX NPT,XPT(625) , YPTI 625 J
COMMON /DAY/ IDW,NH( IDAHO, IHM, IDMT
COMMON /MIX/ ISTA,lYR,IMO,IOA,DAV,IDS,lHO,lST,IOATt,10ATEl,lDATE9
COMMON /PNT/ IC,RT2,IJS«2*),HSI24)IHR(168),CCAL(1200)
COMMON /SFC/ lHR(24),lCHTC2*J.l«DI2<»l,WSI2*JilCLDI2l
COMMON /CAR/ PT12 I2-»J ,PT34 (2*) ,PT6(24) ,PTSAT (24) ,PTSUN(2«1
COMMON /QUE/ UL i 10 ) ,UMJ IOJ ,0121 10) ,Q34< 10) ,061 1UI ,XOUC 101
COMMON /ALLC/ YSO51 J , I YNI 351 1
l,XR(10),YR<10).AR(l(»,EFAC(9ltACR,LR4,LK!>,LR9,LR10(AG,YCl,YC2.XKC,
2YRC,XC2,BY2,BY1,BI2,BI1,CRI,IT,IJ,XG,YG,SA,CA,NN,CR,LRL
COMMON/CNZ*/
IF IN.EQ.l) IC»IC*1
10 CPL=0.0* WSI75=0. 75/1223. 6932*MS(ITI)
IDMT=IUW$ IF ( IDAHO. EQ.O) GO TO 15S IDWT=7
15 DO 120 J'1,9
ITT =IT-IFIXJYRIJ)*HSI75I+IDS
IF (1TT.LT.1J ITT=ITT*24t IF UTT.GT.2*) ITT-1
IF (1DWT-6) 100,105,110
C
C WEEKDAY CONCENTRATION CALCULATIONS.
100 FACT = Q12(J)*PU2(1TTI+(034(J) *UL( Jl I*PT3<»C ITT)
GO TO 120
C
C SATURDAY CONCENTRATION CALCULATIONS.
105 FACT > PTSAT1 ITT)*(Q12(J)»U34(J) *gL(J)J» GU TU 120
C
C SUNDAY AND HOLIDAY CONCENTRATION CALCULATIONS.
110 FACT = PTSUNlITT)»igi2(J>»03
-------
SUBROUTINE STREtTH.N)
THIS bURRUUTlNt COMPUUS THfc STREtT PROFILE OF CO CONCENTRATIONS bASEU
ON THt BACKGROUND VALUES, THE STREET EMISSION VALUES, ANU THE
METEOROLOGICAL CONDITIONS (RESULTS ARE STORED IN ARRAY CCAL).
COMMON
COMMON
COMMUN
COMMON
COMMON
COMMON/CSTM/ VCARIIO) ,Asr ) =UHAS*Wl *H/2»C BG
CCALC IU*12)=UBAS«XLUCBG* CCALt 1D+13) =
CC AL ( I D«-I4 )=OBAS*WI »IUl»CBG* CCALt ID* I
ITIME=1
RFTURNt END
64
-------
SUBROUTINE PPDATA
C THIS SUBROUTINE PRINTS AND PUNCHES OUT THE STATION IU, DATE, MAX-MIN
C TEMPERATURES, ALL SFC DATA. STABILITY INDtX, MIXING DEPTH, AND THE
C CALCULATED CO CONCENTRATIONS.
C
DIMENSION IDVU4)
COMMON/CITY/ ITYP.ISM,IOUT(101
COMMON/OPTS/ NPT,XPT(625),YPTI6251
COMMON /DAY/ IDM,NH,IDAHO,IHW,IDUT
COMMON /MIX/ ISTA,IYR,1MO,IDA,DAY,IDS,IHD,1ST,IOATE,1DATE1,IUATE9
COMMON /PNT/ IC,RT2,lJS(2*ltHS(24l,HRU6B),CCALU2QO»
COMMON /SFC/ 1HR124),ICHT!24),IHO(24),HS(24),ICLD(24J.ITEMPJ2*)
COMMON/CN2W N24
DATA UOV(L),L = 1,141/10H HOUR -,10H CLD-H =,10H CLD-C *.
110H TEMP = ,10H WNU-D *,10H WND-S «,IOH STB-1 «,IOH MIX-
20 =,10H CO-Bb =,10H CU-R1 =|10H CO-R2 -,IOH CO-LI -,
310H CO-L2 =,6H STA/
1 FORMAT (IH ,5X,A10,24I5)
2 FORMAT (IH ,5X,A10,24F5.0I
3 FORMAT (IH ,bX,A10,24F5.1I
4 FORMAT (IH ,bX,A10,24F5.2)
5 FORMAT I/.IH ,*CO COMPUTATION TIME - *F6.3,» StC.*/)
6 FORMAT IA10, 110,4U.F4.0,14,F5.0.5F6.2J
8 FORMAT (IH ,5X,A6,I4)
10 FORMAT (» CITY DATE HR PT CLU OU SPO SI MU CBG
ICR1 CR2 CLl CL2*>
11 FORMAT!* CITY DATE HR CLD DIR SPD SI MD*)
12 FORMAT(8(I4,F&.2)I
C
C LIST OUT OUTPUT DATA IF lUUT(l) GT ZERO.
IF ( IOUT(U .LT.l) GO TO 100
PRINT 1.1DV11),
PRINT 1.IOV13),
PRINT 1,1DV(4), (ITEMP(L),L»1,24J
PRINT 1,IOV<5),
PRINT 3,IDV(6I.
PRINT I.IDV17J,IIJS(LI,L=1,24J
PRINT 2,IDV(8), (HS(LJ,L=l,2'tJ
NCE'5$ IF llSM.tU.O) NCE=1
DO 15 N=1,NPT* ICB*IN-U*N2'»
PRINT 8.IDVI14),N
IC1=ICB*1$ IC2=tCB»N24
DO 15 NCB»l,NCEi N9=8*NCB
PRINT 4,IDV(N9»,(CCAL(LI,L«IC1.IC2I
ICI=IC1*2*0$ IC2=IC2*240
15 CONTINUE
C
C PUNCH OUT OUTPUT UATA IF IOUT(2I GT ZERO.
100 IF (1QUT12J.LT.U GO TO 200
IF (ITYP.GT.2) GO TO 150
IF (ITIM.NE.l) PUNCH 10
DO 105 1=1,24* DO 105 N=l,NPTi ICB»(N-ll*N24»l
105 PUNCH 6,ISTA,1UATE,I,N,ICLO(I),1MD(I),WS(II.IJS(I).HS(II,
1 CCALI1CB),CCALIICB+240J,CCAL(1CB*480»,CCALIICB»720»,CCAL(ICB»V60I
GO TO 200
65
-------
150 I*IHD
IF UTIM.Nfc.l) PUNCH 11
PUNCH 7,lSTA.IUATE,l,lCLD(U,IWOm,WSm,IJSm,HS(l)
PUNCH 12, (N,CCAUN),N=1,NPT)
200 IF (lOUTO).LT.U GO TO 300
IF (ITYP.GT.2) GO TO 250
DO 205 1=1,24$ DO 205 N=1,NPT* ICB= iN-ll*N2**I
205 WHITE
1 ISTA,IDATE,I,N,ICLOm,IWDm,WSll),lJMU,HSm,
2 CCAUICB),CCAUICB«-240),CCAL(ICB+<»801,CCAUICB+72tfJ,CCALtICB*960)
GU TO 300
250 I*IHQ
WRITE (41)
i ISTA.IDATE, i,icLum,iwDm,wi>m,iJsm,Hsm
2 .(N.CCAKNlrNxl.NPT)
300 IHW=IHW+2^
I T I M= 1
IDATE=IDATE1
RETURN i END
66
-------
Appendix B
EPA IBM 360/50 VERSION OF APRAC-1A COMPUTER PROGRAM
C AIR PCLLUTICN DIFFUSION MODEL
C THIS PROGRAM COMPUTES THE CO CONCENTRATIONS AT VARIOUS CITY
C LOCATIONS BASFD ON TRAFFIC AND MFTEOPOLOGICAL INPUT DATA.
C
C SFT ITYP EQUAL TH
C 1 FOR SYNOPTIC MODEL
C 2 FTR CLIMATOIOGICAL MODEL
C 3 FOR GRID PT MODEL.
C
HEAL MRHiMCL
COMMON/CITY/ ITYP , \ SM,IOUTI10)
CCMMCN/CPTS/ NPT,XPT(625),YPT(625)
COMMON/CSTM/ VCAR(IO) ,AST(10),WWST(10),ICST(101,NOCD(10)
COMMON/CHANG/ HANGtACV
COMMON /CCP-D/ XXCtYYC
COMMON /CAY/ IDH,NH,IDAHO,IHH,IDHT
CCMMCN/CSTAB/ SP180.CP180
COMMON/CFUEL/ FUEH 16 I .NHOH10I
COMMON/CPFC/ PF1,PF2,P01,P02,S(,IWD(24),WSI2*), ICLD(24),!TEMP(24)
COMMON /CAR/ PT 12(24), PT3M 24),PT61 24), PTSAT(24),PTSUN(24)
CCMHON /CUE/ CLUOI.QMIIO),012(10),034(10)tQ6UO),XOQ(lO)
COMMON /PNT/ IC,RT2,US(24),HS(24),HRI168),CCALI1200)
COMMON/CSAV/ J5.L7,19,K36,HMI,SXOOI320),S012(9240),3034(3240),
1 SOLO240)
COMMON /CLOQ/ VCH2500) ,VCM(2500) ,CB(3),SB(3)tNRW.NCLtMRW.MCL,CLE,
IGSP.RGS.ZQX
COMMON/CLDAT/IFKI 1200),XI(1200),X2(1200),Yl(1200),Y2(1200),E«1200)
COMMON /ALLC/ YS(35 I),1YN(351)
l,XP(10),YP(10),AP(10),EFAC(9),ACR,LR4,LR5,LR9,LR10,AG,YCl,YC2,XRC,
2YRC,XC2,B>2,BY1,RI2,BI1,CRI,1T,IJ,XG,YG,SA,CA.NN,CR,LRL
I FORMAT (1H1,55X'SYNOPTIC AIR POLLUTION MODEL'/)
2 FCRMAT(1H1,52X,'CLIMATOLOGICAL AIR POLLUTION MODEL'/)
3 FORMAT (1H1.53X,'GRID POINT AIR PCLLUTION MODEL'/)
C
C INPUT THE BASIC PARAMETERS.
IHW=0
CALL INCAT
IF( IABS(ITYP).EQ.l)PfUNT 1
IFl IABS( ITYP).EQ.2)PRINT 2
IFIIABS(ITYP).EC.3)PRINT 3
r
C COMPUTE THE EMISSION FACTOR AND OTHER INITIAL QUANTITIES.
CALL BASIC
C
C INPUT THE MAJOR AND MINOR TRAFFIC EMISSION DATA.
CALL LINKS
C
C PREC.OMPUTE AND STORE ALL THE X/Q AND Q VALUES (IF ITYPE = 2).
IF(IABS(ITYP).EC.2)CALL STORE
C 1
C INPUT THE HEADER CARD AND RAOB OBSERVATIONS AND
C DETERMINE THE MIN-MAX MIXING DEPTH FOR THE FIRST DAY.
CALL RAOBHM
67
-------
c
f INPUT THE FIRST SFC DBS CARD.
CALL SFCPR1
C
C INPUT THE REMAINING SFC DBS FOR THF DAY.
100 CALL SFCDB2
C
C DETERMINE THF FXTRAURBAN CO CONTRIBUTION.
CALL EXTURB (EXTRAO)
C
C BEGIN HOUR LCC1P
DO 200 1=1, 24
C.
C IN GRID POINT MOOFL THE COMPUTATION IS DONE FOR HOUR IHD ONLY.
miABSC ITYP).EQ.3)I=IHD
C
C DETERMINE THE MINIMUM WIND SPFED AND DIRECTION.
CALL MINWIN (I)
IWI=IWD(I)
C
C DETERMINE THE STABILITY INDEX.
CALL STABLE(I,SAS,IJ)
r.
C DETERMINE THE MIXING DEPTH AND AT SUNSFT INPUT THF HEADER
C CARDt THE RACB DATA, AND THE FIRST FOUR HOURS OF SFC DATA
C FOR THE NEXT DAY.
CALL CFPTh (I.SAS.HT)
f
f DETERMINE THE X/0 VALUES USING THE GAUSSIAN AND BOX MODELS.
IF(IABS( ITYP».F0.2)GO TH 140
GC TO 150
C
C SELECTS APPROPRIATE X/0 VALUES FROM STORAGE.
140 CALL LOCXOO(HT,IJ)
GO TP 160
C
C CCMPUTE APPRCPRIATF X/0 VALUES.
150 CALL CALXOOIHT,IJ)
C
C BEGIN POINT LOOP
160 DO 190 N=1,NPT
C
C DETERMINE THE Q VALUES FROM TRAFFIC EMISSION DATA.
IF(IABS( ITYPI.E0.2IGD TO 170
GO TP 180
C
C SELECTS APPROPRIATE Q VALUES FROM STORAGE.
170 CALL LOCQUF(1WI,N)
GO TP 189
f
C m^PLTF APPPCPRIATF 0 VALUES.
180 IF( IABS( ITYP ) .EQ.3)IWI=IWD(IHD)
CALL CALOUE(IWI,N)
C
C DETERMINE THF TOTAL BACKGROUND CO CONCENTRATION.
189 CALL CALCCN (I,N,HT,EXT RAO)
C
68
-------
C DETERMINE THE STREET PROFILE CO CONCENTRATIONS IF IS* IS
C NOT ECUAL TC ZERO.
IFHSM.NE.OJCALL STRFETUtN)
C
190 CONTINUE
IHIA6SUTYP).E0.3)GO TO 210
?OG CONTINUE
210 CONTINUE
C
C PRINT AMP PUNCH THE DESIRtD PARAMETERS.
CALL PPDATA
C
C DAY CF THE HEFK CHECK.
IF (NH.EQ.OI GO TO 300
IF (IDAHO.EC.01 GC TO 295
ICAHO=0
NH=NH-l
295 IF UDATEl.EQ.NHOLCNHM IDAHO«l
300 ICW»IDW*1
IF (IDW.LT.8! GO TO 400
IDW=l
I HUM
IC=0
C
C CHECK FOR THE END DF THE TIME PERIOD.
400 IF (ICATE1.LE.ICATE9) GO TO 100
500 STOP
END
69
-------
SUBROUTINE INCAT
C
C
C THI S
C
SUBPPUTINF RFADS IN BASIC DATA
REAL MRW.MCL
CCfMCN/ClTY/
COMMON/CPTS/
CCMMCN/CSTM/
CCMMfN /CCRC/
ITYP.ISM,I OUT!10)
NPT,XPT(625) ,YPT(625)
VCAR(10),AST(10),WW!>T< 10),ICST( 10),NflCD(10)
XXC,YYC
COMMON/CSTAB/ SP1BO,CP130
CCMMPN /CAY/ IOW, NH, IDAHO, I HW , IDWT
COMMON/CPFC/ PF1 ,PF2,0ni ,P02tSt9) ,KT(24)
COMMON/CFUEl/ FUEU 16),NHOL ( 101
C CM MO N/ Mix/ 1ST At 10) , IYP, I MO, I DA, DAY, IDS, IHD, 1ST, I DATE, I DATE 1,
X IDATE9
CCMMCN /GRAF/ I I(7),IA( 2^), I8« 16,8)iIVI 5,8) ,IGS
CCMMfN /H^M/PCP^,P(25),T(25 ) , MAXT ,M INT , HM AX , HMIN, NSl , SAP
CCMMPN /CAR/ PT12(2
-------
XPT(I)=XPT(II+XXC
YPT( 1)=YPT!I)»YYC
POPS=PnPS*1.0F6
POP4=SORT(SORT(POPS ) I
!GS=IHW
I OIHW- 1
IUAHC=0
SAP=-1.0
ACR=0.017^533
PL18C=ACR*SLAT
SP190=SIN(PLieO)
CP180 = COMPL180 )
PF1=PF1/360C.O
KETURN
ENO
71
-------
SUBROUTINE BASIC
0
C THIS SUBROUTINE COMPUTES BASIC QUANTITIES USED THROUGHPUT PROGRAM
C
COMMON/CHANG/ HANGIACV
CCMMC^/CPFC/ PF1 ,PF2,P01,PO2,5(9),KT(24)
COMMON /CAP/ PT12I24) ,P T34( 24 I , PT6( 24 ) , PTSATJ 24 ) , PTSUM 24 )
COMMON /ALLC/ YS(351 ) ,1YN(151J
1,XR< 10) ,YP(10) , AR (10) ,CFAC(9) ,ACR,LR4,LR5,LR9,LRIO, AG,YCl,Yf.2,XRC,
2YRC,XC2,BY2,BY1,BI2,BI1,C«I,lT,IJ,XG,YGtSA,CA,NN,CR,LRL
ACR = 0.0174533
CR=62.137
HANG=0. 196349625
LRL=5
C
C COMPUTES BASIC FACTORS FTP EMISSION MODEL
C KT(U=PEAK CR NON PEAK TRAFFIC AT ITH HOUR
C PT(I)=FRACTIGN PF OAILY TRAFFIC DURING ITH HOUR
( EFAC( I) = FMISSION FACTOR FOR ITH FACILITY TYPE
( SU) = SPEED FOR ITH FACILITY
C POl AND PF1=INPUT CONSTANTS.
TNU25=TAN(hANG)
RR=7.7671
CRI=10.0/CR
ACV=CRI*CRI*1.0E4
LR4=4
LP5=LR4+1
LR9=9
LR10=LR<>+1
XP=0.85**P01
DC 50 I =1 ,24
IF (KT(IJ.EO. It GO TO 45
GO TO 50
45 PT12(I)=PT12(II*XP
PT34(I)=PT34(I)*XP
PT6( I I = PT6( I)*XP
50 CONTINUE
DO 60 1=1,9
60 EFAC( I)= PF1*SU)**P01
f
C COMPUTES GFOMETRIC QUANTITIES TOR SECTOR SEGMENTS
C AR(I)=APEA CF ITH SEGMENT
C YR(I)=MAX DISTANCE OF ITH SEGMENT FROM RECEPTOR
C XtUI) = MAX HALF WITTH OF SEGMENT
C LR4 AND L"9 = NC. OF SEGMENTS IN SMALLER AND T^TAL SECTO" .
DO 210 L= l,LR10
AG=11.25
IF (L.LT.LR5) AG=?2.5
AG=ACR*AG
TANAG=TAN(AG)
YR(L)=SORT(RR*RR*AG/TANAG)
XR(L )=YR(L)*TANAG
IFfL.EO.LR1)YC1=YP(L)
IF (L.NE.LP4) GO TO 210
YC2=YR(L)
XRC= YC2*TNH25
XC2=XR(L)
YPC=YC2
72
-------
PR=PR+RR
AR(1)=XR(1)*YR(1I*ACV
SUMA=AR(1)
DC 211 1=2,4
ARU ) = XR(I)*YR< 1)*ACV
211 SUMA=SUMA+ARt I I
SUMA=XPC*YR(
-------
SUBROUTINE LINKS
f.
C THIS SUBROUTINE READS IN DATA FOR MAIN ANO LOCAL LINKS.
C
REAL PRW.MCL
DIMENSION V(8)
COMMON/CITY/ITYP, ISM, IOUTC 10)
CCMMCN/CPTS/ KPT,XPT(625),YPT(625)
COMMON /CSTM/ VCAR(10lfAST(10)tWWST(lO),ICST(lO) .NOCOIIO )
COMMON /CORD/ XXC.YYC
CCMNCN/CHANG/ HANGtACV
COMMON /CL DA T/IFM 1300 ), XI ( 12001 ,X?« 1200) ,Y1 U200) ,Y2 ( 1200 ), Et 1200 )
CCMMON/ALLC/ Y1 (351 ) , 1Y Nt 351 )
ItXRI 10) ,YP( 10) ,AR(10» ,FFAC(9» , ACR ,LR4 ,LR5 , LR9 ,LRIO, AG,Y Cl t YC2f XRC ,
2YRC,XC2,BY2,BYl,BI2i ^IltCRl t I TfIJ,XGtYG,SA,CA,NN,CP ,LRL
CCMMCN /CLP.C/ VCL (?500) ,VCM(2500),CB«3), SB( 3),NPW,NCLtMBW,MCI iCLE,
1GSP.RGS.ZOX
1 FORHATI Il,9X,AF5.Cf F6. 0,215)
2 FCRMATU1,PX,3UO)
FORMAT ( IH ,59X, •*<"). OF LINKS =',I5,/)
4 FORMATJ8F10.6 )
'; FORMAT(2X,'GRin INPUT, ',K,« X',m
IF( ITYP.GT.O)GO TO 60
10 LT=KRW*NfL
N = 0
15 READ A, V
DO 20 1=1 ,8
N=N+1
IF( N.GT.LTIGO TO 60
70 VCM(N)=V(I)
GO TO 15
60 M°W = NRW*0.5
MCL=NCL*O.S
LT=NRW*NCL
ARt=l.O/(GSP*GSP**CVl
Rf,S=100.0/GSP
DGS=.05*RGS
ZOX=O.Ol»RGS
CRAD2 = YC1*YC1/0.<56+1000
SVL=0.0
NC5=0
IF( IA8S(ITYP).EQ.3)GO TO 65
NCD1=NOCD(1)
65 XG=XPT(1)
CB( l) = l.O
CP(?)=COS (A75 )
CB(3)=C05(A75)
SB( 1 1 = 0.0
SB(2)=SIM-A75)
SB(3)=-SR(2)
I F( ITYP.LT.OGO Tf 80
00 75 1=1, LT
VCL1 I 1=0.
VCM(I)=0.
74
-------
C "EADS IN DATA FOR MAIK LINKS
C (LINK DATA STORED IN ARRAYS XI, Yl, X2, Y2 , AND E).
IF (LRL.GE.LR9) GC TO 100
GO TO 120
100 READ 1,H1,X1L,Y1L,X2L,Y2L,VV,IF,LL
IF (M1.LT.9) GO Tr 110
GO TO 170
110 VL=VV*LL*0.01
SVL=SVL*VL
EE*EFACUFI*VL
IFK(L)=IF
XP-X1L-X2L
YD=Y1L-Y2L
X1(U) = XU+XXC
YULI-YU+rYC
X2(L)=X2L*XXC
Y2(L) =Y2L»YYC
E(L)-EE*EE/(XD*XD+YD*YO»
L=L»1
r,c TO 100
120 READ l,Ml,XlLtYlL,X2L,Y2LtVV,IF,Ll
IF (Ml.GT.8) GO TO 170
NCD=NCO*l
130 VL*VV*LL*O.C1
SVL=SVL+VL
EE=EFAC(IF)*VL
XD=X2L-X1L
YO«Y2L-Yll
DD=XD*XD*YD*YD
X1L=X1L+XXC
Y1L=Y1L*YYC
C
C THE MAIN LINK DATA ALSO USED TO COMPUTE fNISSION VALUES AT GRID POINTS
C (VALUES STORED IN VCM, NO. OF GRID POINTS = NRWXNCL).
IL2=USORT«DDJ*DGS
ZIL=EE*ZII
YH=YD*2QX*ZII
X»»=XD»ZOX*ZII
YJ=Y1L*ZOX-YM*0.5*MRW
XJ=X1L*ZOX-XM*0.5+MCL
DC 135 1=1, ILZ
YJ=YJ*YM
XJ=XJ+XM
JY=INT1YJ)
JX=INT( XJ )
IF (JY.LT.l .OR. JY.GE.NRWI GOTO 135
IF I JX.LT.l.HR. JX.GE.NCLI GO TO 135
OY=YJ-JY
DX=XJ-JX
IXY={ JX-1>*NRW+JY
DY1=1.0-DY
0X1=1. 0-DX
IXY1=IXY+1
!N=IXY*NRW
IN1=IN+1
VCHHXY » = ZIL*DYI*DX1*VCM(IXYI
75
-------
VCMl IXYl) = ZIL*nY*OXl + VC«M IXY1)
VCM( IN)=ZIL*OY1*DX+VCM{ IN) , .
VCMU N1) = ZI l*DY*r>X+VC»M INI)
CONTINUE
lL=(LL/25)+2
I LM 1=11-1
XCD=XD/ILM1
YDD=Yn/ILMl
IF (NPT.FQ.U GO TO 142
r,0 TO 141
IF (NCDl.fO.NCO) NOCDU) = L
141 E(L )=EE*FE/OD
X1(1)=X1L
YUL)=Y1L
IFKIL )= IF
X2(L)=X2L+XXC
Y21L)=Y2L+YYC
L = L*1
GO TO 120
142 XTL=X1L-XG
YTL=Y1L-YG
OP 160 1=1, IL
IF ( I XTL*XTL+YTL*YTL).LT.CRA02) GO TO 140
150 XTL=XTL*XDn
160 YTL = YTI_ + Ynn
GC TT 120
17C NN=L-1
IZS=0.0
CCN=0.0
GO TO 175
30 PH1NT "i.NRW.NCL
DC 90 1=1, LT
VCL( I )=VCL( I 1*0. C
<30 VCHI I ) = vrM(I )*ARI , - •
RETURN
175 PPTNT 3^'N
C
C RFADS IN LOCAL EMISSITN CATA AMD TRANSFORMES IT ONTO A GRID
C (G^ID VALUES STOBCO IN ARRAY VCD.
250 REAC 2, HI, IX, IY, IZ
IF (M1.GT.8) GO TC 2/5
IZS= IZS+IZ
CX=I I X+XXC)*0.01
C Y=( I Y*YYC) *0.01
Y J=CY*PGS+MRW
JY=INT(YJ )
JX=INT(XJ )
IF ( JY.LT.l.OR. JY.GF.NRW) GO TO 270
IF ( JX.LT.l.OR.JX.Ge.NCL) GO Tn 270
IXY=( JX-1)*NRW+JY
IXYN= IXY + NPW
OYDX=DY*DX
I XYN1=IXYN-H
VCL ( IXY )=I7 *VCLI IXY)
VCL(IXY1)=IZ tVCLUXYl)
76
-------
VCL UXYN)=I Z + VCUUXYN)
VCl ( IXYN1 )=IZ + VCLUXYN1)
271 CGNTIMUE
GO in ?50
"5 !ZS=E«.'>4
IF (! 25* EC, OS GC r" 260
CCN=0.01*CLf*SVL*ARI*EFAC!5)/IZS
280 DO 290 L=1,LT
Z90 VCL
-------
SUBKDUT INF RAHBHM
C
f THIS SUBROUTINE READS IN A HEADER CARD AND A RAOB SOUNDING
C ANPl THE MAX AND M IN VALUES OF THF MIXING DEPTH AKF. CALCULATFO.
C
DIMENSION THFTAI25)
CCM«ON/MIX/ ISTAI 10), IYR , IMO, IDA,r>AY,lDStlHD,lSTtlDATEfI DATE1,
X IDATE9
COMMON /HMM/PnP4,PI 25) , T(25) , MAXT.Ml NTtHMAX,HMI N,NSL , SAP
I FORMAT (/,lh ,'HMIN OENCM IS ZFRO'J
5 FORMAT! ICAl , I 10,315 )
f> FORMAT (22X, 1F4.1, 13X.IF5.1)
C
C READ IN HEADER CARD AND RAGB DATA.
RFAC 5, ISTA, ICAT,MAXT,M INT, IDSC
ICS=IDS*ICSC
I CATF1= IDAT
tVR= ICAT+0.0001
JMri= ITAT-I YR*l 0000
iMn=jMn*n.oi
OAY=30.5*(I Mn-
MAXT=(MAXT-32.0)*r.55'i5+273.7
MINT= (MINT-3?.0)*r.5555*273.7
DO 25 L=ltlCO
^EAC 6,T(L ),P(L )
TlL)=T(L)+273.2
IF (PILJ.EO. 500.0) GO TO 27
25 CCNTINUE
27 NSL=L
C
C DETFRMINF ThF PRFSSUPF LFVEI AT WHICH THE PTTfNTIAL TEMPERATURE
C (THFTA) IS GREATER THAN HP EQUAL TO THE S FC MAX POTENTIAL T FMPER ATU*F .
60 IF ( FLOAT! MAXTJ.LC .T( 1 ) ) GO TO 70
GC TO 80
70 HMAX=C.C
GC Tn 120
PO TETfX=((1000.0/Pll))**0.2«6)*MAXT
OH 10 f I=1,MSL
THFTAII ) = T( I)*(1000.0/P(I )»**0.?86
IF (THCTA1I ).LT.TFTMX) Gn TO 105
GO TO 100
100 IS=I
GO Tn 1 10
115 CTNTINUF
r,c TC us
r
f CALCULATION OF MAX MIXING DEPTH.
110 IS1 = IS-1
PM=P( IS1 )+(P( IS)-P( I SI) )*(TFTMX-THFTA( I SI) ) /( TH^TAt I SI- THFTA (I SI ) )
H"4X=14.7*(MAXT*T ( I S ) )* ALOG ( P ( 1 )/PM )
IF (HMAX.GT.ACCC. D GP n 115
GO TH 120
115 HMAX=4000.0
120 TOP=( T( 2)-T( 1) ) /I P( 2)-P( 1))
DFN = 0.5*TOP*(P(2 )+P ( 1 ) 1-0.143*1 T( 2)*T ( I) )
IF (DFN.NF.C.O) GC TO 130
P^IfJT 1
78
-------
STOP
C
C CALCULATION OF MIN MIXING DEPTH.
UO HMIN»-14.7*(T(2)*T(I))»{I 0.0633-0.298*TOP)/DEN)*PQP4
IF (HMIN.LT.0.0) GO TO 135
IF (HMIN.GT.4000.01 GO TO 135
GO TO 140
135 HHIK=4000.0
140 IF (HMAX.LT.HMN) HMAX=HMIN
RETURN
END
79
-------
SUBROUTINE STORE
C
C THIS SUBROUTINE COMPUTES AND STORES X/Q AND 0 VALUES (IF ITVP FQ 2),
C
CC«*CN/CPTS/ NPT,XPT(6?5),YPT(625)
COMMON/CSAV/ J5,L7,I9,K36,HMI,SXOQ<320),S012«32*D),S03
DO 200 N=1,NPT
N8S=(N-l )*K36»I9
DO 200 K=1,K36
KRS=(K-1)*I
-------
SUBROUTINE SFCTB1
C
C THIS SUBROUTINE REAOS IN THF FIRST SURFACC (IBS CARD WHICH
C CONSISTS CF THE 2400, 0100, 0200, AND 0300 HOURS ORS.
C
CO^MPN /SFC/ IHRJ241,ICHT(24),IWD(241,WS(241,1010(24),1TFMP(2*I
b FORMAT (6X,I
-------
SUBROUTINE SFCOB?
c
C THIS SUBROUTINE REAOS IN THE REST 0* THE SURFACE OBSERVATIONS
C FOR THF DAY.
C
CCMMCK/CPTS/ NPT,XPT(ft25)»VPT(625)
COMMON /PNT/ IC,RT2,I JSI24I,HS(24),HR(168),CCAl(1200)
CCVMCN /SFC/ IHR(2^J, ICHT«2*)t IWD|24),WS(24»t ICLl)(24),
5 FORMAT (6X, I
-------
SUBROUTINE EXTURB (EXTRAQ)
C
C THIS SUBROUTINE DETERMINES THE MAX WIND SPEED FOR THE DAY AND
C OBTAINS THE WIND DIRECTION VECTOR IALLY. IT THEN CALCULATES THE
C EXTRAURBAN CO CONTRIBUTION ON THE BASIS OF GASOLINE GOf;SUrt*VEON IN
C THE SECTOR.
C
C THE FUEL-ARRAY CONTAINS THE RATE OF GAS CONSUMPTION IX THE
C SEGMENTS AS A FUNCTION OF WIND DIRECTION.
C
CCMMCN/CFUEL/ FUEH 16 J, NHOL U0»
COMMON /HMM/POP4,P(25»,T125I,MAXT,MINTtHMAX.HMlN.NSL,SAP
COMMON /SFC/ IHRJ24),ICHT(2*ltIWDC24),WS<24)tICLDC24)tITEMP«24)
1 FORMAT «/,lH ,«MAX WI ND IS CALM.')
C
C DFTFRMINE TH.E 24-HOUR MAX WIND SPEED.
WSMAXxO.O
DO 100 IM,?4
IFIHS1IJ.GT.WSMAXIWSMAX'WSII)
100 CONTINUE
C
C DETERMINE THE COMPONENTS OF THE AVERAGE WIND DIRECTION FDR
C THF MAX WIND SPEED.
NWD=0.0
XWS=0.0
YWS=0.0
00 105 I-It 24
IF (WSm.NE.hSMAX) GC TO 105
WDHAX=IHD(I)*0.0174533
NWD-NWD+1
XWS=XWS+SINIWDMAX)
YWS»YWS*CnS(WDMAX»
105 CONTINUE
C
C DETERMINE THE AVERAGE WIND OPFCTION ACCORDING TO 16 DIRECTIONAL
C POINTS IN RADIANS.
IF (NWD.EO.II GO TO 135
IF (YWS.NE.O.OI GP TO 130
11C IF (XWS) 120,115,125
115 PRINT 1
STOP
120 WDMAX=4.7123889
GO TO 135
125 WDMAX=1.5707963
GO TO 135
130 WD«AX = ATAMXWS/YWSl
IF (YHS.LT.C.O) WDMAX=WDMAX»3.1415927
IF (WOMAX.LT.O.C) WDMAX=WDMAX*6.2831853
135 IOIR=2.54647
-------
SURPDUTINF MINW IN ( IT )
r
f THIS SUBRHUTINP CHECKS THF HfURLY WIND SPEFO ANf) SETS EACH WIN1
C SPEED LESS THAN 1.0 M/S TD 1.0 ( FOR WIND SPEEDS EQUAL TO 0 THE
C 'A! I NO DIRFCTIPNS APE SFT ECUAL TO THAT DF THF PREV IHIIS HOUR).
C
COHMPN /SFC/ 1HRI2*», ICHT«
DATA IOP/27C/
* I=WSI IT)
IF (WI.OE.1.0) GT TO 105
IF (WI.LF.O.n) IWn(IT)=IOP
WS( IT )=1.0
IT5 IDP=IKD(IT)
»
END
84
-------
SUBROUTINE STABLE (I , SAL , IJ)
C
C THIS SUBROUTINE DETERMINES A STABILITY INDEX THROUGH A SERIES
C OF CRITERIAS CONCERNING CLOUD COVER, WIND SPEED, AND SOLAR
C ELEVATION (INSTABILITY INDEX).
C
DIMENSION IX(25),HCOSt25)
CCMMON/CITY/ ITYP, ISM, IOUT( 10»
CCMHON/CSTAE/ SP1BO.CP180
COMMON/M I X/ 1 STA ( 10) , I YR , I MO , I DA ,DAY , I OS , I HO , I ST « I DA TE , I DATE1 ,
X IDATE9
COMMON /SFC/ IHR(24),ICHT(24)tIMDf24i,WS(24)( ICLD(2*I, ITEMPI24I
COMMON /PNT/ lC,RT2,lJS(2
-------
200 IF (WSP.GT.3.C) GH n 210
IWS = 1
GC rr 300
2in IF (WSP.GT.fc.fl Gr TO 2?0
IUS = 2
GT TC 3CO
220 IF (WSP.CT.10.0) GO TO ?30
iw<;=3
00 TO 30C
23? IF (WSP.GT. 12.0 I C-n TO 240
I US = <,
GO TO 300
243 IWS=S
3CO IFX=(rhS-l)*3+IRAn
IJ=IX(IFX)
GC in 310
CALCULATIHN OF NIGHTTIME STABILITY
305 IF (WSP.GT.6.0) Gf TO 310
IF (CC.GF..0.5.ANn.HS".GT.3.0) GO Tn 310
IJ=5
310 IJSII>=IJ
RETURN
ENP
'86
-------
SUBROUTINE DEPTH (IT.SAS.HT)
r
C THIS SUBROUTINE CALCULATES THE ATMOSPHERIC MIXING DEPTH FOR EACH
C HOUR. AT SUNSFT, SUBROUTINES RAOBHMM AND SFCDBS1 ARE CALLED INITIATING
C CALCULATIONS OF NIGHT-TIMF MIXING DEPTHS.
C
CCMMCN/CITY/ ITYP,ISM.IOUT!10)
COMMON /DAY/ lr)V»,NH,lDAHO,lHW,I DWT
CCMMCN/MIX/ISTA(10),IYR,IMO, IDA,DAY,IDS,IHQ,I ST.I DATE,I DATE 1,
X IOATE9
COMMON /HMM/PPPAfP< 25), T I 25) , MA XT.M I N T.HMA X ,HMI N ,NSLf SA »
COHMCN /SFC/ IH«(2A),ICHT(2A),IW0(2A),HSC2A ) ,ICLDI3*),I TEMP!2AJ
CCMMON /PNT/ 1C ,RT2,US(2A) ,HS!2A) ,HM168) ,CCAL(1200)
DATA HSUN.HPSUN/JOOO.O, lfl.0/
DATA ITIM/0/
1 FORMAT!/, 2X,'nATE = •,I X,I 6,2X,•IHHL)•,AX,•NO OF RAOB LEVFLS = ',
1 [3,AX,'SFC PPFSS = • .F7.1.AX,"SFC MAX TEMP =•,14,AX,•SFC MIN TEMP
2 = •,1 A,/)
2 FORMAT!/, ZX.'OATF =',IX,I6,2X,MMON)•,4X,•NO OF RAOB LEVELS =',
1 I3,AX,«SFC PRFSS =•tF7.1,AX,•SFC MAX TEMP = • , fA, AX , • SFC MtN TEMP
2 =',M,/>
3 FORMAT!/, 2X, 'DATE = • , IX, I 6, ?X, •( TUF. J •, 4X ,' NO OF RAO* LEVFLS =• ,
1 I3,4X,'SFC PRESS =• ,F7.1 ,ATE =' , 1 X , 16 ,2 X, • ( SUN) • , AX , ' NO OF RAOB LEVELS =',
1 I3,AX,'SFC PRFSS = • ,F7. 1 , AX,•SFC MAX TFMR =• ,I A,AX ,'SFC MIN TEMP
2 =• ,IA,/)
9 FORMAT! 1H ,63XtlCAl,/)
ITT=ITEMP(IT)
IF (ITT.LT.CINT) ITT=MINT
C
C IS IT DAYTIME I SAS GRE4TFR THAN 0).
IF (SAS.GT.C.O) GC TO 100
C
C IS IT SUNSET
IF ISAP.GT.O.r) GT TO 105
GO TO 120
C
C LINEAR TEfPERATURF I NTf-R POLATI ON IS USEO TO CALCULATE THE
C DAYTIMF MIXING DEPTH.
100 HT=HPIN+IITT-MINT)*(H^AX-H"iN)/(MAXT-MINT)
GO TO 13?
105 HT=hMIN + ( ITT-VINT ) *(HMAX-HMIN)/IMAXT-MJNT)
HSUN=HT
HRSUN=IT
87
-------
C
f.
C
108
AT SUNSFT, PP1NT THE SAS 1C STATION DATA, "FAD THF NEXT SOUNDING
AND CALCULATf THE VAX AN9 MIN FIXING OF<>THS, AND READ THE NFXT
DAYS F1PST FOUR HOURS OF SFC DATA.
IF UrUTlll.LT.l) GO TO 115
NH.l) PRINT 9.ISTA
( ITAHO.FC.Ot GO TO 1 10
i , i OATI ,NSL ,p(i ) ,MAXT ,MINT
115
O, 11, 12, 13, 14, 15, 16), IOW
2,IDATF1,NSL,P< I) ,MAXT,M1NT
115
3,irATFl,NSL, P(l ) ,MAXT,MINT
115
4, ICATF1.NSL, ?( 1 1 ,MAXT,MINT
H5
5, IDATF l,NSl ,p( 1) iMAXT.MJ NT
115
6.IDATF1.NSI ,P( 1) .MAXT.MUT
115
7,irATFl,N*l,P(l) ,MAXT, '4INT
115
«, ICATF1.NSI ,P( II ,MAXT,MINT
SACBHW
CALL SFCPPl
ITIM=1
GO TO 135
1 10
10
11
1 2
1 3
1 '
1 5
16
11 5
IF (I
IF ( I
GO TO
GC TC
PRINT
GC TP
PR I NT
GO TO
PRINT
GO Tn
P3 INT
GC TP
POINT
GC Tn
"RI NT
r,n TT
PRINT
r ALL
C TI^F INTr°PClATirN IS USED TO CALCULATE THE MIXING DEPTH
C 4FTFP SUNSET.
l.'O IF ( IT.GT.15) GO TO 125
GP TT 130
12 b HT = HMIN*(HSLN-HMI M * ( ?^ -I T ) / ( M-HPSUN )
GC TC 132
13C HT=hHIN
13? IP( IAf»S( ITVP).F0.3)Gn TO 1TR
135 IF (HT.LT.50) HT=?0.0
HS( IT ) = HT
END
88
-------
SUBROUTINE CALXOQ (HT.JI)
C
C THIS SUBROUTINE CALCULATES THE X/0 VALUES FOR EACH Of THE
C NINE SECTOR SEGMENTS USING THE GAUSSIAN AND BOX MODELS.
C
DIMENSION R (10) ,A(9,6),B(9, 6),AA(54),B6(54)
COMMON /OUF/ OL(IO) ,QM( 10),012(10),034(10),06(10),XQQ(IO)
COMMON /ALLC/ YS(35 1),IYN(351)
l,XR(10),YP(10),AR(10),EFAC(9),ACR,LR4,LR5fLR9,LP10,AG,YCl,YC2tXRCf
2YRC ,XC2,B>2,'m,BI2,Bll,CRl , I T, I J,XG,YG,SA,CA.NN,CR ,LRL
ECUIVALENCE ( A, AA ), ( B.BB )
DATA R /I.0,125.0,250.0,500.0,1000.0,2000.0, 4000.0,
19000.0,16000.0,3?COO.O/
C
C VALUES USED FOR THF STABILITY PARAMETERS (A AND B).
DATA AA /O.00,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07 ,
2 0.00,O.I2,0. 12,0.12,0.12,0.12,0.12,0.12,0.12,
3 0.00,0.23,0.23,0.23,0.23,0.23,0.23,0.23,0.?3,
4 0.00,0.50,0.50,0.50,0.50,0.50,0.50,0.50,0.50,
5 0.00,1.35,I.35,1.35,1.35,1.35,1.35,1.35,1.35,
6 0.00,3.00,3.00,3.00,3.00,3.00,3.00,3.00,3.OO/
DATA 68 /O.00,1.28,1.28,1.28,1.28,1.28,1.28,1.28,1.28,
2 0.00,1.14,1.14,1.14,1.14,1.14,1.14,1.14,1.14,
3 0.00,0.97,0. <57,0.97, 0.97,0.9 7, 0.97 ,0.97 ,0.97,
4 0.00,0.77,0.77,0.77,0.77,0.77, 0.77,0.77,0.77,
5 0. 00,0.51,0. 51,0.51,0.51,0.5 1,0.51,0.5UO.51,
6 0.00,0.31,0.31,0.31,0.31,0.31,0.31,0.31,0.31/
HP2R=0.7978S5*HT
DO 140 1=1,S
UP2R=0.797885/AR( I )
IJH=HT*fiR(I)
IPl=I + l
BIJ = R ( I,J I)
BIJl=l.O-RIJ
AIJ=A(l,JI)
IF (I .GT.l ) GO TP 25
AIJ=A12,JI)*R(2)***(2»JI)
GO TH 100
25 RT=(HP2R/«IJ)»* (1 .0/8 IJ )
IE (P(IPl)-RT) ICC,100,115
100 IF (BUI) 110,105,110
105 XOOU) = (UP2R/AIJ)*ALOG(R(IP1)/RI I))
GO TO 140
110 XOQII )=(UP2R/ ( A IJ*B IJ 1) )*(B ( IP1 )**B IJ l-R( I > **RI Jl)
GO TO 140
115 IF IRII 1-RT) 125,125,120
120 XOOd )=(P(IP1)-R(I) )/UH
GO TO 140
125 IF (BIJ1) 135,130,135
130 XOOd I = (UP2P/AI J)*ALOG(RT/R( I I ) +(R( IPU-PTI/UH
GO TO 140
135 XCQ(I)=(UP2R/(AIJ+BU1) )*(RT*+BIJ1-R( I)**BI J1)*
-------
SUBROUTINE CALOUF
-------
Bl=BYl
IF (VS2.LT.YC2) B2=BY1
224 IF (ABS(XSD.LT.B1*YS1) GO TO 230
IF (ABS(XS2).LT.B2*YS2) GO TO 232
IF (XS1*XS2.LT.O) GO TO 23*
IF (YS1.GT.YC2.0R.YS2.LT.YC2) GO TO 300
IF |ABS(XS1).LT.XC2.0R.A8SIXS2).LT.XC2) GO TO 234
GO TO 300
23T IF =1
236 IY1=YS1*CRI+1
IY2=YS2*CRI+1
IF ( IY2.GT.3SO) IY2 = 350
IY11)=IYN(IY1)
IF (YSl.GE. YS(I YD) I Y( 1) =1 Y( 1) »1
IY(2)=IYN(IY2)
IF (YS2.GT.YSIIY2)) IY(2)=IY(2)+1
C
C LOCATES THE INTERSECTIONS OF LINK WITH SFCTOR DIVIDING LINFS.
N = 2
IF ( I Y( D.EO.IYI2)) GO TO 250
L=IY(1)
238 YP=YR(L)
XP=AL+BL*YP
IF (ABS(XP),GF.XP(LI) GO TO 240
N=N+1
X IN )=XP
YIN)=YP
JY(N)=L
IY(N)=L
IF (L.EQ.LR4.AND.ABS(XP).GT.XRC) GO T1 240
N=N«-1
X(N )=XP
YIN )=YP
JYJN)=L*1
IY(N )=L+1
240 L=L*1
IF (L.LT.IYI2)) GO TO 238
91
-------
C LOCATFS THE INTERSECTIONS OF LINK WITH SIDES OF SECTOR.
250 K = l
260 IF ( ABSIX(K) ).LT.RY2*Y1K) ) GO TO 280
SR=1.0
IF (X(K).LT.O) SB=-1.0
BI3=SB*BI2
XP=AL/(1.0-eL*Bl3)
YP=flI 3*XP
IF (YP.LT .YPC.OR.YP.GT. YC1) GO TO 270
I YP=YP*CRI+l
IY(N)=IYN(IYP)
IF IYP-YSUYPU 266,262,264
262 IF (K.FQ.l) GO TO 26*
GC TO 266
264 I V(N) = I YCM + 1
266 X(N)=XP
Y(N)=YP
JY(N)=I Y
YP=813*XP
IF I YP.GT.YRC.OR.YP.LT.O) GO TO 280
IYP=YP*CR 1+1
N = N + 1
I Y(N)=I YNII YP)
IF (YP-YSI IYP ) )276, 272, 2Tt
?72 IF (K.EO.l) GO TO 274
GO TH 276
?74 IY(N)= IYJN)+1
276 X(N)=XP
Y(N(=YP
JY(M=IY(N)
?80 K=K+1
IF (K.LT.3) GO TO ?6T
r
f COMPLIES EMISSION CC NTR I BUT I ^N OF LINK TO EACH SEGMENT
f (RESULTS STGRED If Q).
DC 295 Nl=l ,N
IF ( JY(Nl).tO.O) GO TC 295
NS1=IY(NI )
IF (NS1.GT.LRL) GT TO 294
NN1=N1*1
DO 292 N2= NN1.N
IF ( JY( N2I.EC.O) GH TC 292
NS2=IY(N2)
IF (NS1.NE.NS2) GO TT 292
XDN=X(Nl)-X(N2)
YON=Y(N1 I-Y1N2)
OT = SORT ( (XnN*Xl)N + YnN*YD
3SO IF (IFKINO.LT.3) GO TO
GO TO 286
285 Q12(NSI)=C12(MS1H-QT
GO TO 290
290 JY(N2)=O.C
92
-------
GC TO 294
292 CONTINUE
29* JV(N1)=0
295 CONTINUE
C
C tND OF PAIN LOCP
30C CONTINUE
C
C COMPUTES CO EMISSION WITHIN SEGMENTS FROM MINOR LINK DATA.
305 CALL CALLCCIM)
RETURN
END
93.
-------
SIBRPUMNE C AILCC ( Ml
C
C THIS SUBROUTINE USES THE GRID POINT VALUES < VCL AND VCM) TO
f COMPUTE THF AVERAGP CC EMISSIONS WITHIN EACH StCTOP SEGMENT. THE
C RESULTS &RF STHRFD IN Ql AND C|M.
r
. L
Y/ ITYP, ISM, IOUT1 10)
K'PT,XDT(625) ,YPT<625)
QU id ,QMI io) ,oi2( 10) ,Q3M 10) ,Q6(io) ,xno(io)
COMMON /CLOQ/ VCK2500I , VTM( 2500 ) , CB( 3), SB I 3) ,NRW,NCt tMR W.MCL , CLE ,
IGSPtRRSt ZOX
CCMMHN/Attr./ YS( 351 ), IYN( 351)
l,XR(in),YP(10),AR{lO) ,EFAC(9),ACI?,LR'f,LH5,LP9,t1lOtAG,YCliYC2,XRC ,
2YRC ,XC2,B>2,HY1,BI2 ,BI1 ,CRI ,1 T, IJ,XG,YGtSA,CA,NN,CRfLRL
XG=XPT(N)
YG=YPT(N)
nn 30 L = I,C'
YRL=YP(L )
SEVL=C.O
xrs=o.o
KTN=1
LD=1
IF (L.LT.5) Gf TC •>
KTN=3
LD=L-3
WfT=0.25/LO
DO 20 J=l,Lfl
FXOS=1.0-WCT*( 1 .0*2.0*1 LO-J ))
YRR=YRL*FXDS
SSfVL=T.O
TO 10 K=1,KTN
YM=(OIVI*(C**CR(K)-SA*SB(K) )»YG)*ZOX«-MRW
I Xf = I NT ( Xf)
IF ( I YM.LT.l.rR. I VM.GE.NRW) GO Tn 10
IF ( IXM.LT.l .DP.IXM .GE.NCL) GO TO 10
D YM=YM-I Y"
OXM=XM- IXM
IXY=( IXM-1) *NRW+IYM
I XY1=IXY+1
TTFVL= inYMl*VCI ( IXY)*OYM*VCUI XY1)
1 +( CYM1*VC.L ( IXY+ NRW)+OYM*VCL ( IXY1*NRW)
SSFVL=SSEVL+TTF VL
IK( ITYP .LT.O)GO TH 6
IF (I .LF.LRL) GC TO 10
TTHVM= (OYMl*vr>( IXY) »OYM*VCM( IXY1
1 +( CYMl*Vf.M ( IXY t NH W) *OYM*VCM( I XYUNP W) )*nx«*
5 St Vf = SSE W+TTPVN1
n CONTINUE
SFVL=SEVL+SSFVL *FXDS
94
-------
XOS*XDS*FXDS*KTN
20 CONTINUE
OUL) = SEVL*ARtL )/XOS
lF(ITYP.n.O)OLm=SEVM*ARtL)/XDS+QUU
IF(L.GT.LRL)OL(L)=SEVM*AR(L»/XDS+QLIL)
30 CCNTIMJF
RFTURN
END
95
-------
F LOCXCC(HT,J)
c
f THIS SUBROUTINE LOCATES THE PRHPFR X/0 VALUFS FROM STORAGE.
r
COHMON/CSAV/ j5,L7,I9,K36,HMI,SXnQ< 320) ,5012(3240) ,S034
1 SCH3240)
f.CMMPN /OUF/ CL(IO) ,Qfl 10), 012(10), 03*1101, 06(10), XOOdO)
no in L=l,L7
IF (HT.LT.HB) GO TO 12
IT H6=HP+HP
12 IF 1L.GT.LT) L=L7
OP 100 1 = 1, 19
IBS=I+LBS
100 XGQ( I )=SXrlO( IBS)
RETURN
END
96
-------
SUBROUTINE LDCOUE (KD.N)
C
C THIS SUBROUTINE LCCATfS THE PROPER 0 VALUES FROM STORAGE.
C
COMMON/OS AV/ J5 ,L7 , 19 ,K36,HMI ,SXnO( 320) , SO 121 3?40) , 5034(32*0) ,
L SOU3240)
COMMON /OUE/ QL( IC),OM( 101,0121 10), 034(10), 06(10), xnOllO)
IF (K.LT.l) K=K36
K BS = ( K- 1 ) * 1 9* ( N - 1 ) *K 36* I 9
DO 100 1=1,19
IBS=I+KBS
Q12II )-S012f IBS )
034(1 MSC34UBS)
100 QL(I)=SOL(IBS)
RETURN
END
97
-------
SUBROUTINE CALCCN ( I , N, HT ,F XTR AC )
f
C THIS SUBROUTINE COMPUTES THE CO CONCENTRATION AT THE
C RECEPTOR POINT RASED ON THE VARIOUS MODELS. THF DIURNAL TRAFFIC
C PATTERNS FTP KEEKCAYS, SATURDAYS, SUNOAYSt AND HOLIDAYS A3F INCLUDED
C IN THE CO CALCULATIONS.
C
CCMMCN/CITY/ ITYP,ISM,mUT( 10}
CQMMCN/CPTS/ NPT,Xi>T(625),YPT(625)
/PAY/ IOW,NH,IDAHO,IHW,IDWT
/ISTAUOI , IYP, IMO, IDA, DAY, f OS , I HO, 1ST, IDA TE , I DATE 1,
IQATF9
/PNT/ ICtPT?, IJSI24) tHS( 24)tHR(168),CCAL(1200)
/SFC/ IHR(24 I, ICHT(?4», lwr>(24l ,WS<24), ICLDI24), ITFMP124)
COMMON /CAR/ PT12I24I ,PT34( ?4I ,PT6l 24) , PT SA T(24| , PTSUNl 24 )
CTMMON /OUF/ OL(10),QM( 10), 0121 10), 0341 10) , 061 10J , XOQ( I 0)
COMMON /ALLC/ YSOS1I ,IYN(351I
1,XR(10),YR( 10),AP(10),EFAC(9) ,ACR,I R4.LR5 ,LR9 ,LR1 0, AG ,YC1 ,YC2 tXRC ,
?YRC ,XC2,BY2 ,RYl ,11? ,RI1 ,CRI, IT, U , XG, YG.S A, CA.NN, CR ,L RL
CO«MON/CN24/ N?4
IT=I
IF (N.EQ.l) IC=IC+1
in CPL=O.O
WSI75=0.75/(223.6932*WS(IT) I
inv IF( ICV.T-6 1100,105,1 10
r
C WEEKDAY CONCENTRATION CALCULATIONS.
100 IF( ITYP.LT.OIGO TO 104
FACT = 012( J I*PT12 (I TT)t (0341 J) +QLI J ) ) *PT"?4 ( IT T )
GO TO 120
1 04 FACT = OL ( J)*FT12 ( ITT1 + CH J)*PT14< ITT |
GO Tn 120
r
r SATURDAY CCNCENTRAT ION CALCULATIONS.
105 FACT = PTSATt ITT) *[Q12( J)+Q^4( J) *OLIJ))
GO TP 120
r
C S'JNDAY AN" HOLIDAY CONCENTRATION CALCULATIONS.
Ill FACT = PTSUNI ITT) * 1 01 2 ( J) +Q34U ) +01. ( J ) I
1?0 CPL=rPL+FAf T*XOO( J)
CPL^CPL/WSl IT )
IU = N
IF{ I ABS( ITYP) .L T. 3 It 0 = 1 T+(N-I I*N?4
P CCAL ( ID)=(f PL+EXTPAO 1*1001.0
13T
END
98
-------
I2=72C
14=240
30 13=12+240
15=14+240
CCAU ID+I2)=DBAS*XL UCBG
CCAUIO+I3)=DBAS*XL2+COG
CCAUID+14)=DBAS*HI*HZ1+CBG
CCAUID+I5)=DBAS*WI*HZ2»CBG
50 ITIHE = 1
RETURN
END
99
-------
TI WF STRFETII ,M
C
C THIS SUBRCUTINf COPIJTES THE STKEET PPHFILF DF CH CONCENTRATIONS
C ON THF BACKGRniNO VALUFS, THE STRFHT FVISSIPN VALUES, AND THE
C ^FTfTunLOGlCAL CO.'OITirNS ("ESULTS Apr STDRFO IN ARRAY CCAL I .
r
n^PM/CST"/ VCAR(IO) tASTI 10) .WkSTI 10) t ICSTI 11) .NOCnilO )
Cf^MPM /PAY/ IDW.H'H, irnHfl, IHW, IDnT
CCw»t -N/MI X/ISTAdO) ,IYP ,IKl,IDA,l)AY , IDS, 1 HD , 1ST, I DATE, I DATE 1,
X
/CAP/ PT12 124) ,°TT»(24) , PT6 ( 2^ I , PT S AT ( 24 1 , PTSUN ( ?4 )
/PNiT/ 1C ,^T?,1 JS12<») ,HS» ?4> .HRJ1GS) ,CC4LU200)
/ALLC/ YSI351) , IYNI351 I
l.XP ( 10) , YP( 10) , AP(10) ,FFAC(9) ,ACK,LR4,LPS ,1 R?,LP11, AG,YC1,YC2,XMC,
?VKC,yr?,RY?,flYl,P17,HIl,CRI,IT,IJ,XG,VG,SA,CA,NN,CP,LKL
TATA Z3,Z1,2P,X, XLO.CK, ITIME/38.8,3.65,??.3,8.,2.,7.,0/
IF ( I T IMF.FC.l | GC TO 10
XL1=1.0/ISPCT(ZI*ZH-X2)*XLO)
HZl=l Z3-Z1 )/Z3
HZ2=( Z3-Z2) /Z?
n in= i + t\-i)*N?^
CHG^Cf AL( ID)
«=«AST( N)
I L S = I C. S T < M
w ! = i . n / w
1 TT=I+I PS
IF ( ITT.GT. 2M I TT^I
I F ( lO'T-61 12, l^i, 16
i; PT'l=PIh( [ TT)
sn ir 1 3
[TT )
TSiHl I
C »b = Vf »•< ( M
IF (AM).LT.ASHN) ) A'«in= AWD+-360.0
AtV,-A'AL>-AST (N)
IF ( AMG.GT. HC.O) ANG=A»jr,-360.0
IT ( AANG.GT.9C. 0) A ANG= 13 0. 0- AANG
Ir ( AAMG.CT .30.1) OP T'"1 "> 5
CTAl ( I !"l»2Ai) = rP AS*( XL i + WI JttBG
real i ir: + 7?:i i = r^As*i XL i+wi )+r.°G
1^=72 )
IF tAv-,.LT.c.r) r,r
100
-------
SUBROUTINE PPDATA
C
C THIS SUBROUTINE PRINTS AND PUNCHES OUT THE STATION IDt DATE t NAX-MIM
C TEMPEPATURFS, ALL SFC DATA, STABILITY INDEX, MIXING DEPTH, AND TM€
f CALCULATED CO CONCENTRATIONS.
C
COMMON/CITY/ I T YP ,1 SM , I OUT( 10)
COMMON/ CPTS/ NPT,XPT(625),YPT1625)
CCM^riN /CAY/ IDW.NH , I CAHO , I HW , I OWT
CCMMON/MIX/ISTA< 1C) ,IYR , I MO , I DA tDAY , IDS ,1 HD , I ST, I DA TE , I DATE1 ,
X IDATF9
CTMMRN /PNT/ 1C ,RT2 ,IJS(?4) ,HS(24 ),HR(168),CCALU200)
COMMON /SFC/ IHRl 24), ICHTl 241,1 WDI24) ,WS(24) ,ICLni24) ,1 TEMPI?*)
CCMMON/CN24/ N24
DATA ITIM/0/
1 FORMAT (2X, 'HOUR CL D-C TEMP WND-D WND-S STB-I MI X-0 CO-PC
X CO-R1 CC-R2 CO-LI CO-L2 • )
2 FORMAT(2X,«HOUR CLD-C TEMP WNO-0 WNO-S STB-I MIX-D',
X ' CO-HG' I
3 rORMATl2X,I 4,17 ,15,17 ,F T. 1 f I 7 ,F8. 1 ,'>F7. 2)
4 FORM AT (2X, I A, 1 7, I 5, I 7,F 7. I, I 7 ,F 8. 1, F7. 2)
«3 FCRMATI6X,1 STA',14)
10 FORMAT ( A10.I 10,AI<,,FA. 0, 14 ,F 5. 0.5F6.2 I
11 FOPMATIA10, 110,314, F4.0.U.F5.0)
12 FCO^ATC CITY DATE HR PT CLO DIR SPO SI MD CRG
1CP1 CR2 CL1 CL21)
13 FCRMATC CITY DATE HR CL 0 Of SPO SI MD • )
14 FORMAT! R( I4,F6.2) )
IF( IOUT{ 1 ).LT.l 100 TO 100
DC 20 I =1 ,NPT
IC3 = ( I-1)*N'24
PRINT 9,1
IF( lARSf ITYF).EQ.3)GO TO 6
[F( ISM.F3.01GO TO 6
PRINT 1
DC 5 J=l,24
5 PRINT 3,J,ICLD( J) ,ITFMP(J),lWDt J) ,WS( J) , US ( J ) »HS ( J ) , CC AL I I C8+J ) ,
X •CrAL«ICH + (240*J ) I.CCAl. (ICH + I480+J) ) ,CC AL ( I C^H 720* J) ),
X CCAL UCB + ( S6C* J) 1
r,n T° 20
6 PRINT 2
IF( IABS! I TYP).FQ.-5K,0 TO 8
00 7 .1=1,24
7 PTNT 4,J,ICLC(J),1TFMP(J),IWDU),1WD(J),WSU),US(JI,HS(J),
X CCAL (IC,1*J )
Gr TH 20
PRINT 4, IHD, ICLDI I HO) ,1 TE"r>( I HO ) , IWO( IHD) ,WS{ IHD) , US (I HO),
X HS( II-D) ,CCAL ( ICA I
20 CTMTINIJF
r
c PUNCH PUT TUTDUT DATA IF IOUTI?) GT ZFRO.
1C) IF ( I OUT (2) .LT. 1) GO TO 200
IF( IABM ITYP).GT.?)C.n TO 150
IF{ ITIM.NE. 1 (PUNCH It
on ior 1=1,24
DP 1">S N=1,N"T
IC^=( N-I )*N?4+I
101
-------
PUNCH 10, IST'MCATF, I,N, ICLDI t ), IV«nm,HS( I ), tJSI I ) ,HS( I ),
1 C<~AL ( ICH),Cr.AL ( I CO* 240) ,CCAL 1 1C^*'»80 J ,CC AL (ICB + 720) , CC AL (ICR+-960 )
GT in '00
I=inn
I F( ITIM.NC. ] 1PUNCH 13
PUNCH 11, 1ST \.ICATI-, I.ICLDI I), I WO (I ),HSU ), IJSl I) ,HS[ I)
°UNCH 14, (N,r.CAl ( M ,\ = 1 ,NPT)
IF (inum).LT.H r,0 TH 300
IF( IAHSI ITYP) .GT.2IGT TD 250
TO 205 1=1,?^
DC ?0* V=1,NPT
ICP=( N-1)*N24+I
^"ITE I u),iwDm,wsM),iJS
-------
Appendix C
EXAMPLE OF INPUT DATA
BASIC INPUT INFORMATION—CARDS A THROUGH M
1 1
1997.0 2071.0
2001.0 2073.0
38.6
11.02
4.10
43.0
2
1
I
0.0121
0.0659
0.0796
0.017
0.050
0.076
0.0121
0.0659
0.0796
0.015
O.OSd
0.062
0.014
0.028
0.056
710826
71 0906
i Data Kead-i n D
197.0 14692. J
2e)7.0 6083.0
2.36 1931.5 2040.5
15.56
6.19
9.0
2
I
1
0.0080
0.0606
0.0554
0.010
0.045
0.061
0.0080
0.0606
J.0554
0.010
0.062
0.059
0.008
0.036
0.052
710827
15.17
4.78
53.0
2
2
1
0.0046
0.3634
0.03S4
J. 000
0.046
0.053
0.0046
0.0634
0.0354
0.007
0.064
0.051
0.006
0.036
O.J42
4
12.20
4. 78
20.0
Z
2
2
0.0042
0.0628
0.3280
0.005
0.044
0.052
0.0042
0.0628
0.0280
0.058
0.068
J.042
0.005
0.041
0.034
8
y c.' nrout me
24.3
11!. 5
5.0
11.37
6.24
12.0
2
2
2
0.0064
0.063d
0.0270
0,010
0.045
0.043
0.0064
0.063d
0.0270
0.011
0.066
0.034
0.006
0.046
0.030
1
l^Urtl
6
7
700.0
7.o9
3.29
9.0
<:
2
2
0.019
0.0660
0.0268
0.034
0.050
0.038
0.019
0.0660
0.0268
0.020
0.062
0.032
0.010
0.048
0.024
1
2
-0.75
8.59 .
4.19
5.0
1
2
2
0.0514
0.0658
0.0246
0.062
0.058
0.034
0.0314
0.0658
0.0246
0.039
0.062
0.030
0.014
0.053
0.018
4.2o
11.49
1
1
2
0.0726
0.0782
0.01(44
0.064
0.071
0.025
0.0726
0.0782
0.0184
0.054
0.062
0.022
0.018
0.036
0.017
Card*
A
B 1
B'2
C
Dv'l
D 2
E
F 1 i
F 2
F'3
c'l
C 2'
G 3
H.'l'
H(2>
H 3
I.I1
I'2
I 3l
J 1 i
J,2'
J'3'
K'l
K'2:
K 3
L
M 1 I
See Tables 2, 3, and -I.
103
-------
TRAFFIC INPl'T DATA—OAK!^' N
'Data lU-ad-tn bs Subruul i u- • 1
1016
1006
1000
1001
1002
1003
1004
1005
1006
1007
1073
1000
1001
1002
1003
1004
1005
1007
1008
1008
1009
1007
1009
1009
1012
1013
1013
1014
1014
1015
1015
1016
1017
1017
1018
10 ia
1019
1019
1020
1020
1021
1024
1024
1023
1025
1025
1040
1040
1026
1027
102 7
1028
1028
1029
1029
1030
1030
1017
1016
1001
1002
1003
1004
1005
1006
1007
1008
1022
1021
1021
1020
1019
1018
1017
1015
lOit
1009
1011
1010
1013
1012
1013
1014
1032
1015
1031
1016
1030
1029
1018
1023
1019
1027
1020
1026
1021
1025
1024
1023
1025
1022
1041
1040
1026
1042
1027
1028
1039
1029
1038
1030
1037
1031
1036
1998
2004
1981
1996
1998
2QOO
2001
2003
2004
2007
1932
1931
1996
1998
2000
2001
2003
2007
2010
2010
2014
"2007
2014
2014
2010
2005
2005
2003
2003
2001
2001
1998
1996
1996
1995
1995
1993
1993
1991
1991
1990
1978
197-3
1970
1979
1979
1981
1981
198J
1989
1989
1990
1990
1991
1991
1994
1994
2074
2072
2004
2043
2049
2055
2061
2066
2072
2083
2018
2003
2057
2057
2051
2051
2045
2049
2049
2040
2055
2055
2061
2061
2059
2065
2065
2070
207J
2076
2076
2086
2006
1996
1998
1996
1998
2030
2001
2003
2004
2007
2010
19o6
1990
1990
1991
1993
1995
1996
2001
2003
2014
2044
2041
2005
2010
2005
2003
1998
2001
1997
199d
1944
1991
1995
1990
1993
1989
1991
1988
1990
1979
1970
1970
1979
1966
1973
1981
1988
1975
1989
1990
1982
1991
1984
1994
1986
1997
1989
2068
2074
2043
2049
2055
2061
2056
2072
2083
2098
2008
2045
2045
2051
2057
2063
2068
2084
2095
2098
209f
2^-'dO
2103
21 1 J
2103
2095
21 J5
208t
2096
2074
2086
JO 76
2063
2070
2 J5 7
2U65
2051
2059
2045
2055
^049
2040
2^35
2008
2056
2061
2059
2063
2065
2070
2066
20 76
2072
2086
2077
2096
2087
14692
6083
13333
15600
16400
16400
16400
14300
11600
11600
24095
22500
8000
15500
7000
8000
65JO
1300J
20JOJ
1 7500
1 UOJ
8400
16000
86480
15000
16250
16500
1 7500
20000
1750J
1 3)00
6000
20000
6500
2000J
8000
2250J
7000
^2500
1550J
8JOO
8300
8000
9000
15500
8000
7000
7000
4500
4500
8000
4500
6500
4500
6000
6500
15000
AMD 0
LINKS)
2
2
2
2
2
2
2
2
2
4
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
i
2
1 ^
6
43
6
6
5
6
6
11
10
30
43
7
6
7
7
7
7
6
7
29
33
11
15
10
9
7
11
7
11
7
7
6
7
6
7
7
7
7
13
13
11
6
34
6
7
6
6
6
6
6
6
6
10
6
11
6
Card
Ml!
Nf'21
N(3
N( 561
N( 571
The link data (type N cards) given in this example
represents only a fraction of the total link data
lor St. Louis.
104
-------
THAFFIC I\,M
i \TA
'' C on t i n\i!
2096
2096
2105
2105
2iOt>
21 Oh
2106
209b
2098
208?
20d/
20/ /
2077
20/2
20/2
2066
206f>
2056
2056
2061
206}
2063
20 ob
206d
20/*
20?*
2079
20/4
20d8
2088
2099
2 101
2101
2101
2109
2109
2117
2090
2090
2081
20U1
20/5
20/5
2070
20/0
2O6*
206*
2058
2060
2060
2066
2066
2072
2072
2077
2077
1998
1991
1997
1992
1997
1961
1991
1989
1985
1986
1962
198*
1980
1982
1978
1981
197/
1967
1975
1977
1 9 b 1
1 ^69
19/8
19/0
1980
1972
19st
19?*
1985
19/6
19 ?9
1981
19/3
19/6
19/9
19/*
19/*
19 /*
19 /O
1963
19 ?2
1966
19/0
196*
1969
1963
1967
1961
1963
19*7
196*
19*d
1966
1950
1968
1951
2105
2098
2113
2106
2113
21U9
2098
2087
2099
20 7 /
2088
20/2
20/9
2066
207*
2061
2068
2058
2063
2068
2061
206*
20 /*
20/0
20/9
20/5
2088
2081
2099
Z090
2101
2109
2102
2090
2117
2110
2110
2081
2091
2083
2075
2077
20/0
2072
206*
2066
2058
2060
2066
206*
20/2
2069
2077
2075
2083
20dO
8500
170GO
8500
16500
8500
17000
7500
7500
17000
7500
liOOO
7500
6000
6500
7250
6500
6500
15500
6000
6000
7000
7000
6000
8000
6000
8000
6000
6000
6000
1550J
17000
9000
1 7000
9000
9000
18500
9000
5500
16000
16000
5500
8000
5500
dJOO
5500
7000
5500
20000
6500
25000
6500
8750
6500
8000
6500
8000
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
9 N( 58^
6 N(59)
10 Ni'60)
6
10
12
9
11
6
10
6
6
6
6
6
6
6
6
6
6
6
6
6
6
6
6
10
6
11
6
6
8
6
11
10
6
10
9
6
6
6
6
6
o
6
6
6
6
6
15
6
15
6
15
6 N(112)
16 N(ll3)
105
-------
TRAFFIC INPUT DATA—C\Ri>S f, AND 0 (Continued)
Card
1060
1060
1061
1061
1062
1062
1063
1064
1064
1065
1065
1066
1066
1067
1067
1068
1068
1069
1069
1070
1070
1071
1075
1075
1075
1076
1076
1077
1077
1078
1078
1079
1079
1080
1081
1081
1082
1082
1083
1086
1030
1087
1088
1009
1090
1091
1092
1083
1085
1084
1087
1086
1309
9
1061
1067
10b2
1066
1063
1065
1064
1065
1003
1066
10 J2
1067
lOdl
1068
1080
1079
1069
10 78
1070
1077
1071
1076
1074
10 7o
1092
1077
1091
109J
1078
1089
1079
1038
10,-),)
1081
1086
1082
1J83
10o5
1084
Iu8 7
10ci7
1088
10o9
1090
1091
1092
1093
1095
1096
1096
109 7
1097
1354
19,68
1968
1970
1970
1973
1973
19 7-«
1959
1959
1957
1957
1955
19s5
1952
1952
1951
1951
195J
1950
194d
1948
1947
1 v2 J
192f,
192o
1932
1932
193:>
1933
1934
1934
1935
1 93T
1937
1939
1939
1945
1 945
1948
1909
1937
1908
1906
1905
1903
1902
1900
1948
1914
1916
19U8
L909
1600
2083
2083
2091
2091
2 102
2102
21 1 J
2114
21 14
2 Io6
2106
2095
2uV 5
2086
2056
2U80
2 • j a 0
2 J 7 5
2075
2ob9
2069
20059
2u28
2125
2121
2121
2110
21 10
1810
130J
2300
6500
6000
6000
16000
6000
0
20500
21500
13000
21500
14500
25000
14000
29000
6000
8500
29000
8000
29UOO
10500
28JOO
2850J
9000
9500
8000
9500
37000
7000
10750
6500
11000
18000
12000
12000
14500
14000
14000
10000
10000
6000
7 JOO
6000
6000
7000
9000
13000
11500
21500
luooo
10000
6000
6030
26150
2
2
2
2
2
2
^
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
^
2
2
2
2
2
2
2
2
2
2
2
2
L.
2
2
2
2
2
2
2
2
2
2
2
4
2
7
9
15
11
15
8
15
15
9
1^;
11
12
10
16
6
17
15
6
15
6
15
5
15
32
lu
29
5
30
30
6
30
5
30
d
B
30
0
9
31
31
8
3u
7
6
6
5
15
31
14
6
6
4
5
55
N(114)
N(115)
N(116)
N(165)
N(166)
Terminator
0(1)
0(2)
106
-------
TRAFFIC INPUT DATA—CARDS N AND
-------
TRAFFIC INPUT DATA—CARDS ,N ANlJ 0 ' Com lurled
1500
2300
900
1700
2500
1100
700
2500
1300
900
2500
1500
1 100
2500
17uO
1300
27JO
500
1500
2900
700
1 700
2700
900
900
2900
1100
1100
2900
2900
1700
2100
2900
2300
2100
3100
2100
2100
3100
1900
2100
J100
1700
2100
3100
2100
2300
3100
1300
2300
3100
1500
2300
3300
2300
2300
3300
1500
7
3
1
5
6
b
1
8
7
1
7
2
3
2
4
7
1
1
2
1
3
2
1
3
1
<»
7
6
4
O ;i9
o'hiiy
a' 61
0 8(.
0 X7
108
-------
METEOROLOGICAL INPUT DATA—i Uiil.s P, Q, AND R
(Data Read-in by Subroutines RA0IUBW, SR 0BfU and SFC0BS2)
ST LOUIS 7L0326
03879710326120102 -
03879710826120203
03879710826120304
03879710826120405
03879710826120506
03879710326120607
03879710826120708
03879710826120809
03879710826120910
0387971082612101 1
03879710826121112
710826 0 0
710826 1 0
710826 2 J
710826 3 D
710826 4 0
710826 5 0
710826 6 0
710826 7 0
710826 8 0
710826 9 J
710826 10 0
710826 11 0
710826 12 0
710326 13 1
710826 14 1
710826 15 1
710826 16 1
710826 17 0
710826 18 0
710826 19 0
710826 20 0
710826 21 J
710826 22 0
710826 23 0
ST LOUIS 7103
03879710827120102
03879710827120203
0387971082712 J304
03879710827120405
03879710827120500
C3879710827120607
C3879/108271207J8
03879 710827120809
03879710827120910
038797108^712101 I
03879710827121112
03879710827121213
03879710827121il<<
710827 0 0
710827 1 0
710827 2 0
710827 3 0
710827 4 0
710827 5 0
710827 o 0
-
-
-
-
-
-
-
-
-
-
68
66
65
65
64
61
61
65
68
72
76
79
81
82
83
84
82
do
76
71
68
64
61
oO
27
-
-
-
-
-
-
-
-
-
-
-
-
-
53
5o
5H
57
55
54
53
Jt 5o 1
0158095001740
0183094002090
0215048302790
0220032003140
02190290C5940
0154031015110
009 7u310J2000
0076005027900
J055010031230
00^0004039730
-095003058020
63 00 00
62 27 04
59 26 04
uO 33 10
38 1 ! Ot
18 27 04
IB ^7 04
60 iO 05
60 33 10
19 36 05
59 32 05
59 34 05
5o 30 10
59 30 13
56 34 14
54 52 14
55 32 09
56 35 11
57 34 13
56 3& 10
56 01 09
56 J^ J3
56 01 07
r>6 16 07
76 53
01 19J9 70 Jl 740
C144096003230
0132073007750
0106084011920
0104054015100
0113030017190
00680200^6370
0075006027940
0061009031160
0050009036210
J025008041820
-037009049890
-078006038060
55 01 04
56 01 06
56 35 07
5b 3i 04
55 OD 00
54 27 04
53 30 05
09930
09390
09810
09770
09460
08500
07830
07290
07000
06300
05000
09965
0979u
09280
08830
08500
08290
07420
07280
07000
06580
06140
05550
05000
251002
264003
267003
295005
301007
310011
290014
290017
292017
302017
350002
011004
012005
347006
331007
327007
316010
313009
319008
32901t
328019
332018
335020
Card
Q(2)
Q(3)
QUO)
Q(U)
R(2)
R(3)
R(23)
R(24)
P
Q(D
Q(2)
QU2)
Q(13)
R(7)
109
-------
METEOROLOGICAL INPUT DATA—CARDS P, (),
710827 7 0
710827 d 0
710827 9 1
710327 10 7
710827 11 a
710427 12 d
710327 13 6
710327 14 6
710827 15 3
710827 16 2
710827 17 3
710827 Id 2
710327 19 7 00 00
84 55
012 7C94001740
J16S031002 730
J169U73003 730
J112077012 J5 J
0103J640l3JtO
0110J32017530
JOd502U241JO
JObbJOH028l»90
0076004031^30
J079002034060
OJI/60J304 7730
- 3310013058490
30 27 04
09986
Ov«70
u9 75j
OBdlO
08500
u3230
0/05J
07210
07000
06780
05730
05000
030002
031002
0J1003
0110 Jo
002006
010005
CO I 007
350007
359007
358008
35^010
332011
Q(12)
110
-------
(u ru •—• »•**• •»••• m ru ru ru «— ••••• •••••
ru ff> in • r- <* « •* cmom * ru -^ c in »* o o o in in ru ru •* ru •* •* 4) ru r> v
1/1 • x m o in o in ru w o * m
u. rv. < * ••*»• •••••u. ru m — 1/1 ••*«• •••••
tft ru >-4 ru ~* —• •-« ru *-i w) *^mnmfOfumminr)
c. o rri o - ,0 BO r» ^ c o »* o in in *> o * r- * r- ws o in m ^
ru rv**ftj*-~*^**-«flr *-ir **•
c o-
r*) o o * c ru ^ r- c * «* ^-tc-crctr o •» o ^> *^
ru m m ru •*•*• ••••• ru ft in r\t •*»•• »••••
u —* m ru en ru ru*^m««tj -^-<*rvi*fu»*r>j—«*iru
u. u
r'lf*)^ rn ••••• ••*•• (virgin f* ••*
• • ^ r>^u^<£ir. ~- oxomrn t/t m ru "O -jo
irnr^ * •»•••* ••••* N rv -n •* *n ••••• »•••*
^— f^^rt(*>^^ —* —- »v -^ ^^ «v» ^-» in i\* --«r\j ^— ^ ru
m*-ocr^ffin *n^Dtr"^rr ^ ^ * ^ ^ (/• fr*inc»— *v^ TX-**-* ^•x»*»nn
(*irni i*> ••••• •»••• T ru * —•* ••••• •••••
fn _j (/) *-• x — •* f\ rurv^iv ^-— irr^j^n^-r^fvirfu
O t*. X O — 3- « • — iw^-ir^ ^ ^- r- — -y. ^ j.p. ^ ^-jrvT'v ft.'rv^.'-a
*-. j — r\jf*1'\j (V' o *•••• •••••i/* rvjmm o »••••
CL »- O UJ -
5 Q 7 i r>otnox«.c r-crncm ^o-ir^ftj > T — rv c — ~ o xx^aN. r- x ^ -.n o-
^ C UJ '\;"C»'—(t—^
fc-K.^i'»( ^ Tn* 'ft *fc.if».^f\«— **^~-n*#
S1 fyr^m ^> ••«•• •••*• 2: i\' **i »n LP ••••• •••••
TO» r *-/>>X'C Tin-«ir —
in C' (y o —• in ^ nxiT^^* cxru'^'V' •-« »n o j* o — in ^- * —•* T ru •* ""^^r^in
rr>-« ^- 1*1 * TI — »r m "> •* T * T x*^« *n nnrruin in'»*>*
rxj'XjtVi —. ••*•» ••••• u •X'M'V — ••••» *»•»*
.— .\j —« ,» .x ~. fu ru ,v rv. — »i\j^^~aru >«(Vi(N.furu
•* e c o — uri y ^-r-TT"^ p^jsxct *o^n"cln^* —*cc—««T'> j> — -- —1>
rururu — ••••• ••••• fv* rv-i— — ••••• «••%»
H tl II It II H II — II It II Ii II 1^ II H II H n II || II It II II n — H *t U 11 « l\t tl II H H H
'._> o 'O ••— ^ o —" fx ^ ro c —« ru •- (M o o '.^ •- o ID *~ rv- •— fu (_r •«- *\, —* ru
^iniiii — a x -' ' rffai^tj ^ t c. i t i i aDccttTJ-J aifvcij_i
ccjioecox-ii* * i i -J t i i it nn5eto-r---C3c'O.: O^oCOt O 0_Ju.'ZXf-1~-t-OOi*C'O»^OOOC>C
111
-------
m
D
3>
f
o i
? 2
§ ?
5 ?
z «
H O
2. o n>
<° c r-
Jil
jr - - m
- £ fi O
r-
O
NO 3)
,0
0) CD
Q. —
u u.
3 °
Di
-N O
Is
O (B
10 (6
(0 3
ID g
5-D
m
O
z
<
a
Tl
*»
' »2
£ 00
U> i Tl
o
r>
------- |