3SR-844
VOLUME 111
NOVEMBER 1971
SYSTEMS, SCIENCE AND SOFTWARE
P.O. BOX 1620. LA JOLLA, CALIFORNIA 92037, TELEPHONE (714) 453-OO60
A PARTICLE-IN-CELL METHOD FOR NUMERICAL SOLUTION
OF THE ATMOSPHERIC DIFFUSION EQUATION,
AND APPLICATIONS TO AIR POLLUTION PROBLEMS
III: NEXUS/P Operations Manual
^
BY
R. C. SKLAREW, A. J. FABRICK AND J. E. PRAGER
FINAL REPORT
FOR THE
DIVISION OF METEOROLOGY
NATIONAL ENVIRONMENTAL RESEARCH CENTER
RESEARCH TRIANGLE PARK, NORTH CAROLINA 27711
UNDER
EPA CONTRACT NO. 68-02-0006
-------
100R71019
3SR-844
VOLUME III
NOVEMBER 1971
SYSTEMS, SCIENCE AND SOFTWARE
P O BOX 1620, LA JOLLA. CALIFORNIA 92037, TELEPHONE (714) 453 0060
A PARTICLE-IN-CELL METHOD FOR NUMERICAL SOLUTION
OF THE ATMOSPHERIC DIFFUSION EQUATION,
AND APPLICATIONS TO AIR POLLUTION PROBLEMS
III: NEXUS/P Operations Manual
BY
R. C. SKLAREW, A. J. FABRICK AND J. E. PRAGER
FINAL REPORT
FOR THE
DIVISION OF METEOROLOGY
NATIONAL ENVIRONMENTAL RESEARCH CENTER
RESEARCH TRIANGLE PARK, NORTH CAROLINA 27711
UNDER
EPA CONTRACT NO. 68-02-0006
-------
3SR-844
TABLE OF CONTENTS
1. INTRODUCTION
2 . SETUP
2.1
2.2
2.3
2.4
2.5
2.6
3. MAIN
3.1
3.2
3.3
3.4
3.5
3.6
3.7
3,8
3.9
3,10
3.11
3.12
INITAL
WIND
SOURCE
SOURCE Input Routines CARS, MISC, OIL and
POWER '
INVRSN
TIMSTP
PROGRAM NEXUS '
INPUT
SETUP
RESTRT/WFILE
OUTPUT "
KEM and CHEM
ROLLIN/ROLOUT
DIFFUS ' . .
NPART
PARCEL
BORDER
SOURCE
Page
1
8
8
8
9
9
10
10
11
11
12
13
15
15
16
21
21
22
23
25
26
Appendix A - Listings and Flow Charts for SETUP
Appendix B - Listings and Flow Charts for NEXUS/P
-------
3SR-844
1. INTRODUCTION
The Systems, Science and Software (S3) model of photo-
chemical air pollution consists of two computer codes; the
first, SETUP, uses meteorological and source data to create
an input tape for the second, NEXUS/P (Numerical Examination
of U_rban Smog with Photochemistry) . NEXUS/P then moves and
diffuses the pollutants, changes the pollutant concentrations
as a result of photochemical reactions, adds pollutants due
to sources, and stores or retrieves pollutants advected into
or out of the borders of the computational grid. Figure 1
indicates the major subroutines and the flow of data from
SETUP to NEXUS/P,
A flow diagram of SETUP is shown in Figure 2. The
• &
initial concentrations, maximum, minimum and nominal parcel
weights and background levels for all species are read by
subroutine INITAL, which then writes them on the NEXUS/P
input tape. Subroutine WIND reads in the hourly surface winds
as needed, and calculates the wind field for each new time
step by interpolation. Subroutine INVRSN calculates the dif-
fusivity arrays (vertical and horizontal) which vary with the
time of day, due to an assumed linearly increasing inversion
height from 6:00 A.M. to 12 noon. Subroutine SOURCE calls
several subroutines (CARS, OIL, POWER and M1SC) which calcu-
late the contributions of various types of source emissions
for each time step. Winds and diffusivitics (each multiplied
A
At the time this report was written, the Lagrangian mass
particles were called "parcels."
-------
3SR-844
SETUP (separate computer program)
initial
conditions
time step
winds and
diffusion
sources
! r
NEXUS
SETUP
initializes
pro'gram and
sets initial
conditions
L_
KEM
calculates
chemical
reactions
1.
2.
3.
4.
5.
DIFFUS
read in
.winds and
diffusion
NPART
calculates
new
parcels
PARCEL
moves
parcels
BORDER
calculates
mass flux
in borders
SOURCE
reads sources
calculate new
parcels
ROLLIN
reads in
parcel lo-
cations
and concen-
trations
P.
o
o
ROLLOUT
writes parcel
locations and
concentrations
Figure 1 — Macro Flow Diagram for Computing
Photochemical Pollution.
-------
3SR-844
C
START
IN1TAL
READ
INITIAL
CONDITIONS
WRITE
INITIAL
CONCEN-
WIND
WRITE
«t,
WINDS,
DIH'U-
SIVITJES
.AND SOURCES,
NO
READ IN
HOURLY WINDS
AN'D COMPUTE
KINDS 10R CUR-
RENT T IMF
IKVRSN
CALCULATE
DIFFUSIVITIES
SOURCE
CALCULATES
SOURCLS
PRODUCI.D
DURING CUR-
RENT TIMF STFP
TIMSTP
CALCUI ATTS
TIM1 STLP, <5t,
BASLD ON
MAXIMUM WIND
VI.LOCIT1ES
t » t + 6t
CARS
CALCULATES SOURCLS
DU1 TO BOTH FRIt-
WAY AND SURFACE
TRAFFIC
MISC
CALCULATES COURCES
DUF TO HOME HFATERS,
INCINERATORS, I TC
OIL
CALCULA1ES SOURCES
DUE TO OIL REFINERIES
POU'ER
CALCULATES SOURCLS
11UC TO POWER
STATIONS
YIS
hND
INPUT
TAPE
KITH
«t " -6t
END
Figure 2 - Flow Chart for SETUP
-------
3SR-844
by the current time step) and sources for each species are
written on the NEXUS/P input tape. A new time step is calcu-
lated by TIMSTP, and another cycle is begun. This sequence
is repeated until the elapsed simulated time is greater than
a specified final time, at which point the tape is flagged
with a negative time step.
The tape is now ready to be used in a NEXUS/P calcu-
lation. Figure 3 shows the computational flow of NEXUS/P.
First, subroutine INPUT reads cards defining print and dump
options. The initial concentrations and parcel weight speci-
fications for each pollutant species are read by the SETUP
subroutine of NEXUS/P, which then generates the initial par-
cels for the problem, and places their location and weight on
storage files. The OUTPUT subroutine prints the initial con-
ditions, and the main computational loop is entered. First,
the time step is read. Then the KEM routine is called; this
routine handles the reading and writing of concentrations and
fractional changes in concentrations for all species (the
remainder of the code has data for only one specie at a time
in core) , and calls subroutine CHEM (if there is a suffi-
cient amount of pollutant in a cell) to calculate the frac-
tional change in concentration due to photochemical reactions
Now a loop over species is entered; the parcel locations,
weights and the change in concentration arrays for each
specie in turn are read into the main frame from stor-
age files by ROLLIN. Subroutine DIFFUS reads the wind
and diffusion data supplied by SETUP, and creates a pseudo-
velocity described in the discussion of the DIFFUS sub-
routine. Subroutine NPART creates new particles when indi-
cated by a flag set in CHEM. Subroutine PARCEL moves the
parcels using the pseudo-velocity generated by DIFFUS, and
BORDER calculates the mass flux across the outer cells of the
-------
3SR-844
Figure 3 - Flow Chart for NEXUS/P
START
INPUT
SET CONSTANTS
AND PO.AD IN-
PUT CARDS
DATA DFFINING
PRI ,T AND DU IP
OPTIONS
SETUP
)1 HIT, OCLE. PARCtL
VAItlhS. P-Rl^L POSI
TIONS , UElGirS,
CONCENTRATIONS,
SOSDiRS
PFAD DATA NFCtSSARY
TO START PJ03ITM
FROM b IT 8. GENE1
AT)- P«KC LS. >HIT£
COVCLNTRAHONS AND
PARCEL POSITIONS
,
CN LNITS 21-20«i/
PARCEL POSITIONS
AND 'EIGHTS ON \
. 11-10-n '
RUD UVTA FHOM
ILVPOWN ST01-
ACE FOR 0\h
SPLCIE
CrNCrNTRATIO\S
AND F-WCTIONAL
CHANCES hO
-------
Figure 3, contcl
3SR-844
CMCUIATiS MASS
(LUX 1• 1»r
EU1LRIAS CfHS OF
1IIC CR1I) BO'liLRS
CRLATb !»th PARCELS
tut TO severs en
TO VASS nut
SORDtR Cll LS
Ct'-LUIAR
CONCEST^UIC1 1 AT
S.-tCHUD
T!' h AND
kRJTE BI"JP TMFS
PARCEL
A>D IvLlCIITS OS
u.m 11-10.n.
IND 0^
RU.\ Tl'IE
APPROACH-
ING J
(IS TIMS THE
LAbT ClLLf TO
SAVE FCR Kf-
START')
IS
SA\t >
! OD
(IOC.10)
0 T
TIME, crcin,
rCSIT10\S
hFICl'TS. CO\C-ST«
TIO\S. BORERS
COSCI »
TPMIO
A\n/ui
HMSII OIF
DUMP TAilS
-------
3SR-844
grid. Subroutine SOURCE reads the array of sources emitted
during the current time step from the input tape and adds new
parcels as needed. Subroutine OUTPUT prints concentrations,
if requested, followed by ROLOUT, which places the parcel
locations and weights and cellular concentrations on storage
files. This procedure is repeated until either a negative
time step (which flags the end of the input tape) is encount-
ered, the elapsed time exceeds the final time indicated in the
calculation, or maximum run time is approached.
A detailed description of each subroutine in both
SETUP and NEXUS/P, together with a listing of the code and a
computer generated flow chart, follow in the next section.
The computer listings are in IBM 360 compatible FORTRAN IV.
-------
3SR-844
2. SETUP
The main program of the SETUP code defines constants,
calls the major-subroutines, and does all writing on the
NEXUS/P input tape (see Figure 2). The pollutant sources
currently handled by SETUP are carbon monoxide (CO), nitric
oxide and nitrogen dioxide (NO and N0~), and hydrocarbons
(HC). In order to be compatible with NEXUS/P, which uses
five species (NO, N09, HC, 0., and HN09), the CO data are not
£3 £
written on the input tape. Two source records containing
zeroes for the 0., and HN02 are written, corresponding to the
initial concentrations of zero assumed for these species.
2.1 INITAL
Subroutine INITAL reads in the background concen-
tration (in ppm), maximum, minimum and nominal weights for
parcels (in kg) and the surface initial concentrations for
each pollutant specie. The concentrations in cells above
the surface are computed by assuming an exponential falloff
of concentration down to the background at the top level.
Parcel information and the first level concentrations are
printed for diagnostic purposes.
2.2 WIND
Subroutine WIND calculates the wind field at the cur-
rent time. Horizontal ground level winds are obtained by
-------
3SR-844
a linear interpolation of the hourly surface wind data in
meters/second; no additional data are available. Vertical
winds are computed by assuming a divergence-free wind field,
v*u = 0 . In finite difference form, the equation becomes
w - Az|'ui-l
w.
Ax Ay
where u , v , w are the winds in the x , y , z directions
respectively, and Ax , Ay and Az are the cell dimensions
in meters. The wind fields for the remaining horizontal levels
are computed in NEXUS/P using a simple wind flow model; how-
ever, if three-dimensional wind data were available, it would
be appropriate for this subroutine to calculate the entire
wind field. The wind field is edited for diagnostic purposes,
2.3 SOURCE
Subroutine SOURCE calls the source emission subrou-
tines CARS, OIL, POWER and MISC for data to calculate the
amount of pollutants emitted during the current time step.
If the mass of the pollutant in any cell is less than 0.7
times the nominal parcel weight for that specie, that mass is
stored and added to the emissions on the next cycle. This
insures that parcels created in NEXUS/P will not have too
small a mass. The total pollutant emission for the cycle as
well as the_ total cumulative emissions are printed for diag-
nostic purposes.
2.4 SOURCE INPUT ROUTINES CARS, MISC, OIL AND POWER
The data for these routines were obtained from Sys-
tems Applicatons, Inc. Report No. 71 SA1-6, "Contaminant
-------
3SR-844
Emissions in the Los Angeles Basin — Their Sources, Rates and
Distribution." We used tables from the report which provided
data on locations of stationary sources (oil refineries and
power plants), spatial and temporal distributions of freeway
and non-freeway traffic, and automobile emissions factors
(exhaust, blow-by and evaporation). Nitrogen oxide data (NO
and N07) is lumped together, so we assumed 0.95 of the NO
it A.
emissions reported is NO, and 0.05 is N02.
2.5 INVRSN
Subroutine INVRSN calculates the horizontal and
vertical diffusivities based on a linearly increasing inver-
sion height from 6:00 A.M. to 12 noon. The maximum inversion
height for each cell was estimated from the limited soundings
made in the Los Angeles basin. The minimum vertical dif-
fusivity is 0.006 m/sec, which is allowed to increase to
0.6 m/sec in each cell as the inversion layer passes
through, the horizontal diffusivities are twice those of
the vertical.
2.6 TIMSTP
Subroutine TIMSTP calculates the new time step by
allowing a parcel to be carried by the wind field only
a specified fraction of a cell in any direction (in the pres-
ent version of the code, the cell fraction FAC is set to 0.4).
The time step is calculated by finding the minimum value of
— , — , — for all cells, where u , v and w are the x ,
y and z velocities respectively, and Ax , Ay and Az are the
x , y and z dimensions of a cell. At is then set to FAC
tJmes this minimum value.
10
-------
3SR-844
3. MAIN PROGRAM NEXUS
The main program of the NEXUS/P code serves as coordi-
nator to call major subroutines in the proper sequence. A
flow chart of the main program, shoAvn in Figure 3, describes
th'e program flow and all logic associated with the main pro-
gram. In addition, all communications outside of the central
processor are indicated in this figure whether or not they are
completed in the main program. This deviation from a formal
flow chart was done in order to clarify the I/O modes used in
NEXUS/P and to provide a general overall view of the program.
3.1
BLANK COMMON
In order to save space, we have overlayed various
arrays in Blank Common that are not used for computation at
the same time. Since FORTRAN IV doesn't allow PARAMETER state-
ments, the equivalencing of arrays has to be done manually. If
the number of grid cells (NX,NY,NZ), the number of species (NSP),
or the maximum number of parcels (MAX) is changed, all DIMENSION
and EQUIVALENCE statements must be changed.
Schematically, the relevant portion of Blank Common (ARRAY)
is apportioned (not to scale);
ARRAY
C CC
DV
S
1)11 U X Y Z WT
CF CCF
11
-------
3SR-844
This is accomplished by equivalencing the first locations of
indexed quantities to locations in ARRAY as follows:
EQUIVALENCE C(l) with ARRAY(l)
CC(1) with ARRAY(NX-NY-NZ+1)
DV(1) with ARRAY(2-NX-NY'NZ+1)
S(l) with ARRAY(2-NX-NY-NZ+1)
DH(1) with ARRAY(3-NX-NY-NZ+1)
UC1) with ARRAY(4-NX-NY-NZ+1)
X(l) with ARRAY(7-NX-NY-NZ+1)
Y(l) with AR}UY(7-NX-NY-NZ+MAX-f-l)
Z(l) with ARRAY(7-NX.NY-NZ+2-MAX+l)
WT(1) with ARRAY(7-NX.NY.NZ+3-MAX+l)
In KEM, CF(1) with ARRAY(1)
CCF(l) with ARRAY(NSP-NX-NY-NZ+1)
3.2 INPUT
Subroutine INPUT sets up constants and arrays neces-
sary for both computation and output. Default values are sup-
plied for all necessary input parameters. Input to NEXUS is
through the NAMELIST function; variables supplied by input are:
ISP(5) number of pollutants species to be printed
(Default 1, 2, 3, 4, 5)
NCYC frequency of print cycles
(Default = 3)
NOL number of levels (in the Z direction) to
be printed
(Default = NZ)
12
-------
3SR-844
TDELT
CONG
PART
NPSP
RSTRT
SAVE
time interval for printing
(Default = 0)
flag for dumping concentrations on unit 2
(Default = -FALSE-)
flag for dumping particle positions and
weights on unit 28
(Default = -FALSE-)
number of species to be printed
(Default = 5)
restart flag
(Default = -FALSE-)
cycle to be dumped on unit 29 for restarting
(Default = -1)
3.3
SETUP
Subroutine SETUP initializes NEXUS/P at the start of
the problem. The following steps are performed for each
species: the maximum, nominal and minimum parcel weight and
the background concentration for each species are read into
SETUP from the data tape created by the SETUP code; the initial
concentrations (in ppm) are read into the C array; the mass of
pollutant (less the background) in each cell of the computa-
tional grid is then calculated using the following equation:
M = p(C-bk)/106-Mair/Mp-V
(1)
where
M is the pollutant weight (KG)
p is the density of air (KG/M3)
13
-------
3SR-844
C is the pollutant concentration (PPM)
bk is the background pollutant concentration (PPM)
(the background concentration rejiiains constant
throughout the grid and for all times.)
M . /M is the ratio of the molecular weight of air to
3.11 p
that of the pollutant
V is the volume of the cell (M3)
Once the pollutant weight in a cell is calculated, it is
checked to see if the weight is greater than the minimum value
defined for a parcel of that species. If the weight is less
than the minimum parcel weight, no parcels are created and
SETUP investigates the next cell. The number of parcels in
each cell and the mass of each parcel are computed by
n = M/W + 1
on
(2)
VL = M/n
t o
where
n = the number of parcels to be created
M = mass of pollutants in the cell
W = nominal mass of each parcel
W. = parcel mass.
The parcel locations are created with the aid of
URAND, a random number generator (over the interval from 0 to
1), to be distributed uniformly throughout the cell.
14
-------
3SR-844
The subroutine next fills the Eulerian border cells
(which are three times as long as the cells in the main grid)
by assuming an exponential fall-off of pollutant concentration
(less background) from the border. The initial concentra-
tions, parcel location and masses are stored on drum and sub-
routine OUTPUT is called.
3.4 RESTRT/WFILE
Subroutine RESTRT is called from INPUT if the input
variable RSTRT is -TRUE'- The first record read contains the
problem time, cycle number, and arrays of minimum, maximum
and nominal parcel weights and concentration backgrounds from
the previous run for all species. Then, in a loop over species,
the parcel positions and weights and cellular concentrations
are read and stored on drum files by subroutine ROLOUT. OUTPUT
is called to print concentrations. Next, the appropriate num-
ber of records on input unit 8, as well as parcel and concen-
tration dump files if they were written, is skipped over, and
the problem is ready to proceed from where it stopped on the
previous run.
Entry point WFILE is entered from NEXUS/P if a SAVE
cycle is set in INPUT. This performs the inverse operations,
and prepares file 29 for restarting the problem.
/
3.5 OUTPUT
Subroutine OUTPUT converts the cellular mass array to
cellular concentrations (ppm); prints cellular concentrations
(in ppm for HC "and pphm for NO, N02> 03 and HNO-) at requested
intervals; dumps concentrations (ppm) on unit 2, and/or parcel
positions and weights on unit 28 if the appropriate flags are
set.
15
-------
3SR-844
OUTPUT is, of course, optional; however, if the sub-
routine is removed, it is necessary to include the conversion
from cellular mass to cellular concentration elsewhere, as
this is necessary for DIFFUS and KEM.
3.6 KEM and CHEM
Subroutine KEM reads the cellular concentrations for
all the pollutants into the CF array. The routine then loops
over the cells, determining whether the concentration of
pollutants is sufficiently high to undergo significant chemical
reactions, and, if so, transfers to the subroutine CHEM. After
the loop of cells is completed, KEM writes the concentration
array, CF, and the fractional change in concentration array CCF,
for each species on storage files.
Subroutine CHEM solves the chemical mechanism developed
by A. Q. Eschenroeder and J. R. Martinez, ' shown in Table I.
The rate constants for the two photochemical reactions (the first
and last on Table I) vary sinusoidally from 6:30 A.M. to
6:30 P.M., corresponding to the rising and setting of the sun
on September 30 in Los Angeles. Eschenroeder and Martinez spec-
ify half of the hydrocarbons emitted as unreactive, therefore
CHEM uses only half of the HC concentration in the reaction
equations.
The chemical mechanism shown in Table I is translated
into a set of first-order coupled non-linear differential
equations. For example, in a reaction
K
A + B -> C + D (5)
(1) Technical Memorandum No. 1516, Concepts and Applications
of Photochemical Smog Models, General Research Corpora-
tion, June 197TT
16
-------
3SR-844
TABLE I CRef. 1)
RATE COEFFICIENTS FOR EXPANDED MODEL OF THE
HYDROCARBON/NITRIC OXIDE MECHANISM
(Stoichiometry imbalances may occur because of lumped
parameter assumptions.)
Reaction
Model Values from Validation
hu + N02 •»• NO + 0
0 (+ 02) + M -»• 03 + M
0, + NO -> NO, (+ 0,)
0 + HC -*• 2R0
OH + HC -»• 2RO,
R02 + NO -»• N02 + 0.5 OH
PAN
OH + NO •*• HN02*
OH + NO •> UNO
HC
0.4 min"1
1.32 x 10 ppm min
40 ppm" min" ..
6100 ppm" min"
80 ppm min"
1500 ppm" min"
6 ppm min
10 ppm min
ry f\ *" -L • "" X
30 ppm mm
0.0125 ppm min
(H20 +) NO + N02 -*• 2HN02** 0.01
hu + UNO 2 -»• NO + OH
0.001 min"
* Rate constant lumps third body concentration. (M)
**Water vapor lumped into rate coefficient.
17
-------
3SR-844
Species C and D are formed at a rate
F = KAB (6)
where K is the reaction rate constant.
Species A and B are removed at a rate proportional to their
concentrations; the removal rate is usually written
R = RAA = RBB . C7)
where RA = KB RB= KA .
A 1>
The total time rate of change of a specie with con-
centration C is given by
where F and R are the same quantities defined in equations
(6) and (7), but summed over all reactions in which species
C appears.
Differential equations of this type are solved for the
following chemical species; NO, N0~, HC, 0,, and HNO- • The
other chemical species indicated in Table I are assumed to be
in equilibrium at all times, i.e.,
C - F/R
18
-------
3SR-844
The solutions of Equation (8) are approximated for
three cases as described below.
Case I
If the ratio (NO/N02)new/(NO/N02)initial is between 0.9
and 1.1 (i.e. a change of less than 10 percent), and if the
change in 0, concentration is less than 0.01 ppm, the change in
concentrations for all species is assumed small and the solution
for Equation (8) can be approximated by
C1 = C + CF - Rc)At (9)
for all species,
where
C1 = new pollutant concentration (ppm)
C = initial pollutant concentration (ppm)
At = time step (minutes)
Case II
If Case I does not apply but the quantity RAt
is less than 10"' for a given species the solution of
Equation (8) for that species is given by
C1 = C + FAt (10)
Case III
If Cases I and II do not apply the solution of
Equation (8) can be approximated by
-R At
C' « F/R + (C - F/R ) C c (11)
c c
19
-------
3SR-844
This solution is valid only if F and R are constant during
G
the time step, At. Therefore, to assure accuracy, a small
initial time step is used (10 minutes) and the remaining
time steps are calculated using a predictor-corrector method
based on the fractional change in concentration given by
At' = (allowed error) ((C+C')/(C-C'))At (12)
Time is stepped along at these small time steps until the
total elapsed time is equal to the time step used in NEXUS/P.
This procedure is repeated until the time step is completed.
Once the new concentrations are calculated, the ratio of
new to old concentration is computed by
AC = C'/C; (13)
if, however, the original concentration is zero, CHEM returns
AC = -C' . (14)
This value is a flag set for subroutine NPART to create a new
parcel in the cell where AC < 0 and the weight of the parcel
is given by
M = -p(AC+b,)/106-M_,r/M -V (15)
Xv Ct JL I jJ
(the same definitions given in Equation (1) apply).
20
-------
3SR-844
3.7 ROLLIN/ROLOUT
Subroutine ROLLIN reads concentration and change in
concentration array for one specie from unit
N is the specie number, then reads parcel positions and
weights from unit N+10, and rewinds both drums. ROLOUT
writes the concentration array on
N+10, and rewinds the drums.
3.8 DIFFUS
« «
Subroutine DIFFUS calculates the diffusion pseudo-
velocity and adds this to the real wind fields. The routine
reads the lower level horizontal velocities (IL and V,,) and
second level vertical velocities (W-) from the input tape
generated by SETUP. The rest of the velocity field is calcu
lated by
Wl - ° U2 - Ul V2 " Vl
*3 = W2 U3 = 'Ul V3 " -Vl
This procedure constructs a simple divergence-free flow field
based only on the lower level winds. If additional informa-
tion on the wind field is known, this section of DIFFUS could
be removed and SETUP (subroutine WIND) could generate the entire
wind field to be read by DIFFUS. A mathematical derivation of the
parcel technique is given in Appendix A. The finite difference
2.1
-------
3SR-844
form of the diffusion psuedo-velocity , accurate to second
order, (Equation (.8) in Appendix A) is given by
2AX C.
V
d - VCJ+1 - cj,
2AY C.
wd •
2AZ
At the top and bottom of the grid, W, becomes
Wd • .VC1 - C2)
AZ C1
(18)
W, = D (C , - C )
d zv nz-1 _ nzj
AZ C
nz
These psuedo-velocities are then added to the real velocities
and divided by the appropriate cell length to obtain a dis-
placement in meters per time step. The diffusion velocity may
be restricted in order to limit the displacement in each cell
to no greater than 0,8 of a cell per time step.
3.9 NPART
Subroutine NPART reads in the fractional change of
concentration array calculated in KEM. If the fractional
change for a cell is less than zero (no parcels in the cell),
the routine creates a parcel, with a weight given by
Equation (15), in the center of the cell. The fractional
change is then set to one.
22
-------
3SR-844
3.10
PARCEL
Subroutine PARCEL loops over the parcels of each, pol-
lutant specie. The routine multiplies the parcel weight by
the fractional change in concentration array calculated by
KEM for the cell in which the parcel is located. If the par-
cel mass is less than PMIN, the parcel is discarded. If the
parce] mass is greater than PMAX, new parcels are created ac-
cording to;
n = W /PMAX
o
(19)
where n is the number of new parcels
P
W is the weight of the original parcel
W is the weight of the new parcels
P
The original parcel is reweighted to W also. The new parcels
- P
are located randomly about the original parcel. The new par-
cels are placed at the end of the array of parcels.
The average velocity for a parcel is calculated by
"area weighting" the parcel. This procedure assumes the size
of the parcel is the size of a cell; the fractional area of
the parcel in the neighboring cells is the fraction of the
neighboring cells' velocity associated with the parcel.
23
-------
3SR-844
V/77
The procedure is carried out in three dimensions, so it
should properly be called "volume weighting."
An exception to this procedure of "volume weighting"
occurs when a parcel is located in the lowest cell. In this
case, the vertical velocity is linearly interpolated using
half of the second cell velocity at the top edge of the
lowest cell and zero at the lower edge of the cell. Since
the velocities have been multiplied by the time step, the
new parcel positions are calculated by adding the displace-
ments to the old parcel locations. The parcel is now checked
to see if it has moved into the outer half of the outer cells
of the grid; if it has, then it is removed and, if the parcel
is in the lower two levels of the grid (which has four levels) ,
the mass of the parcel is added to the outward flux of the
border cells (see BORDER). If the parcel is not removed, the
cell in which parcel is now located has the parcel mass added
to its concentration.
24
-------
3SR-844
3.11 BORDER
Subroutine BORDER calculates advection in the
Eulerian border cells. The advection scheme is a simple
first-order donor cell method which can be expressed as:
j = Ci + F+ - F'
F = C.+.,*v/3 if v into grid
F+ = C.-v/3 if v out of grid . (20)
F" = C.+,'V/3 if v out of grid
= C. ,«v/3 if v into grid
where C! is the new pollutant mass in border cell i
JL *
C- is the old pollutant mass in border cell i
F is the inward flux of pollutant
F" is the outward flux of pollutant
v is the appropriate velocity (in the x and y
direction) of the border cell on the compu-
tational grid multiplied by the time step.
25
-------
3SR-844
The factor of three divided into the Velocity is due to the
border cells being three times as long as the computational
cells. Fluxes out of the last cell of the border cells are
lost, while fluxes into the computational grid are stored in
the appropriate FLUX arrays and are used by subroutine SOURCE
to create additional parcels in the outer cells of the grid.
3.12 SOURCE
Subroutine SOURCE reads the source array from the
input tape created by the SETUP code. This array consists
of the mass of pollutants emitted in the lower two levels of
the grid during the time step At. The mass in any partic-
ular cell is divided by the nominal parcel weight to deter-
mine the number of new parcels to be created in that cell.
The parcels are then created using URAND in the same manner
as discussed previously for subroutine SETUP. The parcel mass
is then added to the cellular concentration array.
The SOURCE routine next checks the inward fluxes
from the border cells. If necessary, new parcels are created
and put into the two lower levels of the outer cells in the
grid, and their weight is added to the cellular concentration
array. If the parcel mass is less than PMIN, the parcel is
removed and the mass is added back into the border cell array
and removed from the cellular concentration array.
26
-------
3SR-844
Appendix A
LISTINGS § FLOW CHARTS FOR SETUP
-------
?FC3«FLO.FOPFLO,S SETUP/SETUP
FLOrtC"APTED BY FoRFLO /X80B/ ON 13 OCT 71 AT 0?!53:36
COOOOI. C •
cocao:, c SETUP GENERATATES THE INPUT TAPE FOR NEXUS/P
COC-03. C
CC'O^f. C PAPA".ETER NX-25, NY.25. NZ»1, NSP-H
CCCCC5. C COr"CN OX, ov, OZ, PflW(KSP). PMAXf'.'SP), (TWINS'*), BAC<(NSP)»
CC:CC4» C 1 CtNZ.NT>NX,NSP), DVlNZiNY.NX), OH(N7,NY,NX), OT, WINX.NY),
C37oc7t c 2 UICJX.NY?, U2(NX,NY), VIINX.NYJ, V2(NX,NY>, T» TLi TN, FNO, FN02,
C C C C 3 8 « C 3 OTHi NAME(NSP)
cc:c3?. c OIHENSIO.' S(VZ.MT,NX,NSPJ
ccacio. c
ClCCll. COU.^0,^ OX, OY, OZ, PMiN(H), PMAX(H), HNtM), BACK(I),
OC3C'2. 1 CM,25.25,1) , OV(",25,25), OH(M, 25,25), OT , V;t2S,2S),
CC3313. 2 JI(23,2SJ, 02(25,251, Vl(25i25), V2(2S,25|, T. TL,
COCQIt. 3 TN, FNO, Ft«02. OTH, NAMEl")
CCCC15, DIvrNSlO-1 S(U,25,25,H)
CcaclA. EI'JIVAUEMCE (S.C)
Oo:D17« DIVISION XXX(1 .21 ,221
CC3313. DATA OX, OY, OZ, OT, Ti FINAL / 2»3218.6711, 100., 3QO., 21400.,
CICO!?. 125150./
coc^2o. DATA NAME / 'co«, 'NOI, »N02', »HC' /
CCOC2'.. NX-25 . • »
Co:o22. NY=25
or:s26. zz«o»o
c::o27. icvc - o
CC:"23- 00 5 I«l,22
;D:C:?« oo s JM ,21
CCCC33. 00 5 •;»! ,NZ
CCC;3I. 5 XXX(t,J,I)«9«0
co:c32. CALL INITAL
CO:C33. W'lTE (N) PMlNt2),PMAX(2),«N(2),BACK(2> . w
CCCC3--. KFITE (V) t ( (C(K.J,I,2) ,K = 1 ,NZj ,J«H,21) ,I»3,2S) W
CjC335. *RITE (N) PMJNO) ,PMAX(3) ,WN(3) ,8ACM3) JO
cc;c:4- WRITE INI 11 ICIK.J. i .3) .•c-i .NZI ,j«H,2M> ,i»3,2t) i
CCCC37. WRITE (N) PM!N(M) ,PMAX(1) ,nN(H) ,BACK(1) C»
:;:c33. *RITE INI t c 2. *?ITE(MXXX
C;jC"3. JO 1CYC • !CYC + I
c;c:-u- CALL '.re
PPflT 4013, ICYC, OT, T
6310 FCPt-AT I'lCYCLE ', [1, JCX, 'TIME STEP " ', F10.3, • SECONDS, TlK£
j - •, F12.3)
-------
ODOC17. CALL SOURCE
C033SO. WRITE (N) OT
C"352. 1 ( (IV( ; ,J) , 1-3,21) ,»;»•»,21? i
003C33. 2 ItIDV!K,J,1),K-t,1),J»1,21>,1-3,21) , (((OH(XtJt!>iK»l.4 It
C3C3S1. 3 J*1,21>,1-3,21)
ccross. WRITE IN) ! (ISIK.J.I.2',K»I,2)•j-i.zi),i-3,zi)
C3C054. I-RITE If.) I IUltI>J),I«3,21>,J-1,21), ( ( V 1 ( I , J) 11-3 ,21) , J-1, 21) «
3C:037. I t,j»i,21),1-3,21), < I IOHIK.J,I),*•!,i)•
C'CO^1?. 3 J = u > 21) i 1 "3 ,2D
CC:c6a. WPIJE IN) ! I (SIIC, J, I ,3) .K-I ,2) i J'1,21) , 1-3,21)
CCC04I. WRITE IN) ( I I I3CC. J) I ,1> ,K-1 ,2) , J"1 ,21) .1-3,21)
CG:044. *R(TE (N) 1 (UII I,J) , 1-3,21) ,J=1,2«) , tI V1( I,J) ,I«3,21),J«1,21),
CCCC67. 1 CCut I.J),1-3,21) ,Jof,21), ."
C3CC63. 2 (< (OV(K,J,II ,K"t,1),J«H,21),1-3,21), (( I OH(K,J,t).K-I,1),
CC034'. 3 J=H.21).I"3,21)
COC370- ''RITE (N) ( ( UXX 1 K , J , I I ,K-I ,2 ) , J" I . 2 I ) »I - 1 , 22 I
02:371. WRITE (N) ItU|(I,J),1.3,21),J-"|21), ((V|(I,J)i1-3,21),J-1,21)i
CICC72. 1 (!u\I.J),1=3,21),J-1,21),
OC:-73. 2 ((IOV(K,J,I) ,K«I ,1),J«1,Z1),1-3,Z1), (((OH I X,J,I),K«I.11.
C"07". 3 J = u>21> , I =3,21)
CC:C75. VyRITE (N) ( I ( XXX I K , J , I ) ,K- t , 2 ) , J- I , 2 t ) , I -1 ,22 )
o:-::74. CALL TIMSTP
CCTC77. IF (DT .LT, 120.) RETURN 0
CDC073. T = T » OT
OCC079. Ir (T «LE. FINAL) GO TO JO
cccc11:, OT <= -OT . .
COCCS!. ARITE (N) DT
cc:cs2. STOP
C3CG83. END
CO
73
I
CO
-p..
-------
SETUP
FLOWCHARTED BY FORFLO /xeoe/ ON is OCT 71 AT o?:sj;39
BEGIN
1 C SETUP GENERATATES THE INPUT
CTAPE FOR NEXUS/P
—c PARAMETER NX-ZS, NY-2S, NZ-S,
3
CNSp»1 3
C COMMON ox, PY, DZ, PHIN(NSP). 3
CP^AXIVSP), WN(NSP), OACK(NSP), 3
C ! C 3
COMMON OXf DY, DZ, PMJNM), PHAXM)> riN!
1). BACKCtl, CO. 25, 25, 'f ) , OV(1,2S,25),
OHt">25,25!. 07, Y;t2S,25), Ul(25,25), U2(
25, 2S), VJ(25,25), V2(25,25), T, TL, TN,
FK'O, FN02, DTH, NA^EIH)
DIMENSION SCf^S^S,1!}
EQUIVALENCE (S,C)
3ATA DX, DY, DZ, oT, T, FKIAL / 2«32I8,
69M1, 100., 3CC., 21400. ,25150. /
DATA NAME / 'CO', 'NO', 'N02', 'HC« /
NZ-1
,vSF = f
N - I]
zz*o.o
ICYC " 0
A
A
A
A
A
A
A
A
A
A
A
3
3
3
3
0
8
0
C
^
C
00 5 1=1,22
DO 5 J-l,21
I
I
00 5 S=J.NZ
00
C 53
-------
. 1 •'>•
: : CALL INITAL : :
i
i
t WRITE (N) PMlNlZ) ,PKAX(Z) ,WN(2) ,BACX(2/
I
I
/ "RITE (N) ( ( 1C (K • Ji I .Z) •*•! iNZ) t J-t (21/
I
I
/ \IRITE (N) PMINI3) ,P1AX?3> i«N(3) .BACKI3/
I
I
/ WRITE (N) U(C(K,J,I ,3> rK»| tNZ) tj»1,zo/
I i
I
/ "RITE (H) PHIN(I) (PMAXC<) tHNII) iBACKf/
I
I
/ WPITE IN) 1 ((C(K,J,1 tM> tK-1 tNZ> tJ'1,2H/
I
I
/ WRITE PKtN(3> >PKAX(3) ,nN'<3) >BACK(3/
I
I
/ WRITE (N) XXX /
• f
I
/ WRITE (N) PMINI3) tPMAX(3> i*N(3) tBACKO/
I
I
/ V.RITE (N1 XXX /
I
C 103 I
: ICYC * icrc * i :
,
•
i *
,
.
z
1
I
z
z
I
oo
-------
: : . CA.LL WIND : : '
I
i
: : CALL INVKSN : :
1
I
,,,,,,,, , ,,,,,, ,,,,,,,,,,,,,,,, ,,,,,,,,,, ,/
/ PRINT 6010! tcrc, OT, T /
/ , /
I C 4010 FOSMAT (MCYCLE '• 111 10X3
I C 'TIME STEP » ', FJ0.3. • 3
i CSECONDS, TIME i « •• FIZ.JI 3
I
5
: : CALL SOURCE : :
t
i
/,i,t,,,,,,,,,,,.,,,,,,, ,.i ,,,,,,,•,,,,,,,/
/ WRITE tN) DT /
/ , , '. ,.,..../
I
I
/... t t ., t .,,.,.,, T •.' 1 '.• f !•« t .... I . t ....!/
KSjTE ''UltI-3>21) ,J»1,71) , (( /
VI(I,J>,I"3,2M) ,J"H,2M) , «l»U,J),t"3, /
.J,tJ,K*1,1),J»1, /
iD»K»l .1> r v"1. /
) /
, ,.,/
I
I
,,..,,,,,,,,,,,,,.,,,,,.,,,,,,,, t/
(N) ( ( (S(K,J,I ,2) .K-I ,2> ,J-1 ,21I/
/ 1-3,21) /
. , /
I
I
/..t.».,...,f,,i.,,t•,t.,,«.t,,,,,,,, t ,,, v
WRITE (N) ( (UK I ,J) ,I«3.2H) , ,)•"(,2 , I»3, /
2t) ,J=",24) , l«DV, <(«(I,J),i"3, /
211 ,J»",2M) , ( ( ICVIK, J, I ) ,l"M ,") ,J«a, /
2" ) > I"3, z«M , ( ( (OH 11 , Ji I ) <\" 1 ,1 ) ? J"M. /
^
Z
Z
Z
Z
Z
Z
Z
Z
I
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
•>z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
CT
CO
-------
J
I
/ WP:TE 'N> u J)>l-3, /
21) ,j»H ,2") , I JIDV(K,J,I ) iK«l ,1) ,J-H, /
21) , 'O,2H) , < ( (DHCCiJt I) tK»l >1) i J»M, /
2") , JB3,2") /
J
1
/ WRITE (N) < ( (XXXtK.J, I) ,r.»l ,21 , J-l f21 )/
f : = [ ,27i ' /•
I
I
WRITE (N> (IUl(IiJ>2<))t t( /
VI ( I , J) , i-.3,2l) ,J'M,2«I , ( Ut I ,J) ,!»3, /
2") , J»3i2'O i tl'"! ,•») » J-Ht /
2tJ . I»3,2f) /
;
i
/• WRITE: /
/ I«I,22) /
:
i
: : CALL TIMSTP : :
j
i
! TRUE t
I
: RETURN o :
i
: T • T + OT
I
I
1 FALSE
I
z
z
z
z
z
z
z
I
z
z
z
z
2
z
z
t.
z
z
2
Z
z
z
z
z
z
2
Z
z
z
z
z
%z
z
j
Z
z
z
z
2
2
Z
z
z
z
z
z
z
z
z
z
-Pi
-------
OT « -OT
I
I
WRITE CM or
j
/»tt«ttti,i,t.iii.t,N
: STOP :
\..... /
I I
en
00
-------
FLCnCHARTEO BT Fo"FI-0 /*808/ ON 13 OCT 7t AT 07153:18
CCCC01. S'JSICUTINE INJTAL
coc:c2- cc:"^ ox, or, oz, PMJNII), PMAX 1
GCC3'2. 00 30 N n t, NSP
COCO.3- PCAD 6COO, °MAX(N)t PMIN(N>» WN(N), BACK(N)
^> COC;!J« 6CGOF-f"1"T(6£lZ«6)
OO CCS015. 00 13 L = I, NT
CCCC'6. J--T-1-L
c-::i7. ^^A.5 so;o, (Ci I ,J, I ,N> > I«I,NX>
CCColS. 5010 FC'^AT (25F2.0)
DOCCl'- DO 10 I » 1, NX
CC3323. IF (,S ,£3. 2 ,OR, N ,E9, 3) C(l,J,t»N> " . 01 • C t ! »J, 1 . N)
CC~32!« C(2,J,I,N) » (C(1,J,I,N) - 8ACK(N))».S «
SZ'tZZ- C!1,J,!,S') » tCt:>J,I,N) - BACMN>>»«25
c:c:^3. ic C(«,J,I,N> » e*ct
-------
SNITAL
FLOWCHARTEO 8Y FORFtO /X808/ ON 13 OCT 71 AT 0?:S3:H9
a
A 3
A 5
A 3
A s
A 3
A 0
A3
f 3
A e
A B
A S
A 3
A 3
A 3
A 3
A 3 c
A 3 c
A P C
A 3 C
s 3 c
A 3 C
A 9 C
SUBROUTINE INITAL
COMMON OX, OT, OZ, P«INM)» P1AXI1), <1NI
t), OACICCH, CtH,25,25,1) , DV(1,Z5,25),
CKt",25,25), OT, ,v<25,25), U1(25,Z5), UZ!
23,251, VI(2S,25), 72(25,25), T, TL, TN,
F,%0, FN02, OTH. NAfE(M)
OIXESSIOW Sit,25,25,M)
EQUIVALENCE (S.C)
NX^ZS
NZ't
N S P = t
-?:
oo 30 N
READ 6000, PMAXCM, P«IN(N), T/NtN), /
8ACK(N) /
...., ..... . ................. ,.,......./
I
I --- C 6000 FORMATtiElZriJ
I
I
_>•
00 1C L » I, NT
J=NY*I-U
RCAO soio,
r-i,NX>
,,,/
/•
./
I
I --- C 5010 FORMAT tZSFZ.O)
I
1
CO
f*~*
I
c»
•p^-
->:
oo 10 I » it NX
FALSE
IF IN
.£0, 3)
)
-------
,
A 3 C t / «
A 3 C :cU»J»l,N>"«OI»C(l»JiI,N> t ,'/ I
A c c :............,....V,.... f...: I
A 3 C i I
ARC 0< — ....................
A S C I
A 3 C • • ••
HOC Cf2,J,I,N> ' (C(tiJilfN) - 3ACK|N))».S » '
A c c SACriv;
A 3 C CI3,J,I,N) « (CU.J.IiN) - BACIUN»«.2S *
ABC BACK(N)
ABC « « '
A 3 C I
ADC c IC3 i
ABC .,.,.....,... ,
A . ... C(t,J,I,N) » BACK(N) :
A : , :
A !
A I
A /.'.» t t . t f t . t ,.,,,!.., t ,.. t , t ,,.,. I,,!,, t ./ 1
A / PRINT 4010. NAME(N) / ||
A /.. /
A I
A I C 6010 FORMAT (M1NITIAV. 3
A I CONCENTRATION OF 't A6) 3
A I
A t
A /,,,.,,,.,,,,,,,,,,,,,,,,,,,,,,.,,,,,,,,,,/
A / PRINT 6011, PHAX(N), PKIN(N), rtN(N), /
A / BACK(N) /
A /.,... , ,/
A
f, ---£ 601: FORMAT I // IX, 'PHAX', 8X3
A . C 'PMIV, !CX, '«N«, 8X, 'BACK* /3
A CHF12.3 I /// ) 3
A
A
A ...,«.«....,...,...,...,,...............,,...
A 8 • ->; oo 20 u « i. NT :
A 3 :.. :
A 8 t
A 5 I
A 3 ,
A 3 : J«NY»I-U ;
A 3 :. :
A 3 I
A3 C 203 I
A 3 /t»t»it«t,ttitftt,tn ,tt?,«,«.,»! titt •»»«!/
A . . ._ / PRIVT 602C, tCM.JiJt"), I-IfNX) /
A / /
A I
A I—-C 6020 FORMAT (2SF5.O 3
A I
A C 303 I
A ...«
—..—..—.. —; CONTINUE :
i
!
/'"••••"" """"N
: RETI/PN :
\ /
en
JO
I
c»
-------
.S WIND
BT F0P,FLO /X803/ ON 13 OCT 71 AT 0?:53:55
CCOC01. SU3ROjT;N£ YIJND
OC3CC2' CC"1CN OX, OT, OZ, PMIN(H), PMAXCDi YrN(f>» 3ACK(t)t
CCC:03. 1 Cf1,25,25,1), OV(1,2S,2S>, OH 11125t25)f DT, W(25,25j,
CC;r:i« 2 Ul (25i2!>) t U2(25>2S), Vl(25,25), V2(25,25), T» TLt
CaCCC'S. 3 T'., F'-'C. FS02, OTH, NA^EIM)
Ccr:C4. DICE'S.0^ S<<",2S,25,H )
oo:?c7. EQUIVALENCE (S.c)
0:333". M*=2S
CCCC'Ot N Y = 2 3
credit NZ = "
COCCI2. !.S? = t
cocc:--r
C C C C 2 C t
CCCC2! .
OC"22.
CZ'Z'I .
CCC03D.
003335.
CDCC36.
C07C37.
I-
I sv
! i
B
READ
I
400!
&OCC
2
10
3
REA
00
00
U 1 (
v : ;
PEA
=?CA
D
1
1
! ,
1 ,
I
0
S>< ,
1
6001
6000
J »
I -
J ) »
J) B
6201
6COO
VE. OJ GO TO 10
t
t
J
I
p
L-
1
r
I.
X,Ul ,V1
, NY
t N J(
UK i «j)»ci
vi i i i j>»ci
Q
L
X ,U2, VJ
FOPVAT ( 15 )
FCT
00
DO
U2 (
'2 (
72
nYM
NXf
IF
REA
DO
oo
U2(
V2(
VAT <6E 12 t4)
)
2
; ,
i ,
IS
J
1
(T
0
3
3
I >
1 >
J =
! =
J) •
J) =
U*36
<• NY
« NX
• LT
0}
J "
I «
J) •>
Jl B
i
i
u
,
L
1
1
f NY
> VX
U2( I , J)»C1
V2( I , J)»C1
0
- 1
- t
T2> GO TO 20
> U2, V2
• NY
, NX
U2I 1 ,J)'C1
V2(I,JI«CI
CX5
20 IF (T -UT. Tl) GO TO 50
CO 30 J * 1 , Nt
DO 30 I « 1 , NX
UMI..I) * - Ul{1*1,J)>/OX » (VldtJ-1) - V1(J,J»1>>
I /f>r)».5
-------
OOC052.
END
tvJ
oo
-------
FLOYICHARTEO at FORFLO /xsos/ ON ta OCT ?» AT 07:53:57
A
A 9
A 3
A 3
A 3
A 3
A 3
. . . ...
N SSJBROuTIvc; r/lNO /
I
I
CO"MON OX, OY, DZ, PMINCDf PPAXCO, WN (
"), 3ACKCO, C(1, 25,25,1) , OVCf.25,25),
OH(M,25,2S), OT, \, (2S,2S>, Ul(25,25), U2(
25i25,, V|(25,25), V2(25,25!, T, TU , TN ,
"NO, FN'02, DTw, r.'AfEC")
EOUI VAl.£f«CE 1 S .C>
OATA C'/ .1H70H 1 /
NZ't
I
I
I FALSE
:
i
: is" • i ;
* • * '
i
i
/ PEAo 4001 , L f
I
/ READ 6000, X,U1,V1 /
I
1
•»----• •>" 00 1 J « 1 , NT :
i
I
...... ->• no i lot. NX l
I
1
UK i , j> • ui c 1 1 J) *ci :
!
C 13 I
Z
Z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
c
z
z
z
OO
-------
A
A
A
A 3
A 3
A 8
A 0
A 3
A 3
A
A
I /'
: TI - u«340o :
j
i
/ READ *oot, L /
i
i
/ READ 4000. x,u2,vz /
i
j C 4001 FORMAT(IS) 3
1 C 40CO FORHAT(6E:i2»4) 3
1
1
I
I
I
: wti.j, - uzn.jj.cj
i
C ?3 1
;
T2 = L'3600
Nt Ml " NY - 1
NXM1 » NX - 1
I
C 103 I
I FALSE
I
/ READ 19) Lt U2. V2 /
I
I
I
I
.
Z
Z
z
z
z
z
z
z
z
z
z
z
2
z
z
z
z
I
. 2
z
z
2
z
z
z
z
z
z
z
z
z
z
z
.
•
z
z
z
z
z
z
z
z
z
z
z
7
CO
?a
i
oo
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
K
A
A
A
3
3
3
9
8
rj
B
"
B
B
g
8
C
g
B
3
1 //
I
: izn ij> « uzu ,j)»ci :
i
t 33 I
1
I
I 72 * L • 36QO I
1
C 2C3 I
1 FALSE
I
I
I
I
I
1
: ui = {U2(I,J)«OT •» ui
-------
A3 .1
A3 C 103 I
A 3 •' •
—; w!l,J> = DZ*( (Ul U-l > J> - UM 1*1 , J) )/DX * I
: (Vic ,J-j) - viu ,j*i» /DYJ«.S :
I
I
/t 11111 f 1111111 > t, t
-------
SOURCL/SEIUP
TEO si FO'FLO /xsos/ ON- 13 OCT 7j AT o?:5f;o3
OC3:22.
c:?c:;,
C~:o24.
CC7C27.
CC"0?S.
Cr:c2'-
c;.c;c.
C2C03I.
CC"032.
C3C3J3t
NE SOURCE
OCJ02. CO^-CN ox, or, oz, PMJNCO, PKAXII), WN,
OCQC3. 1 Ct M, 25, 23,1) , DV ( M , 25 t 25 ) , OH < <) , 25 , 25 > , OT, W(2S,2S>,
ICaC". 2 Ut!25,2S|, 02(25,25), VHZSiZSJt V2(25,25), T» TL i
CCCCCJ- 3 T"-. r"0i 'N02, OT", NAKE(M)
CT?C«.. OI«C'.'SlC« S(M,J5,2S,M)
C"007. Ei'JiV«LEMCE 13, C)
::c~r. c or^sicti TOTINS"), CT(.MSP). SC(NSP), STOREtsZtMir
Cc:c,:°. OI"E»3'0'l TOTf), CTO), NP(M), S TORE 1 1 , 25 , Z5 , H )
cccoiO. DATA N»MC / 'co'i 'NO', 'N02', «HC' /
c c o c . : , N r, = 2 s
OCICN»I,NSP
00 10 I = 1, NX
oo 10 J = i , NY
ooio.'=-i,\z
si>;,j,:,N)oo«
OTH » or/3600.
C»LL CAPS
CALL -ISC
C..LL OIL
CAL', PO/.Ert
00 3C N =
C T ( N ) » 0 .
30 30 I =»
30 3C J =
00 30 X •
P « SCC, J,
IF (P/i(NlN
STCRCfK.J,
1 i NSP
1 , NX
i , NT
It 2
It"? * STOR£flC,Jt I ,,V'
) -GE« «7) GO TO 20
I ,N| « P
S(K , J, I ,N)cQ.O
GO TO 30
20 S|K, J, I ,N)
C1(") " CT
• P
(N) * P
p
STCPEIX , j, I ,N] • o.
30 CON'T.N'JE
F=!!NT 6030
33 F00",,- CC', 10X, 'NO. PARCELS', 2IIOX, 'POLLUTANT EMITTED1) /
i i3x, 'CPEATEO', isx, 'THIS CYCLE*. iax, 'THUS FAR* // »
DO 10 H * i , NSP
N'Ui " CT»M)/WN(M) • .5
"0 "PI'^T 60")C, NA^E(V), NPtNl, CT(N), TOTIN)
cto FC^^-AT IA?, no, iex, I=>EIC.H, i?x, EIC.M)
-------
FLOflCHARTED BY FORFLO /X808/1 ON J3 OCT 71 AT 09:S1!Oq
oo
\ SUOROUTlNE SOURCE /
\ • « /
I
1
COMMON OX, 0", DZ, PKIMH), PPAXMJ, WNl
•)), OACKCMIi CO,25,25.•)) , 0 V (1, 25 , 75 ) ,
C"Ct,25,25) » Of, «(2S,2S), 1)1(25,25), U2(
2SiZ5>, V1(2S,25), V2IZ5.25), T, TU, TH,
FNO, FN02, OTH, .SAKE ("I
DIMENSION S<<4 ,25,25, H)
EQUIVALENCE (S,C)
I
i---c OIKENSION TOTINS?), CTINSP>, NPJ
I CNSP), S70REtNZ,NY,NX,NSP) 3
I
I
DIMENSION TOTO), CTItli NP < H ) , STOREOi
25i25,'i)
A,"E / 'co», 'NO', «NOr», IHC' /
NX»25
NY.25
NZ»"
NSP»1
I
I
A._ . _.>; 00 10 N o 1, K'SP
A :,, , ,
A 1
A I
A • ,,, t ,,.,.....,,,..
* B—————_">; 00 10 I » 1, N*
A 3 :........
A e :
A 3 I
A S
A D C >: 00 10 J ' 1, NT
A 3 c :,
A 3 C I
A 3 C I
A 3 C
A 3 c c———-—->: DO :o K • u NZ
A 3 c c : ,
A 3 C 0 I
A B C 0 C 1C3 I
A E C 0
—: stK,j,i,Nj r, o-
i
i
OTH - OT/3/.00*
I
• <~ ,1 i t C' P p T
C/)
50
I
oo
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
8
e
a
3
g
s
0
8
3
0
p
3
I '
I /x
: : CALL CARS : '•
i
i
: : CALL «isc : :
i
i
: : CALL OIL : :
i
i
: : CALL POKER ; :
I ;,
I
I
I
: CT IN) • o. :
i
i
i
i
C I
c i
C 5 I
C D I
CO : P « S(ICfJtI,K) * STORTlK.J, I ,N> :
C 0 I
C 0 I
C 3 i FALSE
C 0 !
C 3 1
r- 2 : STORE
-------
to
o
A S C 0 N / •/ X Z
A B C D / X Z
A S C D • " X Z
A 5 C 0 o<-«--~- —- — - — -•-"•"••-«- — — ••-•«•«•••• — •••••-•«•••— — ••'••••••••-• —•- — .0
A 3 C 0 C 203 I X
ABC; , t x
A S C 3 5CIC, J, I ,N! • P X
A B C 0 CTtN) » CTIN) + P X
A 3 C 0 . TOT(N) » TOTtN) » P X
A S C D '. STOREiK, J, I ,N) • 0. ! X
A s c o •••• .t..: x
A 3 C 0 I X
A 3 C 0 o<----------------------------•-----•--•-••-----•----• — ••A.-T —.. — --• —• — ---. —0
A 8 C D C 303 I
J 3 C :> •
—• CONTlf>.UE
I
I
/tf ... tt ,, f(t ,,,,,...,t.,,tii,.,t t•.t,.f1,/
/ PRINT 6020 /
/,... / "
I
j c 6030 FOPMAT ('0', 10X. 'NOt 3
I CPARCELS', ZfiOX, 'POLLUTANT 3
I CEKITTEO'J / I I3X? 'CREATED'/ 3
I C15X, 'THIS CYCLE'i I8X, 'THUS 3
I CfA1!' // ) 3
I
I
->:
DO «!0 N
1 , N'SP
NPtN) « CT(N)/WN(N) » .5
"03
N)
6010, NACE(N)i NPtN), CT(N), TOT/
/
I
I --- C 6010 FORMAT (AT, [10, I8X,
I Cl?E10.tf 17X, ElOtH)
I
I
CO
JO
I
oo
-pi
-------
SFCR'FUO'FCRFLO.S CARS
BY Ft 04t1i25t25), OT, VJ125.25),
2 Ul(25.:5), U2I25.251, v!(2S,2S), V2(2S,25), T» TLi
3 TAi FNO, F\'02i OTH, NAHE(H)
Ox^f.TlON 3(«,25,2S,S)
E301/A'-ENCc (3,C)
Ol-r^stO.1 Fy tr,x ,,'Y) , S«(NXiWT) » TFft21)» TFS(21) ,
Oi-E'3!0;i F!M23,23>, 5^(25,25)1 TFF(2M), TFS(2t)
DATA SHC. 3N07, OCO, 5MCB3, 5HCEV, FNO, FN02 / 7.7t H.t 51.2. |.J,
1 ,3>-3, -7, . ! /
OATA T-e t 5«,0078, ,0173, .05"1, .0768, ,0418. .0536, 1«.0"»?S,
1 .0547. 3».07"6, ,0578. 5»,0302 /
OA'A TFS / A».OC68, .QZ73, 2«-.065l, 2».OS02, S».060?t Z>.0819,
1 2«.03tl, 1..C308 /
cc:; '
c o '.. o r
COC012.
IF ('F «NE» 0) GO To 30
IF =• 1
00 10 J • 1, Nt
10 fi£AD 5010, (FM(I,J),J.I,N*>
5010 Fo°HAT (25F3.0)
00 23 J= I , NY
20 PEAO 50;s, (s«(i,j),1-1,NXJ
33 73 = T
T!"£ » TS/34CO"
IT . TIKE • 1.
EVAP = 0.
i*- (TI.UE .GE. 7, .AND. TIKE .LE» i7.) EVAP • QHCEV
TT « 36oo«i IT - n
IF (TS - OT .(.T. TT) SO TO 10
Tf * OTH'TFFIIT)
FS •> 3"P«TFS( jT)
GO TO So
"0 FF • ((TT - TS * OT)»TFF * 9N'0!(»FNO«F
6C S(I,J,I,1| » Stl.J.t.H) » (QHC « QHC33)»F * £VAP»FS»SM«I,J)
Co
•;:-'• * •
c 3:317 1
-------
CARS
9Y FORFLO /X808/ ON 13 OCT 71 AT
SUBROUTINE CARS
COMMON OX,
OZ, PfIN(1|» PMAXC1I, (VN(
D"!M, 25,75), OT, 'A(25,2S), J1(2S,Z5), U2(
25,25), v:<2S,2SJ, V2(2S,;5;, r, TL , TN,
Ft,0, TN02, DTH, NAfEtH)
0 I M 2 N S j 0 N S(H125,25i't)
CE (S,C>
I
t—c DIMENSION rn
-------
I /' Z
, / 2
A—. . ->: oo'20 j»i, NY ; z
A : : • 2
A I 2
A C 203 I 2
A ,,.,.,.,,,.,. ,,,,,,i,.,,,,t.t,,,,,,/ Z
—.- ....... / Ri:AD 5010, {SMtl, JJ.I.l ,NXJ / Z
<- t 2
t 2
0<---------------«----«-«-«---«-------------------------------- — --.-«----«.--..---.-,o
C 303 I
TS = T
TIME - TS/3AOO*
IT » T I.vS + 1 •
EVAP » 0.
I
x,,t,.,t,,,.,,,,,,,,,,,,,,,,,,,,, i,,,,,,,.N PAUSE;
< IF (TIME .GE« 7. .ANO. TIME ,LE. 17.) > •
\..«. ,./
i TRUE
I
; EVAP " QHCEV ;
i
o< . ............
I
: TT » 3600MJT - 1! - :
I " •
/tttntt,...,.,,, , r. ,t,.tt,iitt<•»»«>•»t»i\ TRUE »
< IF (TS - OT .LT. TT) 60 TO 1Q >- — -----....,---. — .----- — - — ---- — ------—....... .. — -. 0
\.» .........,../ Z
I FALSE ' Z
I Z
I Z
t.. Z
FF > OTH'TFF(IT) 5 Z
...t.... ........,..,: • z
t z
1 Z
/,,,,.,,.,-,,.f,.,,.N I
: GO TO so :-••— ....................—........—...............—.—..o z
\ / X Z
" X Z
X Z
o< .... .—.......—....__.,._........_.-_...._._.—.........—. .—_0
C "(03 I X
FF « ttTT - TS » OT)'TFFIlT-1) * (TS - X
TT)»TFF(IT)1/3AOO. X
FS « ((TT - TS + DTi»TFS(lT-I) » (TS - X
T:i»TFS(:T>)/3*cc. x
I * X
0< ,. . ._...._................_.-_........-.-..-.--.-..-..-...-_........0
c so: i
oo
-------
*•««•-"------ — —-«-»-x • »»v OW W - II rt'
A :.,....................................
A I
A ' I
A ,
A n-----—---.-_->: oo 60 i • 1 > <•>*•
A s :....
A 3 I
A 3 I
A S t
A B : F • FF«FM(I,J) * FS'SKUiJ)
A 3 ' S(l,J,I.l) o StliJiIil> » QCO'F
A 3 ! SCtJ.IiZ) » SIJ.J-I,?) * ONOX»rNO»P
A 3 : S(ItJi!i3) = S(l,J.Ii3) * QNOX»rN02»F
A s :............, ,
A 3 1
AS C 603 1
A S
—; sii.j,:,") " sii,j.i,i) * IQHC *
J F « £VAP«F5«S>'< I , J)
I
I
/•"••••••"•""" '\
: RETURN :
\,................../
CO
x>
I
oo
-------
B'C?!«rLC.FOpn.O,S
r LO,,C"AR7£0 ET FQPFLO /X308/ ON 13 OCT 71 AT 0?:51:11
CC3301* SUO'OJTINE INVRSN
-CCOC2. C0f1"0"l OX, 3Y, DZ,
OCCCC3. i CO,23,25,1J, OV(1,25,25), On(1,2S,2Sl. DT» W(25,25),
•D"-:1-. 2 Ul(25'2S), U2(25,25/, Vl(25»25), V2I2S.2S), T. TL.
CSCCTJ. 3 T'., FTO, FN02. DTn, NAMEO)
c:::ci. DIMENSION 30,25.25,1)
CCC737. EQUIVALENCE (S,C)
C Oi'E^S.ON H!NV(NX,\Y)| OVHfNZ,NT.VX), DHM(NZ,NY,NX)
O.vjNSlCV HI\V<23,25), OVH(1,25,25), OHH(1,25,25)
IF (17 .GT. 0) GO TO 20
17*1
^ X = 25
c;:ci3. NT = 25
CCCOIS. N3p . 1
>• cc:r' 4. oo ic L » i, NT
^ CC53.7. J - hv . , - L
^ CI-; 3. "£AD 5COC, (HJNVtl.J),J.J,25]
5000 FOPSAT125I3)
CC 1C t ° 1t NX
vi'.Vl.fiJ' » HPIV( I i J>«3. 0180061/21600.
co 10 r « i, Nz
DVHIr,J,u « ,006
10 o u H t
0" 1 >
- £ •
"••).
*" M - f
c *• s •
3-6.
OVH
0"i*
! F
OVH
DM"
30 COW
10 00
00
00
ov<
50 C-1(
c. ,0
(
;
(
(
(
T
5
S
5
Y
"
J
*\ ,
K ,
•C
<
^ t
K,
IN
L>
Q
C
1 J
• J
J ,
J,
• £
Ji
j,
VE
i
j
K
it !
L.
I )
J )
0.
1
; i
I )
3
e
a
)
'
= ,591'FRAC
a 1 , J 8B«FRAC
I I GO TO 30
• .6
= 1-2
I, NX
! , NT
i, NZ
« DVHtK.J, 1 ) •
" OHKCt.J, I )•
* .006
* .012
•
OT
OT
en
oo
-------
1NY?SN
FUOr»CHAPTEt> BY FORFLO /XB08/ ON 13 OCT 71 AT 0?:St:iS
to
\ SUBROUTINE iNVRSN /
\.. t ./
1
1
CO^KO!". DX, 0" , OZ, PMIN(I). PHAX(I), YIN<
"), B»CUH>, Clt.zS.ZS,") , OV(H,25,251 ,
CH<<<>ZS,25) , OT, ft<25,25), U1(25,25J, U2<
?S,2SJ, VK?S,25), V212S.2S), T, TL, TN,
FNO, FN07, DTH, NAKEIH)
- OI^lcTf SiON SLEMCE (S,C)
C DIMENSION KINV(NX,NYIt OVHINZ, 3
NTiW)() 3
I DIMENSION HINV(2S,25)| DVH(1125,25 I, OHH( :
: 1,25.25) :
I
I
,,.,..,,,.,,,,.,,,,,,,,, . , 11 > i . t • t i . iv TRUE
< IF CIT ,GT« 0) GO TO 70 > .-.„». —_.-..-—.—-. ... — -.-..--_..- . ... — ...5
\ / z
I FALSE Z
I Z
1 I
. z
IT = I Z
NX - 25 Z
NY • Z5 ' Z
NZ = 1 Z
NIP « " . Z
: z
i z
I Z
, ,. z
A-. .->: ao 10 L « i, NT : . z
A : •• : z
A I Z
A I Z
A Z
A : j « NY • i - L : z
A : ....: z
A 1 I
A 1 Z
A /. t,i...i,,,«,,,,,,,,,,ittt..i.tittf.tt.i,/ Z
A / READ scooi (HINVIi,J),i»i,25) / z
A / , / Z
A 1 Z
I C 5000 FORf!ATI25I3l 3 Z
* I Z
A I Z
A . Z
j p...............>; DO iG i • it NX. : z
7
CM
en
oo
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
f
fl
A
A
A
A
A
A
A
A
A
A
A
3
e
3
3 C
6 C
2 C
3 C
S C
B C
8
E
3
B
g
3
B
B
3
3
3
3
3
3
S
D
5
e
3
i
: HiNvdiJ) » HIVYU ij)o.o«iS506i/2j6oo« :
i
i
i
: OVH«K,J,I> - -CC6 :
C 1CD I
t '
C 203 I
< TO MO f >
I FALSE '
I
I
I
I
I
1
X1NV « HINV(!,J)»17 - 21600. >/OZ :
: K ' XJNV :
i
i
I FALSE
1
FRAC - XJNV - 1?
K » K * 1
DVVIKtJiP » .S?H«FRAC * .006
OHH(K,J,|) B 1.188'FPAC + .012
I
1
1 FALSE
•
.
A
X
X
X
X
X
X
X
x
X
X
z
z
z
2
z
2
Z
Z
Z
2
Z
Z
2
z
z
z
z
z
z
z
2
2
Z
Z
2
Z
Z
z
z
z
2
Z
z
z
2
Z
z
z
z
2
Z
Z
z
z
z
CO
?3
i
oo
-pi
-pa.
-------
(S)
oo
A
A
A
A
A
A
A
A
A
A
A
;,
A
A
A
8
3
3
3
S
3
3
8 C
3 C
B C
^ C
<5 C
3 C
3 C
3 C
r C
K .* K ' I '
DVMl'iC.Jil) « .4 '
DHHIK, J, I ) » I -2
I
C 303 I
I
t 103 I
1
I
1
I
I
I
: ov(Ktj,n * OVHJK.JI! I»DT :
i
c so: i
i
i
: RETURN :
x z
X Z
X Z
X Z
X 2
I
_
z
z
•
.
%
•
CD
oo
-pi
-pi
-------
FLCAC«-ARTEO BY FC^FUO /X808/ ON 13 OCT 7j AT 0?!5M:i?
02=301. SUBROUTINE KISC
C3CIC2. CGVKCN OX, Ov, OZi PHlN(«f>t PMX(t),
CC::CJ. 1 C( S,25,25,•)), DVO,2Si2S) i DHtl,25,251, OT, M25,2S|,
CCCCC1-. 2 UM23.25)r U2(2S,ZS>, Vl(2St25), V2C25,25), Ti TL,
3 IN, rno, FNOZ, OTH, NAME(HI
OjfCf.SiO^ 3(1,25,25,1)
(S.CI
A\OX (NX ,N'Y> , AHRCVX,NY>, ALR(NX»NT) !
5), AHR(25,25), ALR(?S,25)
NSP=S
>F •= I
00 2 L • 1t NY
J»NT»1-L
c::c's. SEAD 5000, UNOXIt,j>tJ»J»NX>
CCCC.'!. Z COMIVjE
cc:523. oo H u » i« NT
Cr:c22. PEA3 5COO, (AHR(JiJ)il-liNX)
cc:r2J. oo 6 L = i» NT
=£,,0 scoo, ULRC I fj> , t«i ,NXI
6 CCNTiN'JE
C;_32S. SOCO FO^"AT (2SF3.0J
C:C^2=. 1C 00 23 I « 1 , NX
CC:r J3. 00 20 J • 1 i NT
CCC-33,. 3tl,J,i.2) - SU.J.1,2) + ANOX < J . J) »OTH«FNO
OCCC32. 3<:,J,I,3> = S(IiJ,I,3) * ANOX t I , J J »DTrl«FN02
CCC:33. 20 5(1, J, I, M) a SUiJtltH) * OTH« ( AH1^ ( I , J I
CO
JO
I
c»
-Pi.
-ps.
-------
F(.OIVC"ARTE!> BT FORFLO SXS09/ ON 13 OCT 71 AT
A
A
A
A
A
A
A
A
\ SUBROUTINE MlSC /
i
COUPON ox, OT , oz, p^iNCn, PMAXtt), WNI :
H), OACK(H). C(M,25,2S,H), OV(1,25,25) ,
DM(«,25,2S> . DT, v,(2S,2j), V1(2S,2S), U2<
25.25), V[I2S,2S), V2(2S,25), T, TL, TM, •*
FNO, FN07, OTH, NA«EI1)
DIMENSIOf S (1 t2S,2StNY) J
I
1
DI."E'2S1
NX.25
I
I
I FALSE
1
T
: JF - i :
i
i
i
i
; J.NY-M-L :
i
:
/ READ SOOOt UNOXC I ,J) i I"l ,NX) /
I
C 23 I
i
%
2
Z
Z
z
2
Z
z
2
z
Z
Z
2
Z
2
Z
Z
Z
z
z
z
z
-
z
7
OJ
cn
?=>
OS
-------
A
A
A
A
A
A
M
A
/t
A
f>\
A
A
A
A
A
A
A
A
A
A
r*
^
A
A
g
S
3
o
S
3
=5
S
3
....... .7; DW i u • it n' . x' '
I
I
: j-f-r-M-1. J
i
:
/ READ sooo, IAHTI i ,j> , i«i ,N*) /
i
',. * *
1 <-0 T .
I
I
t
!
: J»NY-M-U ' :
i
i
/ R£A!J 5000, f A(.P< I , J) ,1-1 ,NX) /
I
C 63 1
1
1 c 5000 FORMAT (ZSF3.0) 3
I
C 103 I
1
t
*' **„'*. ,v »•««•
1
I
S!l,w',l,Z) = S(liJitiZ) * ANOX ( 1 . J)»OTH»
S(J,J,!,3> •• S(],J,I,3) • ANOX 1 1 ,J)«OTH«
FACZ
I
c zo: i
: A L o i : i j > j :
4
z
z
z
2
i
z
z
z
z
z
2
z
z
*•
z
z
z
z
z
z
z
z
z
z
z
•
,
OJ
CO
?3
I
oo
-------
RETURN
OO
-------
SFO=«FLO.FORFLO,S OIL
A'TEO BY Fo"fUO /X800/ ON 13 OCT 71 AT 09:51:23
CDC301' SU3POJTI.VE OIL
coccc:. co-"o\ ox, DY, oz, P«IN, PMAMD, TNCO, BACKIMI.
CCC3C3. 1 C<" ,25,25,1) , OV(1,25,25), OH<1,25,25), OT, YM25.2S),
CCC-0". 2 UH25t7S), U2(25,25», VM25.25), V2(25I25)« Ti TL,
COC3C5. 3 T'i, F'-Q, FNC2, 07H, NAKEtl)
O-CCCA. DI-iNSlC^ SIM ,25,25,M)
CCCIC7. EOU1V«LCVCE (S,CI
CCljr?. C OIUC>5!ON UOCX(NCR), ICCYINOS), AKOXINOR), AHt(NOR), AURtNOR)
CDCO:J. 3ATA LOCX / 12, 11, n, 12, 9, 11, 12, is, 7, is, 2*17, it /
CDC012. DATA LOGY / 2«8, 2»9, 3«10, 11» 3»12, 13, 11 /
CCOCiJ. OATA A.> C:C01S. DA'A AHS / 31., 2., 3., 11,, 2*.» 2«, 9., 3,, 13., 5., 10., 6.,
CM C?CT 6. i :. /
CM Co:?.7. OAT,', ALP / 333., 20., 26,, 132«, 2*0., 20., 87., 32., 1J9., 52.,
cochin. i 1C"., 59,, ic. /
ccc:;2. NJ=H
CCCC2J. SS'sl
C^C72 + OTH«(AHR(N) » ALRfN))
c;a;;a. REYU^N
oc:a2i. END
OJ
en
oo
-------
FLOBCHARTED 8* FOfFUO /X808/ ON 13 OCT 7J AT 07:51:21
O-l
\ SUBROUTINE OIL /
\.,f « /
1
• COMMON ox, or, oz, pMisciji PHAXCO, wm
: S), BACIC(t), CO,25,25,1 ), 0V(«,23,25 ) ,
; OKI",2S,25>, OT, r,(25,25), Ut(2S,25>, U21
: 25,23), VJ(25,23), V2(25,25), T, TL, TN,
: F\0, FN02, OTh, NA'"ttu>
S(4,25,25,HI
NC£ (S,CI '
1
1-—C PARAMETER NOR-13 3
1 C DIMENSION LOCXISOP), LOCY?NOR)i3
I CA'-OXfOR) t AIRING"!), ALR(NOR) 3
I
I
DIMENSION LOCXI13) ,LOCt(l3',ANOX(13) ,AHR(
13),AUP<13)
OATA LOCX / 12, 11, 11, 12, 7, 11, 12,
!S, 7, 15, 2»17, J1 /
OATA UCCr / 2«8, 2*7, 3«10> 111 3«12, J3,
11 /
D/TA AN'OX / 115., Z7., 31., 573., 311,,
27., 116., S2,, 556., 67., 135,, 77., 7.
DATA AHR / 31., 2., 3., HI.. 26., 2,, 7.,
3'» 13., 5., IO-, *•> I« t
DATA AU"? / 335., 20,, 26., H32., 260.,
20., 67., 32', flf., 52., 101., 57., 10.
NX-25
I
I
00 10 N • 1, NOR
I
1
t « LOCXtN)
J «
OTH
S(l,J,ti3) « SI],J,I,3) » FN02«ANOX(N)
3T."
C 1C3
I
oo
-------
RETURN
en
I
oo
-------
o«Ft..t., -i
FS.07C"AR7EO BY Fo»F>.0 /X803/ ON 13 OCT 71 AT a?:si:Z6
sy.acuf INE POWER
CO--ON ox, DY, or, pMjNt«o> PKAXJ»I» «NtC
CC~:~7, £ffJ!VALJ!N.'Ct (S.C)
ceres'1' c "A'A.'C'CK N-'P-.! ,'
c:u3C°. c o:cuN5ios I,OCX>.JPP), S.OCYINPPJI UOCHNPP). AMTINPP)
CCCC.3. DI-I«SiON LOCXU I ) ,LOCY| 11) .UOCZtl 1 ) f AMTI II J
ccv?1'. OAT\ u3c~ / u, 7, 3, :?, nt 1*1 4, at i5i 10, 11 /
c::o:2. DATA LOCV / 7, 12, 10, 3, o, 7, it, 21, 20. 22. 21 /
c;csi3. OATA LCCZ / o«2, 3*1 /
CCCC!". ?frA « T / ZOP7,, 990. t ZOOS.., 1225., 333., 60?. i 78,, 166.,
c::? 3. i 3; . , 21 1 . r AC. /
•<•<.-??
oo 10 N « i,
I <= LOCXlN")
J - LOCY(N)
C2C026. JO S(K,J,I,3) = SUtJil»3> * FNOZ^OTH^AMT ( N )
c::az7. RETOKN
occcrs. END
CO
S3
i
00
-------
FLOKCHARTED BY FORFLO /XSCB/ ON 13 OCT 71 AT 0?:5i:27
SUBROUTINE POWER
CC1KO«J OX, OY, OZ, PXINItlf PMAXCO, WN(
2S),
C--tu,I5,25l , DT, «t23,?S>, Ul(25,25), U2 (
25,25), vi(2S,25>, vZ{ZS.2$>, T, TU, 7N,
-NO, Ff.02, OT^, NA^IECI)
DIi'E'lSlON SO,25,25,H)
.CE: (S,C)
..-c PARAMETER NPP-IJ 3
C DIMENSION LOCXIVPP), LOCY(NPP),3
CLOCZ(N°P) > A.-T(NPP) 3
CA!
A
A
A
A
A
A
A
A
A
01 HENS I ON LOCX(ll),1.0CY
-------
rr, Sr
/xacs/ ON 13 OCT 7j AT 09:51:2?
CCCrcj.
CC""-Orj OX, OY, cZ, PMiNl1)f P^AXO), «N(<1>,
1 Cl M ,JS,75,<.) , DV(H,?Si?5!l DHltiZD.ZSI) DT,
2 UK25.25J. U2(25,25), Vl(25,25), V2(25,25), 7i TL ,
3 I" F'tO, F«J02, OTH, fAMEtl)
OIIEN:!OM so ,;siZ5ii)
ESU, VALEMCt (StO
D*TAFAC/tH/
N.--2S
00
.', S F * M
CT o 100CO.
oo 2 u « ), NY
- - f Y * J - L
Pfil,,7 4COC» (Ul ( I r J) i I»l i
2 CC-J7I, UE
DO ", L = i. NY
J - J Y » l - L
00 6 L
1 l NY
P'I','7 4002, (H| I ,J) , l«t ,NX)
& CO"7,NJ£
A030 FCf^AT <*1U' / |25F5.2>J
6001 FCP.-HT t'!V» / (2SF5t2)>
60C2 C0"-A7 (MWt / (23FS.Z»
DO 1C J « 1, NY
OC 10 I « 1, NX
IF (UiU.J) .NE. 0«) OT t> HINtOTt OX/ASS (UHI i J)I 1
IF (V>(1,J) .ME. 0.) DT " KINtOT, 0V/ABS(VH1,J)})
IF («(I,J) ,N£. 0.) OT " MINtOT, DZ/ASSIVf( : iJ) ) )
OT = 0"«FAC
CO
oo
-pu
-fc.
-------
FLOWCHARTER QY FORFLO /X808/ ON 13 OCT 71 AT 0?:Si:30
SODROU7INE
OX, 3T, 02, PMJNCM, PMAXCO. WMC
S), 8ACICCH1, Ctlr^S.ZSf 1) > OV[U,25.23),
OM(" i25,2S> , DT, Vf(?5i25), Ui(:5,2S), U2(
25.2S)i Vl(25i2S), V7(25,2St, T, TL. TN(
F\0, FN02, DTH, NAfl'CM)
-3UIVAL-f'C." (S,C)
DATA F/C /."/
DT n 10000.
A
A
A
A
A
A
A
A
^
A
A
A
A
A
A
A
A
A
I
I
J«NY»l-L
I
I
/ PRINT ACOOi IU1 ( I §J) t 1«I iNX)
I
C 23 I
I
I
I
: J=NY»I-'-
t
i
/ P?1NT 6001. 1 VI ( 1 . J) • 1»1 .MX )
I
C "3 I
,
.
•
.
I
.
*
*
I
CO
-------
I
/ x
A , .->; 00 6 L « 1» NTf S
A : » »•« :
A I
A I
A t .. t t ..,.....<.......,.
A :,,.......... ................;
A I
A I
A /( ....I,.,... ....t .,.. ...,t... t I».t......,/
A / PRINT 6C02, (Hilij).1-1.NX) /
A /,.... , • .,....,./
A I
A C 63 I
A « ,.,.,,,...,.,
—. _--: CONTINUE :
i , ;
I C 4000 FORMAT (MU* / (2SF5.2IJ 3
I C 6001 FORKAT ; DO 10 J = it NT '
A ;.,.....«.. ....................<....,..;
« ;
A I
A ..........,.,.,,...
A 3 . .—.>; oo 10 i « i, NX :
A s :,...,..,, , ,,...:
A 3 I
A3 I
A S /' «»•»•»••••• » • f • > .<» • » « •»!«•« ( • • •• «t •«. ««\ FAI.SE
A B < IF ai(IiJ) .NE. 0.) > . — — ——.
A 5 \.. ., ,,,,./•
f, B 1 TPUE
^ 3 !
A 3 ....t .(.i... ............... >.
A= : OT»KINiOTiOX/ABS(Ul(I.J))) 5
A B : :
A3 t I
A 5 0< " • — ----------
A 3 1
AO /" «•«»"«» • « •• • ttff. « • i t <»«• fnt»t»f»i ,,iv FALSE
A3 < IF !YJ!I»J) ,NE. 00 > — - — •
A 3 \ t
A 3 t TRl/E
A 3 I
A 3 ,1
A3 t OT = MINiOTi 0»/ABS(VI(IiJ)I) !
A 3 :. ,,,.,...,..,,,,......., ,.;
; 5 i
A 3 o< —--------------- - -- - --- - ----
A 3 I
A ? /t.tt.tttt....•••.it,.,...,...,...i....,,.v FALSE
A : < ir (nn>j) .NE. o.) >..-.-...._.,.
r. Z \ /
A 3 I 7?UE
r a 1
/• 3
OJ
00
oo
-------
A
A
V
B
3
' I
t 103 :
i
i
i
f
t
OT « DT«FAC
: RETURN
\
I
CO
-------
3SR-844
Appendix B
LISTINGS § FLOW CHARTS FOR NEXUS/P
-------
eroa«rt.O'FOPFLO»s NEXUS
FLO«C"ARTED BY FoRFLO /X808/ ON 13 OCT 71 AT 09:55:29
033COI .
C3CCC2 •
Code", .
030003-
C"CC6 •
Ccc;c7 >
c;c::a .
C2C3r°.
CCCC 3 •
C:D: i i »
C2" I Z»
c:c: . 3.
CG'Ol "i
C2C2 ' = •
C C 3 c 1 6 ,
OCCC i 7 »
CC3C13-
c;c: 1 9 »
crc j23 •
circa i .
ccccJj.
CCT023 .
zzzzn •
CCC-J5.
CCCC.'i.
c;r;? 7 .
cccc; 3 .
CCCC2'.
C ~ Z C 3 2 •
OCC333.
ceo::-'.
C ~ CC 35 •
C33C36 «
o : - c 2 7 «
CC3I30 •
CC303'.
cccc-o.
p ••» i r u i
. . - _ I •
c:::i2.
C23C u 3 »
CCCDi".
c::oM5.
ccrr H «, .
c;ccu7 .
C t-C "c •
C3C3-9 .
cccc=c«
CCCC5I .
C
C
c
c
c
c
c
£
c
c
c
c
c
c
c
c
c
r
V
c
c
c
;
C
TEXU3/F, KUf'E'UCAL EXA" I NAT I OS' OF URBAN S^OG WITH PHOTOCHEM 1 STKY
IS A 3-Olf ENSiONAL PAPCEL- I N-cELL METHOD CODE
PARAMETER NX»2Z,N'Y=2i,"Z"i,Nsp«s,MAX»750Q
PAPA"CTEP NX~1»)X-I, NYMl"f.'Y-l, N Z .'. 1 » N Z - I
PAI:A'/ETER NA = N'X»«Y»NZ, \TOT»1»1AX»7*NA
PARA^^'E" NA1«NA*1» NA2^NA 1 ».'.A , NAl = NA2*NA, NAl»NA3*NA,
J \f,:"7 "IA+'. , MA Y'NAX *-»,*X i «AZ«^A*»"AX , HA •!* VA Z «1AX
CO-xO' OX, OY, OZ, TOX, TOY, T ? Z . PMIN, PMAX, v,N, BACK, DT,
1 Pf )3T"(NSP) , PMXS*P(;,SP) , «>|STOR(N5P) , BCxSTR(NSP), NPSTOR(NSP)
2 FACTlNSP), F'fX, TNT, F^•Z, V, "V, NP, T, ICYC,
3 A = r>AY f STOT) , FLUXl(NX), FLUX2INX), FLUYl(NY),
S FLL-T25NT), Qxl diNXiNSP) , B X 2 ( " , N X , N3P ) , 0 Y 1 ( 1 , NY , NSP ) ,
6 BV2(H .NY.I.SP)
DI-;NSIO. CU,Z,NY,NX), cc ,
1 v. ) > ( U , A RR A Y ( N AS ) 1 , ( X , ARR A Y ( NAX ) )
2 (V,AP1AY (NAY) ) , (Z,A"?PAY(NHZ) > i ( »T , ARRAY 1 NA» > J
CO",MON OX,OYfoZtT9XlTDY,T!>Z,Pr'.tN,PHAX,MNt3ACKlOTl
|Py'STP(5), PMXSTR(S), Y-NSTOPtSl, 3CKSTRI5), NPSTOR15),
2FACTIS). FNX, FNY, FNZ, V, RV, KP , T, ICTC, AR? A Y t 5093* ) ,
3F'_jj;t22),FL'JX2(22),FL'-IYl(2:),FLUY2(21), BX 1(1,22,5), 6X2(1 ,22,5
"3rl (1,21 ,5) ,8Y2(".Z1 ,S) ,
5i*Jx.NY,^-?p,^/y,KXKi ,NYM! ,NZMI ,NA,NTOT,NAX,NAY,NAZ,NA^,NZ,
6HAI N t "* \'A3 F'AU DUf^NY
OJ-EJS'ON C(1,2!,22), CC(1,2I,22), OV(f,21,22>, BH(1,21,2Z(,
1 U( 1 > 2 1 . 27 , 3) , 5(1,21,22), X(9SCO), Y(95oa), Z19S03), "TI9500)
ECU]V;LENCE (C, ARRAY), ( CC , A* .NEX02700
NEX0200C
N£XO?*3QO
NE/030CO
NE xoo i c o
NEX0320Q
NEX0330C
I (S, ARRAY ( 3477) ), ( OH , A«R A Y ( 551* 5 ) ) i I U , ARK A Y ( 7 39 3 ) ) , ( X . ARR A Y < 1 2737ME X03 <( OC
C
C
C
C
c
k.
c
c
c
c
c
21 1 1 X'jA^KATlZ^TJ/M, \^,h«KflT*JiVJr/l» \n(,AKnttTlT|-IJ/IJ
II-TE^E"* DUMMY
CO">-ON /OUT/ isp(S), f.CYC, NOL, TDELT, PTIME, CONC, PA«T,
1 S?i5), !PSP, RSTRT, SAVE
INTEOE" s?> SAVE
LOGICAL CONC, PART, RSTRT
I,J,' ARE USE3 AS IVOICIES IN THE X,Y,Z DIRECTIONS,
"ESPECTIV!." T"RCUG1CL'T THE CODE, AND N IS THE
POL'.UTA''D SPECIES NU"SER
•
-------
c
r
0
c
;
C
C
c
;
c
c
c
0
c
£
c
£
r
c
c
c
r*
r
C
C
~
c
~
c
z
c
c
c
r
c
c
r
3
c
c
c
c
£
c
c
c
^»
c
c
7
~
~
~
c
c
„
;
-
-
c
-
2
r
n
£
-
:
-
;
5
5
-
5
-
s
2 .
^ .
s .
4 .
7 ,
zccss.
2
3
c
-
-
*
™
2
J
o
£
1
-
-
-
7
^
c
w
;
~
^
3
Q
n
£
o
£
o
2
2
2
*
5
A
i
i
o
6
J
6
4
6-
6
9 ,
- ,
1 ,
J ,
3.
u .
J .
<> .
7 ,
S '
7 .
c DC-C.
3
3
-
c
^
-
r
-
~
z
-
~
o
?
~
2
2
7
5
c
-
3
c
c
^
r
7
c
~
2
c
-_
~
^
-
J
^
-
;
™
~
c
"
"
2
~
c
~
^
z
^
^
*
c
2
3
2
f\
Q
c
0
s.
-
;
-
;
2
"
f
~
"
-
-
c
2
;
c
J
;
"
;
7
7
7
7
. .
3 ,
*• t
75.
7
7
'
7
3
s
c
3
J
o
3
5
s
c
?
»
c
?
9
'
'
c
•- *- -1 O
c
c
7
"
~
•
' ~
£
2
c
~
~
~
z
~
'2
I
'
I
I
•
I
^
£
^
~
"
^
;
6 •
7.
C •
? .
•f •
1 •
3 ,
u ,
S.
6.
7 •
a .
0 .
-.
i .
Z'
: .
" .
5«
s .
7 ,
= .
^ •
;.
i •
2.
2 .
" .
S.
/• .
C7 •
r^ f. f, , -. Q
-< - -• . J 3 •
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
(•
^
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
r
c
c
c
c
".AX
ox
DY
oz
TOX
TDY
TCZ
PH'N
"•I AX
>'N
BACK
07
PC'.STP
PFXSTR
AN5TOR
BC'STR
N'STO?
FACT
FI-X
FNV
rf-Z
V
RV
NP
T
ICYC
C
OH
u
cc
X
V
Z
«T
FLUX:
" L U •! 2
FLUYl
FLUY2
SX!
BS?
EY1
BY2
S
C/
IS?
"^CYC
I.OL
TDCLT
PTI-E
co> c
PA'T
SP
9
e
=
a
=
a
«
e
n
=
9
a
«
=
»
a
»
«
0
•
a
n
•
»
e
w
e
»
•
a
e
e
K
e
"
e
=
a
a
a
9
K
a
a
a
Q
D
B
ti
a
•
MAXIMUM NUMJEP OF PARCELS FOR ANY SPECIE
CELL SIZE IN THE X DIRECTION (METERS)
CELL SIZE IN THE Y DIRECTION (METERS)
CELL SIZE IN THE Z DIRECTION (DETERS)
2*OX
2»DT
2-OZ
MIMCuM PARCEL HEIGHT ALLOWED (KG)
NAXIfi/f PARCEL WEIGHT ALLOWED (KG)
AVERACE PARCEL WEIGHT (KG)
SAC^-C^UND CONCENTRATION (PPM)
TI..E STEP (SECOi.DS)
STORAGE, BY Sr'CCItS, FOR P,«I\
S*0"ACE, 3Y SPECIES, FOR pfAX .j
STORAGE, BY SI'ECIES, FOR *N
STORAGE, 9T SFlCiES, FOR SACK
STO'AGE, BY SPECIES, FOR NP
F,,CTOR FOR CONVEPSION FROM WEIGHT TO PPM
FLCATT'G NX
FLOA'tt.'O NY
F^OAT.IIG VZ
VCLUHE FACTOR FOR CONVERTING WEIGHTS TO PPM
1 . / V
NUMBER OF PARCELS
ELA'SEO TIKE (SECONDS)
cuHtEr T CYCLE NUMBER
CC (CEN'PAT ION AP^AY ! >< G UNTIL CONVERTED TO PPM pOR
UOPIZOST«U DIFf'JSION APPAY (M/SEC)
VELOCITY ARRAY- U(K,J,I,1) o VELOCITY IN X
U(K,J. i ,2) « VELOCITY TN y
u(Kij,i,3) ° VELOCITY IN z
FRACTIONAL CHANGE IN CONCENTP AT I ONS
PARCEL X COORDINATE A"°AY (CELL FRACTIONS)
PAPCFL Y COORDINATE ARRAY (CELL FRACTIONS)
PAFCEL Z COCHDI 4ATE A"RAY (CELL FRACTIONS)
PARCEL -VEiGHT ARRAY «G)
FLUX OUT DOTTED Qf VERTICAL GRTO (KG)
FLUX OUT TCP OF VERTICAL GPID (CG)
FLUX OUT LErT OF VERTICAL GRID (rG)
FLUX OUT RIGHT OF VERTICAL GPID (KG)
8CROER CELLS THAT STORE POLLUTANTS GOI^G OUT TOP OF
BOTOE1? CELLS THAT STOPS POLLUTANTS GOIVO OUT BOTTOM
NEXCSZOO
NEX05300
NEXC5HOC
NEX05SOO
NEX05630
NEXOS700
NEXOSUOO
NEXOS700
NEXC4C.J
N E X 3 4 To 0
NEXC42CO
MEXOA300
N£X06SCO
NEX04SOO
NEX04600
NEX047CQ
N'EAC6300
NE*04?00
f'E*C7C90
II E A 0 7 I c 0
"E1C7Z03
N£*07300
NE'-07-iOO
NEX07SCO
NEX07400
NEro77EO
OUTPUT JNCX07800
NEXQ79CC
OIKECTION
n 1 0 c r T T riN
y 1 KC, c i lui.
NE<063CO
NEX08HCO
NEXOBSDO
NEXOB600
NEXC3700
NCX06803
NEXC8900
NEX09000
NEX09IOO
GRID (J-NV)
OF GRID (J-l)
SOJ5E» CELLS THAT S70FE POLLUTANTS GOING OUT LEFT OF GKIO (1=1)
S3^DER CELLS THAT STORE POLLUTANTS GOi' OWING PARTICLE POSITIONS ON UNIT 2s
(DEFAULT » ,FAL3F«)
I'AI-ES OF POLLUTANT SPECIES
OF GRID (I = ,NX)
NEXO
-------
C30107.
03C in.
C33. i7.
C ~ C 4 ' J «
C33: IS.
CCCi i 7.
CD"1 i 1 S .
C 33 1 ' ° •
oc: ' 23.
C C 3 ' 2 i .
03C 22 ,
CC3.23.
C33 1 2J .
333 I 25.
033 • 26 .
C3-3' 27,
c::- 23.
ccciz?»
033 31 .
CC3132.
C.C3I 33.
033 1 3" .
C33 . 35.
C33-; 36.
C33 . 37.
033-39.
33313?.
C C 3 ! 'i 3 .
C33 ."! .
C 3 3 i M 2 .
C33 ' "3 .
C3C i HH .
CC3!"5-
C33 ' S6.
C331"5.
£ « o t (J O
C33 i S3-.
033 .51 .
033 I S2«
1331=3.
C3-C.S-.
C33135.
C33 15.7.
w - - 1 3 O «
C33-S7.
C33 : 6C >
33CI6' .
CC3 , 62.
CC3 <0.
C3-C1 6" .
033165-
C NPSP « NUCBER OF SPECIES TO BE PRINTED
c PSTRT « RESTART FLAG (DEFAULT • .FALSE.)
C SAVE » C'CLE TO BE DUMPED ON UNIT 27 FOR RESTARTING
DATA DX,Ov,DZ,T,FINAL/2"32I8.6'?f",100.,21600.,77ZOO./
C fCLECULAR "EIGHT OF POLLUTANTS
C NO N02 HC 03 HN02
DATA FACT / 3o.» •**«» s*«i *
-------
CC3I66.
COC:67«
C3314S.
CC3167.
tc
C2"73.
ccs, 7«.
OC:.73.
CZ: 76.
CC:'77-
C3r SC.
cc:- a i .
cc:: 52«
cc'.sj.
C"!S1.
C32.SS.
ca;;e6.
INPUT NEW PAPCELS
CALL SOuPCE(N)
CALL OUTPUTCN)
CALL R0|
NEX16700
20 CO iT.'Jl-E
c c-ccs TO SEE IF MAXIMUM RUN TIME is APPROACHING. IF so, GO TO
C C4LL S3iiAKN(200, 1, S30)
ir i!c*c .jo. SAVE! GO TO 30
!c (SAVE «G7- 0 -AND. «OD(tCYC»10) «EQ« 0> CALU WFILE
IF (T ,L7. F!\ALI GO TO 10
30 IF ISAVC ,LE. 0) GO TO 10
DUMPED ON DRUM FOR RESTARTING')
6CM3i
iO-^O FC°,-/.T t'CCYCLE'i I*.
uo T = -T
co so N » I i MSP
c TH:S MLU FIMISH OFF ANY DUMP TAPES
SO CALL OUTPUT(-(N»IOO) )
5TCP
END
NEXI4900
NEX170CO
NEXJ7100
NEX172CO
NCX 1 7530
NEX1 74C(T
NEX] 7700
IIEX17S3C
WCX J 7?00
NEX1800C
NEX1 8200
SEXI B3CO
HEX18MCO
NEX18600
NEX18700
cn
»
oo
-------
NEXUS ' FUOWCHARTEO BY FO&FI.O /xsos/ ON 13 OCT 7i AT o?:ss:32
BEGIN
I
I C'NtXOOrDO 3
I C NEVJS/P, NUMERICAL EXAMINATION 3
I COf t'"3AN SMOG WITH PHOT OCHEM I STR 3
I C, NEXOOJIG 3
I C IS A 3-OlMENSIONAt. PA°C£L-JN- 3
I CCELL METHOD CODE NEXCOMQO 3
I
I --- C PA'SA^ETE'? NX.22.NY-2) tN2»H,NSP-3
I CS,NAX«9500 NEX00700 3
I C PAP/^ETER NX11«='X-1, NYMl»NY-l,3
I CNZlI*NZ-t NEXCC300 3
I c PARAMETER NAoNx«r,Y«NZ, NToT-1)* 3
I CMAX«7<"IA NEX0070C 3
1 C PA^CZTE1! NA,=\A»li NA2-NA1»NA,3
I C .1 NAX*7«NA+l, NAY»NAX»MAX, NAZ-3
, I CSAY*,I/>X , NA,,->MAZ»MAX NEXOHOO 3
I C COUPON OX, OY, OZ, TOX, TOY, 3
I CTOZi P1IN» P"jAX, AN, BACK, DT, 3'
I C.'EXCIZOO 3
I C 1 P1NSTPINSP)i PMXSTR(NSP), 3
I C'VNSTC1? U.S?) , BCKSTR(NSP), NPSTOR3
I CN3P) , NEX01300 3
I C 2 FACTU.S?), FSX, FNY. FNZ, V, 3
I CPV, NP, T, ICYC, «£XOJSOO 3
I C 3 AtRAYINTOT), FLUXl(NX), FLUX23
. I CNX), FuUYl(NY). NEXOISOO 3
I C 5 FLUY2(NY), 8X,(H.NX.NSP), BX23
I CliNX.NSP). BYM1 ,MY,NSP> , 3
I CNEX01600 3
I C 4 3Y2[t,HY,NSP) NEXOI700 3
J C DIMENSION CI^Z.NY.VX), CC(NZ,NY3
I C'JX), C/IVZ ,NT ,MX I , OH(NZ T'YiNX) > 3
I CIlEXOIBOO 3
I C I U(NZ,NY,NX,3), S(NZ,NY,NX), X3
I CflAXI, Y(HAX), Z(MAXJ, IVT{MAX| 3
I CNEX01900 3
i c EQUIVALENCE ic,ARRAY), ccc, 3
1 CA'""»Y (^nl ) ) , (OV,ARKAY(NA2) ) , 3
I C 1 IS, A-(9AY (NAZI ) , ( OH , ARR A Y ( N A3 3
I C), tU, ARRAYINAS1}, (X , ARRAY(NAX)3
1 C, NEXC2100 3
I C 2 (Y.ASSAYINAY)), (Z,ASRAY
-------
CO
H.ZltSl, NX,NY,NSP,MAX,r,X«l,NYHl,NZKl,NA,
NTCT,NAX,NAY ,NAZ ,NAW,NZ , NAl, NA2 , NA3,
OIW:KSICN en, 21, 721, cc, XC?500), Y(7SOO), 2(7500), N7<
9SCC)
(D^,APPAY(36?7) > , (S,AT«AY(3i97> I , (OH,
AP=i-Av(S3')5)),(U, ARRAY (7393') i I X, ARRAY!
12?37>), ( Y ,AFRAY(22M37> ) , (Z,A?SAY(
3 1 V 3 7 ) ) , (V. T,/1RtiAT(HlM37))
I'.'TECCR DOI'MY
COM"ON /OUT/ 15=15), NCYC, NOL, TDtLT,
PTI.-E, COSC, PART, SP(S), t.'PSP, RSTKT,
SAVE
INTEGER SP, SAVE
LOGICAL CONC, P^RT, RSTRT
1
--- C
I.J,< ARE USED AS INOlCIES IN 3
X,Y,Z DIRECTIONS, NEXO"500 3
C- RES"ECT!VLY THROUGHOUT THE COOE3
c ANO M is THE NExo<(60o
c POLLUTAVD SPCCIES NUMBER
C NX - NUr3Ei? or CELLS IN THE X
CDlr 2»OY NE705700
C 'TOZ » 2«OZ NEX05630
c P"IN » MI HI -UK PARCEL HEIGHT
CALLOAEO I^GI NE/CSVOO
c PMAX » fAxiMjs PARCEL AEIGHT
CALLO^EO CGJ NE.'caooo 3
C «.l = AVERAGE PARCEL '.EIGHT (KG)3
CSEX0610C 3
C BACK " BACKGROUND CONCENTRATIONS
C(P°M NEX06200
C CT = TIlE STEP (SECONDS)
CNEXC6300
C "i-STR= STORAGE, BY SPECIES,
C PKXSTR= STCRACE, BY SPECIES,
NExCoSCC
ST3SA&C, BY SPECIES,
liH .VEX06SCC
I
CO
-------
to
CFOR OACK NEX067CO
C VPSTOP« STORAGE, BY SPECIES,
EFOR NP NE*C6000
C TACT » FACTO1! FOR CONVERSION
CFR?« .rEIGHT TO PPM NEX06700
FUOATiSG NX NEX070CO
FLOATING NT NEX07100
FLOATING NZ NEXC7200
C FNX
C FvY
C FNZ
C v =
VOLU"E FACTOR FOR C Of. VE9T 1 N3
C "EIGHTS To PP1 NEX07300 3
C Rv m 1 • /V NEXC7HOO 3
C 'IP « NL'fSER OF PARCCLS NEX07500J
C T « ELAPSED TI1E (SECONDS) 3
CS'EX076CO 3
C 1CTC « CURRENT CYCLE NUMBER 3
CNEXC77CO 3
C C » CONCENTRATION AR^AT (KG 3
CUiUlL CONVERTED TO PPM FOR 3
COUTP'JT)N£X070DO 3
C OH = KO«!Iro»."AL DIFFUSION ARRAT3
C(^/SEC) SE707TOO 3
C U - VELOCITY ARRAY. U 1 1C , J , I , 1 ) 3
c VELOCITY IN x oiREcTNExcaoco 3
C U I IN THE 3
OJ
en
/a
I
oo
-------
CNEXJOOCO ,
C (DEFAULT "
-------
w
A
A
A
A
A
A
I
I
C SET UP CONSTANTS NECESSARY TO 1
CSTART TCE PUN AND READ INPUT I
CNEJI36CO I
C CAPOS OEFINIM& "RINT AND DUMP I
COPTION3 NEX1J700 I
I I
I
: : CALL IN=UT : :
i
I C INPUT INITIAL CELLULAR 3
I CcONCENTPATIONS AND SET ARRAYS 3 "
i c IF FSTRT FLAS is SET, SUBROUTINJ
I C PEST°T IS CALLED FROM IVPUT 3
i
I
1 TRUE
I
: : CALL SETUP : :
i
I
I
I
Q » 1,/fACTIN)
I
C S3 I
/ «NS701|N> , ECKSTRjN), FACTIS),
-------
to
I-1
o
I
!' T - 7 + OT
I
I
: : . CALL KEM
i
i
A.-.--- ->: DO 20 N " 1) NS°
A :.« ..,....,
A I
A I
A .
A NP . NPSTORIN)
A V.N « (.NSTCRtN)
A PMIN-P1NSTFMN)
A P1AX - PPXSTR(N)
A BACK • BC
-------
A
A
A
n
A
A
A
A
A
A
A
A
: : CALL SOURCE(M) : :
i
: : CALL OUTPUT IN) ; ;
I
I
: : CALL ROLOUT(N) : :
i
i
: NDSTO«MN> » HP :
i
C 203 I
I
i — c cfEc* T° SEE IF MAXIMUM RUN 3
CTIKE IS APPROACHING' IF SOi GO 3
CTO 3Q.VEX17300 3
c CALL SOWARMZOO, 1, $30) 3
CNEX17"00 3
I FALSE
I
I
< 0! CAvL .VFILE >
I
:: 0) CALL "FILE : :
1 I
I FALSE
I
C 303 I
I FALSE
1
: . CALL '-FILE . :
i
X
x
X
X
X
X
X
X
X
X
X
X
x
X
X
X
X
X
X
X
X
X
X
X
X
x
X
X
X
X
X
x
X
X
•
z
z
z
z
z
z
z
z
z
z
z
z
2
Z
z
z
z
z
z
z
z
z
z
z
z
2
z
z
2
Z
2
Z
Z
z
z
z
z
z
z
T
OJ
I
CO
-------
t f . t . , t ..... r . f f • t f ... f . 1 t 1 f t , / / Z
PRINT 6010, ICYC / I
- ...................... / -x ^
i z
I --- c 60HO FORMAT COcTCLE'i 14, • 3 Z
1 CDUMPE5 ON DRUM FOR RESTARTING*! 3 Z
I Z
C fOD
A----—--—---..->: oo so N ° i« NSP :
A :...>.... ....,:
A 1
A J C THIS »I!-L FINISH OFF ANT DUMP 3
A ! CTAPES NEXI8SOO 3
A I
A £ 5C3 I
A
.. —; ; CAUU OUTPUT(-tN»ioc>) ; :
1
w -
/' ..... ....... m .t,N
is) : STOP :
en
oo
4s.
•fc.
-------
GFCp»FLO«F07FLO,S
FLO*CUA?TED BT FQPFLO /xsoe/ on 13 OCT 7i AT o?!55:i6
COC03'.' SUBROUTINE DORDER(N) 80R00100
CCC:32. C THIS ROUTINE UPDATES THE BORDER CELLS, EACH OF RHICH IS 3 TJME5 AS "IDE
33CCC3- C T.i'cE AS HIGH AS THC CENTRAL GRID CELLS. PASS REMOVED FROM THE OUTER
CCCrC''' C CELLS IS REMOVED FROM FURTHER CALCULATIONS, BOROOtOO
Cr;;35« COVMCN ox ,OY ,5Z,*ox .TO* , TOZ,PUIN,PMAX • «N,BACI<,DT , SOROO&OO
C C C C 0 6 • i P M " 3 T =5 ( 5 ) , p i' X S T rt ( 5 ) , v, F^, FJV, FNZ, V, RV, UP, *, ICYC, A'RAY(509361, OOSC070C
3FLUxl(22),FLU)(2(22),FLUYl(211,rLUY2(Zl), DXI(1>22,5), 3X2(1,22,5),30 SOCbOC
CCC309. «Brl ( ",-21 ,S) ,OY2(1,21 ,S) , BOR00900
CC^r I • 6.-.A1, NA2, NA3, NA", DUl'MY UO^OllOO
C2C312. D!u:N5'.C>I C(°,21,22), CC (11 2 1 , 22 ) , 0V (1 , ? i , 22 ) , OH (1, 2 1 , 22 I , BOO12CO
CC0313. I 0(1,21,22,3), 5(1,21,22), XI9500), *(9SOO>, Z(?500), ftT(9SOO) 80^01333
CCCC'I' EQUIVALENCE (CiAR"AY), ( CC , ARR A Y ( 1 El 9 ) 1 , ( 0 V , A^ R A Y ( 3697 ) ) , B0t?0i"03
CC:C15- 1 IS, ACSA" (3697) ) , ( OH , AP« A Y ( 5 5 1 S ) ) , ( U , A>*t! A Y ( 7 3 9 3 ) ) , ( X , A RS A Y I 1 2737 B ORO ! SOO
C::C!6. 21), (Y.A'^AT(22M37)1 , (Zi ARSAYI31937)), (,IT , ARRAY(11137)) BOR01630
CCCCi7« INTESEP 00".MY BOROI703
C C 3 ^ '. 3 t C •t««««»»»»»»«»»»«»»»»»»»»i«»»»»««»«»»»«»»«»»«««»»»«»««»»»»»»«»»»«»BOH01flOO
CC3C ?. DO 100 I • 1, NX UOPOI7CO
croc:c. ^EL » uti,i,i,21/3. BOKOZOCO
C:73Ci. If C/EL .GE. 0>) GO TO 25 BO«02lOO
C:Z??Z- F-,1 = FL'JXl(I) BCKC2200
2;CC23« 00 2C L » 1, 1 BOfCZjSS
c:cc:i. Fr « -BXIIL, 1 ,N> «VEI. BOK02ioo
C:C';25. BX!(L,I,N) o BX1(L,I,NI - FX * FX1 80P02530
C3CC26. 20 ri(l = rx 80r02600
C02CC7. rwL'sltJ) = 0. BOS32700
co3:rs. GC TO 35 oo-iozsoo
CC:D2?« 25 FX1 •> 0X1 f 1 , I ,N)»VEL 30^02900
cr:::3. FLUXKD •= FXI eo?C3ooo
c;r:3 . oo 30 L •= i, 3 BOSOSICO
co:c32. FX >= sxi (L*I ,: ,N)»veu 80^03200
CDCC33. BX](L,1,N) = DXJIL.I.N) - FX1 * FX 80KC3330
CCC33". 30 FXl = FX 80^03100
C0:33i. bXl|u,l,N) « BXl(1,I,Nl - FXl BOS035QO
c;::36. 35 cor"NUE BOR03&30 ;C
C:CC37, VEL=Ull,SY,I,2)/3, DOS03730 ^
cc:;3e. IF iV£L .LE. o«) GO TO is B0'2(I,I,")«VEL
50t,QL='4,3 BC°0"bOO
FS«-TX2 (L*l , 1 ,N) »VEL BOKO''90Q
5X7(L,!,'M"3X2(L,1,>11-FXI-'FX BOKObOCO
5C FXIn'X
-------
CCC052.
CCOC53.
CCCC75.
C3 ce=.
cc
G'-7.
j:°.
:c 070.
ICO CONTINUE
00 230 J • 1t NT
IF |V£L -GE. 0.) GO To 225
FXI « FUUYI(j)
DO 220 L » It 1
cx = -BY:IL.J,N>»VEL
3TML,JfM « BYItL.JiN) - FX » FXI
220 FXI •> FX
F.UYI(J) » 0.
CC TO 335
225 CX I = BY 1(1,J,N)«VEL
FLUY I ( *>) * FXI
00 330 L • 1. 3
FX » 3*1(L*ItJtNl'VEL
SYKL.J.M » BTI(LiJiV) - FXI » FX
3CO FXI = FX
3y!|H,J,N) = BV!(«,J,N> - FXJ
335 C C f TIN U E
•'EL « u ( I i J t Nr 11) /3.
iF (VEL «LE. G«) GO TO 215
CII"FU-T2CJ)
00 2HO L=ltM
FX==Y2(LtJtNI.VEL
2HO FX I«Fx
FLU>2(J)=3.
GO TO 200
Fl_Uv/( J) "FXI
DO 2^0 L=li3
rX=-BY2(L«l,J,N)«VEL
OY2(*.,J,N)»3Y2
-------
FO"FLO /X3C8/ ON 13 OCT 7J AT 0?:SS:S3
ts
H-'
en
SUBPOUTINE BO'DERIN)
1
1 c TKIS ROUTJSE UPDATES THE BCRDER3
I Cc^LLSi EACH OF t.HlCH IS 3 TtPES 3
1 C T.Y1CE AS HIGH AS THE CENTRAL
J CcRIC) CELLS. MASS REMOVED FROM
3
3
3
C CEuLS !S "E.'OVEO FROM FURTHER 3
CCALCULATIO'IS. SOsOO-fOO 3
COMMON OX,OY,DZ,TOX,TDY,ToZ,PMlN,PttAXiWN,
SACftST, P>!N5Tff(Sl» PMXSTR{5I, rtv5TOR(5li
BC>CSTR!S), NPSTO'itS). FACT(S), FNX, FNY ,
r^Z, V, IV, HP, T» ICYC, A1?? A V ( S09 36 | ,
FLUX!<22>,FL'JX2<22>,FH.Yl(2l),FH,'Y2t21>,
3X1 (t,22,5) , OX2Ci.22,S) ,0^! O,ri ,SI ,S"2(
NAi,
NA3,
NTC7 ,'-AX ,f,AY,'.AZ ,MA,V, ARPAY( I 8"?J ) i
, AR1?; Y ( 3697) ) » IDM.
1 73?J) ) i (XiA"FAY(
I2C37)), (Y,«f
INTEtER
• >:
:—c
i c •«
I C»BOS01803
I
I
00 ICO I - li NX
YEL * U( 1 ,1 , i i2)/3.
i
'•3
••3
3
IF (VEL .6E- 0.) GO TO 25
! FALSE
FX- « rs uxi f 1 1
*\ TRUE
>..-...
./
en
:»
00
-------
A
A
A 5
A 3
A D
A 5
A e
/, 0
A
A
A
A
A
A
A
A
A
A
A
A 6
A S
A 3
A 3
A 9
A 3
A 8
A
A
A
t
A
A
A
I /,
1 //
1
: . FX = -3x1 » o. :
1
I
C 253 I
: FXJ . sxit 1 1 1 iN)«V£L :
: FLUXI ti j - FXI :
t
i
i
i
: FX » sxi ('.»! , I ,"<; *v£t.
: 3XS OXl(M|IiN) - FX1 :
I
C 3S3 I
CONTINUE :
: VEL » u( i ,f»r,i ,21/3. :
I
I
X
X
X
x
X
X
X
X
X
X
X
X
x
X
X
^
X
X
x
x
x
X
.
z
Z
7
7
7
Z
Z
7
Z
z
z
z
7
7
7
Z
Z
7
Z
7
Z
Z
Z
Ul
CO
oo
•f*
-------
A
A
A
A
A
A
A
A
A
A
>
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
1
A
A
A
fl
t,
e
,\
A
A
A
3
C
5
5
3
2
3
e
5
3
3
a
: FALSE
i
i
: Fxi-Fuuxzu) :
I
I
i
i
: F7,»BX2 (L , I iM) «VEL :
: 3X2.')»SX2(L) I .Ni-FXI+FX :
I
C 5C3 I
I
I
: 3X21 1 , i ,N! • 8X21" ,: ,N)-FX!
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
i
z
z
z
2
z
jr
z
z
2
7
z
z
2
z
z
T
z
z
z
z
z
z
L
z
z
C/3
*3
i
co
-------
t—'
oo
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
'
A
A
0
3
5
3
0
a
3
0
r
I
I
1 J
I
: VEL « u< i , Ji i , i )/3. :
t
i FALSE
I
I
: rxi - FLUYICJ, :
I
i
i
i
: FX • -CYI IL, J.NJ.VEL :
: 3vl
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
t
A
A
A
A
A
A
A
A
A
A
A
3
3
3
O
9
9
s
0
s
: BYl(LiJ»M • 3YKU,J,N> - FXl * FX : /
i
c acc3 • i
i
: OYId.J.N) . OYId.J.N) - FXl I
1
C 3353 I
: CONTINUE : -
: VEL • ui i , J.NX, i j/3. :
I
i
i FALSE
i
i
: Fxi.FuoT2(j, :
i
i
i
i
: FX«3Y2iL, J.N) »vcu :
i
C 2HC3 t
I
I
: rurrxJi-o. :
i
C 2H53 1
; FXl^ur2:i;';xrVEL ;
i
X
X
X
X
X
X
,
X
X
X
x
X
X
X
X
z
z
z
2
z
2
z
z
z
s Z
z
z
z
z
z
Z
z
z
z
z
z
2
z
z
CM
CO
X3
I
oo
-------
to
tsj
O
A
A
A
A
A
A
A
A
A
A
A
5
3
3
B
3
0
J
i
; ^a^^uXo-pM.Fx •:
i
C 2S33 !
I
t OY2(M i J,N)=DY2 (1 , J,N) -FXl t
I
C 2003 I
:
i
X
X
X
X
X
X
X
X
*
X
X
or>
X!
i
CO
-------
SFO*«-LO.rORFLO,S CHEK
ri.Or(CHA=>TED BT FflRCLO /XB08/ ON J3 OCT 71 AT 07556:10
w
02:00: .
cco:c2«
CSCCC3.
co;oot .
sroocs.
033C1 ! .
C3"3!2.
c:coi 3.
0:0334.
CC C2 3 7 •
SU3POUTINE CHEK (NO, N02, HC, 03, HN02 t ONO, ON02, OHC» 003,
1 DH-.02iOO,T!H)
C THIS ROUTINE SOLVES THE CHEMICAL RATE CHANGE EQUATIONS
C THE FRACTIONAL CHANGE IN CONCENTRATION OF EACH SPECIE.
DATA CO,C2,C3,CH,C5,C6.C7,C3,C',C10,C11,C20/.M,
1 2.4tj370cij£f'icAL REACTION RATE VARJES SINUSOJDALLT FROM 4:
C17-.C.O
Tl-ii«7lM/3600.
;F(T'V-H.GT.18.S.O'?.TI'!H.LT.6.5) GO TO 5
SI,>.C = s:N!ITI^H-6«5>», 16"01502)
C:2=C23«SINC
5 CO.-Ti.^'JE
D7 i l.E-9
"(1! * SO
Y12) « .02
YI31-HC
"13|«HC» 0,5
Y(") » 03
Y(5I ' ".J02
10 CONTINUE
EO » CI»Y(Z)/(C2 » Ct«Y(3))
Y(A! > EO
DI=C6«Y(11 * C7-Y(2)
D2=CS«Y(3) « C8*Y!1) * C?*Yt7)
F2=(2.«CH.Y(4)-Y(3)+C10»Y(S)»Y[3))/01
FU=C12*T(51/D2
012=1.0-F1«F3
EP02= (F2*F"i»F 1 1/012
v ( 7 ! = E. •) 0 2
EC" =
F-"C2 = C S • Y I 8 ) • Y t 1 )
• FHN02
2.«C11
CHM005QO
001 / CT-006CO
CM1C0700
Ch«OCHOO
TO 6
CH-OIOOO
C M -I C 1 100
CHf.OI200
I30PM
CH101HOO
CHfOiSCO
CH.-01 400
CH-OI70C
CHMC1dOO
CHM01900
CM-02000
CHf-02 100
CH-022CO
CH-02HCO
CHrC2300
CH-024DO
CHK027CO
CHM02flOO
CH1030CO
CH^OJ I 00
Chf03200
CH«03303
CHM03SOO
CHf03400
CH''037CO
CH103800
CH"t03aOO
CH"OHOCO
CH.-OM 100
CH-CM600
CHM3S700
CO
^3
I
CO
-pi
CHiOSlCO
-------
tsj
OCOC52«
CC3C33.
C3CCS1.
CC3333.
C3C356"
C"CS7.
cc::ss.
c c c ' " ? •
c ;3;;i3 «
CDCC41 .
cctc^z «
c-:-63.
c;:; j1- «
c::3*i-
c:f:47^
c ; : : 6 s .
c c ;• 3 6 ° •
cc:3"3 •
Co;c7 i .
C3"372.
ore:'-'.
£33 : 73 •
cr:c76.
C33-77 .
uC - w 7 *• »
cj:c 7° t
ccccsc-
c;;:o i •
CC"g 2 .
c c 3 : a 3 .
cccsay.
c;:cc3.
CC -CC6 •
C:33S7.
33CZ33 .
co:;:?.
C3CC9 1 .
cc;;r2 >
o::s* z •
c;";'0 .
cc:cc3.
C33C?6 .
C33C"7 «
CCCD93 .
c;;:?? .
C"3 ' 3O
3?t C ' .
C3 . 7 J •
~~~- ~? •
•3- -.-'
:::!cT.
C;C i C7 .
: : c . : ? •
RNO • C3«Y(t) » C6»Y(7) * C8»Y(8) + Ctl*T(2)
3 ; 1 ) • R.'- 0
RN02 • Cl * C7»YI7) • C7«Yt8) » ClMYtl)
8(2) « PN02
PMC " C1»Y<6> * C5»YI8) * ClO»Y(1)
3(3) » P"C
P03 • C3«Y( I > * CIO»Y<3>
E ' " ) » R03
Si T) « C12
c TEST TO S:E IF SIGNIFICANT PHOTOCHEIICAU REACTION WILL OCCUR THIS
C TIME STEP, AND, IF NOTi USE SIMPLE RATE EOUATJON AND RETURN
IF ( ITEST .St. 01 GO TO 1?
1TE3T •> 1
P'=PATE( 1 )
F3--Z (31 «Y(3)»TT
f " - 3 A T C ( 1 1
K3=^nTE \S )
XTFST • A93((Ytll » Rl 1/IYI2) * R2 )) / ABS { Y (1 ) /Y ( 2 ) )
•F(Pi.£C.C..ANO.R2.£9.0.l XTEST»1.0
IF(Ai5S GO TO 2Q
Y3( i ! = A(i )»OT * Y( I 1
GO TO 30
20 '511) = AID/BiD » (YC) - At t)/B(I»*EXPI-B( I )»OT)
30 CONTINUE
£0 » Cl«YS(2)/(C2 * Cf«YSC3n
TSI6) - EO
01=C6«YS( 1 ) » C7»YS(2)
D2=ci'TS(3) » C8«Y3tl) » C'*YSl2)
F ;=>2- 'C3>VS (31/01
F2'-(2."C'"YS(6)«YS(3* + CJO'YSt1l«YS(3))/Ol
F3=3.5'C*»YS( 1 1/32
F1=C12*YS(5)/02
C.Z= 1 .C-FI «F3
£-02= (F2-rH»Fl 1/D12
Y3 !7 1 » £"-32
£0J •> (F2»F3 + F1 1 /OI 2
TS(O) •= ECH
' .0 • C1-YS(2) » C12«YS<5)
t- ( : ; = F c
F .52 « C3"»S(tll «YS( 1 1 * C6«YS< 1 )»YSI7)
/• ! 2 1 = F 32
F03 « C2-YSI61
A(M) = F03
rt"02 » C9'vS(ai»YS( I 1 » 2.»CI I«YSf H»YS(2»
CHH05200
CMM03300
CHM05«00
CHrosseo
CHMC56CO
CH-05700
CHJ'C'iSOO
CH!-039CO
C H C C 4 3"J 0
CMf.CAlCO
CHM06230
CPf063CO
CHP06MOC
CHI106S03
CH1C6700
CIJM06SOO
CKf'36?00
CH-07000
CKf'07 I 00
CHM07200
CHf 07HOO
C""07 303
Ct'K076cO
CH"077CO
CH-0 '800
CH»'079CO
CH^OSCOO
CH"08IOO
CHf,08200
CHMOOJOO
c^-oe-oo
CHf 05b30
CH*!086C3
CM r, 00730
CHM08800
CM1C87CO
CM109030
CH-OVi 00
CH"0?2CO
CMHOV3CO
CHMCV'OO
CH"0?^00
CH1376CO
CHMOV700
ClM09«00
CHMOV930
CH11C033
C"iJ 0 1 CO
C M f 1 C 2 C 0
CH,"10JCO
CHr 101CO
C H M I 0 S 0 Q
CMMlCiCO
CHMI0700
OfftlOUOC
OJ
I
oo
-------
COO 07.
c:o.io-
CCC .!.
CCC!IS.
c:r,is.
err;i4.
CCO, .7.
cc;-'s.
CC~ 17-
co:.22-
C70.23.
* C8«YS(8)
c:3*35.
o:ci36.
roi 3s.
ccci°z>
c;o "3.
(* f- r* 14 y
ccoi-'S.
C3C116.
OOC1MS.
- C3»TS(t) • C<
= F ^0
• Cl * C?»YS(7) »C9»YS<8> * C11»Y3(M
RNO
B! 1 )
R .02
612)
^f-C » C1'YS(6) « C5«YS(8) • CIO*YS(1)
0(3} » RHC
••03 = C3»YS(1) * C10*YS(3)
0 (S ) = P03
B 13.) = C 1 2
DO 30 I = 1 , S
IF Oi:>»D* .oT. l«E-7> GO TO MO
f ( 1) «• A! ! ) «OT * Yt I )
GO TO 30
MO T(l) * A(1)/B(I) * (Y(l) - At!l/BtI))»EXPC-Btl)»OT>
50 CO'JT'N'J;
T |>.r B T-7-E * DT
IF (T J1E .GE. TT) GO TO 900
C * 0.
DO AC I " 1,5
iF (T(I) .LT. 1 .E-1) GO TO 60
C a MAX(C, ASS(tYd) - YSCIJ/tYtl) * YS(II|M
60 COV-MVJC
IF (C -LE. .1»ER) C « •1»ER
07 » f I JIE"?/C«DTi TT - TIUE)
CO 70 10
900 C3f,71\tC
DN'T » Y(1)/NO
,r MO .£0. 0. ) ONO - -Yt 1 )
0 10? " Y ( 2 ) /f'02
»F (N02 .£3. 0-) ON02 » -Y(2>
DMC » Y(3)/KC
DH'C * Y (3) /"C * 0<5
IF (riC .C"« C.I DHC « -Y(3)
003 - T!»;/03
I- (03 .E5. 0. ) 003 •» -Yd)
C.)
-YtS>
CHM10700
CHK1 1003
CH-I 1 100
Cf«l 1200
CH11 13CO
1SOO
l 1500
CHM i : 400
C^-M 1 i!3
Curl lb-50
CHM1 1900
CHM 2030
CM«122CO
CHMI2300
CH«12HOO
CH"1 2300
CH-J260C
CMM 2700
CH^I 2300
CH'112900
CHM 33CO
CH-13IOO
CHH13200
CH"13300
CH"U 3"00
CHMI 3600
CHMI 3700
CH.-I 3JCO
CHrt| HOOO
CH111300
CH11 ^600
CH«|H700
END
co
-e»
-Pa.
-------
BY FoRFlO /X803/ ON 13 OCT 71 AT 0?:S.*:i?
\ SUBROUTINE CHEM (NO, N02i KC, 03, /
S Pt.QZt DUO. DNQ2, DHC, DO3, OHN02,DD, /
\ Tin /
\ /
I
I 1 THIS ROUTINE SOLVES THE 3
I CCHcllCAL "AT: CHANGE EQUATIONS 3
I C THE FRACTIONAL CHANGE IN 3
I CCONCESTRAT ION Or C.ACH SPECIE. 3
I CCHfOOHOO 3
I
I
DATA co,C2»C3,c-<,cs,c6,c7,C8,c?,cio,ci i, :
C20/.", 2-6SE6, "0., 61CO«, 80., 1SOC.,
•4., 10., 30.. .0125, .Cl, «COl /
Olw£NS|0'J Y(8)i r5(8), A(6), B(3)
REAL '-'Of W02
DATA ER /.C5/.
RATE(I) ' tA(I I - BCIl'Yt1))»TT
1TCST - 0
TIME = 0.
I
i—c PMOTOCHE"ICAL REACTION RATE 3
I CVARIES SINUSCIDALLY FROM 6;30AH 3
I CTO 4:3C"M01300 3
I
I
Cl-0.
Cl2=0.0
TJMH-TIM/36QO.
I
I
/t,,,,,,,,,,,,.,,,,,,,,,,,,.,,,,.,.,,,,,,,% TRUE
< IF(TirH.GT.18.S.OR'TIHH.LT'4'S) GO To 5 >..---......—........-.........._-.....-..-...._...._..._..._,,0
N / Z
I FALSE Z
I I
I Z
, Z
SlNC«SiNUTlKH*6.S>». 1*981582) Z
CI=CO»SI ,'C Z
Z
Z
I Z
0 <-.---..--.--.._-..-....._-...........-._....-.._...-.........•...-__-_..._...-....».. .-0
C 53 I
COST1UJC
CT • 1.£-?
TT = DO/60.
V ! 1 ) ' '»3
-------
I C T(3)-HC CHMC2&00
I
YI3)»HC» 0.5
Yls) • 03
Y(5) • HN02
1C3
I
0<-
I
CO'JT INUE
£0 » Cl«r(2)/(C2 « C1
Y<6> •• EO
DJ»C6»Y(1) » C7«Y(2)
02-C!i»Y(3) + CS-YU) * C9«Y(2)
NJ
F2«<2.»CH«Y!6)»Y(3<+C'0«'T("t)«Y(3l)/DI
F3^C.S-Ct "t ( 1 ) /02
F«»C1:*1< S)/D2
0'.2'1'0-F1«F3
ER02» (C2»FM»FI)/012
Y(7) .
EOH - (F2-F3 •>
Y (8) => EOH
FNO » C1->Y<2) « C12»Y!5)
A (!) " F MD
FN02 » C3»Y(H)-ril> * C6»r(j)«Y(7)
A ( 2 ) = F t. 0 2
F03 •= C2»r,'i)
AIM) = F03
FHNOZ « C8»Y(8)»Y(1) » 2,«C1l.Ytl)»Y(2)
A<5) « FKVOZ
>?NO • C3«Y(") » CA*Y(7) * C8»Y(8) + Cll»
Y(2)
3(1) = PMO
RN02 - Cl • C7«Y(7> « C?»Y<8) * Cll'YIl)
SI2) « RN02
PHC t. Ct»Y(6) » CS*r(3) » CIO«Y(S>
8(3) • RHC
R03 » C3«Y(I ) * C10*Y(3)
3
-------
to
o\
A
A
RS-PATEI1*) "
XTCST « AOSMY(l) * *l )/ ,
I FALSE
I
I
Yd) " Y { 1 ) » P.t
YIZ) = Y(2) » R2
Yd J = T(M ) * R1
YtS) m Y(S) * RS
I
i
I
C 103 I
1 TRUE I . '
I I
I
C 193 I
: COH.TINUE :
i
-
•
.
y
V
U
U
y
U
y
U
U
y
U
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
J(
X
X
X
X
X
X
X
X
X
X
X
X
I
z
z
z
z
z
z
Z
^
^
i
z
2
Z
Z
Z
z
z
z
z
z
I
^
z
I
Z
z
z
z
z
2
Z
z
z
z
7
Z
2
2
en
oo
-------
w
A
A
A
A
A
A
A
A
A
A
A
A
A
f
A
A
A
A
A
A
1
C IF t5(tJ*07 • GT • I • E—7 ) GO TO 20
I FALSE
I
1
TSC I 1 « AC « )«OT -Yd)
I
I
C 203 I
• YSCI) » ACII/BCI) * CY(I) - A(M/B(IM*
: EXP i-oc i )«DT>
i
C 303 j
CONTINUE
I
I
EO « C1«YS(2)/CC2 » C"*»YS(3M
YS<6> » EO
Ol«C6»YS( i ) * C7»YS(2)
02«CS«YS(3) * C3'YS(l! - C?«YS(2>
Fl»2,»c5«YS(3)/Dl
F2«(2»»CHoYSlA)»Ys!3)*C!0«YSt')>»YS(3>J/
01
FS^C.S'C^'YSt 1 )/02
FM-C12*YS(S) /B2
Oil" .Q-F' «F3
ER02* tF2*F't«F J ) /D]2
YS(7 ) = E.TC2
EOH " (rZ»F3 » r1)/012
YS ( 3) « EOH
FNO = CJ«YS(2) » C12'YS(5)
All) = FNO
FN02 « C3«YSCn «vsl 1 ) * CA'/5 C 1) »YS ( 7 )
A (2) = FN02
F03 » C2»Y5(6)
A (u ) = F 03
FHN02 • C8»YStS) *Y3( 1 I « 2 • *C 1 1 »YS ( 1 ) »YS (
2)
A(5) » FHN02
R^^ « C3«YS1M) • C4»Y5(7) » CO'YStB) +
Cl !»YS(2)
Bl 1 ) «= R^O
1 '02 » Cl » C7»YS(7) -c?»YS(8) * Cll'YSI
; i
0 (?) = R N o 2
FHC B C"»YS(6) * CS«Y3(°) + CIO'YStH)
SI3) = -'HC
•?03 « C 3 • r 3 ( 1 ) * C10«YS(3>
SCI) a «;03
TRUE
;
5
S
S
S
. 5
: s
: s
; 5
S
.
.
•
•
U
u
u
u
u
y
u
\J
u
u
u
u
u
u
u
IJ
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
_u
"u
u
u
u
u
u
(J
u
u
u
u
u
u
y
u
u
u
x
X
r.
»
X
X
X
X
X
X
X
%
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
C/3
i
oo
-------
A
A
A
A
A
A
A
A
A
*
5 *
00 A
A
A
-
A
A
A
A
I
I '
....... _>• 00 50 I « I i S !
I
I
1 FALSE
I
1
: Y! i > • At i >"OT » Yin :
I
I
C fG3 I
: Yd) « Atn/od) » (Yd) - AU)/Btii)« :
: EXPC-DI i »-»OT> : ,
I
C 503 I
I
I
; TIME • TIME * o*
•
I FALSE
1
: c «= o. :
i
i
I
1
1 FALSE
I
5
5
S
S
s
s
s
.
•
y
U
U
y
y
U
U
y
U
(J
U
U
U
U
U
U
U
•J
U
U
U
U
U
U
U
y
U
U
U
•V
U
y
U
U
U
U
y
y
y
y
U
y
U
U
U
It
X
X
X
X
X
X
X
X
X
x
X
X
X
?
Z
z
-
I
z
z
z
z
2
Z
Z
i
Z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
.
z
z
z
z
z
z
z
7
CO
CO
-Pa.
-fa.
-------
N)
A
A
A
A
A
C « MAXtC, ABSHYtl) - YS(I»/(Y(U » YS 1 : ''.
I
C 403 I
I
I
I TRUE
I
C » .I-ER :
I
I
; OT = «i.V(ER/C»OT, TT - TJKE)
I
Z
C 9003 I
CONTINUE :
: o»o • YU ) /NO s
i
i
! TRUE
i
: ONO - -Y( i ) :
I I
I
: ONC2 . Y(2)/N02 :
I
1
I
: 0-102 « -Y(2j :
I
u
u
u
u
tj
u
u
u
u
u
u
(J
u
u
u
u
u
u
u
y
u
u
« ,
X Z
X Z
y> z
x z
z
2
Z
Z
z
z
z
z
z
z
z
z
2
CO
S3
I
00
•fe.
-------
I
w
w
o
: OHC « Y t3!/MC * o.s :
i
I
/t , t ,.,,,,,,,,,,,,,,,.,,,,, i ,,,.,,,.,,,,, ,N FALSE
< " IF (HC .EC. 0»l > ---- - —
\ ..................................... ..../
1 TFUE
I
: OHC « -T(3> :
i
o<- ---».---•-•------------—-••
I
: D03 o Yt«)/03 :
i
i
/,,>,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, , ,,,, ,^ rALSE
< IF (03 »ES. 0') > -------
\ ......... . ............................... /
; TRUE
1
: coo • -Y(H) :
t
o< ---- - -------------------- ...
I
: DHNOZ * Yisi/H'402 :
i
!
/t ... t ,,.,,,,,,,,,,,,,,,,,,,,,,,,.,.,,,,. ,N FALSE
< IF JHN02 .CO. 0.) >" ----
\ ......................................... /
I TRUE
I
t OHN02 » -Y(3)
1
o<
1
/•""••'•••'•"•• "\
: RETURN ;
\ ................... /
o-i
CO
00
-------
SFOs«FLO.roRFLO,5 OlFfJS
rLCrCHA'TED BY FQ?F'LO /X808/ ON 13 OCT 7J AT 0?!56:s3
cccois.
C
c
c
c
c
c
c
c
c
SVMOUTJNE DIFFUSIN)
T"!S "OUT.NE CALCULATES T".E DIFFUSION PSEUDO
CPI'ON D/>,3Y|DZ.TDX,?0YtTO!'tPKlNfP>'AXtftN,8AC
1 P *; N 3 7 ? < S ) , ?MxST''(S)t *, BCKSTR(S),
2F,.CT(S>, F.-X, Ff.v, FNZi V, RV, NP , T, ICYC,
3FL';Xi(22>,"LUX2(22),c'LVJYlI21),FLUY2(2!), 3X1
"3Y ' ' « ,2. i3) .SY2(ui2l ,'jl ,
3 1* , r,.sSt',M/.v,SX^l, 1Y.'.! ,HZK1 ,N'V,NTOT,N'AX,NAY
6 S A i i »i A 2 > . f 3 , '•! A J > 0 J v> M *
Oi"E'.SIO.I CM, 2,, 22), CC(4, 21,22), OV(0,2i,2
1 0(4,21,72,3), S(t,21,22l, X(9500), Y(7500),
f;t/IVALE>C- 1C, ARRAY), !CC,A«9AYIiat9)), (DV
1 (S,«="UY( 3A77 ) ) , (OH, APHAY ISGtS) ) , (U, ARRAY (
2)). (Y, A"P,iY(72'i37) ) , 1 Z ,ARhA" 131937) ), (rtT,
INTEGER DL'I 'Y
CP" '0,4 /OUT/ JSPI5), NCYC. NOL> TOELT, PTJ^E
i s <• i s ' , N P s p , RSTRT, SAVE
I.-TECES SP, s,,'.c
LOG,CA,L coHCi PA^T, RSTHT
01F00100
VELOCITY OIF00200
>C,OT, OIFC030C
NPSTOR«5) , DiFOCSOO
ARRAY 15093*) , 01F005DO
(S.22,5), eX2<'),22,5> .DIFCQCCC
CIFCC700
,VAZ (NAf'tNZ , OIFOS-'CO
OIF30900
2), DH1",21 ,22) , 01FC1300
ZI7500), mT(7500l OIF011CO
,,,RPAY(3677| ) , CIFOUCO
7373)), (X.ARRAYII27370IF01300
ARRAY 1 11 137 )) OIF01H30
D1F01S03
, CONCr PARTi , 01FOI600
CIF31700
o ; FO i acc
0! FO 1 VOO
OIFC2100
SX, Sv, SZ AFE THE LASCEST FRACTION OF A CELL IN THE X, Y, Z 0 IRECO I f C2200
"ESPECT1VCUY, THAT A PAPCEL WIuL BE ALLOWEO
OATA sx ,S" , sz/3«.3/
R:AD tsi { tut i , j, i , i j , :»! ,NXI ,0=1 ,NYJ , uuu
l J=:,;,Y), t iu(2 , j, : ,3) , 1-1 ,vx) , j=i ,NYJ , DV,
IF IDT .LT. o.) RETURN
00 5 I • 1 i NX
CO 5 J » 1 , \Y
00 2 L •• 1 • 2
U< I ,v', ! ,L I « U! 1 i J, I ,L'«DT
U(J,J,I,L! = U(1,J,1,L)
1/13, J, I ,L) « -Ut 1 ,J,I ,L)
2 u ( •• , J , I , i. ) ' U(3iJ,I,L/ .,
U ( 1 , J , I , 3 ) a C •
U ( .•> , w , 1 , 3 1 •> 2«*OT-U(2.J,t.3>
U I 3 , J , I , 3 ) » 0 '• 2 , J , I , 3 )
J I " , J , 1 , 3 ) a 0 .
S CO- TIN'JE
; o » c ' > 2 , f; « -; i
DO " G J = 2 , •, Y ,« 1
CO "C K = ' > ' I
IF (CC ,J, I ) .57. BACr ) GO TO 8
L". " , J, i , I 1 = L-CC , J, I , 1 )/OX
U(X>Jtii2) « L'lSi Ji i i2>/0v
L 1 1' , J , • , 3 ) = UC. ,J.!,3)/OZ
OC *0 "0
TO MOVE IN ONE TIME STOIFCZOOO
D I F 0 2 '1 0 0
, J, I ,21 , 1»1 ,NX) , 01FC270C
OH 01F02ECO
01F02700
; OIF03COO
D1F031CC
OIF0320Q
01F03300
01F33M30
01F03SOQ
01F036DQ
OIF03700
01 FQ3BC3
01FOJ700
01F3"OCO
DIFOH ICO
DIFOSSCO
o: Fc<<6oo
0 I rOM 700
OlFOtoCO
01 FOH'30
0 ! c 3S3CO
OIFOblCO
CO
?=>
I
CO
-------
CCCC52'
!F !< •£«
GO TO 10
13
c::c6i.
c :'. c s J •
FRO" M/SEC TO FRACTION OF CELL/DT
» (f(K,J,l.l) » 3YKI/OX
J.I ,7) + OYYl/OY
( ,,33 lUKiJiI'!'1
I '21 '
( ; 3 S ! 'J ! < . J . I i 3 ) )
T i' '- E
43 K •> I . NZ
DO SO 1 * 1 > N!t
•J ( K , 1 , 1 . I ! • U < K • 1 • I
U ( ' , 1 , I • 2 ) = V, ( K ' 1 . I
0 I R , i i *• . 3 ) = U ( K i '• i J
IF
IF
Ir
CO
CO
J!S;JI:5:;:5:!:5!:
,I I/OX
i2'/OT
U (>C , I- t i I i 2 )
o-*,'.?«: -3) • Jf-.NY
oo so J ° 1. NY
50 C!K,J,!> • 0.
00 60 J » 1. NT
y i K , J, '• . 1' • u U t J11
^(t,J.l)2I « UlKiJ.l
•J ( " , J . 1 . 3 1 " U ( ',. J i I
U ( v> , J , M X • I ) « U I " i J i
U!X,v!,'X,2) » U«iJ.
U(K,J, -X,3) » UK.J,
60 CONTINUE
•5ETUPN
E'O
u (x|N*,i.2)/DY
Jtn.Ny. i >3)/OZ
,1I/DX
,2'/DY
,3)/OZ
.xti)/ox
OIF05200
OIF0530C
OIF05-400
OJF05500
OIFC56CO
01F057CO
OiF05SCO
0 I F 0 6 0 CIO
D1FC6IC3
OIF0620C
01F06M30
OIF06500
01 F06600
01F06700
OIFC60CO
OIFL'6?00
OIFC70CO
01F07100
OIF07200
OIF073CO
01F07MOO
D1F07jCO
DIF076CO
01FO/QOC
OJF07900
D1FOUOOO
01F03100
DIFCB2CO
DIFC83CO
01FC3«00
otFcasoo
D1FOU(,CO
OIF08700
01F08000
OIF08900
01F07000
O-l
en
CO
-------
OIFFL-S
FO"FLO /X808/ ON 13 OCT 71 AT 07:56:16
SUBROUTINE OIFFl'SlN)
C THIS POUTINE CALCULATES THE
CDIFFUSION PSEUOO VELOCITY
C01F30Z03
DX,DY,OZ,TDX,TDY,TDZ,PMIN,PMAX,TfN,
B/lCltiOT, PS'.STP (5) , PPXSTRIS), «NSTOR(S)j
OC'CST'MSji NP3TORI5), FACTtS>, F"X, FVY,
r>Z, v, rv, UP, ;, ICYC, ARRAY(S0f3&),
",21 ,5)
,S>, 3X2(1i2Z.5),SYH1,2)i5),BY2(
NX ,NY ,NSPiKAX>NXMl ,NY"H >H2M1 ,NA,
,NAY,N/ Z,MAif ,NZ, NA1, NAZ, NA3,
OJ ".ENS: ON C ATPAYCt t"37 ) )
INTEGER DJ!"MY
COMMON /OUT/ IS°(S), tCTC, NOLi TOELT,
PTllE, CONC. PAST, SF(S), NPSP, RSTRT,
SAVE
INTEGER SP, SAVE
LOGICAL CONC, PART, RSTRT
DATA SX.SY
t:«0!F02DOO
C-OIF07SCO
[.•DIFC2600
C SX, SY, SZ ARE THE LARGEST 3
[FRACTION OF A CELL IN THE X, Y« 3
CZ OIRECOIF02200 3
c RESPECTIVELY, THAT A PARCEL 3
BE ALLO-TED TO MOVE IN ONE 3
STOIF023CO 3
00
»
I
CO
-------
/ READ (8) («J( 1 ,J,I ill .1-1 .NX! iJ-I ,NT) . / '/
/ ((U(I,J,I,2),1«1.NX), JM.NYlt «U(2, /
/ j>: .31, i = i ,HXI ,j»i ,NYj, ov, OH /
/....,.,.,.... ,,,.......,.,...,...../
I
I
,.,.,,, ...,,,,,,,,,,,t,,,.,.,,,,,,,,,,,,, ,x FALSE
<: IT (OT .LT« 0') > ———
\, ,.,.,...,,..,,, ..,......../
I TRUE
I
/"">"•""""" «N
: RETURN :
N... /
0<------------------------•----.
I
..............I.........................t.... I
A--- — -->: - oo s i • i • NX :
A I..,........,...,..,.,.........,........,..,:
A I
A I
A
4 ., . _>; DO 5 J » 11 NY :
A B ! ;«•• ...I
A S t
A 8 I
A 3
A 3 C -->: 00 2 L « I, 2 !
ADC : :
A 3 C 1
A 0 C 1
A 3 C ,......,, . . . .
A S C UUiJitiU « 0(1 ,J.|iLJ»3T
A S C U(2.J,I,L) <= UtliJiltL)
ABC UOtJtitD'-UdtJtliU)
A 5 C ;....,,
A 3 C ' 1
A S C C 23 J
A S C •••• •
A E ...—; U(H,J,!,t_) ' UOfJiI.L)
A 3 :
A 3 1
A 3 I
A 3
A3 0(1,J,I,31 * 0-
A 0 U(2,J.I,3) » 2.«DT«U(2,J,I,3)
A3 U(3,J,I,3)-U(2>JiIi3)
AS UO,J,1 ,3) - 0.
A 3 • •
A S I
A3 C 53 1
A 3 •• •
. .. —; CONTINUE :
i
i-.-c ...3
I C •»••« 3
I C'D]FO"20C 3
I c «••.»•«•••• ;,
I C-»« »••• ••••« 3
I C«OIFOU300 ]
en
oo
-------
w
w
on
A
A
A
A
A
*\
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
*
A
„
A
^
A
A
A
*
3
3
6
3
P
3
0
3
D
f
n
0
3
3
3
3
C
2
«
3
D
B
a
3
a
c
o
2
G
n
S
2
3
5
3
8
3
B
3
C
0
8
3
S
3
a
a
3
3
3
3
0
R
£
3
3
S
C
C
c
c
£
c
c
c
c
c
r
c
c
<-
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
r
c
c
c
c
r
c
^
c
;
c
c
c
00 SO I « 2f NXKJ
7
1
1
I
• r\n ^ t » » •
i
i
I FALSE
I
I
U(K,J,!,lJ •> UIKrJtl >1 >/DX :
UU,J, I ,2) • im»J,t .ZJ/OY :
1
I
C 83 I
CONTINUE
CVX o DH(KtJ - C(KfJ,l*m
/(CtK t J» J J»TOX)
OVY . OH<«:,J, : i.tcir., j-i ,i> - c
-------
A 3 C
A 5 C
A 3 C
A 3 C
t 5 C
A 3 C
» = C
A 2 C
A 3 C
A 2 C
A 3 C
A 3 C
A 3 C
A 9 C
A 5 C
A 8 C
A 3 C
A E C
A 3 C
A 3 C
A 3 C
A S C
A 3 C
A ; c
A 3 C
A 3
A
A 3 C
A 3 C
A t
A 5 C
A 3 C
A
A
A 3
A 3
A E C
A 3 C
ABC
A 3 C
A 3
A S
A B
e c
A
A 3 C
A a c
A 2 C
A 2 C
A 3 C
A 3 C
A 3 C
A 3 C
A 3 C
A 3 C
A 3 C
» 2 C
A 3 C
A 3 C
A 3 C
A 3 C
A 3 C
. 2 C
i n '
c
: ovz
: i ,J
i
: ovz
: /! c
c
:
: U :
i
3C3 1
CONTI.\UE I
t
i — c CONVERSION FROM M/SEC TO 3
I CF"ACTIO'. OF CEH.L/OT OIF06300 3
I
I
.J.1,1) "
-------
w
A S C
A 3 C
A 3 C
ABC
A S C
k 0 C
A
A
A
A
A
» 3'
A 3
A 3
A 3
A 3
A 3
A 3
A 3
A 3
A 9
\ 1
I D
.. D
„ 0
A 3
ADC-
ABC
A 3 C
ABC
A 3 C
A ----
A
A
A
8
A
A
A 3
A 3
A B
A 6
A 3
A S
A B
t. 3
A 3
A 3
» 3
A 5
uUiJ.1,3) • si6msz,u{t,j,i ,3> j
103
I
0<-
I
CONTINUE
•>:
00 60 K • i, NZ
-> •
00 50 1 » It NX
IHK.l.i.l) - UU, 1 , ! i 1 )/OX
<-'Uilt!tZ) <• U(K>lit>2)/DY
U(Kt 1 i t ,3) - UUiI.I «3)/Or
Ui'Ci'lY.I , I ) » uUfNY,! ,1 1/DX
V U i N Y , I , 2 ) •> UUiN'YiIf2)/DY
•>:
DO SO J
i , NT
503
CU i Ji 1 I • 0.
•>:
DO 60 J • 1•
V(K.,J,l,l) . (J(K,J, I
U ( K i J > 1 , 2 ) o UU»Jili2)/DY
U(X.J,1 , 3 > « U(KiJ.li3)/OZ
U I K , J , N X , 1 ) » U ( K i J i N!'
U ! K , J , N X , 2 ) = U U i J i N i
t'UiJiNX,3> » U(K>J|NX,3)/OZ
603
CONTINUE
/
en
CO
-------
'>FLO,S TNPUT
ST FoRFLO /XOQO/ ON 13 OCT 71 AT 09554:55
oo
C30CO 1 •
C33332 •
C33333 •
C3333" •
C3333B «
C3.3334 •
CC33-37-
C 3 - C - E •
3' C7 •
c 3 : 3 '. I '
C333 i 2 •
C333 I 3.
C? 3 C » 4 .
C3C3* 5 ,
333. ' 6.
C3-. i 7.
33:3.3»
C333.9-
333321.
C 3 "0 22 •
C 3 '- 3 C 3 .
333^7J »
3333J3,
C33335 ,
C333J 7 .
C 3 3 3 2 S .
C-33r 2' •
333333.
c
c
c
c
c
c
c
f>
c
c
c
c
c
c
c
c
c
c
SUBROUTINE INPUT
THIS ROUTINE SETS UP CONSTANTS AND ARRAYS NECESSARY FOR BOTH
CO-PJTI'G AND INPUT/OUTPUT. DEFAULT VALUES ARE SUPPLIED FOR
ALL NECESSA1?" PA^A1ETE"5. VARIABLES TO BE READ INTO NAMELIST
iN"i;7 A-iE A3 FOLLOWS.
,3=" « \A'ES OF POLLUTANT SPEC]E3 TO BE PRINTED
f-CTC » FPE'vUCVCY CF P'tlNT CYCLES (DEFAULT » i)
NOL = UUf-ER Oc LEVELS ( I >1 THE z DIRECTION) To BE PRINTED
(DErAL'LT 3 MI
TDELT = T1"E INTERVAL FOR PRINTING (DEFAULT » 0.)
CO\C " FLAG FOR DUMPING CONCENTRATIONS ON UNjT'2
iDE-AULT o .FALSE.)
PA'T = FLAG FSi? OU"i=i,\G PA°*lCLE POSITIONS ON UNIT 28
COC^AULT » .FALSE.)
NPSP * Nl/fBE" OF SPECIES TO BE PRINTED
(CEFALLT . NSP)
RSTKT - RESTART FLAG (DEFAULT « .FALSE.)
SAVE • CVC'.E 70 BC DoKPEO CN UM 7 29 FOR RESTARTING
(DEFAULT » - I )
COI',"CN DX,OY,DZ,TDr,TOY,TOZ,Pl!N,F.u.AX,»N,3ACK,DT,
jPTSTPlS), P^X3TS(S), ''VSTO"(5), 3CX3TRI5), N'PSTO')(S),
2rAC7(5), CNX, 7\", F(!7, V, (?V, NP, 7, ICYC, . <* " A Y ( SC936 ) ,
jFLvM (22' ,«-Ll/n2(22) ,FLJYl(21),FLUY2(2I), 3X1 (".22, 51, 3x2(1,22,5)
"5Y1 (S ,21 .5) ,3*2"! ,21,5),
S% x , t. Y , NSP , SAX , Nx,i i ,N»^ i ,NZMI ,NA,NTOT,NAX,NAY,NAZ,NA*',NZ,
6 -A!, i.A2, ^-AJ, NAU, D'Jf""l
OI"E 3IC' C!M, 21,72), CC(H,21,22), 0 V ( 1 , Z 1 , 22 ) , OH(H,21,22),
1 U(« ,21 ,22,31 , S(",21,27I, X(9503), Y(9500), 2(9500), «T(9SOO)
E"j;v/LEvCE (C,A^RAY), (CC,ARPAY(18M9)), (OV,ARRAY(3697)I,
1NP00100
INP002CO
INP033.--0
INP2CMC3
INPCOS30
KIP00630
INPC07CO
1NPC3SOO
I NP0090C
I N D 0 1 0 3 0
INPC1 ICO
INP012CC
ISP01300
1NF31 "CO
INP01SOO
INPOUOO
INPD1 700
1 N P C I 3 C 0
INPOI703
INPC2100
INPC2230
IN°023CO
, INP02100
l'.P02503
IN"02tOO
INP32700
INF02030
U-PC29CO
1N"033CO
C333:!!. 1 IS,/13'. Y ( 3697) ) , tljH.AR^AYlSStS) ) , ( U • ARIA Y { 7393 > ) , 1 X , ARRAY ( I 2937 I FCJ'CO
333332.
333333.
33333J.
CC333 7 .
3 3 3 " "* ' *
C "* 3" 3 r »
C 3 3 3 • 3 .
C t 3 3 *• I •
-I-:!'*'
C;:^s*.
CC3C-5.
33 33 J 4 .
::3:-.s-
3:::-^.
»» „ - j 5 „ •
c
c
c
c
2)), I '. A83/."! 22H37) ) , ( Z, A^RAT 131 737) ) , ( t»T , ARS A T ( 1 1 u 3 7 ) )
if IEGER OL""".Y
CO-,',0 /3UT/ ISPI5), NCYC, NOL, TDELT, PTIME, CONC, PART,
1 SP(5), NPSP, °STRT, SAVE
I..UCE3 SP, SAVE
L C .', ; C A L C3'C, PA'JT, RSTf?T
NfELiST /II'PUT/ is?, NCYC, NOL, TDELT, CONC, PART, NPSP,
I P'T'-T, SA^E
Dfl'A IS" / 1 , ?, 3, <), 5 /
D,,7A r,C YC.NCL.NPSP, SAVE, CONC, P ART, RSTRT/1 , H , S , - 1 ,3«. FALSE./
"V = I .23£-6«OK»3T»or
T3X » 2. 'TX
7DY » 2,'CY
"or = 2.«DZ
F , y «• Jx
F 1, » « < Y
r-, •• c .'z
"C'D (5, INPUT)
"FMTt. 16, I'JPuT)
! ,v P 0 3 2 0 C
INP03300
INP03MC3
1 NP035CO
Ifj°03600
I VPC3703
I N P 0 3 £ 3 0
1KP03900
IMP01C3C
INP3M100
IN"Ct2CO
INPOHS03
If'OtbOO
I'fPOf 600
I'|POH703
I NPOS 703
INPOSCOO
ItiPOSlOO
CO
^
i
oo
-------
CC3C53*
IFtRSTPT) CALL RESTRT
RETU
END
te
INP05200
JNPOS300
INPOSMOO
CO
CO
-pi
-------
IsPUT
FLOfiCHARTED BY FOPFLO /X60S/ ON 13 OCT 7l AT
INP'JT
CC
8A
I
I C THIS POUTJNE SETS UP CONSTANTS 3
I CA*D ARPAYS NECESSAPY FOR BOTH 3
I CINPC0200 3
i c COMPUTING AND IN°UT/OUTPUT. 3
i COCFAULT VALUES A-JE SUPPLIED FOR 3
I C I JFCC3CC 3
i c ALL 'cCESSA"?" PA°AKETERS. 3
I CVARUDLES TO BE rEAO INTO 3
I C;STS!5), NPSTCJitS), FACT(5>> FNX, FNY,
'Z, /, RV, s'p, T, ICVC, ARRAY (53936) ,
v.'Xt(72),FLUX?!22),FLUYl(2i),FLUY2(2!),
It",22,5), 3X2(".22 5),QY,(', , 2 1 . 5 ) , 3 Y Z I
21,'i'1) • X , Jv,>S',f'AX,',)"'l,Nrv'liN'Z"l.NA,
CT, AX,NAY,"AZ,\/.*(,NZ, ^AI, NA2, NA3,
-E'.S'O.. C(1,2|,22), CCO,21 ,27) , OV(S,
.22', D«M«,2I ,Z2> t u(M,2l ,22,3) , S(H,
)) , Z(?50C) , WT(
OJ
CO
oo
-------
,..,-• ..,.,.,...,.,. . v ,«MA* • > , < VV •-••<••• I •->/> I ,
')), INPUT) INP05100 3
I
I
< IF1R3TRT) > .......j
\ / I
! TPUE
I
: : CALL RESTRT ' !
0^*w«~*>*w«>«v>**
-------
e"-C7«FLO.FORrLO,S
F_cf,C"A"7EO ST
/*ecs/ ON 13 OCT 7i AT o?:37:oH
NJ
OC2C15.
Cc;r!6.
C::C17.
c::c;s.
CCC02S.
C::=2».
CCCC34-
c::c37.
Cc:c--2.
C:C3'-3.
KEM KEM00100
C 7h,s POJ7JNE DETERMINES WETHE" THE CONCENTRATION OF POLLUTANTS IN EACH
c CELL is SUFFICIENTLY HIGH TO WAVE A SIGNIFICANT PHOTOCHEMICAL KE«ccrco
C EF-£C". A'.O, .F SO, CALLS C"£M- KClOaHOO
CO--CS BX.OY^Z.TCX^DY.TOZ ,PKiNiPMAX,YiN,3ACK,DT, KEfOOSOO
IPX 3TPIS), PSXSTIIS), *NSTOP(S), BCKSTRtS), NPSTOR(S), r.EMC0600
2FACT15), F*X, F'lT, F."Zt V, PV, NP, T, ICYC, APRA Y ( S0?34 ) , XE''OC700
SFLUXI (zz) .FLUXZC'Z) ,FLJYI i OXMH.ZZ.S), ex2tf|22i5) .KEMCOSOC
21 ,si , KE
.^YKI ,;JZMI ,NA,NTOTINAX,NAY,NAZ,NA«INZ, KC.
6.'Ai
IA?,
HA".
KECOI200
CCt,21,72), CCd.ZI ,22) , OV(H,?| ,22) , OH ( M , Z I , 22 ) ,
I '-! " , 21 ,27, 3) , S(ui71j?Z)i X(?300), Y» «7!?500)
£0'J!VA'-:,JCC (C.A'PAYl, (CC .A^RAY! 18M9) ) , ( D V , Af° A Y ( 3 6 ?7 ! ) , KEMOi»00
1 (Si A7RAYI3677) ) i I OM , AR«i AY < SStS) ) , < U , ARRAY ( 7393 I ) , ( X i ARR A Y ( I 2?37K Er 0 I SCO
2M. (Y > A"!0AY:22H37 ) ) , I Z , ARRAY ( 3 I 937 ) ) i ( WT , ARRAY ( 1 J 137 )) KEfOltOO
INTE&E" OU11Y KE101700
...t«. ...... ........... ...... ...... ............... ..... •••*••••••*. xEKCidoo
PA"t*E"E'< WAS « N'A«NSP»! I
ESUIV/LENCE IcF.A^RAYJi ( CCF , A^RAY ( 'IAS > )
Oi^£NSio-- CFO ,2: ,22,5) , ccFts.21 ,22 15>
£:„! V^LE.-CE (CF.ARFAYJ. (CCF^PRAYITZI i j i
DATA EPSI/UE-S/
/
00 1 -I =• 1, N3?
i' » 2 C » N
READ (M) ( t (CF«i Ji I fN) ,K«1 ,NZl ,J=1 ,NY) , I
°E> I .0 1
1 CO'T.'OE
oc ;c I = 1, NX
CC iC J • 1, SY
00 1C t » t, NZ
SUM » o.
DO 7 M « i > NSP
,NX)
iF ISU' «ST. EPSU GO To ?
oo s K o i . us?
e ccFf , j, i ,.\i , i.
GO TO '0
9 CONTIVUI
CALL CHEKtCF(K,J,I, I ) ,CF(K, J , I , 2 ) ,CF ( K , J, I ,3) ,CF(K,J,I,1) ,
i CF«,J,I,5), CCF(K,J, I , I ) ,CCF(K,J, I ,21 ,CCF(K,J,I ,31 ,
2 CC'i'i Jt i ,1) .CCFK.J.I ,5) ,oT,T|
13 caiTiNUE
DO 20 >• = 1, USP
KE"C:COO
KEf-o?ioo
KEf02«00
KEr02&00
KE>027CO
KES02800
KEHC320Q
KEV.03&00
KErcmos
Kcnoaeoo
KE1C3700
-Z) i
U(CF(K,J,I ,1)
«l ,CY! i I«] ,NX)
I ,NT)
,NX) ,
I (CCF(K,J,I ,
i o
20
100
-------
END
CO
KEHOS200
O-J
CO
i
CO
-------
FLOflCHARTED BY FORFLO /X80S/ ON 13 OCT 71 AT 0?:S7.:C6
SUBROUTINE
C THIS ROUTINE DETERMINES TIHETHSR3
CTHE CONCENTRATION OF POLLUTANTS 3
CIN-;E1C02CO 3
c CULL is SJFFICIENTLY K:GH TO 3
CHA'.E A SIGNIFICANT PHOTOCKE11 C AL 3
C>-^t':C300 3
C EFFECT, AND, IF SO, CALLS CHEM.3
3
BACK.OT, PMNSTR(S), P!1 X 5 T H ( 5 ) , .'MSTOR(S),
SC'~ST?(S), NPSTO->t5), FACT<3>, FNX, FNY,
rv-Z, V, °,j, tp, T. JCYCi AR^Ar (SC'36) ,
(22), FLUX2(22),FLJYl(ZI),FLUy3!21),
H , 2 ' , 5 ) , t» X , N * , «, S P , K A X , •'• * M I . " Y K | , N Z M 1 , N A ,
NTCT,NAX,\AY,N,,Z,NA;, ,(Z, K A 1 , [IA2, NA3,
DIMENSION C(l,21,22), CC<1 .2! ,27) , OV(M,
21,22), Oh("i?l >22) • 00,21,22,3), SIM,
21-22), X(4S031, YfVSOO), ZCPSCO), *T<
EQUIVALENCE (C.ASPAY1, (CC), (WT,A^-'AV(M1''37)!
IN7FSER OU'-f'Y
•--C » 3
r.,,.,...,....., •••..*.....]
E PASAI'ETER NAS " NA«NSP»I 3
E<£M01?00 3
C EQUIVALENCE (CF,ARRAY), (CCF, 3
CARR«Y(NAS)) KEMC2COC 3
I
CF(1,21 ,22,5) , CCF<«t,21 ,22,5)
;E (CF^R^AY), (CCF,APRAYI?2M1 ))
DATA EPS1/1.E-3/
I
t
A------------—_->. oo i ,'. » ;, MSP
A :
A 1
CO
»
i
00
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
e
A
A
A
A
r
i
f.
i>
3
3
n
9
3
P
j
•}
3
5
^
R
3
O
3
O
^
O
3
q
j
0
c
c
c
c
c
c
c
C 0
C 0
C 0
c
c
c
c
c
c
C D
C D
C C
^
r
: K • 20 * N ; * /
v;
i
/ READ (K) ( ( (CF(K| J.: .N) iK"l >M) iJ"l ,NT/
I
J
/ REWIND P /
C 13 I
*
1
I
1
1
i
I
: SUM » o. :
i
i
i
C 73 I
I
I
I FALSE
1
1
C 83 !
1
•
•
%
•
C/}
i
Z -P-
z
Z
2
z
z
z
J
jr
Z
7
-------
A
A
A
A
A
A
A
A
i
A
A
A
A
A
H
A
A
A
A
A
A
A
A
V
3
3
Z
3
E
3
a
e
c
2
3
8
3
8
e
5
5
s
2
3
g
i»
^
C
;
C
c
c
c
c
c
c
r
c
r
C
c
c
c
c
c
c
c
c
c
C 73 I
: . CONTINUE :
i
i
: CAUL CHENtcf iic.j.if n iCrtx,j,i,2» ,CF{
: K,J,I ,3) ,CF(K,j, t >*)) , C^U , J. I ,5) ,
: CCMK .J, I , 1 J ,CCr(K . Jr I .2) rCCTIKt J, J ,
: 3), CCF
I
C 103 1
----_--: CONTIWJ: :
................... — ,0
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
Z
Z
Z
Z
I
A
A
A
A
A
A
f.
A
A
A
A
A
A
I
I
I
1
: M » 20 * N :
i
i
/ nxnt (H) I UCF(K> Jt I iN> >X»1 .fvZ) t J"l •
/ NY),l«l,NX)f ( ( ( CCF(K i J, 1 iN) , K>lfNZ)f /
I
/ REWIND M /
c 20: i
. — . — • CONTINUE :
«/
C '.003
OJ
CO
-pa.
-------
or FflTLO /SBC'S/ ON1 13 OCT 71 AT 0?:S7:16
occ;;: •
C3CCC2.
CD;c33-
cic^C" •
3;":cs«
Co-";..
OC'CC 7 •
r-:~- - ,
SUBROUTINE NPAPT(N)
C Tu;s PCUTI'.'E 3ETES»!1NES WHETHER NEW PARCELS HAVE BEEN CREATED BY
c A.DI if so. PUTS THEM INTO THE CENTER
COI^O'I DX,DY,OZ,TOX,TOY,TOZ,PMlN,PMAX,
jP>" 3T°I5> , "(XSTSIS), >.:J370R(5> , EC^ST
2F,,CT:51, F^X, FN", FNZ, V, RV, HP, T,
3'\. J7\ 121! ,Fl_uy2(2?) .FLUY1 (2, ) ,FLUY2(21
fSrif'S'M'SJ.e^c'.^i.s),
OF THE APPROPRIATE CELL.
WN.BACK.OT,
" t 5) , KP3701? (5) ,
1CYC. A015AY t 50?36 1 ,
', EX1 (M ,22.51 , BX21M, 22rS)
|
MPT00100
CIIEM
NPT03:CO
NPTOCMOO
NPT00533
NP*CC600
,NPTC07C3
NPIOOSOO
CO:rC9. 5 X ,NY , •'SP.'IAX ,NXri ,I,'Y,',1 ,NZM1 tNA ,NTOT ,NAX ,NAY ,NAZ .NAW.NZ, NPT00?30
c:rc 3.
c?:r' ..
or:;1 2.
CCO- i j.
C"G 1 " <
c : '."> i = •
c::?u.
c::oio.
CC ^0 1 7 >
0:0323.
CC -22 1 ,
Crr?22 •
- : -T 2 z •
c::~:-<.
C _ j _ ' j «
::;:f4'
C3CD2S.
C CC 7 2 7 «
c:c::i.
cc':J2.
C 3 0 C -1 3 •
C jC;3"-.
ANA', 'JA2, NA3, NAU, OU"I.;Y
0 "E'iSiOf. C(«,21,2?), CC(1,21 ,2?) , DV(1, 21,22), OH|S,21,22)«
I UM<21,2;,3), SM, 21,22), X(?SCO)P Yt
E^'JIV^tEsC: (CiAHcl A^ 'AYISS'tS) ) , (U,
2)!i I YiARPAY(?2H37) ) , ( Z i A1PAY<31937 ) }
UITESCT 3L' ^Y
XX = -1,/FACTtN)
00 : 3 ! = I . NX
DO ic J = 1 1 vr
DO 10 K«l ,NJ
ir tec:1? , J, 1 1 1 s.i.io
« CCf TIMt'E
:=• (ciK,j,i) ,E3, o.) cc;, j. ! ) = i .
X ( ; => ) •> ;
Y t • o \ a 1
1 t ' I / 3 w
Z > N 9 ) » X
10 co- TI.'.UE
R £ T I, "',
E«D
'SCO). Z(9500), *T(9SOO)
1 i (0^ .ARRAY ( 3697 ) ) ,
tIPTOlOCO
NPT01 IOC
NPT01 200
UPTO. 3CO
ARRAY(73?3>), (X, ARRAY I12937NP7DM03
, (WT.AR^AT IM1M37 1)
tJP'O 1 b30
NPT3 1 A CO
N°To:eoc
NPTC17CO
NPIC2030
NPT02 ICO
NPI02200
NP702303
NPT02"33
NPIC25BO
N»T32403
NP~02t)00
NP702930
NPT03QCO
WP70J200
NPT033CO
NPT03-J03
en
oo
-------
BY FORFLO /xscs/ ON ID OCT 7t AT 09:57:17
63
.p»
00
C THIS ROUTINE DETERMINES "METHER3
CNErt PARCELS HAVE BEEN CREATED 8Y3
CCNPTCC20C 3
C AND, IF SO, PUTS THEM INTO THE 3
CC-'.TCP OF THE APPROPRIATE CELL. 3
CNPTC0300 3
: CO"«ON CX,DY,CZ,TOX,TOY,7DZ,P1IH,PKAX,tN.
3ACr.,OT, P.u'.37PC5!t P^XS7R(3), wrSTCRCSIt
OC^STFIIS), '.PSTO1?! 5) , FACTPAY(50936),
FLUM(22>,FLU X2(22>.FLU Yl(21),FLl!Y 2(21),
3X.(U,22,5)) 5X2M,22,5>,3T1CU,2I,!>),BT2(
1,2' ,5.) , 'JX,!,'" tNSP ,!• AX, NX MI ,NYM1 , NZMI ,NA ,
NTC7 , '..' X ,VAY ,"AZ ,t Aft ,NZ , NA1, NA2, NA3,
^ Af, 3tMfT
Oi.'.E! SION C(",2l,22), CC(",21 ,271 , OVC«,
2',22), OH(" ,21 ,22) , U(1,21,22 ,3) , S11 ,
21,221, X(9300)i YC9303), Z(°300)i HT<
75CC)
E:JiVtLENCE (CtAR^AY)t (CCtAKRAYIIOM9)) ,
COV,A'?PAY|3697))I (S,AKPAY136971), (DH,
ARRAY(7393)1> (X,APRAY(
YJ22M37)), CZ,ARRAY<
317371), ( V, T , I
A
A
A
A
A
A
A
A
3
a
5
a
3
3
3
C
C
C
C •••••••••
c
C'NP7C17CO
XX"-J./FACT(N>
00 10 I • 1, NX
oo 10 J = i, NY
oc 10 K = I,NZ
v
i
CO
?a
00
-------
to
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
rt
A
£.
A
A
J
3
5
r
c
0
J
B
?
n
a
a
-
3
8
s
3
0
o
?
B
2
C
^
C
r-
c
c
r
C
*»
C
r
*»
C
^
/•
c
c
£
/•
c
c
C 13
: CONTINUE :
i
i
I TRUE
I
: cc: it
i
I
C S3 I
Np o NP * I !
WT(NP) » (CC(«»Ji!> * 8AC»:)«XX
CCtKi Ji 1 ) » 1 . !
X(M?) » i :
Y ( 'i P i a j :
Z = K :
i
c 10: :
:
: RETU°N :
\..... /
X
*
X
x
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
x
X
.
L
2
z
2
J
2
z
z
z
2
z
z
2
z
J
z
z
z
z
CO
"73
CO
-------
SFC»«FLD«rO?ri.O,S OUTPUT
FLOrCM-STED 6Y Fo°FLO /X308/ ON 13 OCT 71 AT 0?.'S7:30
c 3 : c j 7 .
Cd
On
O
C C - » I 2 *
c : : c 3 3 •
OJTPUT(K)
C THIS ROUTINE PRINTS CELLULAR C^NCCNTRAT J ONS, IN PP«»
C AND vE'GHTS ON L-1T 28 IF TME APPROPRIATE FLAGS ARE SET.
REQUESTED
AND v E
C 01M C • C X
J"»H13TR'5)
2FAC7,'5J, FNX, FNY, f
3FLUXI(22),FLUn2(Z2),
« 3 Y J t 1 , 2 '. ,5) , a Y 2 ( " , 2 1
4- t ; , * A 2 , "A3 , f.AI ,
Bii'E'is IC'i CI M ,21 ,22)
V.NSTORI5), OC ) , DU7DMOO
1 (S,ARPAY(3697! ) > (DMiArtRAY1S515) 1 ,(U,ARRAYI 7393 ) ) , (X.ARRAYI 127370UTCIiOO
2)>i C-,A = FAY(22137) ) , (Z,ARRAY(31937) ) , (AT,ASRAY(1 I 137) ) OUT01600
1"TEG;R Ol'N.'Y OU701700
1
10
20
3C
35
"0
CC","ON /OUT/ ISPI5), NC^C, MOL, TOELT, PTir"E, CONC, PAPT,
s->(5), N?SP, PSTPT, SAVE
L.-EGE* s°, SAVE
LOGICAL CO C, PA'JT, RST«T
OI-E. SIC CLIM22) , IT<2)
LOGICAL ,tOT,
D\TA Sp / 'MO1, 'N02', 'HC', '03', 'HM02* /
,•1 •= A 3 S ( M 1
- "(>»• 57. ICO GO TO 35
IF SJ,!) « CU, J,I)«FACT(N) » BACX
IF <•' • E 0 • 1) ICYC = 1 C v C * 1 •
•or? = . T»(.E»
IF ( IC"C .£•!. 1 ) GO TO "0
'.F I^CYC .»E. 0 .A'JO. M00( ICYCiNCYC) «EO. 0) NOPR « tFAUSE,
IF (T -57. PTjME) GO TO 20
PIleE " P T 1 -. C * TOEuT
JCP" = .FALSE.
• F ( -0°R) GO TO ICO
00 30 I = 1 , MSP
iF ( 'S-'l 1 ) .EC. M GO TO HO
CO ,7 ' J;
C2 TO 'CO
•• = . - : c o
:TT = T » .5
! T t 1 ' = i TT/3600
I T (?) « f OCI .TT/40. 401
OI.TC1900
CU702030
OU7C2IOO
OUT C2«!QO
OUT02HOO
OUT02SOO
OU7024CC
OU7027CO
OUTC2800
O'J 702900
• OUTC3QGO
OU703 1 00
CUT03200
OJT03300
OUTC3MCO
OU703SCO
OU7C36CO
CU'03700
00703800
OU703900
OUTOHOOO
OUTOH 100
OUTC1200
CU70M303
OU7L'"HCC
OU7d"300
OU7CH6CO
OU701700
uuToiaco
OU701700
OU70SOGO
OUTOi' 00
CO
i
CO
-------
CC2C53.
CC3Ci°-
c:c:'5.
DO 60 K • 1, NOL
"PINT 4Q.1G, SPIN), K, IT,
60!0 rC"'"^T! ' IDENS'TY PLOT OF »
i ICXi 'CTCLE1, I1, 10X, I*
00 50 L « 1, Nv
j « NY - i « i
CO »? I^lt'-X
CLI <'.i i'C(K,j,n
u« ..ME«3I CL I » t I ) "CL I N < I ) • 1 00 •
50 F?!,7 60SO, CLIN
6C53 FORMAT (22F6.M
60 CC.Tt.f'JE
1FI-.3T.1C3) PETUPN
ioc i~ ; j? .EC. o GO TO 123
IF (.,JOT. PART) GO TO 120
NT=2»
*P;TH: INTJ NP, T, cxt t ) ,i«l ,NP>
; ( »T( I > i !•*! ,NS1
IZc'lF (.NOT. CONC) GO TO 130
-*"[T£: <2I T, ICTC, C .
130 P£T^°N "
£NO
A6i ', LEVEU'i 12, IOX.
' PARCELS' //>
I2
12.
OUTOS200
OUTOS300
OUTOSMOO
OUT05500
OUT056CO
OUT01.700
ourosaco
OUT05900
CUTC6000
OUT06IOO
OUT06703
UUTC6303
ouT06uco
out06bao
OUT046SO
OUT06900
OyTO?OCO
OU707100
OUT07200
CO
XJ
I
oo
-------
OUTPUT
FLOKCHARTEO 3Y FORFLO /KZO&f ON 13 OCT 71 AT 07:57:2,1
SUBROUTINE OUTPUT
I
I --- C THIS ROUTINE PRINTS CELLULAR 3
I CCONCENT»AT IONS, [N PPM, AT 3
I CSEllcSTED OUTC023C 3
C INTERVALS, AND QUfPS CONCENTRAT3
C0">'5 0V UNIT 2, PARTICLE
CPOSITIC'^S CUTQ030C
C A"0 ^EIGHTS ON UNIT 23 IF THE
CA^'O^I ATE FLAGS ARE SET.
On
NJ
COMMON PX,OY.DZ,ToX,TOY,TOZ,Ph'IN,PMAX,1N,
SAC<,OT, P«*J3TP(5). P.1XSTRI5), VISTORIS),
3C"3TK(S)> NfSTORCS), FACT(S), F"X, FNY,
F'.'Z'. V, P\ , 'JP, T, ICYC, A"P»Y(50?35) ,
FLUXl(22),!rLUX2(22',FLUYl(2l),FLUY2(21),
tif I M,22,!i) , DX2('(,22,S) ,BV1 <<),2! ,S) ,3Y2(
NTOT ,\AV ,'JAY ,SAZ il.AV ,NZ , NAJ, NA2, HA3,
f.;H, D'JI'MY
OlME-SIOV CO,21,22), CC(1,21 ,27) , OV{1,
21,72), D^C<,21,22), Ul",21,22,3), SO,
21,;2)i XI7SOOI, Yl'SOO), ZC?500), «T (
(0' ,'.?"AV(36
/>?c/iY(55-)3>)
3l?37)!, (1.
tCjA^oAY), ( CC , AR R A Y ( 1 3 M9 ) ) ,
7) ! , (S,APHAY(36'7 ) ) , (OH,
,IU,ARRAYI739J)), II, ASSAY!
'OT) ) , (Z,ARRAY(
INTEGER DUXMY
I
I C
1 C'
I C'
I
I
CO«VON /OUT/ ISP(b), NCYC, NOL, TOELTi
"TIKE, CONC, PART, SP15), NPSP, RSTRT,
SAVE
INTECE1? SP, SAVE
LOGICAL CONC, PAFT, RSTRT
•••3
'••3
3
CO
SO
I
CO
-p-
C«OuT323CO
DATA SP /
LOGICAL
jo'
'Me1, '031
-------
AOS(«)
CM
A
A
A
^
A 8
A S
A 3 C
A 3 C
-. 2 C
* S C
A 3 C
A r C
t /
I
I FALSE
I
I
1 FALSE
I
1
I FALSE
I
1
I
I
I
I
I Z CONVERT KG TO PpM OUT03MOO 3
I
C 103 I
1
1
I TRUE
I
:cvc « zcrc + i
i
: MOP"? = .TRUE. :
i
!
i FALSE
I
i
y
U
U
U
\j
\j
U
U
u
g
U
U
y
U
U
U
U
« U
u
U
U
u
U
• y
y
U
J
U
U
u
U
u
u
u
I)
u
u
IJ
x
A
X
X
X
^
X
X
X
x
x
X
X
X
X
X
X
X
x
x
X
X
x
X
X
X
X
x
X
X
»
X
x
X
X
x
X
X
z
z
z
z
J
z
z
z
z
2
z
z
z
2
2
Z
Z
z
2
Z
Z
2
2
Z
Z
z
z
z
z
z
z
2
Z
Z
2
Z
Z
Z
z
z
z
z
z
z
2
z
I
z
C/J
OO
-------
Cn
A
A
A
A
A
A
„
< 0) NO°R » .FALSE. >
I TRUE
1
0) NO?R o .FALSE- !
i
I
I FALSE
I
1
: pTiMc = PT;«E » TOELT :
; NOPP •= .FALSE. : ' t
i
C 2CJ I
I FALSC
I
----•»>* D030Inlt*»SP •
I
I
I FAI.SE
I
C 301 1
I
I
C 353 I
: N.N-IOO :
i
c sc: :
: :TT « T » .5
: ! T ( 2 ) - S00< ! TT/60i 60 )
S
s
s
s
s
s
s
s
-
.-•
u
u
u
u
u
II
0
y
U
U
y
U
U
u
u
u
y
u
u
u
u
y
U
U
y
U
U
U
y
y
U
U
y
U
U
U
y
U
U
U
X
X
X
X
X
X
X
X
X
f
X
X
X
X
X
X
X
X
X
X
X
X
X
X
J(
X
X
X
X
X
X
X
X
X
X
X
X
X
X
T
7
z
z
7
z
7
z
2
z
z
z
7
z
z
2
z
z
z
z
z
2
,
T
z
z
7
z
z
z
7
7
Z
Z
7
Z
Z
OJ
cn
50
i
oo
-------
->:
DO
Cn
Cn
,,
A / FPJNT 4010, S°(»;, K, IT, IC*C, «P /
A / .............. . ....................... .../
A I
A 1 --- C *0!C FOPMATt MDENStTY PLOT OF '3
A I C A4, ', L£VrLi, 12, 10X, 12, 'S'3
A I C 12, 1 !CX, 'CYCLE', 'H. JOXt 3
A 1 Cl&, ' PARCELS' //> 3
A I
A !
A ...... ............................ .. .........
* - --------- - ---- ->: DO so L • : , NT
A n : ...... ..... ........... ..... ...... ..... ..... :
A 5 1
A 3 I
A B . ....... . .................... . .......... .....
AD : J»NY-L*I :
A o : ........ . ......................... . ........ :
A 3 I
A 3 I
A 3 ........ ...... . .................. ..... ...... ,
A 3 c ------------ ->: oo "? 1 = 1, ax
A 3 c : ........................................ ...:
A S C I
ABC I
A ° C ............................ > ................
A 3 C t CLIM 1 )«CCK, J, 1 ) !
A 3 c • ........................ .. ............. ....:
/ 2 C I
A ; c c t?3 I
A 3 - /..., t ,,,,,,,,,,,,,,.,,,,,,,,,,,,,,,,,,,, ,N
A s ------------- -_< tr
S .................................. • ...... /
I TRUE
I
X
X
X
X
X
X
X
x
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
r.
x
x
X
X
X
X
X
X
X
X
X »
x»
X
X
X
X
X
X
X
X
X
X
K
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
C/)
I
00
-------
Cd
en
ON
\ .../ / t
I
I
C 1003 1
I FA..SE
I
I
N /
I FnLSE
I
I
: ,..7.?a
I
I '
/ >' ^ 1 T £ (NT) N», T, /
/ I=1,N°), (Z( I ) , 1«1 |NP) , (WT( I ) ,1«1 ,MP> /
/ /
I
C 1203 I
I FA'.SE
1
I
t 4.RI7E (21 T, IC^C. C /
1
C 1303 I
X
X
X
-
z
z
z
z
^
z
z
z
z
7
z
z
z
z
z
z
z
z
I
T
z
z
z
z
z
%
O-l
CO
CO
-------
£0 ST Fopc">-0 /X008/
13 OCT 7J AT 09:57:31
w
cocoo: r
"ICCC2.
cc:cc«.
?:::-~*I
- , C ~ 7 *
c c : c c n. .
ItOCC" •
c c ; ; i z '
cccc: : »
c::c ; 2 •
::r_-' 3.
c;:; • •> .
o;:_ 5.
Sl/ = ROUTl»E "AoCELlN)
C THtS POUTJNE CALCULATES AN INTERPOLATED (DIFFUSION ANQ AOVECTION)
C VFuOCiTY FOP EACH PARCEL, AND "HEN MOVES THE PARCEL TO ITS NEW
C POSi7IC\. IF 17 MOVCS OUT OF THE GRID, THE BOEDER CELL FLUXES
c A^r JPD/, TED*
CO-'-C'J DX.OY.OZ.TDX.TOY.TDZ.PKlN.P^AX.rtN.CAC'. ,DT,
1 "• ' 5 T " I 5 ) , Pf^STfl(5}, «NST01(5'» BCKSTR(5), N"STOS(5),
2FACTIS), F^X, F.-Y, FM7, V, "V, NP , 7, 1CYC, ASSA Y ( 50936 ),
3FLU't!!2^),FLU<2(27),FLU^:(2J),!rLUY2(2:), 0X1(1,22,5), 8X2(1,22,5)
"SYIC<,Z1I5)IDY2(H>21I5)I
sr j ^.Y ,N3p,.iAx ,NXM ,NY,I! ,NZKI ,SA,NTOT,NAX,NAY,NAZ,NA",NZ,
6"A1, NA2, N,\3, NA«, DUMMY
OIv£.3iCN C(H,2I,22), CCI",21 ,22) , D/(u,2],22), DH|M,2J.22)f
1 'uC-,21 ,22 ,3) , SIM,?:, 22). XC95CO), Y(95CO), Z(9500), ^T(9SOO)
ESJ.VALE..CE IC.AR^AV), i cc, A«fAY 1 1 819) i , (ov ,APRAY (3697 ) ) ,
PAR001CO
PAP003M
PAR3C100
PACCC600
A P 00 7 C 0
pAKCoaoa
,PAVOO?CO
P A r\ 0 ' 000
PAROI ico
P ;. « a I i c o
PA"0 13D3
PA^OIMOO
PAKOISPO
Cr:?!6. I (3 , « = = AV(3697) ) T ( OH , A = R »' I 55 M 5 ) ) , ( U , ARP ;, Y 1 7 3 * 3 ) ) , t X , APR A Y tl 2937 ° A KO 1 6 30
c:::< 7.
C3331B.
cere* 9 .
orrrz • .
crccz; .
c:-c23-
;^:;-- '
c:.c24.
cc:o23.
CDC;;?.
r~ ~ - •> -« ,
c:::3 1 .
c;:c: 2'
c : r r 3 3 .
c;:: 3".
csrr3S •
c;";*.
;c:337.
3;3c; s •
:ccn39.
c:;o-3.
::::-:.
C7;;U3.
;;; ;- j .
: : : c •• = •
Z~'Z-'t •
:?::'••'.
C3;--'i
2>). ( viA1"*,:Y!2?'j3-> ) , tZ,/,R,'?AY(31937 ) ) , , WT , AR'iA Y I 1 1 137 ) J
IN7EO:K CU'IY
DIMENSION VEU(3)
ir ( ." .El. 0) FETURN
SJO a ^?
>J,5 '
JJ « "(L) - ,S
US •= ZIL) - .5
VT(L> * «T(L)»CC(KK,JJ( i I )
C SUBTRACT PARCELS
i- (.7(L! .LT. PM1N) GO TO 2CO
l * X (L !
J « Y(L)
r « Z
j 1 « I •» 1
Jt = J • 1
K 1 « \ » 1
FSI • X(L) - I
FY1 « TIL! - J
.r (- .r;. C) K • 1
FXe 1 .-FXl
Fv= 1 .-" Y 1
FZ=1 .--Zl
F =FZ-rI *FY
-^"'r!!1"'!
r^-z.'-xirM1
PAKOI700
PA-<03500
PAR0360!3
PAH03/00
PAfQSDOO
PAKC3VOO
PAKO"OCO
PA*\'0U 1 CO
PA-fOIJOO
PAKCIMOO
PASOMSOO
PAf OH6CO
PAf OM7CO
PA-tO"900
CO
CO
F6-FZ.
1«FY
PAKOblCC
-------
CO
F7»FZI*FX1«F«1 PA"?05200
F8=F:i'FX«FYl PAKOS300
IF (L -GT. N0») GO TO 20
A:O PARCELS
lf It'TtLI .LE. PMAX) GO TO 20 PAROS&30
NC" = "T(U)«PP^AX PAKOS7CO
/.-^ « 1 » Ul^iJ'tiMI'Fl » U(KiJiIl,M)*F2 * UtK>JIitltM)*F3 » PARC7CCO
; Vtf, Jl i 1 f.'U TM » U'KJ ij, 1 rM> »FS » U 1 K 1 , J i I I 11) »F4 » PAK07103
2J(I,3)'rl * UlKiJ»It,3)'F2 + UfX,Jltll»3)»F3 » PAR07HOO
1 V/(X i Jl . I |3)'FM « U|K 1 > Jt I ,3) »F5 » UtKI , J,I 1 .3) »F6 » P«K075CO
2 t'lc I , Jl , I 1 ,3 | »F7 * U(Kl , Jl , I ,3) »FO PAKO/400
GO TO 50 PAK07700
HO V£t_(3) - (Z(L) - .5) •.5»(U(2, J, I ,3)'(FUFS) + U ( 2 , J , 11 , 3 ) • ( F2*F4 ) PAK07B33
; » '^12, Jl , > 1 ,3) • !F3*r 7) » UI7, Jl , I ,3) • 1FM + F6) P.\H07vOO
3 » UtI>J. t!i3)•(F1*F3)J PAPOtIOC ,
MOVE PARCEL PAMC02CO
50 XtL) = X(L) » VEL(I) PAROS300
Y(L) " rIL) * V£L(2) PA^OetOO
ZtL) » Z(L> * VEU3I • PARC3500
IF (XtL! ,LT, 1.) GO TO 160 PA- * .5
C/)
CCC'tJJt 1 I ) » CIKKfJJ.Il) » WT'L) . - •- ^u ,
100 ">.>=' « ''»" » 1 PAK09600 QQ
;('•=•) = X(L) PA"097CO 4^
Y(K,P) = TIL) PAKO'.COO *>.
Z ( ' ,° ) •= ZtL) PAKCV70C
ATI . .Pi » »T(|.) PARIOOOO
GO TO 200 PAS101CO
160 Ir tZIU) .LE. 2.5) FLUYl(J) • FLUYltJ) * «T(L) PAP10200
170 1- :ZIL) -Lt. 2.5) FLUY2IJ) • FLUY2IJ) » ?'TtL) PA*1C"03
C: 70 2CO PAKICSOO
180 lc IZ'L) .LE. 2.5) FL'jXJtJ) » FLUXltJ) » WTtL) PA"
-------
0001!Q.
200 t o L * 1
IF II. 'LE» S'P) GO TO 1
IF INP .GE. MAX) PETURN 0
E.NO
W
PAR10700
PA^llCOO
PAR1 1 100
PAP; lies
PAR1 1HOO
C/2
s*~>
\
oo
-------
FLOHCHA3TED 8* FCTLO /X803/ ON ]3 OCT 7)
09:57:35
PARCEL'*)
I
I C THIS ROUTINE CALCJLATES AN 3
J C INTERPOLATED (OI'FL'SION AND 3
I CASVECT1CN) PARC0200 3
I C 'ELOCITY FOR EACH PA1CEL. AND 3
I CTUE" "CVES THE PAPCEL TO ITS NEA3
C POSITION. IF IT MOVES OUT OF 3
CTUE C
I TRUE
RETt'P.l
CO
CO
-------
c
''
•
NOP « NP /
NNP » 0
RPKAX » l./PHA*
L » 1
I
13 J
CONTINUE
It « X(L) » .S
JJ » TIL? * .D
K< - Z (U * tS
WT ZtL)
I I » I » 1
Jl * J + 1
Kl • K * I
rxj " xcu - i
FY1 » TtL) - J
fZ) « Z(U - K
J
I
I TRUE
I
< " 1
[
1
FXol .-FX1
FY»1 .-FT!
f z=i .-rzi
F1='Z»FXTY
F2=FZ«FX!«FT
F3-FZ°F^1 "FT1
F t » F I • F X • F Y 1
F5 = F^! T)"rT
F6-=F? 1 «rX I «FY
F7=Ftl •FXI'FVI
rC-=FZi 'FX»FY:
»
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
x
X
X
X
X
X
X
X
X
X
X
X
x
• ' K
X
I
z
z
z
z
z
z
z
z
z
z
z
z
z
z
2
z
z
z
z
z
z
z
z
z
J
z
z
2
Z
Z
Z
j
Z
Z
z
z
z
z
z
z
z
z
z
z
z
z
7
z
z
7
oo
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
I FALSE
I
i c ADO PARCELS PAROS300 3
I
I
1 FALSE
I
1
' NGPnwT{L)»npPAx I
: WTN = HTCD/INGP * t) :
I
I
I
1
NP « NP + 1
•YT(NP) •> nTN
XCNP) « X(L) « UTANDlbUMMY) »«01
YtS'Pt •> Y(L) » URANS [D'JFKY ) ».01 ,
j
I TRUE
I
ii» « N? - i :
i
C 103 I
I
1
: I/TIL) = WTN :
i
I C CALCUI-ATC THE AVERAGE PARCEL 3
1 COlhECT'^A' 06300 . 3
I
C 203 !
I
. — ...„ — ._...._....„
U
U
U
U
u
u
u
u
u
•J
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
J
u
u
u
u
u
u
u
u
t'
u
u
u
•
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
x
X
X
X
X
X
X
x
x
X
X
X
z
I
z
z
z
z
z
z
z
t
z
z
z
2
2
2
z
z
z
z
z
z
z
2
z
z
z
z
z
z
2
z
z
z
z
z
2
z
z
z
z
I
z
2
2
2
z
z
2
CO
I
CO
-------
tMK, JJ ,1 1 ,« )«F3 • OlKtJl »I fMJ»F4 « U(
Kl tjil t")»F5 + U!<1 , J,I1 ,K)«F6 » U » UIK tJt 1 ,3>*F1 » U(K,J, I ! ,3>»F2 •»
: UltiJl ill ,I)»F3 t U((C,JI il ,3)«F1 * U(
: K: .j,: ,3>»F5 » UIM , j,u ,3'-F6 » uui, •
: *.lrllr3>*r7 • U(M ,JI , 1 ,3) *F8
I
!
C HO: 1
V E L( 3 ) • ( Z ( L ) - .S)».S*(U<2tJiIt3>*(Fl«
FS> * U(2,J, I! ,3)«(F2+FM « U(2,Jl,tli3;
• ( F 3 » F 7 ) » UtriJi,Ii3)*(FS*F8) * Uli.J,
Ii3)»tFl*FSI * U( i t Jt (I i3)*(F2»F6) » Ulli
J1»!I.3)«(F3»F7) » UlliJl»Ii3!«(F1»FO»
I
[ c MOVE PARCEL PAR08200 3
I
C 503 I
X(UI • T (L) + V£L ( I )
Ytt) « T(L) » VEL<2)
Z(L> « ZILJ * VEU(3)
I
I FALSE
I
I
I TALSE
[
i
I FALSE
1
S
s
S
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
c s
c s
5 S
« s
(I 5
-n n «•
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
0
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
1]
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
1
z
z
4.
z
z
z
z
z
z
z
z
z
z
I
I
I
z
z
z
z
z
z
z
z
z
z
z
Z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
f
z
z
z
z
z
T
z
z
7
CO
oo
-PI.
-------
i FALSE ,
1
I
1 FALSE
I
I
1 TRUE
1
: Z(L> » .501 :
i
i
! I n X(L» « -S
JJ = Y(L) * «S
KK n ZCLJ + -5
C(K * CCKfJJ.II » * WT(t>
1
c 103: i '
X(NNP) « X(L)
Y ( N N P > = Y | L )
•• ft T t N • P > " « T ( L )
I
I
C 1AC3 1
i
: FLUYnjj = FLUYKJ) * van.) :
i
I
C I7C3 1
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
M
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
0
Q
Q
0
Q
Q
Q
Q
Q
Q
3
Q
9
0
0
Q
Q
Q
Q
9
8
Q
0
Q
Q
9
Q
Q
<5
Q
Q
e
Q
8
S
0
Q
5
Q
Q
0
Q
Q
G
G
a
q
0
Q
Q
0
0
a
Q
Q
Q
Q
0
f]
S
5
S
S
s
s
s
s
s
s
s
s
s
s
s
s
s
5
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
5
s
s
s
s
s. .
s
s
s
s
5
s
s
s
s
s
s
s
u
u
u
0
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
y
u
u
u
u
u
u
u
y
u
u
u
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
* ,
X
X
X
X
X
X
f.
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
n
t
z
2
z
z
z
t.
I
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z *
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
7
00
00
.£>.
-------
to
o
tn
:
c
'•
c
5
C
:
1 TRUE • / I
I , ' I
FLLY2CJ) o FLUTZ(J) * WT'L) : I
1 I
I
1803 1
I TRUE
1
FLUXt(J) - FLUXMJ) * nT(L) :
1
I
1703 I
1 TRUE
1
FLUX2U) « FLUXZIsM » flT(L) ;
1
2003 I
L - L * 1 :
I
1
I FALSE
I
1
I
: RETURN c •'
N
N
N
S
N
N
N
N
N
N
N
N
N
N
N
.
Q
0
Q
Q
9
Q
Q
Q
Q
••
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
z
z
z
I
z
z
z
z
z
2
z
z
z
2
z
z
z
z
z
I
z
z
z
I
z
z
z
2
z
2
z
z
O4
on
do
-Ca.
-------
CO
»
i
co
-------
-.o.s RESTRT
ST FQ'FLO /xsoas ON 13 OCT 7i AT o?:s7:i?
OCOCC ! •
o c .• " A .
c::::7.
ccc: :.
c c i : ' 2 .
c:::. 3.
ccc;i i.
t C *~ - i / •
L. _ _ „ >- _ »
C C - - 2 3 •
•->-' •(,
CCCC20-
' £"' j , >
''.:-,:>>•
P - -i - ^ - ^
c r ~ n 3 ** ,
C ~ " " *• " •
** - - — u
t •. c u r •
w v - , S ' •
SUIoOjTjNE RESTRT RST00100
C Tb.1: POUT NE 'F.AOS ALL I NFCPrAT 1 ON NECESSARY TO RESTART A RSTC0200
C pPCpLrv F*?Cf UnJ7 27, RS70C3C?
lP* S*rl5)| <"» «7I9SOO) R5T01200
E"'J ". ALC »CC iCiAPR^Y), ICC, ARRAY ( !3H'l )» ( [)V , A CR A Y ( 3697 ) 1 , RiTC1300
1 (5, AT %AV(36"7) I , (C",ARRAY(55.1&) ) , ( U , ARR AT ( 7 3?3 I ) , ( X , AHR A Y I 1 2737 R 5 T 0 11 00
2)1, (Y.APRAYI22H371 ) , ( Z , ARK A Y ( 3 1 9 j7 ) ) , ( ,fT , A"R A Y C 11 1 37 ) ) RS701SIO
INTEC£'J DU-IY RS70U30
C3i«c\ /ou7/ :spis), SCYC, NOL, TCELT, PTIJ-E, CONC, PAKT, ,!"! ,!.P) RS7029CO
S -iF/1,3 '29) C "RS703000
C-L' SCLO'J7(M KS703IOO
f-v- »i -Q R5T032CO
CA' ' D'TrOTf") RiT03300
10 CCI 711 Or RS703MOO
'FA: <2?) 8X1, o*2» SYl, 8Y2 VS703S03
N <= '-S = "2 * (Z.NS° + 1)«JCYC "5T03700
00 1 5 i = 1 , N PS7C 3aOO
JS "E B (8) RS70J90Q
IF (.'07. CONC) GO TO 18 RS7CMCCO
M = ,lt;s.jCYc RS70M100
00 17 I •= I, N PS701ZOO
is IF ,. cr. PAPT? RETURN PS70H'i?o
-r-j^« HSIC4130
E^TRY'fFILE SSTQ^vOS
C 7-5 POL'7,NE DU^PS ALL INFORMATION NECESSARY TO RESTART A PROBLEM RSTObOOO
CO
oo
-------
OCCCS4.
OCCCS7.
CCGCS9.
c;c 367
'/RITE (2?) T, ICTCi PMNSTR, PMXSTR, liNSTOR, BC^STR, NPSTOR
00 2C N « 1. N3P
N? « NP3TOR!H)
M = ? 0 • *J "
PEAOIni c
IF (NP .EO. OJ GO TO 12
«="-IO
FCAO'.") (X( ! ) , I«l ,Np) , (Yl I) i I"l iNp) t (Zt I I , t-l »NP) »
! Ift"(i ) i>»:«NP)
MI7E (2°) (X(I)il"t,NP) , tY(UtI-tiNP)t (2(I),I"t|NP).
1 ( '• T( ! I , i = I |!'D)
12 ' "?' TE 1291 C
2C CO' 7 INUC
r,oiTr (2?, SXJf B)[2l BYI, 8T2
ENBFILE2?
K£ po 2?
RST05200
RST05300
RST03-'00
RSTC550C
KS105600
RSTC5700
RSTC3300
RST06000
RSTCAICO
RSTC6/00
RSTU6MQO
KST06SOO
RST36400
RST06900
RST06700
W
Ov
O3
CO
?0
I
OO
-------
?£STRT
FOTLO /xscs/ ON 12 OCT 71 AT o?:57:so
to
\ SU3°OUTI.SE PESTRT t
\ • /
I
I C THIS ROUTINE PEADS ALL 3
I C If-FORKAT ION NECESSARY TO RESTARTS
I CA 1IS,PKAX,«iN,
3C",ST0(5!, NPSTOnIS). CACT- I ,»ZH1 .NA,
NTCT , N«X , MAY,,"AZ ,HA«, NZ , MAI, NA2, NA3,
i-AJ, O'JI'IY
D;",;Nr!C^ CIS,21,22), CCIS.21,22), OVd,
7;i22), OH(1i?l ,22) , Utt.Zl ,22.3) , SIH,
Zi,ZZl, X(95CC!, Y(?500;, Z(?3C3), «T(
e^l/i'-.iLENCE (C.A'^AY), (CC. A^RAYC 1 8 M ? ) ) ,
AR»AT(55t5)),(U,ARKAY(73?3))> (X.ARPAYC
12737)), ( Y,A-?RA?;22°37) ) , (Z,ARRAY!
INTEGER D'JICY
I
I C» ••« • *.,...3
I C«^5T01700 • 3
I
I
COMMON /OUT/ !SP(S), NCYC, NOL, TOELT,
PTiME, CONC, PART, SPtS), NFS", RSTRT,
SAVE
INTEGER SP, SAVE
LOGICAL CCNC, PAfU,.<*STRT
I
I---C • • 3
1 C ..........................3
I
I
/ PEAD (2?) T, ICYC. P^USTR, PMXSTR, /
/ ''NSTOH, SCKST", NPSTO1? /
i
I
->• oo in N - i, MS" :
OJ
I
OO
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
W A
0 A
A
A
A
A
A
A
A
A
A
I
; FACT(N) . '28.H.V/FACTIN) :
!
1 FALSE
I
1
f READ tz?> n< n • :°! t">?> • tvtn(i»j .NPJ, /
i
C 53 1
/ RSAOIZ91C /
I
J
: . CALL ROLOUTINI : :
i
I
i
: : CALL oUTpuuMMi : :
i
C 103 1
1
/ READ tz?> 3x1, BXZ, ot\, BYZ /
i
i
<2 / REV, IND 29 /
\
I
: N « Nsp«2 » (2»''SP 4. i)«icrc :
i
j
J
z
z
z
z
z
z
z
'
0-J
en
CO
-p»
.
-------
te
A
A
A
A
A
--------- / RCflU l»J /
I
1 F „ L 3 E
I
I
: M.NSP.ICYC :
i
i
C 173 1
I
c IB: !
I TRUE
I
: RETURN :
i
: N.NSP.ICTC :
i
i
i
C 173 1
--------- / RrAO (28) /
1
i
: RETURN :
z
z
z
z
2
Z
z
z
z
1
z
•
CO
oo
1 c THIS "OUTlNE SUS°S ALL 3
I C I'lFO'cAT JCI f'ECEjSARr TO PESTART3
-------
I C ONTO UNIT 27. RSTOS100
I
I
r-o
A
A
A
A
A
A
A
A
A
*
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
/ H°JTE (29) T, JCYC, PMNSTR, PKXSTK, /
/ v.:iS7CK( 3CKSTR, NPST01 /
I
I
i
M»2C»N : *
!
1
I
i FALSE
i
i
: -J..M-.O :
t
/ REAOIK' Ui I ) , I«! ,NP) f (Y( I) ,I»| ,NP) , { /
/ZlIJ.Iol.Nrli (VTIIlil-l.NP) f
I
1
/ • (HI) i 1*1 |N?| , UT( 1 ) i fl iMP) /
I
c 12: i
/ «BlT£ (29) C /
I
C 203 I
I
I
/ '.?IT£ (29) BXl, B*Z» 3Y1, STZ /•
1
I
Z
Z
z
z
z
z
» z
z
i
2
Z
Z
z
z
z
z
,
en
1
-pfc
-------
1
I
REWIND 29
• «**•**•*•<
1
1
04
i
CO
-pi
-------
FLO.S POLL
EY Fo"rLO /X308/ ON 13 OCT 71 AT C?:57:5?
fed
CCC32.
tN) ROLCOlOO
CO-cO>, Or ,OY,OZ,TOX,TOT,TDZ,PMtN|PKAX,«N,SACK,CT, ROLC0200
IffNSTRIS; , PMX5TPIS), ANSTOtlSlt 3CKSTR(5)i NFSTC^IS), ROLC0100
2F«CTI3>, FMX, r-NT, FNZ. V, PV , MP , T, 1CTC, A" P A T I 50934) , ROLCO*t03
3-LUXI !Z7) ,FLUXZ(72) .FLO"1(21) ,FLUTZ(21 ) , DX[(1,22,5), 0X2(1,22,5),ROL00500
5 -'J ,NV , -S" ,^.',v. , NX"1 ,:i»", 1 ,NZ«1 ,NA,NTOT ,NAX ,NAY ,J,'AZ ,NAW|NZ , RCL00700
6NA!, r .';, Sk3, fAI, CuI'fY ROLCCjOC
D,>EN3I3" CM,21,77), CCM,J1 ,22) , DV(" ,2) ,22), OHtt,21.22>, ^OLOQuCO
1 UIH ,21 ,22,3) , S(l,21,22)i X17SOO), Y(?500), Zt?SCO), »iT(?500) ROL0100C
ESL'l VALC'ICE (CtA.TiTAT), ( CC . APRA Y I 1 81? I I i ( 0V , A ^R A Y ( 36?7 1 ) , ROL01100
I 15,A°"., v!36'71), (0--,A'J^\Y(55'J5)),(U,APRA1'(7373):, (X,AhRAY(12?37ROLOI200
2)), ( Y . A~JAY ( 22137 )), ( Z , AR ° A Y I 3 1 7 3 7 ! 1 , ( V(T , AR1? A Y ( 1 1 1 3 7 ) ) RCLOI30C
i.i'rc*:' SU>'»,Y RCLCISOO
• ••••«• ••«• ».«.«».4».«»••««••.•«•«• •«••»••••.ROLO i soo
r o 20 * N
"E:/D («) c, cc
°£» i '0 M
IF (N? .to, 0) RETURN
" • 1 - 10
''CAD ("•) ixi i j , i«i ,NPJ ,
1 (f T( I ) • I-l tNP)
> ,1 = 1 ,NPJ ,
EI'T'T POUOUT(N)
M . -I • 10
!F ( .? .E3, 0) GO To 10
(•SITE (•<» (XI I ) ,I«1 ,NP) , (Yt I ) .1-1 .NP) i (Z(!) , 1-1 »NP) ,
1 I i T ( ' ) , i ° 1 t N p )
pt» !' 0 K
10 y ' ' * 10
1 R \ T E ! H ) C
ENO
POL01600
ROLCI700
HOLO)bOO
ROL01903
ROLO^CCO
ROL22IOO
ROLC2200
ROLOJ300
KOLOr"00
POL02530
POL02503
ROLO/7CO
PO(.C2tOO
ROL02700
ROL03COO
P 0'_ 0 3 I 0 0
RCt.03700
POL03300
ROL03SOO
ROL03500
O4
CO
5=3
i
CO
-pi
-------
Ft,CrtCuART£D sY FOTLO /XGC8/ ON ]3 OCT 71 AT 07158:00
--0
Cn
C ROLLININ)
N OX,3Yi02iTO'<.TOYiTDZt<'PINiP'''A>r , .(Nt
OTi °MlSTF(5)i PhXST S ( 3 ! » Vf3TO<"5)i
TP(S), r?STOR!5!. FACT36)i
i),Fi.'Jvlt2: >iFLl/Y2(2l)»
M , 2 ' , 5 I , ', X , •< Y , '-1 3 P , « A X , „ X M 1 , t! Y " 1 , '. Z !! 1 , t, A ,
• - 0 7 , ' ., X , > A T , V A 7 , - A i, i V Z , ,-. * ! i N A 2 , N ,\ 3 ,
I. / t , C ^ v "
C • " ; N S i C N C I 1 . 7 ! , Z 2 ) , CCtM,71,22), 0 / ( 4 f
2'i22), OMC-tZS.ZZ), 'Jvu,2l,72,3), S(",
2 .Z2!i XI73CO), '(95CO), ZI950C), I^TI
SSCO)
£3Ul\A_E\CE (C-APflAY), (CC • *R"AY ( |Et9) ) ,
(0V ,«peAY ( 3fr97 I ) , ( S .ARHAY (3i?7 ) ) i OHi
3 i ' : 7 I ) . ( W T i A r P A Y t "i ! « 3 7 ) )
I. TEGS R 0'- MY
1 —C
I C
I C
I
I
H « 2C » N
PEAO (") Ci cc
REV/»NO
»POL01SOO
3
I
I
I
r
i
Oo
IF UP .EO. 0>
I TRUE
I
'\ FALSE
-------
:
/ READ
/ (X ! 1 ) t 1=! ,NP) , /
I
I
I
I
; RETURN :
i
ENTNP) t (Y( 1 ) ,1-1 ,NP) , /
i
REV. 1 NO K /
I
M a K •» 10 I
I
1
ViRlTE (M> C /
i
I
I
I
'
z
z
z
z
z
z
z
z
z
z
2
z
•
,
w
C/3
CO
-------
ShTUP/NhXUSP
FLO .cn-'TED BY FQ^FLO /xsoe/ ON 13 OCT ?i AT 07:53:05
S'JCRO'jTlME 5ETUP STPOC100
TM 3 °CL'f!^E -SACS INFOPM/ITION NECESSARY TO START A PROBLEM STPOCZOO
F?OM UNIT 0. STP00300
•Of S7PIS), PIXSTR(S), ViNST01(S)» SCSSTRIS), NPSTOR(S), STP00500
2F,'C7!r.>, t^.x, FNV. F'U, V, RV, N'P , T, ICYCt A;! ,NZM1 , \M , NT OT , VAX , VA Y , NA Z , NA "> » NZ , STPCC°00
6.,1, '..i2i ••i/'3, N,M, OW""Y STP'JIOOO
C;:c... 0!~=>3IC.I CM,21 ,:, CCO,21 ,22) , DVtu,21,2?), DHt1,21,22), STPOMOO
C2:C12. 1 U( <• .21 , 27,3) , S(<(,21,22), M9SCO), Yt?500), Z(9500), IVT(9500) STPC1200
C:.313. :C,l-'j V^LE'-CE >C rwS'l'.Y) , (CC.rtRRAYI Iflt9) ) , (DV,«P°Ar <34?7) ) , S7PC1300
CrC".1-. 1 «',A»?AV( 3497) ) , (OH.ARf AYtSSl";) ) i (Ut Af!RAY(73?3t ) , ( X i ASRA T < 1293 7S TPO 11 CO
60 C::C15. 21), I " ,AK'~AYI22137) ) , ( Z , AKR A Y I 3 1 937 ) ) i ( «T , ARRAY ••....<...»»«•.•••..«<.•..«»...STPOI700
DO 30 •' = 1 , ,JSP ' STP CleOO
o::E2c« RE-o '81 PMK;, PPAX, AN, DACK STPOZOOO
CCC;2>. READ (8) C S7-02103
CC.cn^. \P * C. STP02^00
c;:c/"3. 'NS70=!(N) » «n STP02joo
cc;:2--. 3C"sT-!(iS) = SACK srrozfco
CCr-25. PJ"I37"('1)"DKIM STP07DCO
CCC":*. °,',v.37n(M » P.",AX S7PC2600
cr;-27. oo 2C ; » i, NX ' STPo;70c
c:;.. r 3 • C'r:j=i,NY S7P02ooo
CO 2C < « 1, M S7P02»00
C ST/P71NG ..EIGHTS FOR EACH PARCEL ARE DETERMINED BY THE CONCENTRATIONS
C CM 7"E DATA FILE AND THE AVERAGE HEIGHT _ S7PC3100
C *> IS 70TAL HEIGHT OF POLLUTANT IN EACH CELL " STP03200
« = (CCCiJit) - BACKI/FACTtN) STP03300
IF (<• •(.£• P«;M CO T0 20 STF03-OO
c 3'°E^c '. AMONG ..'0 PARCELS, EACH '-JTH AEIGHT PT» S7P035oo
^0 . A/.,.. » 1 S7P03&00
°r, r ,v/\'0 STP037CO - ^
00 10 L » I, NO ' S7P03JOO -g
'IP » \P * 1 S7"OJV30 i
MCI") « "ft S7PO"000 CO
C STARTING P'RCEL POSITIONS A"?E DERIVED WITH THE AJO OF URAND, A UNIFORMLY -P"
c c:s7--1°u7Eo RANDOM NUVBE^ GENERATOR COVER THE INTERVAL FROM o. TO i ) •*"•
xi.'.p! « : - .5 » u»Aiiotcuv*;Y) S7PO'-303
• C") = J - -5 * V,»Af D(OJ^.1<') ST°0'(103
:t.«P) » K - ,5 » URAf,D(OLv>1Y) S7POS3CO
C lr "APCELS FALL OUTSIDE TME GRlo LIMITS THEY A^E THROWN A»AT. STP016CO
ir'Zt P).CE.F"Z,Or>.XI',P).GE.Fl,X.OR.Y(NP).GE.FNY.O^.X(NP).LE.l««OK.STfON7oo
1C CO 7, •.'-'£ S7P01VOO
23 CO Ii' JE S7P03000
-------
w
^J
oo
CC-77S.
CC777S.
E..CH SORCER CELL IS 6 T[MES THE NORMAL CELL SIZE! THREE TIMES
A3 LO'.G AND TV.1CE AS HIGH.
!>J!TUL COT CENTRA! I ON ASSUMES EX°OS'ENTIAL DECAY TO DACKGROUND
LS/E1- - 1/2 I'l 8 CELLS
CC7377.
DC 7 U=lil
00 7 i = 1 > '• X
D^ (L , I i • ) «
7 DX2IL,I,N) -
CO 3 J«1 > U Y
9-1(L,J,N) "
8 3 v i ( L , J , M =
? F-5-f-C-FAC
f «• ;Q » (.
.." • T E ! M > C
«= P - 10
fr^!Tc f) (X( I ) i
! (> T( I » i I-! »NPl
PC' I' - M
CAL1- OUTPUT(-N)
E'JO
(Cll.l.II - EACK)«FUG'&./FACT(N>
(ClltNY.I) - BACK>'FUG»*«/FACT(H)
(C(l,J,i) - SAC1?) *FUG«6./FACT(N)
(CtliJ.MX) - 6AC<<>«F'JG-6./FACT(M
(Y( I ) i I
IZ11 ) ,1-1. NP> ,
STPC5300
STPOSMOO
STPC&500
STP05600
STPC5/00
STP05800
STPC2900
STPOAOOO
STPOA100
STPC62CO
S7P06300
STP06"-00
STP06SDO
STP064CO
S7PJ6/CO
57P063CO
STP047GO
STP073CO
STf-07 100
ST°07?00
STPC7JOO
STP07M03
STPU7SOO
STI'07400
STP077CO
in
oo
-------
SETjp
FLOWCHARTS gf FORFLO /xsoe/ ON 13 OCT 71 AT o?:s8:o7
'E SETUP
tc
-J
i
j c THIS ROUTINE READS INFORMATION 3
CNECCSSAKY TO STAcjT A PROOLEM 3
CSTPOC203 3
C F), YINSTop(S),
BCKSTPISIi ,PSTO,->!5>. r/,cT(S), TNX, ^MY,
FNZ, V, RV, N?, T» ICYCi AH " A Y I C 0' 3 6 ) ,
FL'-'/l [22! ,FLUX2(22) ,FLt'Yl (21 ) ,FLUY2(21 ) ,
N-OT ,,-AX ,,'AY.I AZ.'JA,, ,I.Z i N « 1 ,
'-\H, OU"MV
o : '• E N s • c N c 1 1 1 z 1 1 ? : i > cctiiz
ZliC?). D«( 1 .21 ,22 ) , l/("-,7l
: DO 30 N • I i NSP !
A :.. ^ >.......... ....... ........,.,..,.., ,....:
t. i
A :
A t.,t........... .......................... ••••
A : FACT(N) » ?8."«V/FA'CTtV)
A : ............. . ...... . ..... . ...... . ....... ..:
A I
A 1
A /• • t • t t . i i i. i . f , t n t . m .. f t i i t m t i i,, t,,/
i / READ (8) PMIN, FMAXl VN, BACK /
A / ......................................... /
A I
' i
/.•,.,,,.,,,,,,..,,.•,,,,, ..... ,.,,.,,, ,,,,
A / READ 13) C /
* / ......................................... /
A I
s I
A .. .............................. .... .........
••IP = o ;
C/)
JO
i
CO
-------
to
00
o
A
A
A
A
A
A
A
A
A
A
I
t
A
A
A
A
A
A
A
A
i
A
A
A
A
A
A
A
A
£
I
A
A
A
A
A
A
A
A
A
A
A
A
A
A
5
3
5
0
D
3
3
3
3
a
3
n
0
5
3
e
3
S
B
a
a
5
c
r
3
0
S
3
B
5
D
3
a
3
a
c
a
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
r
c
c
c
c
r
c
5
c
0
3
0
c
r
3
3
0
0
C
0
3
D
3
0
*%
D
D
0
0
D E
D £
t* ~
0 t.
0 r
D E
^ r
^ r
0 E
; c
PSSSTIMM
PWXSTKIM
I
I
1
I
I
I
I
i
I
I
I
I
I
I
: x - tcix.J.t i -
i
i
< IF ffl .LE. P,1I
: NO » W/.Y
: PYI •>
-------
Co
A
A
A
t
A
A
A
A
A
A
A
A
A
A
A
A
t
>
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
i
A
A
A
A
A
&
C
e
3
s
r,
s
3
r
r
3
e
3
3
Q
3
E
e
e
c
3
3
e
n
c
c
c
c
c
c
r
c
c
r.
C
c
c
c
r
C
c
c
c
0 E
3 E
^ r
•* r
1 E
^ F
D E
D £
3 E
S E
^ r
C £
2 r
r i
0
D
I
Xt'JP) » I - .5 * UPANO !DU!"ir ) , '
VC1P) » 0 - .5 » U7AN-' (3UU1V>
--
< CE • «•':», OR, X(NP ) .LC« 1 • tO
I
1
: t.— i
C 103 i
I
C 103 I
I
I
: NPSTORlM »
I
I
I
I
I
I
I
1
I
: FAC-.A
FUG'FAC
1
1
I
-c IF PAPCELS FALL OUTSIDE THE 3
CG^IO LMTTS TKET APE THKO»N A«AT3
C ST"OM620 3
,T.Y<"P> .LE. i . ) NP« > 1
TRl-E I
I
: i
i
:
•
•
•
HP '•
-t SETUP BORDER ARRAYS STP05ZCO 3
C EACH BOEDER CELL !S 6 TJ«£S THE3
CNORMAL CELL SIZE. THREE TIMES 3
CSTPOS300 3
C AS LONG A.VO TWICE AS HIGH. 3
C iNiT AL COMCENTRAT ION ASSUMES 3
CC*Ptf'.E'IT [AL OECAr TO BACKGROUND ]
C LEV-;t - 1/2 IN 0 CELLS STP0560C3
•
s :
j x r
Z
Z
z
z
z
z
z
z
z
z
z
T
z
z
z
t,
z
T
t,
z
z
z
z
z
CO
>3
oo
-------
A
A
A
M
A
A
A
A
A
A
A
A
A
A
A
W A
oo t
•^ A
A
A
A
A
A
A
A
A
A
A
A
;
A
A
A
A
,\
A
A
3
**
T
a
3
D
3
3
e
D
n
n
n
3
r
c
c
£
C
c
c
c
c
f
c
c
I
: 3x:«.,:,N') « icd.i.n - oAcxi«FUG>6./
: F A c 7 1 N ) :
i
C 7] I
._: arz(L,i,N) - ten, NY, n - 3ACKi«FUG»6./ :
: Ffcrt n :
i
i
->• 00 8 J-> 1 ,'!T :
i
i
: SYllL.J.NJ - (Ctl.J.l) - 9AC^J»FUG»6./ :
: FACT (r ) :
I
C 03 I
•- _--• 3Y2(LiJ»N) " (C(1|J|"X) - BACK) «FUG«6./ :
: r A c T ! H i -_ :
f
c ?: i
1
1
: » " 23 » M i
I
;
/ ARlTC (n) C /
1
I
/ pE/'INOH /
I
I
: M « c. - 10 :
1
I
/ i^lTE (M) (X( I) > 1 = 1 »NP) • (Y( P f 1 = 1 .'.'?>. /
en
OO
V
;p,1 = 1 ,MP> , ,
izi'!.:=>:,%?), {^TI;),' = I^
-------
* : : CALL OUTPUT(-N)
A :.:
A I
A C 333 I
CONTINUE
O3
CO
-------
zr<-- -FLC.FOSFLO.S SOURCE/NLXUSP
FLO^C^A^TED ar FQRFLO /xeoe/ ON 13 OCT ~i AT 09:58:15
to
n r " - n 4 ,
cc-c:u»
C 3 D ™ ' u »
i- r •" — •» ,
CC:c22«
c r z o *. u •
==:-z?:
-:=-:^:
ceres-.
>•----'•
cc:--:.
C * ~ " "* »
*- •• — - •
SU3FOUTIN-; SO'^RCEtN)
C "HIS SL'DRCL'T i f E PCAD3 THE SOVCES FRCH UNIT 0, AND CONVERTS THEM
c INTO "ATCELS, JS:JG THE AVERAGE WEIGHT V.N. PA"CEL POSITIONS
c APE ASSIGN u3i..a T"E RANDOM NUMQEP GENERATOR.
C fRfC IS >;SEO TO DETERMINE THE VERTICAL FRACTION OF THE CELL THE
C ME' P A R C E 1. 1 S ) ..ILL D E PUT I N T 0 t
]P1 \5TRlS), PI'-'STMSJt *SSTO^(51i SCKSTRfS), NPSTO'flS),
2rACTI!:!, -<>r, F.'Y, F'.Z. V, "\ , HP, T. 1CYC, ATRAYISOV3&) .
2FL'Y- ' (22! ,F'_'JX2(/;) ,r.UvI(ZIJ,FLCY2(2ll, ;XI(i BXZ<1,ZZi5)
HLY 1 C- ,2 . ,i) .SYZCi.Z! ,3' ,
6 ' ,1 . ' ^2 , 'S3, "AM , C'j.^MT
D'«t,s C'- C'",Z!,ZZ', CC< ° .2I-.Z2 J , 0V ( 1 , 2 1 , 22 ) . OH(1,Zt,2Z)f
I 0(1 , 2 i , 'Z , 3) , S ('< ,21 |Z7) i .Xt'^SOO), Y(?500I, ZIVSOO). «Tt?500>
E-l.'IVALC^CE (CiAc"AY), (CC.APSAYIIR I'll), (OV,A = RAY(3697])» <
1 ) , tZ,ARRAY(31»3D ) , ( '> T , A^RA Y ( 1 1 137 ) )
J TA F = -,C /. 33333333 /
READ (0) 11 (S 2, NYM'
so rc : i » ?, Nxfi
00 20 "' = . . 2
c NO ;s "r-s -ij"r,r^ Or Nen PARCELS CREATED
IF'S (!"., JJ, ' ! ) .LE. 1 .OF.-20) GO To 20
NO * sttr , jj, i II/V.N » .5
1- ( 0 .Ea. 0) GO TO 20
P,I - nr,*;, jj, , i )/no
00 1C L i 1 , NO
\" = ,;F » '
•i T ( • P ) •> r /
X^r) « i; - .5 * L'RA-,D(D'JMVY)
ZI.H ) « « - .5 " UPAnOOUMPY 1'
1C CC' TIf.'.'E
20 CC'i-!M.'E
C S01.°CES rOR 30ROER
C
C Y = 1 . EDGE OF GiID
C IF T"E f.I'D 13 BLOWJN& INTO TME GRID, AND THERE ARE POLLUTANTS
C I.I T"E CO-DER CEuLS, !,E»; PARCELS APE CREATED A«0 PUT INTO THE Z
C LO.iC VCPT'CAL CFuLS. THC'R , 'EIGHT IS ADDED TO THE EriSTING
C CCLLv.LAp CO'CE-^AT'p-jS.
C IF TiC ..EiOr'r OF T"E '!£" PARCEL IS "00 SMALL. IT IS AOOEO SACK
^0 : 30 ! 1*1 ,NX
I-IFLL'X! I ' ; ) .LE.3. ) GC TO 115
; r _ x - F L '. " i : ' i
SRC00100
SRC00200
SRC00300
SKCCOMOO
SKCOOSOO
S9C00600
SPCCC300
.SPCOIOOO
5f
-------
oo
o'-:c'j<;'
cr-c 7 o•
1C?
lie
1")
115
: ! ?
12C
; 30
rL'JXl ( 1 I 1-0.
NO = XrL */•••'». 5
Jr (f,j ,E5. 0) GO TO 111
f i * r F i_ x / ' c
00 1 i 3 '.»! ,S'0
. o - .. c . 1
v 7 ( '. P ) « D .'
X i;,r I •- I I-.S»UPAND(DUO!0'J"1/Y ) » • j
2c o) -.5* (U";,'.0(0l;»t '.T ) ) 'Z'
.r(r' f) •GE.FMf.CT.xUP) «LE« 1 -'GO TO 10?
' r. = z i ' ~ ) * s
C !"K t ' • 1 1 ) = C IKK, 1 , I I I * PB
03 TO i 10 -
t ?=.--!
c ; T i . -' E
i 0 TO J 1 S
BX, ( ' i I i ,r'J » 8X1(1, 1I.N1 * xrLX
v = .» EtrE or GRID - PARCELS APE TREATED IN A BANNER SIMILAR TO
T"-T CF 7KE OOT70," CELLS.
!c IFLJ72 I ! 1 ) .LE»0« )GO To 120
X-t-«FLk S21 I I !
F •_ ' X Z ' i ' - 0 .
> ' -xrLX/A )»,.&
IF ( .0
00 '20 L=l ,VO
I, o c . o.i
X" »)«• ' I-.5«L!»AN3<3'JKKlr>
Yi i-'/="'--'>A<;!0'Jt"")«'5
n r > = • 3 " ( u ~ « ,•> o ( o i. " K T j ) « z «
[r i x (• P J ,G£.F J/.OR.XINP) .LE. l.)GO TO 119
rr = z ( ^° ) » .S
C'-X,' f , I J ) - C(K<,NY,II ) * P"
00 TO 120
i "•= ,?- s
CO1 T ' •'-'£
00 TO !J3
s;2(i,'i,M • a/Z(i(!!i.M) *
co' T U.JE
» = 1. E05E OF GRID - PARCELS
T" T Or Tr-E BOTTOM CELLS'
00 23C JJ=1 ,'.Y
If (FLl-'f. j jj> .LE.OJGO TO 215
XrL/ = c-'JTi (JJ)
F'.O'i ( J-' -C-
..O-i-L''/*. »«5 v
;=• CO .Er.. 0) GO TO 2M
P, =?-.X/I.O
20 210 L=l T'O
i 7 ' =• ) = P i
t , • ' I = i . - ij •; t j o ( D 'J v M " ! • • 5
» ? 1 =.'j-.i-o-A '0 [DUMMY )
Z!10.--5-(l.:j.SE(0'/"<'v))«2.
' f • - ' I," ) . ' ~ .r T .S1-'' Y I • " ) «LE' 1 • 'GO TO 7C?
SRC35200
SPC0530?
SRCC563'3
SFCOi7C3
SKCO = l!CO
srtccbvco
SSC04C30
S'KOi i 30
S i C '-' ft 2 3 0
SRC06600
5NC06700
SPCC6300
S3C07COO
SKC071CO
S^C07203
S1C07100
5Kt07600
SPCD7700
5^:07200
S"C07"CO
SPCCOODO
3PCOU1C3
SSC08100
sncoesco
SRCC3700
s P c J a i o o
S'-COdoCO
S-00
5 '( c 1 C / 1 3
b<
-------
corjc?.
CC?)1C-
00
f< •= ZIN-*! » .5
C("*tJJ,l) - CUK(JJil> » P«
CO TO 210 •
2C? •' " = .' p~l
210 CC«TI"J£
CO TO >iS
Ztt 8-1 ( i .JJ,',') » BTHtiJjtN) + XFLX
c y • %.x EOCE or G^jO - PARCELS A^
215
TREA'EO IN A MANNER SIMILAR TO
SRC10730
SRC! 1000
SRC1 I 100
s"d 1200
SPC11300
SRC I I HOO
S^CMSOO
. )GO TO 230
o co TO 22?
00 T2C L-l.'JO
1F( T C.P I .GE.f |'Y.O".Y«NP) • U E • I . ) GO TO 219
KX » z't.B) * .5
GO TO 220
2 7 ,p = r,p-i
220 CC'.T. ,1-E
CO TO 230
229 BV7(1,>J,N
OY2(I,JJ,H) * XFLX
S"CI10CO
SPCI1TCO
SSc'2000
SSC12100
s^c-zzco
SHCI2300
SKCI2"00
SRC127CO
SKC12UOO
SRC12900
3RCI30CD
s\ci3a?o
SPC13MCO
"5RC13SOO
S°C13600
cr
SRC13900
SRCISOOO
CO
oo
-------
8? FOTLO /X308/ ON' J3 OCT 71 AT 0».'53:i6
50Ui~CEJ CF TMC CELU 3
CT«E 3RCC05CO 3
: ••.£¥• Pt^cELisi WILL BE PUT INTO.3
CSKCCC6CC 3
Co
CO I'M ON OX,DY,?Z,TOX,TOYiT5Z.r.-lIN,P«AX,WN,
pfNSTt(5j, P.MXSr.- (51 , f.fJSTOR(5) ,
NFSTCVS). rAc">S), FHX, FNT,
r Vi RV, 'JP, T, '.CYC.
FL'- 11(22), rLU?2(22),FH/Yl
SrlJ", 27,31. 3X2(S,77,5I,!)T1(«,2I,5),OY7(
•4,21 ,S) , ,X,' Y,f 3P,M,,X.NX"1 ,f'YM| , N Z Ml ,NA,
N70T ,'JAX ,MA» t%AZ iNAid^Z , !<'>!, •, A 2 , HA3,
N ' H , 0 U " r v
C'"£ 5(0. C!«,21,?7), CCO.ZI ,22) , OV(t,
,2! ,22)
UtH,7l ,27,3) •
Z(7SCC>«
SIH,
l.Ft,CE (C,Ar
-------
A
A
t
k
f.
A
f.
A
A
A
W A
00 A
A
r\
A
A
A
A
A
A
A
A
A
J
tt
3
3
3
3
3
o
3
~
r
r
s
J
3
p
D
D
3
n
5
3
3
n
"
2
•s
C
C
r
C
"
^
C
c
c
r
c
-
~
c
c
c
c
c
c
C 0
C D
C 0
C 0
C 3
C 2
C 0
C D
C D
C D
c
f
1
I
I
• fs »/ 9 •
I
j, — C MO IS THE NUK3EP Or NE t PARCELS3
i CCQEATED SRCOJAOC 3
i
i
i FALSE
i
: NO • SIKK, jj, : i )/«N » .5 :
i
i
I FALSE
Pn o S"'1
CIKKiJJ,!!) = CIKKiJJ.tt! f P«V
1
C J03 I
1
C 203 I
z
z
z
2
Z
z
z
z
2
Z
z
z
z .
z
z
z
z
z
I <"
2 CO
z f
2 CO
•Ti
z
z
I
! C SOURCES FOK BChOCR SRCOtOOO
i
! C T = ]. E?CE OF
-------
CG'MO, AND THESE ARE POLLUTANTS 3
CSCCC"3C3 3
C IN TsE 3CRDEP CELLS. NEW 3
CPA»CE'.S AJE CHEATED AS3 PUT INT03
CTHT 2 SSCOSMOO 3
C L?v,£:r> VE^T.CAL CELLS. THEIR 1
CrtEiCuT !3 AUOED 73 THE EXISTING 3
csPcd-isoo 3
C CELLULAR CONCENTRATIONS. 3
C3°COS6C'0 3
C IF TJE HEIGHT OF THE NE« PARCEL!
CIS TOO 5"ALL. IT IS ADOEO SAC1? 3
CS"CCS7C" 3
C INTO TME SORBET CELL ARRAY. 3
CS*CO'tOC?0 3
CO
A
A
A
A
A
A
A
A
f
'•
A
A
A
^«
A
A
C
A
A
\
A
A
*"
A
A
3
3
q
n
3
a
3
u
i
i
I FALSE
i
: XFLX--LMXI i n )
FLUX: ( i ' ;=o.
; f<3»XFLv«/.''N».S
:
I
I FALSE
I
I
: PA«=XFL*/\O :
i
t
•---••-- •> * OC 1 'C L«l »NO "
1
^Vp»^,'+ |
fi~ • NP)=?,-.
XC.PJ'I I-.i + O-A'.'DiOUKHY)
YC P ) = 1 . »J?A». ")(DUM11')»,S
Ztfl»l*.5«(i;pANO(OUi'«''))»2.
I
I
I FALSE
I
X
X
X
X
x
X
x
X
X
x
X
X
X
X
x
X
X
x
X
x
x
X
X
u x
u x
U X
2
z
z
z
z
4.
Z
z
I
z
2
Z
z
z
I
i
I
I
z
£
2
2
Z
Z
2
Z
Z
z
z
z
2
Z
z
z
z
z
CO
70
I
CO
-------
to
+ P« :
I
C 1C°3 I
: NP-NP-I ;
i
C I 103 I
I
I
c in: i
: BXIII.It.N) . BMC1,iI,K> * XFUX :
5
S
5
S
s
S
s
5
U
u
u
u
u
y
U
V
x
X
X
X
X
X
X
x
x
X
X
x
X
X
I
i—c Y » >JY E?GE OF Gdo - PARCELS 3
I CA°E 7^LATEO IN A 1AHNCR SIMILAR 3
i CTO S^C07000 3
I c THAT or THE BOTTOM CELLS. 3
I CS-IC37IOO 3
o<
C 1153 I
,,,,,«,,,,,,.,,,,,,,,,.,,,.,,,,,,,
< !F(FLUX2(I I).L£*C«)GO TO 130
\ ,...,
1 FALSE
I
•\ TRUE
>__-.-.
./
FLUX2< 1
IF (NO .EG, 01 GO TO 129
i FALSE
i
i
•\ TRUE
>......
./
-0
Z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
i
CM
CO
JO
I
CO
-pi
-------
w
A
"
A
<,
t.
t,
k
A
A
A
n.
A
A
A
•*
A
A
•
A
A
A
t
A
*
8
0
r
3
0
3
-
a
S
a
B
2
£
3
I ,/
1
ZtSP) -i,5*{t'9AND tOU ^Y)j»2.
,'< T ( N =" J a P 1
I
t
t
: CC"!»NY|It) - ClKliNY,! I > * PW
1
I
c si?: :
; N p » •, P - 1 i
i
C I2C3 1
- * * .
I
C I2'3 I
: cxzu.n.N) = BXzdiiiiN; * xrux :
i
C 1303 I
X
^
^
X
X
•s
X
X
X
x
a
X
u x
t A
U X
U X
u x
u x
U X
•J x
U X
s u x
s u x
S U X
s x
s x
s x
s x
S X
X
x
X
X
x
X
X
z
z
z
t.
2
Z
Z
z
z
z
2
z
Z
Z
z
z
z
z
2
Z
Z
z
z
z
„ z
z
z
z
z
2
Z
z
z
z
2
Z
,
t/J
en
50
CO
•*
CO.MS..JE
I
.—c * » i. EJGC OF GPIO - PARCELS 3
I CA>C TREATED It A MANNER SIMILAR 3
1 r TO 5pc i1? no 3
-------
I C3RC09SOO
vo
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
f.
A
A
/*
A
-
A
A
A
A
A
A
A
**
*
3
a
S
3
3
3
S
3
3
5
3
3
-
n
3
3
I
I
I
i
i FALSE
i
i
! NoHIrLX/'.H^B
i
I
! FALSE
i
I
I
I
1
I
.* 7 ( N P ) « P 'V
X(NP)"l»+URAr 0 < D J M « Y ) « , 5
Y(NP)«JJ-.S*'J'«'i .0(OU,*VY>
I
I FALSE
I
: ciKrfjjii)«c<>*CiJJ»i)+pa !
I
I
C 20*>3 I
u
V
0
u
u
u
u
y
s u
s u
s
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
I
2
z
z
z
2
Z
Z
Z
z
z
z
z
z
z
z
z
•>z
z
z
z
z
z
z
z
z
z
z
z
z
2
z
z
z
z
z
T
z
z
7
en
X"
i
oo
-------
w
VD
* 3
A "
A 5
A
A
A
I
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A 3
A 2
A 3
A r
, ?
A 3
Np»NP-i : /
I
c 210: :
i
i
C 2!f3 !
3fi . DTI : i , JJ,N) - xrtx :
i
, C X -, NX EDGE OF flSIO - PARCELS 3
I CARE TREATED IN A MANNER SIMILAR 3
I CT9 SPCI 1 60C 3
i c THAT OF THE BOTTOM CELLS. 3
I CS'-Cll/OO 3
i
C 2153 I
I FALSE
I
I
: FLUY2I JJI-O.
I
J
I FALSE
i
: PYI-XFLX/NO •
J
I
i
i
: fi'"--p», :
, 7 ( •>' p ) = P /(
X (N'° I "NX-u^A •: (BUM"* ) »«S
: T! «p )nJJ-.5»Ufi/* D I O'JMKT)
Z ( • «• 1 3 . S - ( 'J .J A >• 3 ( C J 1 •• 1 ) I •> 2 .
-
S X
S X
S X
X
x
x
x
X
X
x
X
X
x
X
x
X
X
x
X
x
X
X
X
X
X
I
Z
z
Z
z
z
z
2
t
I
z
z
z
2
z
z
z
z
z
2
z
z
z
z
z
z
z
z
2
z
z
Z CfJ
2 (/)
z ^
2 I
z »
2 *
z
z
z
z
z
2
Z
z
2
7
-------
(33
A
*
A
A
rt
"
A
A
A
\
A
A
A
A
A
A
3
2
J
5
3
3
J
fl
5
o
•
i /'
I FALSE
1
I
; KK » 7
-------
fi.Iv,C"A?7EO 3" FcRrl-0 /X808/ ON 13 OCT 7\ AT 07!58;N6
c mANO ;s A MACUINE INDEPENDENT RANSOM Nurse*? GENERATOR
c it-: np5T VALUE or NN SHOCLO BE EITHER i CR s ro^ BEST RESULTS
C "HE CYCLE UTNGTH IS S1?Z NU^SERS OEFC^E REPEATING
C ISTEGE'''1'1 i-'N
N •( » lOo !' '•'•4359,32768 )
UROOOI03
URDSCZOO
UROOOJOO
UP000500
UH0006CO
URCCC7CC
UHOOOOOO
URU009CO
tn
CO
CX3
-------
L'PAsO FUOi*CHART£0 8* FORFUO /X8C8/ ON 13 OCT 71 AT 0?:58:i7
\ FUNCTION
\
U1AND IS A KACHJNE INDEPENDENT 3
I CRANOOf, NUMBER GENERATOR URD002003
I C THE rIRST VALUE OF NN SHOULD BE3
I CEITHE? I OR 5 FQ1 SEST RESULTS 3
I CU«D;3300 3
I C TV,; CTCLE LENGTH IS 8!?2 3
! CnU"3CRS OEFOSE REPEATING 3
I CUTOC'-CC 3
t C !NTECEP«H NN UROOOSOO 3
1
I
NN » MODtHN«60S9.32768)
U?ANO «• FLOAT(NN) • .C0003Q5176
w
OJ
00
po
OO
------- |