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)- PKC LS. >HIT
COVCLNTRAHONS AND
PARCEL POSITIONS
                                         ,
                          CN LNITS 21-20i/
                          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 1r
                                                    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

-------
?FC3FLO.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.  NZ1,  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,  Ft02.  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  /  23218.6711, 100., 3QO., 21400.,
CICO!?.          125150./
coc^2o.           DATA  NAME  /  'co,  'NOI,  N02',  HC' /
CCOC2'..           NX-25                                                   .                                                       
Co:o22.           NY=25



or:s26.           zzoo
c::o27.           icvc  - o
CC:"23-           00  5  Il,22
;D:C:?           oo  s  JM ,21
CCCC33.           00  5  ;! ,NZ
CCC;3I.          5 XXX(t,J,I)90
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 ,JH,21) ,I3,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 ,jH,2M> ,i3,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),J1,21>,1-3,21) ,  (((OH(XtJt!>iKl.4 It
C3C3S1.           3  J*1,21>,1-3,21)
ccross.            WRITE IN)  ! (ISIK.J.I.2',KI,2)j-i.zi),i-3,zi)
C3C054.            I-RITE If.)  I IUltI>J),I3,21>,J-1,21),  ( ( V 1 ( I , J) 11-3 ,21) , J-1, 21) 
3C:037.           I  t,ji,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) ,I3,21),J1,21),
CCCC67.           1  CCut I.J),1-3,21) ,Jof,21), ."
C3CC63.           2   (< (OV(K,J,II ,K"t,1),JH,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) ,KI ,1),J1,Z1),1-3,Z1),  (((OH I X,J,I),KI.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
                                               CNSp1                            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  /  232I8,
                      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> iN(3) .BACKI3/
I
I
/ WRITE (N) U(C(K,J,I ,3> rK| tNZ) tj1,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) ,J1,71) ,  ((      /
    VI(I,J>,I"3,2M) ,J"H,2M) ,  lU,J),t"3,       /
                       .J,tJ,K*1,1),J1,        /
                         iDKl .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) ,I3.2H) , ,)"(,2 , I3,       /
   2t) ,J=",24) , lDV,  <((I,J),i"3,       /
   211 ,J",2M) , ( ( ICVIK, J, I ) ,l"M ,") ,Ja,        /
  2" ) > I"3, zM ,  ( ( (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) ,jH ,2") , I JIDV(K,J,I ) iKl ,1) ,J-H, /
21) , 'O,2H) , < ( (DHCCiJt I) tKl >1) i JM, /
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,2I , ( Ut I ,J) ,!3, /
2") , J3i2'O i  tl'"! ,)  J-Ht /
2tJ . I3,2f) /
;
i
/ WRITE: /
/ II,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
 /ttttti,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(6lZ6)
OO      CCS015.            00 13 L = I, NT
        CCCC'6.            J--T-1-L
        c-::i7.            ^^A.5  so;o,  (Ci I ,J, I ,N> > II,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,tN> "  . 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,  PINM)  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,  PIN(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              :cUJl,N>"OIC(lJiI,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                :                  JNYI-U                   ;
A 3                :.	:
A 3                                       I
A3                   C    203              I
A 3                   /ttitt,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,JICI
                                                                                                                                             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,J1>>
                  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.fCE 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 - u340o :
j
i
/ READ *oot, L /
i
i
/ READ 4000. x,u2,vz /
i
j 	 C 4001 FORMAT(IS) 3
1 C 40CO FORHAT(6E:i24) 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..            OIC'.'SlC  S(M,J5,2S,M)
C"007.            Ei'JiVLEMCE  13, C)
 ::c~r.     c      or^sicti  TOTINS"), CT(.MSP). SC(NSP), STOREtsZtMir
Cc:c,:.            OI"E3'0'l  TOTf), CTO), NP(M),  S TORE 1 1 , 25 , Z5 , H )
cccoiO.            DATA  NMC  /  'co'i 'NO',  'N02', HC' /
c c o c . : ,            N r, = 2 s
                   OCICNI,NSP
                   00  10  I  =  1,  NX
                   oo  10  J  =  i ,  NY
                   ooio.'=-i,\z
                   si>;,j,:,N)oo
                   OTH    or/3600.
CLL 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"? * STORflC,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  "  CTM)/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' /
                                       NX25
                                       NY.25
                                        NZ"
                                       NSP1

                                           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 fiAD 5010,  (FM(I,J),J.I,N*>
5010 FoHAT (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  36ooi IT -  n
     IF (TS -  OT .(.T. TT)  SO  TO  10
     Tf * OTH'TFFIIT)
     FS > 3"PTFS( jT)
     GO TO So
  "0 FF   ((TT  - TS * OT)TFF  * 9N'0!(FNOF

  6C S(I,J,I,1|   Stl.J.t.H)    (QHC    QHC33)F * VAPFSSMI,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
                                            tc  DIMENSION rn
-------
                                          I                           /'                                                           Z
                   	,	     /                                                            2
A.	.	->:              oo'20  ji,  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<>ti\   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 +  DTiTFS(lT-I)    (TS  -                                                                  X
                     T:iTFS(: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  FFFM(I,J)  * FS'SKUiJ)
A 3                '      S(l,J,I.l) o StliJiIil>   QCO'F
A 3                !  SCtJ.IiZ)   SIJ.J-I,?)  * ONOXrNOP
A 3                :  S(ItJi!i3)  = S(l,J.Ii3)  * QNOXrN02F
A s                :............,	,	
A 3                                      1
AS                  C   603             1
A S
	;  sii.j,:,")  " sii,j.i,i)  * IQHC *
                   J  F  VAPF5S>'< I , J)

                                         I
                                         I
                               /""""" '\
                              :          RETURN      :
                               \,................../
                                                                                                                                       CO
                                                                                                                                       x>
                                                                                                                                        I
                                                                                                                                       oo

-------
         B'C?!rLC.FOpn.O,S
         r LO,,C"AR70 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(2525),  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 8BFRAC
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
                                       FUOrCHAPTEt> 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<
                      "),  BCUH>, 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),ii,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.oiS506i/2j6oo :
i


i
i
: OVHK,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?HFRAC * .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! IDT :
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(NXNT)         !
                                      5), AHR(25,25), ALR(?S,25)
                   NSP=S

                   >F  =  I
                   00  2  L    1t  NY
                   JNT1-L
c::c's.            SEAD  5000,   UNOXIt,j>tJJNX>
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> , ti ,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) OTHFNO
OCCC32.           3<:,J,I,3>  =  S(IiJ,I,3)  *  ANOX t I , J J DTrlFN02
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, NAEI1)
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> , ii ,N*) /
i
',. * *
1 <-0 T .
I
I


t
!
: JNY-M-U ' :
i
i
/ RA!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,  PIN,  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.            EOU1VLCVCE  (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 /  28,  29,  310,  11  312,  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, 217, J1  /
  OATA UCCr / 28, 2*7, 310>  111  312, 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)  FN02ANOX(N)
                             3T."
                             C   1C3
                                                                                                                                                  I
                                                                                                                                                  oo

-------
                                                 RETURN
en
                                                                                                                                                     I
                                                                                                                                                     oo

-------
   oFt..t.,      -i
FS.07C"AR7EO BY  FoF>.0  /X803/  ON  13 OCT 71  AT a?:si:Z6
                  sy.acuf INE  POWER
                  CO--ON  ox,  DY,  or, pMjNto> PKAXJI 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-ISiON  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  /  o2,  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)  =  SUtJil3>  *  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
                               CC1KOJ 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(NP) > 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, FJ02,  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 Il 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) , lt ,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
JNYl-L
I
I
/ PRINT ACOOi IU1 ( I J) t 1I iNX)
I
C 23 I


I


I
I
: J=NYI-'-
t
i
/ P?1NT 6001. 1 VI ( 1 . J)  11 .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=                 :   OTKINiOTiOX/ABS(Ul(I.J)))             5
A B                 :	:
A3                                       t                                     I
A 5                                       0<	"	  ----------
A 3                                       1
AO                  /" "     ttff.   i t < fnttfi ,,iv  FALSE
A3                 <            IF  !YJ!IJ)  ,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  DTFAC
:          RETURN
 \
                                                                                                             I

                                                                                                            CO

-------
                                      3SR-844
            Appendix B



LISTINGS  FLOW CHARTS FOR NEXUS/P

-------
eroart.O'FOPFLOs    NEXUS
FLOC"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 NX2Z,N'Y=2i,"Z"i,Nsps,MAX750Q
PAPA"CTEP NX~1)X-I, NYMl"f.'Y-l, N Z .'. 1  N Z - I
PAI:A'/ETER NA = N'XYNZ, \TOT11AX7*NA
PARA^^'E" NA1NA*1 NA2^NA 1 .'.A , NAl = NA2*NA, NAlNA3*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
NXO?*3QO
NE/030CO
NE xoo i c o
NEX0320Q
NEX0330C
I (S, ARRAY ( 3477) ), ( OH , AR 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, \^,hKflT*JiVJr/l \n(,AKnttTlT|-IJ/IJ
II-TE^E"* DUMMY
CO">-ON /OUT/ isp(S), f.CYC, NOL, TDELT, PTIME, CONC, PAT,
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
2DT
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
UOPIZOSTU 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
NX06SCO
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(tCYC10)  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(-(NIOO) )
      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. PACL-JN-  3
                                                    I    CCELL METHOD  CODE  NEXCOMQO       3
                                                    I
                                                    I --- C PA'SA^ETE'? NX.22.NY-2) tN2H,NSP-3
                                                    I    CS,NAX9500 NEX00700             3
                                                    I    C PAP/^ETER NX11='X-1,  NYMlNY-l,3
                                                    I    CNZlI*NZ-t  NEXCC300              3
                                                    I    c PARAMETER NAoNxr,YNZ, NToT-1)* 3
                                                    I    CMAX7<"IA NEX0070C                3
                                                    1    C PA^CZTE1! NA,=\Ali  NA2-NA1NA,3
                                                    I    C .1 NAX*7NA+l,  NAYNAXMAX, NAZ-3
,                                                    I    CSAY*,I/>X , NA,,->MAZMAX 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,Xl,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,ATAY(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 2OY  NE705700
                                                        C 'TOZ  2OZ  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(PM 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)NX070DO                  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 PESTT 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
: NDSTOMN>  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(-tNioc>)           ;  :
                                                   1
w                              -
                                        /' ..... ....... m .t,N
is)                                     :          STOP         :
                                                                                                                                                en
                                                                                                                                                oo
                                                                                                                                                4s.
                                                                                                                                                fc.

-------
GFCpFLOF07FLO,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      tiBOH01flOO
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                                          BO02lOO
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  iVL  .LE. o) GO TO is                                          B0'2(I,I,")VEL
                   50t,QL='4,3                                                        BC0"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 |VL  -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 IFx
    FLU>2(J)=3.
    GO TO 200

    Fl_Uv/( J) "FXI
    DO 2^0 L=li3
    rX=-BY2(Ll,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    CBOS01803
                                                     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)VL :
: FLUXI ti j - FXI :
t
i


i
i
: FX  sxi ('.! , I ,"<; *vt.
: 3XS OXl(M|IiN) - FX1 :
I
C 3S3 I
CONTINUE :
: VEL  u( i ,fr,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(LiJM  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
: FX3Y2iL, 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.4tj370cijf'icAL  REACTION RATE VARJES SINUSOJDALLT  FROM 4:

   C17-.C.O
   Tl-ii7lM/3600.
   ;F(T'V-H.GT.18.S.O'?.TI'!H.LT.6.5) GO TO 5
   SI,>.C = s:N!ITI^H-65>, 16"01502)

   C:2=C23SINC
 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   CIY(Z)/(C2    CtY(3))
   Y(A! >  EO
   DI=C6Y(11 * C7-Y(2)
   D2=CSY(3)  C8*Y!1) *  C?*Yt7)
F2=(2.CH.Y(4)-Y(3)+C10Y(S)Y[3))/01

FU=C12*T(51/D2
012=1.0-F1F3
EP02= (F2*F"iF 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
                                                                       ChOCHOO
                                                                                         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
                                                                                               CH03303
                                                                                               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  C3Y(t)  C6Y(7) * C8Y(8) + Ctl*T(2)
3 ; 1 )  R.'- 0
RN02  Cl * C7YI7)  C7Yt8)  ClMYtl)
8(2)  PN02
PMC " C1Y<6> * C5YI8) * ClOY(1)
3(3)  P"C
P03  C3Y( I > * CIOY<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 XTEST1.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  ClYS(2)/(C2 * CfYSC3n
TSI6) - EO
01=C6YS( 1 )  C7YS(2)
D2=ci'TS(3)  C8Y3tl)  C'*YSl2)
F ;=>2- 'C3>VS (31/01
F2'-(2."C'"YS(6)YS(3* + CJO'YSt1lYS(3))/Ol
F3=3.5'C*YS( 1 1/32
F1=C12*YS(5)/02
C.Z= 1 .C-FI F3
-02= (F2-rHFl 1/D12
Y3 !7 1  "-32
0J > (F2F3 + F1 1 /OI 2
TS(O) = ECH
' .0  C1-YS(2)  C12YS<5)
t- ( : ; = F c
F .52  C3"S(tll YS( 1 1 * C6YS< 1 )YSI7)
/ ! 2 1 = F 32
F03  C2-YSI61
A(M) = F03
rt"02  C9'vS(aiYS( I 1  2.CI IYSf HYS(2
CHH05200
CMM03300
CHM0500
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
ClM0900
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.
                                                *  C8YS(8)
c:3*35.
o:ci36.
 roi 3s.
ccciz>
c;o "3.
(* f- r* 14 y
ccoi-'S.
C3C116.

OOC1MS.
                         - C3TS(t)  C<
                         = F ^0
                          Cl  * C?YS(7) C9YS<8>  * C11Y3(M
    RNO
    B! 1 )
    R .02
    612)
    ^f-C    C1'YS(6)    C5YS(8)  CIO*YS(1)
    0(3}   RHC
    03   = C3YS(1)  *  C10*YS(3)
    0 (S )  = P03
    B 13.)  = C 1 2
    DO 30 I =  1 ,  S
    IF Oi:>D*  .oT.  lE-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.  .1ER)  C   1ER
    07  f I JIE"?/CDTi  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
                                                                                       Cfl 1200
                                                                                       CH11 13CO
                                                                                           1SOO
                                                                                          l 1500
                                                                                       CHM i : 400
                                                                                       C^-M 1 i!3
                                                                                       Curl lb-50
                                                                                       CHM1 1900
                                                                                       CHM 2030
                                                                                       CM122CO
                                                                                       CHMI2300
                                                                                       CH12HOO
                                                                                       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,C2C3,c-<,cs,c6,c7,C8,c?,cio,ci i,  :
  C20/.",  2-6SE6, "0., 61CO, 80.,  1SOC.,
  4.,  10., 30.. .0125, .Cl, COl /
  OlwNS|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
                       ic 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
      SlNCSiNUTlKH*6.S>. 1*981582)                                                                           Z
                CI=COSI ,'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    Clr(2)/(C2  C1

                                             Y<6>   EO

                                       DJC6Y(1)   C7Y(2)

                                  02-C!iY(3)  +  CS-YU) * C9Y(2)
NJ
F2<2.CHY!6)Y(3<+C'0'T("t)Y(3l)/DI

           F3^C.S-Ct "t ( 1 ) /02

            FC1:*1< S)/D2
             0'.2'1'0-F1F3

         ER02  (C2FMFI)/012

              Y(7) .

        EOH - (F2-F3  >

              Y (8) =>  EOH

       FNO  C1->Y<2)    C12Y!5)

              A (!) "  F MD

FN02  C3Y(H)-ril>  *  C6r(j)Y(7)

              A ( 2 ) =  F t. 0 2

             F03  = C2r,'i)

              AIM) =  F03

FHNOZ  C8Y(8)Y(1)    2,C1l.Ytl)Y(2)

             A<5)   FKVOZ

>?NO  C3Y(")  CA*Y(7)  * C8Y(8)  +  Cll
Y(2)

              3(1) =  PMO

RN02 - Cl  C7Y(7>    C?Y<8)  *  Cll'YIl)


              SI2)   RN02

PHC t. CtY(6)  CS*r(3)   CIOY(S>

              8(3)   RHC

       R03  C3Y(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  E7 ) 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  C1YS(2)/CC2  C"*YS(3M
YS<6>  EO
OlC6YS( i ) * C7YS(2)
02CSYS(3) * C3'YS(l! - C?YS(2>
Fl2,c5YS(3)/Dl
F2(2CHoYSlA)Ys!3)*C!0YSt')>YS(3>J/
01
FS^C.S'C^'YSt 1 )/02
FM-C12*YS(S) /B2
Oil" .Q-F' F3
ER02* tF2*F'tF J ) /D]2
YS(7 ) = E.TC2
EOH " (rZF3  r1)/012
YS ( 3)  EOH
FNO = CJYS(2)  C12'YS(5)
All) = FNO
FN02  C3YSCn vsl 1 ) * CA'/5 C 1) YS ( 7 )
A (2) = FN02
F03  C2Y5(6)
A (u ) = F 03
FHN02  C8YStS) *Y3( 1 I  2  *C 1 1 YS ( 1 ) YS (
2)
A(5)  FHN02
R^^  C3YS1M)  C4Y5(7)  CO'YStB) +
Cl !YS(2)
Bl 1 ) = R^O
1 '02  Cl  C7YS(7) -c?YS(8) * Cll'YSI
; i
0 (?) = R N o 2
FHC B C"YS(6) * CSY3() + CIO'YStH)
SI3) = -'HC
?03  C 3  r 3 ( 1 ) * C10YS(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/COT, TT - TJKE)
I
Z



C 9003 I
CONTINUE :
: oo  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. 0l               > ---- - 
 \ ..................................... ..../
                      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

-------
SFOsFLO.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,A9AYIiat9)), (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
NPSTOR5) , 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 , 11 ,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)   UiJ.
                    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
                                                                     01FC300
                                                                     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),11.NX),   JM.NYlt  U(2,  /
                    / j>: .31, i = i ,HXI ,ji ,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.|iLJ3T
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.DTU(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   COIFOU300                        ]
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)  imJ,t .ZJ/OY :
1
I



C 83 I
CONTINUE
CVX o DH(KtJ - C(KfJ,l*m
/(CtK t J J JTOX)
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 UUJili2)/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
IN023CO
, 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) ) , ( tT , 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-6OK3Tor
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
Ifj03600
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 INUT/OUTPUT.     3
                                           i   COCFAULT VALUES A-JE SUPPLIED FOR  3
                                           I   C I JFCC3CC                          3
                                           i   c ALL 'cCESSA"?" PAAKETERS.       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', DM,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"-C7FLO.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      KEccrco
            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,NAINZ,           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 ! ) ,         KEMOi00
                  1 (Si A7RAYI3677) ) i  I OM , ARi 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'ANSP!             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 (CFi Ji I fN) ,K1 ,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
                                                                                      KEf0200
                                                                                      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  " NANSPI        3
                                               E.              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> >X1 .fvZ) t J"l 
/ NY),ll,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>" 3TI5> , "(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 Kl ,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, "',
ED
'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
NTo:eoc
NPTC17CO
NPIC2030
NPT02 ICO
NPI02200
NP702303
NPT02"33
NPIC25BO
NT32403
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: :


:
: RETUN :
\..... 	 /

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

-------
        SFCFLDrO?ri.O,S   OUTPUT
        FLOrCM-STED 6Y  FoFLO /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, RSTT
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 ( -0R) 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 ..ME3I 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 ) ,il ,NP>
    ; ( T( I > i !*! ,NS1
 IZc'lF (.NOT.  CONC)  GO TO  130
     -*"[T: <2I  T, ICTC, C         .
 130 PT^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    CCONCENTAT 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"PY(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
                                COVON /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-
                                                         COuT323CO
                                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) NOR  .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,  ',  LVrLi,  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                :              JNY-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 ) , 11 |NP) , (WT( I ) ,11 ,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/ = ROUTlE "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, AfAY 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.
                                     1FY
                                                                                                PAKOblCC

-------
CO
                           F7FZI*FX1F1                                                       PA"?05200
                           F8=F:i'FXFYl                                                        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 * UlKiJIt,3)'F2  +  UfX,Jltll3)F3         PAR07HOO
  1 V/(X i Jl . I |3)'FM   U|K 1 > Jt I ,3) F5  UtKI , J,I 1 .3) F6                PK075CO
  2 t'lc I , Jl , I 1 ,3 | F7  *  U(Kl , Jl , I ,3) FO                                 PAKO/400
   GO  TO 50                                                             PAK07700
HO Vt_(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)  *  VL(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
FY1 .-FT!
f z=i .-rzi
F1='ZFXTY
F2=FZFX!FT
F3-FZF^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 'FXFY:








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 fMJF4  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)*(F2F6)  Ulli
J1!I.3)(F3F7)  UlliJlIi3!(F1FO
I
[ 	 c MOVE PARCEL PAR08200 3
I
C 503 I
X(UI  T (L) + VL ( 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
C3ic\ /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( ! ) , Il ,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
  \             SU3OUTI.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 , NX , 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 ? ) ) ,

    ARAT(55t5)),(U,ARKAY(73?3))>  (X.ARPAYC
    12737)), ( Y,A-?RA?;2237) ) ,  (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(ij .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  Nsp2  (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 lJ /
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 SUSS  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





/ HJTE (29) T, JCYC, PMNSTR, PKXSTK, /
/ v.:iS7CK( 3CKSTR, NPST01 /
I
I


i
M2CN : *
!

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
                           2FCTI3>,  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 ,NZ1 ,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 , ii ,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 ) ,I1 ,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,CrtCuARTD  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 ,           STPCC00
                          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,PAr <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                                                       S7P0200
                     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--1u7Eo RANDOM NUVBE^ GENERATOR COVER THE  INTERVAL FROM o. TO  i )                                                   *"
                           xi.'.p!    :  - .5     uAiiotcuv*;Y)                                       S7PO'-303
                            C")  =  J  - -5   *  V,Af D(OJ^.1<')                                       ST0'(103
                           :t.P)    K  - ,5     URAf,D(OLv>1Y)                                       S7POS3CO
                     C     lr "APCELS  FALL  OUTSIDE  TME  GRlo  LIMITS  THEY  A^E  THROWN AAT.      STP016CO
                           ir'Zt  P).CE.F"Z,Or>.XI',P).GE.Fl,X.OR.Y(NP).GE.FNY.O^.X(NP).LE.lOK.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 EXOS'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 J1 > 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?) *FUG6./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
ST07?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 THKON AAT3
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 TJS 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 - oAcxiFUG>6./
: F A c 7 1 N ) :
i
C 7] I
	 ._: arz(L,i,N) - ten, NY, n - 3ACKiFUG6./ :
: Ffcrt n :
i
i
	 -> 00 8 J-> 1 ,'!T :
i
i
: SYllL.J.NJ - (Ctl.J.l) - 9AC^JFUG6./ :
: FACT (r ) :
I
C 03 I
- 	 _-- 3Y2(LiJN) " (C(1|J|"X) - BACK) FUG6./ :
: 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(313D ) , ( '> 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 A0 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-.SUPAND(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 ) .LE0 )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-.5L!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 CCTI"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".YNP)  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
                                                                                               SC13600
                           cr
                                                                                               SRC13900
                                                                                               SRCISOOO
                                                                                                                                                  CO
                                                                                                                                                  oo

-------
                                                      8?   FOTLO /X308/  ON'  J3 OCT 71 AT  0.'53:i6
                                                      50Ui~CEJ CF TMC CELU    3
                                                           CTE  3RCC05CO                      3
                                                           :  .  Pt^cELisi  WILL BE PUT INTO.3
                                                           CSKCCC6CC                          3
Co
CO I'M ON  OX,DY,?Z,TOX,TOYiT5Z.r.-lIN,PAX,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 PV
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
                                                        CPACE'.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
                                                        C3COS6C'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<3XFLv/.''N.S
:
I


I FALSE
I
I
: PA=XFL*/\O :
i
t
----- > * OC 1 'C Ll NO "

1
^Vp^,'+ |
fi~  NP)=?,-.
XC.PJ'I I-.i + O-A'.'DiOUKHY)
YC P ) = 1 . J?A. ")(DUM11'),S
Ztfll*.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 1C3 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
                                                    ic Y  >JY E?GE OF Gdo  -  PARCELS   3
                                                    I    CAE 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)j2.
,'< 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<>*CiJJi)+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

NpNP-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-.5Ufi/* 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*CHART0  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    CUD;3300                          3
                                                   I    C  TV,; CTCLE LENGTH IS 8!?2        3
                                                   !    CnU"3CRS OEFOSE REPEATING         3
                                                   I    CUTOC'-CC                          3
                                                   t    C  !NTECEPH NN UROOOSOO           3
                                                   1
                                                   I
                                      NN   MODtHN60S9.32768)
                                  U?ANO  FLOAT(NN)    .C0003Q5176
w
                                                                                                                                                OJ
                                                                                                                                                00
                                                                                                                                                po

                                                                                                                                                OO

-------