NBSIR 80-2157

?30713

A Computer Program for Analysis

of Pressurized Stairwells and

Pressurized Elevator Shafts

John H. Klote

Center for Fire Research

National Engineering Laboratory

U.S. Department of Commerce

National Bureau of Standards

Washington, DC 20234

January 1981

Final Report

f-QG

100

I.U56

80-2157

1981

U.S. DEPARTMENT OF COMMERCE

NATIONAL BUREAU OF STANDARDS

NBSIR 80-2157

» * > A COMPUTER PROGRAM FOR ANALYSIS

OF PRESSURIZED STAIRWELLS AND

PRESSURIZED ELEVATOR SHAFTS

National Bureau of Standards

Library, E-01 Admin. Bldg,

FEB 2 7 1981

f\Oi &Ct- - Clooo

S'o fQSt

John H. Klote

Center for Fire Research

National Engineering Laboratory

U.S. Department of Commerce

National Bureau of Standards

Washington, DC 20234

January 1981

Final Report

Sponsored in part by:

The Department of Health and Human Services

Washington, DC 20201

U.S. DEPARTMENT OF COMMERCE, Philip M. Klutznick, Secretary

Jordan J. Baruch, Assistant Secretary for Productivity. Technology, and Innovation

NATIONAL BUREAU OF STANDARDS, Ernest Ambler. Director

PREFACE

This report is an interim product of a joint effort of the Department

of Health and Human Services and the National Bureau of Standards (NBS)

,

Center for Fire Research. The program is a multi-year activity initiated in

1975. It consists of projects in the areas of: decision analysis, fire and

smoke detection, smoke movement and control, automatic extinguishment, and

behavior of institutional populations in fire situations.

This report describes a computer program which analyzes pressurized

stairwells and pressurized elevators. The program was initially intended as

a research tool to investigate the feasibility of specific systems. However, this program may be of interest to design engineers responsible for

pressurized stairwells or pressurized elevators.

iii
TABLE OF CONTENTS

Page

Preface iii

Abstract 1

1. INTRODUCTION 1

2. PROGRAM CONCEPT 2

3. ASSUMPTIONS AND LIMITATIONS 3

4. EQUATIONS 3

5. PROGRAM DESCRIPTION 7

5.1 Main Program 7

5 . 2 INPUT Subroutine 8

5.3 CORR Subroutine * 10

5.4 INIT Subroutine 10

5.5 BLDGP Subroutine 10

5.6 SHAFTP Subroutine 10

5.7 PZAD Subroutine 11

5.8 OUT Subroutine 11

5.9 INPUT2 Subroutine 11

6. FUTURE DIRECTION 11

7. REFERENCES 11

APPENDIX A. DATA INPUT DESCRIPTION A1

APPENDIX B. INPUT EXAMPLES B1

APPENDIX C. EXAMPLE OUTPUT Cl

APPENDIX D. PROGRAM LISTING Dl

iv
A COMPUTER PROGRAM FOR ANALYSIS OF PRESSURIZED

STAIRWELLS AND PRESSURIZED ELEVATOR SHAFTS

John H. Klote

Abstract

Pressurized stairwells and pressurized elevators

can be used as a means of providing a smoke free exit

route during fire situations. This paper describes a computer program which analyzes systems intended to

pressurize stairwells or elevator shafts. The basic

assumptions and limitations of the program are also

discussed. The appendices contain a program listing

and examples.

Key words: Air movement; computer programs; egress;

elevator shafts; escape means; modeling; pressurization;

simulation; smoke control; stairwells.

1. INTRODUCTION

Pressurized stairwells are being used as a means of providing smoke

free exit routes during fire situations. A discussion of several designs

for pressurized stairwell systems is provided by Benjamin and Klote [1] . In addition, the concept of pressurized elevator shafts has gained considerable interest as a means of fire escape for the handicapped. The purpose of

the program described in this paper is to calculate the pressure differentials

produced by such systems. Because a stairwell or an elevator shaft is

connected to a building, the program calculates air flows and pressure

differentials throughout the building in order to obtain the differential

pressures across the stairwell or elevator.

A number of computer programs have been developed which are applicable

to smoke control. Some of these programs calculate steady state air flow

and pressures throughout a building [2,3]. Other programs go beyond this to

calculate smoke concentrations throughout a building that would be produced

in the event of a real fire [4-9]

.

In general, most of these programs are

capable of analyzing pressurized stairwells and pressurized elevators.

However, the program described in this paper has been specifically written

for analysis of pressurized stairwells and elevators. The data input has

^"Numbers in brackets refer to the literature references listed at the end of

this paper.

1
been designed to minimize the quantity of required data and still maintain a high level of generality in the model. The output consists of the differential pressures across all of the building shafts in addition to the flows

and pressures throughout the building. In addition, the assumptions and

limitations of the program discussed in the next section were chosen specifically to allow modeling of pressurized stairwells and pressurized elevators.

This program was intended primarily as a research tool to investigate

the feasibility of specific systems and to determine the interaction between

these systems and the rest of the building. This program has already been

used to analyze pressurized stairwells without vestibules and to evaluate

factors which affect the performance of these systems [10] . It is also

possible that the program may be used as a design tool directly or to generate

quantities of design data. This paper is not intended to be a design guide

for either pressurized stairwells or pressurized elevators. The state-ofthe-art of these systems is still developing and designers of these systems

should seek the most current data available.

2. PROGRAM CONCEPT

In this computer program a building is represented by a network of

spaces or nodes each at a specific pressure and temperature. The stairwells

and other shafts are modeled by a vertical series of spaces, one for each

floor. Air flows through leakage paths from regions of high pressure to

regions of low pressure. These leakage paths are doors and windows which

may be opened or closed. Leakage can also occur through partitions, floors,

exterior walls and roofs. The air flow through a leakage path is a function

of the differential pressure across the leakage path.

In this computer model air from outside the building can be introduced

by a pressurization system into any level of a shaft or even into other

building spaces. This allows simulation of stairwell pressurization, elevator shaft pressurization, stairwell vestibule pressurization and even elevator lobby pressurization. The pressures throughout the building and flow

rates through all the flow paths are obtained by solving the air flow network including the driving forces such as the wind, the pressurization

system or an inside outside temperature difference.

2
3.

ASSUMPTIONS AND LIMITATIONS

1. Each space is considered to be at one specific pressure and one specific

temperature

.

2. The flows and leakage paths are assumed to occur at mid-height of each

level

.

3. The net air supplied by the air handling system or by the pressurization

system is assumed to be constant and independent of building pressure.

4. The outside air temperature is assumed to be constant.

5. The barometer pressure at ground level is assumed to be standard

atmospheric pressure (101325 P ) . The results of the program are not

3. very sensitive to changes in atmospheric pressure. For altitudes

considerably different from sea level the more accurate value can be

substituted by changing an assign statement in the subroutine INPUT.

4.

EQUATIONS

A. Flow equation

F = C (AP) x (3.1)

where

F = mass flow rate

C = flow coefficient

X = flow exponent

AP = differential pressure across flow path

The following equation has been used by Sander and Tamura [3] for

2 corrections for flows not at standard conditions .

C - C

o

<£->

X 1 —2 X

o where

C = flow coefficient at standard conditions

o

p = density of entering air

2 This relationship was originally derived by Tamura and Wilson [11]

.

3
P Q = density of air at standard conditions

U = dynamic viscosity of entering air

UQ= dynamic viscosity of air at standard conditions

This author has observed that the dynamic viscosity is proportional to

absolute temperature to the 3/4 power in the range of -40°C to 60°C (-40°F

to 140°F) with a maximum error of 1.2 percent. Using this and the fact that

density is inversely proportional to absolute temperature the above relation

reduces to a function of temperature only.

where

:

T = absolute temperature of entry air

T = standard absolute temperature

B . Mass Balance Equations

3 . For building compartment l

N N

E F d,j)

+ E

j=i k=l

F ,. . . + F - ... =0 o(i,k) f(i) (3.3)

and for shafts

i=N.

N N

S F

(i,j)

+ ]C Fo (i,k)

+ F

f(i]

j=l k=l

= 0

where

:

(3.4)

F

(i / j ) mass flow rate from space j to space i. For building

compartments this flow can be either horizontal or vertical,

however for shafts this flow can only be horizontal.

o (i ,k)

mass flow rate from direction k outside of the building to

space i.

3 In this paper the term building compartment refers to a space in a building

other than in a shaft.

4
Ff ( i) = net mass flow rate of air due to the air handling system or

due to a pressurization system.

N

c = number of building spaces connected to space i.

N = number of connections to the outside from space i. o

is the space number at bottom level of shaft and the spaces in the

shaft are numbered consecutively up to ^ which is the space number at

the top of the shaft.

C. Shaft Pressures

The following relation is used to calculate the gage pressure, P

(i)' at

level i of a shaft in terms of P^_^ at level i - 1.

P

(i)

P (i-D P

z - P, ( 3 . 5 )

where

P = hydrostatic pressure difference

P^ = pressure loss due to friction

The following equation is used to calculate the hydrostatic pressure.

= 2Z

RT

h d)

" h d

( 3 . 6 )

where

h ^

= height of point i

h^_u = height of point i-1

g = gravitational acceleration

R = gas constant - _ T

(i) + T(i-1)

_

P

(i) + P (i-1)

+ P, 2 b

P^ is a constant used to convert an average gage pressure to the average

absolute pressure P.

The following equation is used to calculate the pressure loss due to friction,

5
( 3 . 7 )

where

:

= upward flow from i -1 to i in shaft

C = shaft flow coefficient

s S = sign of F ^ u

D . Outside Pressures

Outside pressures can either be entered by the user or can be calculated by

the following method.

P / ' \ ~

P-U /

• \ + C P , . » 0(1) h(i) w v(i)

where

:

( 3 . 8 )

PQ

^

= outside gage pressure at height h(i) above absolute pressure at

ground level.

Ph(i)

= hYdrostatic Pressure difference between h(i) and ground level

Pv ^j = velocity pressure due to the wind at height h(i)

C = pressure coefficient

w

Because the outside temperature is constant

where

:

P atm

= ab s °lute barometric pressure at ground level

T Qut = outside absolute temperature

When the outside pressures are calculated by the computer the wind

velocities are assumed to be described by the power law.

V = V {£-) n

o h

o where

:

V = wind velocity at height h

o o

n = wind exponent

6
This relationship has been extensively used to describe the boundarylayer velocity profile of the wind near the surface of the earth. This

equation assumes that the terrain surrounding the building is homogeneous.

That is, that there are no large obstructions near the building which could

produce local wind effects. A value of 0.16 for the wind exponent is appro

priate for flat terrain. The wind exponent increases with rougher terrain.

For very rough terrain such as urban areas a value of 0.40 would be

appropriate

.

The equation for the velocity pressure at height h(i) is obtained by

substituting the velocity from the power law into the usual relation for

velocity pressure (P v = jpV 2

)

.

where p is the outside air density.

The pressure coefficients are in the range of -0.8 to 0.8 where positive

values are for windward walls and negative values are for leeward walls.

The pressure coefficient depends upon building geometry and varies locally

over the wall surface. Numerical values for and n as well as practical

engineering information is available from a number of sources [12-15] .

5. PROGRAM DESCRIPTION

This program is written in the computer language FORTRAN V specifically

for use on the UNIVAC 1108. A program listing is provided in appendix D.

The following is a detailed description of the main program and the major

subroutines

.

5.1 Main Program

This program was written to handle four different run types as follows:

1. RUN: In this run type subroutine INPUT is used to read data and then

the pressures and flows are calculated and printed. After this, the

subroutine INPUT2 is used to read modifications to existing data or to

determine if execution should be terminated. Under this run type, if

modifications were read then the pressures and flows are again calculated and printed. Any number of modified data sets can be determined

and solved in this manner.

7
2 . RUN 1: This run type is similar to the one above except that no

solutions for pressures and flows are produced for modified data read

under the subroutine INPUT2 . A solution is obtained for the original

data read under subroutine INPUT.

3. TEST: In this run type data are read and tested but no solutions are

produced. Run types 2 and 3 can be used to test input data which are

to be used later under another run type and under a low cost priority.

4. FIND F: In this run type, the rate of pressurization air necessary to

produce a specific pressure differential across a shaft is determined.

Figure 1 is a flow chart for the main program logic under the conditions

of a run type RUN.

5 . 2 INPUT Subroutine

This routine reads the data that are necessary for a flow analysis of

the stairwell or elevator, including an analysis of the rest of the building.

These data consist of the following:

1. Outside temperature.

2. Temperature throughout the building.

3. Outside pressures. These can be entered or calculated as described

earlier.

4. Description of the flow network including flow coefficients and flow

exponents for all connections and the net air flows to each space due

to the air conditioning system or due to a pressurization system.

The data above can be entered in either SI units or in engineering

units. There are two modes for the description of the flow network. The

first is a general mode which allows any building geometry and any connections between building spaces that the user chooses to define. The second

description mode is simplified in that each floor is represented by one

building space and the flow coefficients and exponents are all the same in

the vertical direction unless they are defined by exception to be different.

This simplified data input greatly reduces the quantity of data required to

define a large building. Appendix A contains a detailed description of both

of these methods of input.

8
START INPUT Subroutine

Read and Print Data

Initialize P Terms

z

CORR Subroutine

Calculate Flow Coefficients

Corrected for Temperature

Determine if Another

Solution is Desired

With Modified Data

Figure J . Flow chart for main program logic under a run type RUN.

9
In addition to reading data, this subroutine provides temperature and

pressure data as well as a complete description of the flow network. This

routine also calculates initial estimates of the hydrostatic pressure

differences. When data is entered in engineering units the subroutine UNITS

is called which converts all units to the SI system.

5.3 CORR Subroutine

This routine calculates flow coefficients corrected for the temperature

of the entering air using eg. (3.2). Two sets of these corrected coefficients

are calculated for each flow path to allow for flow in either direction.

5.4 INIT Subroutine

This routine calculates initial estimates of the building pressures by

a technique used by Sander [2]

.

In this technique, mass flows are considered

linear functions of differential pressure and therefore the flow equations

can be expressed and solved in matrix form. In this estimate, shaft pressures

are considered hydrostatic. The resulting pressures form a starting point

for the iterative solution which follows.

5 . 5 BLDGP Subroutine

The iterative solution for the building pressures and flows consists of

the three subroutines BLDGP, SHAFTP and PZAD. The subroutine BLDGP operates

on the building compartments sequentially. The sum of all the mass flows

into compartment i is calculated. If the absolute value of this sum is less

than a convergence limit then eq. (3.3) is considered satisfied and the

computer proceeds to the next compartment or returns to the main program.

However, if the absolute value of the sum is greater than the convergence

limit, then an improved estimate of the pressure at compartment i is obtained

by the regula falsi method [16] . When none of the pressures need to be

modified this routine passes a convergence signal to the main program.

5 .

6

SHAFTP Subroutine

The structure of this routine is very similar to that of BLDGP except

that it operates on shafts sequentially. The sum of all the mass flows into

shaft i is calculated. If the absolute value of this sum is less than the

convergence limit then eq. (3.4) is also considered satisfied and the computer proceeds to the next shaft or returns to the main program. However,

if the absolute value of the sum is greater than the convergence limit, then

improved estimates of the shaft pressure are calculated. This is done by

10
changing the pressures at the bottom of the shaft and then recalculating the

shaft pressure by eq. (3.5) . Again the regula falsi method is used, and if

none of the shaft pressures need to be modified a convergence signal is

passed to the main program. It can be seen from figure 1 that if convergence

is achieved in both BLDGP and SHAFTP, then the subroutine OUT will print the

solution. Otherwise, the hydrostatic pressure differences are adjusted in

the subroutine PZAD.

5.7

PZAD Subroutine

This routine calculates hydrostatic pressure differences by eq. (3.6)

using the most recent pressure estimates.

5.8

OUT Subroutine

This routine outputs mass flows and pressures for the flow network as well as the differential pressures across each shaft. If the data input was

in engineering units then the subroutine UNITS is called to convert variables

to the engineering system before output.

5.9

INPUT2 Subroutine

This routine reads modifications to the flow network. The modifications

can consist of changes in the net flow to a stairwell or elevator pressurization

system or of changes to flow coefficients between either between building

spaces or to the outside. After the modifications are read, a solution to

the new flow network is obtained.

6. FUTURE DIRECTION

In the future it is planned to analyze pressurized stairwell systems

with pressurized vestibules as has already been done for systems without

vestibules [10] . In addition, the use of this program is a part of a large

project to study the feasibility of using pressurized elevator shafts for

fire evacuation of the handicapped. The further development of this program

so that it can easily be used as a design tool is being considered. However,

this development would depend upon the determination that a real need exists

for such a design tool.

7. REFERENCES

[1] Benjamin, I. A. and Klote, J. H., Stairwell pressurization systems,

Nat. Bur. Stand. (U.S.), NBSIR 79-1749 (June 1979).

11
[2] Sander, D. M. , FORTRAN IV program to calculate air infiltration in

buildings. National Research Council Canada, DBR Computer Program No.

37 (May 1974) .

[3] Sander, D. M. and Tamura, G. T. , FORTRAN IV program to simulate air

movement in multi- story buildings, National Research Council Canada,

DBR Computer Program No. 35 (March 1973)

.

[4] Yoshida, H., Shaw, C. Y. and Tamura, G. T., A FORTRAN IV program to

calculate smoke concentrations in a multi-story building. National

Research Council Canada, DBR Computer Program No. 45 (June 1979)

.

[5] Butcher, E. G., Fardell, P. J. and Jackman, P. J., Prediction of the

behaviour of smoke in a building using a computer. Symposium on movement

of smoke on escape routes in buildings, Watford College of Technology,

Watford, Herts, England (1969) pp 70-75.

[6] Barrett, R. E. and Locklin, D. W. , A computer technique for predicting

smoke movement in tall buildings, Symposium on movement of smoke on

escape routes in buildings, Watford College of Technology, Watford,

Herts, England (1969) pp 78-87.

[7] Rilling, J., Etude Sur Le Desenfumage - Methode De Calcul Du Mouvement

Des Fumees Entre Volumes d'un Batiment, CSTB, France (1978)

.

[8] Evers, E. and Waterhouse, A., A computer model for analyzing smoke

movement in buildings, BRE Fire Research Station Borehamwood, Hertsfordshire, England (Nov. 1978)

.

[9] Wakamatsu, T., Calculation methods for predicting smoke movement in

buildings fires and designing smoke control systems, Fire Standards and

Safety, ASTM STP 614, A. F. Robertson, Ed., American Society for Testing

and Materials (1977) pp 168-193.

[10] Klote, J. H. , Smoke control by stairwell pressurization. Engineering

Applications of Fire Technology, SFPE and Nat. Bur. Stand., (April

1980) .

[11] Tamura, G. T. and Wilson, A. G. , Pressure differences caused by chimney

effect in three high buildings, Appendix A, National Research Council

of Canada, Research Paper No. 346 (1968)

.

[12] Sachs, P., Wind forces in engineering, Pergamon Press, New York (1972).

[13] Houghton, E. L. and Carruther, N. B., Wind forces on buildings and

structures, John Wiley & Sons, New York (1976).

12
[14] Simiu, E. and Scanlan, R. H., Wind effects on structures: an introduction

to wind engineering, John Wiley & Sons, New York (1978).

[15] MacDonald, A. J. , Wind loading on buildings, John Wiley & Sons, New

York (1975)

.

[16] Carnahan, B., Luther, H. A. and Wilkes, J. 0., Applied Numerical

Methods, John Wiley & Sons, Inc-, New York (1969).

13
'

'

j

!
APPENDIX A DATA INPUT DESCRIPTION

APPENDIX A. DATA INPUT DESCRIPTION

Data input consists of the following elements:

1. Initial data

2. Building heights

3. Temperature profiles

4 . Outside pressure profiles

5. General flow network input

6. Abbreviated flow network input

7. FIND F data

8. INPUT2 data

Each of these input elements is described in detail in the following

sections. Elements 1 through 4 are always required. The flow network can

either be defined by a general method of input using element 5 or by an

abbreviated method using element 6. The virtue of the abbreviated method is

that less data are required. In the abbreviated method each floor is represented by one space and flow coefficients and net flows are all the same

vertically. However, the abbreviated flow network input element also allows

a number of specific exceptions to this simple model.

Element 7, FIND F data, is used for a FIND F run type. Element 8,

INPUT2 data, is used to modify the existing data for another analysis. In

the following sections the input required for each of the eight data elements

is described in detail. Each block or group of blocks below represent an input

card. Unless otherwise stated these cards are unformatted, that is the

numbers do not have to be placed in specific columns of the card. However,

separate pieces of numerical data must be separated by one or more spaces.

Examples of input data are provided in Appendix B.

1 . Initial data

run type^ either RUN, RUN 1, TEST, or FIND F (col. 1-6 must start in

col. 1)

project title (col. 1-72)

A1
outside unit indication summary output

temperature (°C,°F) (1 for SI, 2 for Eng) (0 for none, or file number)"^

2 . Building heights

N. , no. of building

levels

input parameter

(either 0 or 1)

If input parameter = 0, then heights for each building level are to be

individually inputed as follows:

( 1 ) h

(2)

h

( 3)

h

(i) (V

o o o

where is the height of the center of level i above the ground (m, ft)

If input parameter = 1, then the following card must be entered.

h,^. distance

between floors (m, ft)

3 . Temperature profiles

no. of temperature profiles

^All net flows are inputed at standard conditions of 21°C (70°F) and 1

atmosphere

.

2 The user must assign this file before program execution.

A2
For each temperature profile the following data must be supplied.

no. of temp,

points

level temperature level temperature

no. (°C, °F) no. (°C, °F)

level temperature

no. (°C, °F)

O

4 . Outside pressure profiles

N

po

no. of outside

pressure profiles

input parameter

(either 0 or 1)

If the input parameter = 0, each outside pressure profile is entered as follows

:

P P P

o(l) o(2) o(3) Po(i) o(Nh

)

where p

Q (i)

i- s the outside pressure at level i.

If the input parameter = 1, the outside pressures are calculated and the

following data are required.

V

wind

velocity

h

height at which

velocity is measured

n

wind

exponent

pressure coefficients for each pressure profile

CW(1) CW(2) CW(Npo)

o o o o

A3
5 . General flow network input

5 . 1 General building data

N

no. of building input parameter

compartments (0 signals non- abbreviated input)

For each compartment the following data must be supplied.

N

eC no. of connections to

other spaces

N

j-O no. of connections to

the outside

F

f net flow

( l/s

,

cfm)

temperature

profile

number

compartment

level

For each connection between this compartment and other spaces in the building

the following connection data are required.

number of CQ 3 X

connected flow coefficient"3 flow exponent

space

For each connection to the outside the following connection data are required.

outside pressure

profile number

C

flow coefficient

X

flow exponent

5.2 General shaft data

no of shafts

3 All flow coefficients are entered for standard conditions of 21°C (70°F) and

1 atmosphere. The units for all flow coefficients are Us -1 (Pa) or cfm

(”H 2 0)

1

" X

. A4
For each shaft the following data are required.

shaft title (col 1-18)

C

shaft flow

coefficient

bottom

level of shaft

top level

of shaft

temperature

profile

number

For each level of this shaft the following data are required.

no . of connections

to other spaces

no . of connections

to the outside

net mass

flow (£/ s, cfm)

For each connection between this level of the shaft and other building

spaces the following data are required.

no. of other C

Q

X

spaces flow coefficient flow exponent

For each connection between this level of the shaft and the outside the

following data are required.

outside pressure

profile number flow coe¥ficient

X

flow exponent

A5
6 . Abbreviated flow network input

6 . 1 Specific building data

N

no. of floors input parameter

in building (1 signals abbreviated input)

The following data are entered only once because it is the same for all

floors of the building except at locations where specific exceptions are

defined

N

eO

no. of connections to

the outside

no. of the

temperature

profile

flow c8efficient

through floors

X

flow

exponent

through floors

no. of

exceptions

For each connection to the outside the following data is required.

outside pressure C X

profile number flow coefficient flow exponent

For each exception the following data is required.

exception type

(1, 2 or 3)

building level

for the exception

The next card depends upon the exception type. For exception type = 1, a net flow is defined for this building level.

F

f net flow for

building level (£/s, cfm)

A6
For exception type = 2 , an exception to an outside connection is defined.

no. of outside C X

pressure profile flow coefficient flow exponent

For exception type = 3, a new compartment is added at the building level.

In this case the following cards are needed to define this new compartment.

C

flow coefficient

between new compartment

and rest of the floor

X

flow

exponent

no . of connections

to the outside

Each of the outside connections is defined as follows.

outside pressure C X

profile number flow coefficient flow

exponent

6 . 2 Specific shaft data

no. of shafts

For each shaft the following data are required,

shaft title (col 1-18)

C

g

bottom level top level temperature

shaft flow of shaft of shaft profile number

coefficient

A7
The following data are entered only once because it is the same for each

level of the shaft except at locations where specific exceptions are defined.

C flow no. of

net mass flow coefficient exponent exceptions

flow to building

( l/s

,

cfm) »

For each outside connection the following data are required.

outside pressure C X

profile number flow coefficient flow

exponent

N

rO no. of connections to

the outside

For each exception the following data are required.

exception type

(1, 2 or 3)

level of

shaft

The next card depends upon the exception type. For exception type =1, an exception is defined to the net flow into the level of the shaft.

net mass^flow

( H/s

,

cfm)

For exception type = 2 , an exception is defined to the outside connection

for this shaft.

outside pressure C X

profile number flow coe¥ficient flow

exponent

A8
For exception type = 3 , an exception is defined to the flow coefficient

between the shaft and the building.

C

flow coefficient

7 • FIND F input .

For a FIND F run type the following data is required.

shaft level

number of shaft

desired differential

pressure across

shaft

convergence

limit

8. INPUT2 data

This routine is used to modify the existing data for a rerun.

N

no. of

modifications

If N

m

= 900, program

execution is stopped.

For each modification the following data are entered.

modification

type

IS

Integer

B

real

number

For modification type = 1, the net flow (F^ = B) at each level of shaft, IS,

is changed.

For modification type = 2, the flow coefficient (C = B) between each level

of shaft IS and the building is modified.

For modification type = 3, the flow coefficient (C Q

= B) to the outside is

changed for each outside connection from shaft IS. However, if IS = 0, then

the flow coefficient for each outside connection from all building compartments is modified.

A9
For modification type = 4, the flow coefficient (C

q = B) between two spaces

(IS and JS) is changed. In addition to the modification type, IS and B the

second space number is entered on the next line.

JS

second space

number

A10
APPENDIX B. INPUT EXAMPLES

APPENDIX B. INPUT EXAMPLES

1. Example 1

A ten story building with a pressurized stairwell and no vertical leakage

within the building is heated to 70°F when the outside temperature is -20°F.

The stairwell temperature is 60°F at the tenth floor and 50°F at the bottom

floor. The stairwell is pressurized by a net 400 cfm* per floor. The wind

is 30 mph at a height of 30 ft and the wind exponent is 0.14. This building

has connections to the outside in two directions. The flow coefficients are

the same vertically and are listed in table Bl and the flow exponents are

expected to be 0.5.

Connection Cq/ flow coefficients*

location cfm per /"^O

Between stairwell & building

900

Between building & outside into the wind

800

Between building & outside away from

the wind

800

Table Bl. Flow coefficients for example 1.

1.1 General Input

In the general flow network input the building compartments can be located on

any floor and there can be any number of compartments per floor. In this

case there is only one compartment per floor so the floor levels are chosen

to agree with the compartment numbers.

Therefore, space 1 through 10 refer to floors 1 through 10 and spaces 11

through 20 refer to the spaces in the stairshaft. Connections between two

spaces only need to be defined for one of the spaces. In this example the

connections between the building and the shaft are all defined for the shaft

spaces

* At standard conditions

Bl
RUN

Initial EXAMPLE OF GENERAL FLOW NETWORK INPUT

data -20 2 0

building 10 1

height* 5 10

te.mpzAatu.sie 2 pAofille* 1 1 70

2 1 60 1 50

out*lde 2 1

pAe**uAe 30 30 .15

pAofille* 0. 7

- 0.7

General

building

data 10 0

Ut 0 2 0 1 1

{^loOA 1 800 .5

2 800 .5

1 nd 1 0 2 0 1 2

1\looA 1 800 .5

2 800 .5

3 Ad 1 0 2 0 1 3 filoOA 1 800 .5

2 800 .5

4th 0 2 0 1 4 filOOA 1 800 .5

2 800 .5

5th [ 0 2 0 1 5 filooA 1 800 .5

l 2 800 .5

6th 1 0 2 0 1 6

t\l00A , 1 800 .5

I 2 800 .5

7th 0 2 0 1 7

{\l00A 1 800 .5

2 800 .5

Sth 0 2 0 1 8

t\l00A 1 800 .5

2 800 .5

9th 0 2 0 1 9

faloOA i 1 800 .5

2 800 .5

1 Oth 0 2 0 1 10

ilOQA 1 800 .5

8 800 .5

B2
1

I

Gumnal

ihCL^t

data

STAIRWELL

| 72000 1 10 2 u t 1 0 400

llooh. 1 900 .5

2 nd 1 0 400

hloon 1 900 .5

3 A.d 1 0 400

hloon \ 1 900 .5

4th. 1 0 400

6 Iooa 1 900 .5

5th 1 filoot

1 x 0 400

1 900 .5

6th 1 0 400

hloon 1 900 .5

7th 1 0 400

iloofi

l 1 900 .5

8th 1 0 400

faloon. 1 900 .5

9th 1 0 400

filooh 1 900 .5

10 th 1 0 400

faloon 1 900 .5

znd

data

j 900

1.2 Abbreviated Input

This problem can also be solved by use of the abbreviated flow network input.

This reduces the quantity of data required and to that extent reduces the

chance of input error. The abbreviated data is as follows:

RUN

Initial EXAMPLE OF ABBREVIATED FLOW NETWORK DATA

data -20 2 0

building ( 10 1

height* { 5 10

te.mp<inatune.

pn.o kll<i.t>

21 1 70

2 1 60 1 50

B3
outAtdz

pHZt>bUHZ

psiofitZeA

ab bazvtatzd

bui.ldi.ng

data

abba 2.ui.atzd

Akafit

data

2 1

30 30

0.7 -0.7

10 1

I 2 1

i 1 800

1 2 800

.15

0 .5

.5

.5

!

1

STAIRWELL

72000 1 10

0 400 900 .5

end { 900

data

0

2

0

1.3 Example 1 Output

The output for this example is the same for both these forms of input. This

example output is given in appendix C.

2 . EXAMPLE 2

This is a 30 story building which is 70°F inside. The outside the air

temperature is -5°F and there is no wind. This building has a stairwell and

an elevator. The flow coefficients which are generally the same vertically

are listed in table B2 and the flow exponents are taken to be 0.5.

Connection location C flow coefficient*

ocfm per /"H ^0

Between stairwell & 900

building

Between building & 1600

outside

Between elevator & 1600

building

Table B2 . Flow coefficients for example 2.

* At standard conditions

B4
In this example a total of four separate cases are analyzed. These cases

are:

1. First floor outside door open.

2. First floor outside door closed.

3. First floor door between the elevator and the building open.

4. Stairwell pressurized at a rate of 900 cfm* per floor and all doors

closed.

2.1 Example 2 - Input

The last three cases are handled by the use of the INPUT2 routine. The data

are as follows:

RUN

tntttal | EXAMPLE INPUT - - 30 STORY BUILDING

data -5 2 0

butldtng 1 30 1

kztghtA 1 5 10

tnmp 2.f1 atan.11

j pn-ofitlH 1 1 1 70

outAtde. 1 1

pn.HAAun.2. 0 1 1

data 1

abbn.Hvtatnd | 30 1

butldtng 1 1 0 0 1

data 1 1600 .5

2 1

[ 1 25000 .5

2 Abbn. 2.vtatnd Ahafit 1 STAIRWELL

Shafit Vata 72000 1 30 1

o 0 900 .5 0

Ahafit 2 ELEVATOR

80000 1

i

—

CO o

0 0 1600 .5 0

Cai> 2. 2 : mo dtfitcatto n 1

type. 3 3 0 1600

* At standard conditions

B5
Col6<l 3:

modtfitcatton 1

typo, 4 4 1 25000

31

Ca-ae 4:

modtfitcatto n 2 typ<U> 1 and 1 2 1 900

1 1 400

end ( 900

data

2 . 2 Example 2 Output

The output for example 2 case 1 (the data above not including modifications

for Case 2, 3 and 4) is given in appendix C.

B6
APPENDIX C. EXAMPLE OUTPUT

APPENDIX

C.

EXAMPLE

OUTPUT I X £ eg in in eg O'

O

N eg <0 O' ffi eg ffi eg eg

N ffi eg O'

© in eg o O' o eg

O O' ffi eg O' o - eg ffi in o ffi ® O' in in O' in o O' © s N O' s N N ffi N ffi o ffi m in ffi ©

U. 10 O'

N

®

N o eg to s in ffi ffiN O' O' O'

O o «* ffi eg o ffi CM©

m

H in ffi eg in ffi eg

n eg ffi eg <o ffi eg ffi o ffiNO ffiN * ffi

N o fn

K Z * ID O Z J o

u.

ax

UJ

in tn in in in in in in in in in in in in in in in in in in in in in in in in in in in in O o O • • • o o o o o o

O' 03 03

O

O

O • • • O © O O o o

O' CO

® o o o • • • o o o o © ©

O' CO® o o © • • • o o o 0 O O01 CO 00

© o O • • • o o o o o o

O'

®

a o o © • • • o o o o o oO' « ® o © o • • • o o o o © o

O'

®

® o o o • • • o o o o o oO'

®

® o o © • • • o o o o o o

O'

®® o o o • • • o o o o o o

O' ffi

® - - eg eg - eg ffi •< <g «• eg in - eg <o - eg n -• N « - eg » - eg o — eg

eg o

UJzzou

UJ UJ a o h N H z in in

-> iK 0 3 3 a o o

i- »• « z in tn ** b- b- 0 3 3 a o o

UJ UJ a o h M H z in v) ~ i- 0 3 3 a o o

UJ UJ o o H * w Z U) W

i- tH 0 3 3 a o o

UJ UJ a o

I—

w « Z V) V) -> b~ b- 0 3 3 a o o

UJ UJ o o z in to 0 3 3 a o o

UJ UJ o ©

|*MH z in w "hh 0 3 3 a o o

UJ UJ o o z w V)

IH 0 3 3 a o o

UJ UJ o o z in tn

*-*

HK 0 3 3 a o o

UJ UJ o o h M H z in in «-(-- 0 3 3 a o o o

UJ

X

X o H J

IL u.

iszH UJ a a J 3 UJ

*-«

h -t 3 <

ID

a IL UJ o

>-

a

a

a z a o UJ h t- ©O

M*

ao 1 o 1 j

u.

_J <

UJ b- a < o O »- z z

*<4 UJ

X z o UJ J- J _l QC IL a < z a a < z M X o < 111 u

Cl
CONNECTED

TO

COFFICIENT

EXPONENT

SCFM

a

a

ujz N o \0 m O' ro o 0> N .

• • • • • • • . . • 1 n in o N s <o <0

O' N N CD CO O' O —

•

«\j n

ro m ro n ro m * 4 <r

I I I I I I I I I I

lAtninintninwinwin a

o

uj x X o « _l

u. u.

a

a

o

-i

IL

oooooooooo . ........a oooooooooo

9 9 l> 0> O1 O' O'OO'

«NI<l<HliMH>0

zzzzzzzzzz

oooooooooo aaaaaaaaaa

oooooooooo .......... oooooooooo

•-CVItA^IA'ONOO'O

UJa

UJI

ClIAo o

CMz

zM

o W

z

IX (A o H

IL M < 3 z

a> V) > IA a IA o o

*- a w o

< . o a

•4 UJ o

a t- 3 3

i in O

1 ia

UJ -1 a _J < a UJs a < a. a < « O M < o

1- *- z

z <A «

UJ UJ _j a

a UJ a

z u. o

< u. o

X H j

UJ o u.

NIDO«l04O0'N0> OCOOONSNNOxO oooooooo-* —

SMHDONKIOOS W-OCMMDSSSO

<\J04<\|- — —

O<M0S'0IDtfnNM w >4 O H O

C2

2*674

2

2*545

3

2.397

4

2*244

5

2.090

6

1.933

7

1.776

8

1.618

1.460

10

1.301

STAIRWELL

2.862

12

2.715

13

2.567

14

2.419

IS

2.271

16

2.123

17

1.975

18

1.826

1.677

20

1.527
AIR

FLOW

MODEL

EXAMPLE

2

DATA

—

30

STORY

BUILDING

! L

h- h h h t- K

UJ UJ UJ Ui UJ UJ UJ UJ UJ z z z Z z z z z z

* y. o n O' <M in (VI n O' <0 © © O' © PH P* CM © CVJ o © © (VJ © o © © © •h p* © pH © © © © ©

o fJL J u © in o o N © o O' © O' © © <\J pH in © N pH © CVJ © © N N o © N

u. in o pH w N (0 © O' * © pH © ft © w <VJ O o © O' N © N N © PH N

in O' ft w in IN (VI N OJ ft N (VI <t 'O (VJ © (vi ft © pH © © pH © © ©

! i pH 1 1 1 1 1 T i 1 1 1 1 1 i i i 1 1 1

h

Z

* UJ O Z

_i o ininto © © © © © © © © ©

u. a ••• • • • •••

x

Ui

© © © © © © © to in in in in in in in in in in

o

iux o

44 _i

u. u.

U o o o o © © o o © o o o o © o o o o O o o o o o O o o o o o 44 • • • • • • • a » • • • • • • • • • • • • • • • • • • • • •

u. o o o o © o o O © o o o o o o o o o o o o o o o o o o o o o

u. o o o o o o o O o o o o © o © o © o © o o o o o o o o o o o

o O' © o O' © © © © © © <0 © © © © © © © © © © © © © © © © © © ©

u pH ©

(VJ

44 pH pH pH 44 44 pH 44 44 pH 44 44 44 44 pH pH

01 pH 44 (VJ (VJ 44 © © pH * 44 © © pH © © 44 N N © © *4 © © 44 © o pH z

© © © « © © © © © © © © © © © © © © p N

U

44H UJ UJ UJ UJ UJ UJ UJ UJ UJ UJ U a o o o o o o a o a

UJ i- t- w 1- J- 44 H J- 44 H 44 H K 44 N N 44 K N 44 K »- 44 1- 1- 44 H H 44 z z z © z Z © Z z in Z z © Z z © Z Z w Z Z © Z z © z z © Z Z ©

2- 44 PH J- 44 44 »- t-4 44 H 44 44 H 44 44 H 44 44 1- 44 44 1- 44 44 H 44 44 »- 44 44 H

o o o 3 o o 3 o o 3 a o 3 o o 3 o 3 o o 3 o o 3 o o 3 o 3

u a a O a a O a a O a a o a a O a a. o a a O a a O a a. a a a O

o O © o o o o o o o

X • • • • • • • • •

C3
POINT

42

900.0

.5

-122

POINT

72

1600.0

.5

-243

OUTSIDE

I

1600.0

.5

366

Ul2

N <0 <t O

IDZ

CVJ T- CVJ CVJ

UJ Ui Z Z

cm*) in n co © in —

IDz

in w

K

UJz

— © ro o> o to — — o cvj in co o cvj © — © — o in.ni n —

n n in

O' O' O'

I — CVJ

I

in in —

O PI o

, -CVJ © ^ 14

CVJ «(• hI o n n

N <t —

- CVJ

I

o co co o O' O' — — CVJ

I

cvj o cvj

CVJ 10 — cvj ro

I

o in to ^ N — — CVJ <t

I

>0 N I*) in o o

— m <*

l

m in in in in in in in in in in in in in in in in in in in in in in in

o o o o o o © o o o © © o o o o © o © o o o © o • • • • 0 0 0 • • • 0 0 • 0 0 0 0 0 0 • 0 • 0 0 o o o o © © o © © © © © o o © o o o © o o © o o

o o o © © © © o © © o o o o o o © o © o o o o ©

O' © © O' <o © O' © © © © © O' © © © © © © © © © © ©

n m -

N

* * —

N

in in —

4 N

© © —

* N

N N —

* N

© 00 —

* N

© © —

* N

co- rn ©

UJo - H -

Z Z W

-> - J- 0 0 3

a a o

UJo

i- i- «

z z tn H M h

0 0 3

a a o

UJo

hi- —

z z tn — — h

0 0 3

a a o

z z ©

0 0 3

a a o

UJ©

h" H —

Z Z © — — J— 0 0 3 a a. o

UJo

H H -

Z Z «

— — H 0 0 3

a a o

UJo

j- »- —

z z ©

- - I- 0 0 3

a o. o

UJa h h -

z z «

— — H

0 0 3 a a o

C4

POINT

51

900.0

.5

172

POINT

81

1600.0

.S

335

OUTSIDE

1

1600.0

.5

-507
IdZ

bUiz

Idz

idZ

-

idz

idZ

ft id idz 4 n n * © © ro kocom eg eg -i eg © m o< eg o> « ro -< s s 4 »* eg © 4 © n ui « n

** N O' © 4 W CO lO <t g\ CO © *-• O

© © — © © eg 4

i i

<X 4 O' © © o 4

eg eg © © 4 ©

© eg 4 © <VJ 4 ©

i i i

© © ©

4 © -

eg 4 K

I

K © 4

© © 4

N < K

© 4 eg ©Of*

eg © t*

i

(0 — 0 SN» eg © n

l

© © © © © © © © © © © © © © © © © © © © © © © © © ©

© o o o o o o o o o o o o o o o o o o o o o o o o o

X z o u. o © <x © o r> ©

• 1* » • • • • • • • • • • • • • • • • • • • • • • • _1 u • • • • • • • • o o o e o o o o o o © © © o o o o © o o o o o o o o 11 </> © o © © *4 s eg o o o o o o o o © o © o o o o © o o o © © o o o o o o N in 4 n eg O o © © o © © o © © o © © o © © ex © © o © © o © © o © © © eg eg eg eg eg eg ©4

z

x id eg •4 © © 4 * in © «4 © © *4 N © © ©4 o cx «4 o © •* o z © © © © © © © © © © © © © © © © o _l o

n a

X

© in © in © © © ©

id id id id id id id id

id o o o o o ea o o o

h* 4 i- 4 h i- i- 4 »- b- fr4 b- 4 1- i- 4 *- b

M »- 1- *4 z in z z © z z w z z V) z Z w Z Z © z z © z Z w z z ©

4 H *•© 4 i- M «- 1- 4 4 b- 4 *«4 1- *4 t- *4 *4 l- 4 *4 i- *4 4 t- O 3 o o 3 O O 3 o O 3 o O 3 o O 3 O O 3 o o 3 o o 3

a O a CL O a a o a CL O a a O a a O a a. O a a O a a O

b

•

Z

id o o © o o o © ©

x — • • • • • • • • o o o o o o o o o u o o o o o o o o

• • • • • • • • J « o o o o o © c o

11 11 o o o o o o o o

u.o

u

©

eg 4

eg ©

eg © seg eg ©

eg o o

eg ©

o

»-Q

idH

U

tdz

zou

M«in4lA<MI

h- t~ b- b- z z z z 4 44 *4 1-4 o o o o

a a cl a

b- 1- H t- z z z z 4 »-4 M M

a o o

a. a a a

o© x x o oooooooo «JU- u.

© 4 ©

eg eg eg © n ©eg eg eg o> oeg ©

-eg©<t©©N©

C5
H

IUZ

OMNMflJiTONO

'0'3CUO''OCUr-OCU4 - - - I I — - -

I I I I I I I I I I 1 I

* s H CU co h- (\J ro • O U- N 4 <0 O' CM n o *4 •4 in <0 N ro 'O in >0 o - co m 4 in o CU 10 4 CO

CU in >0 N CO 00 il in in o o O' co in 4* in s vC (\J ro N in 4 ro CO o in N in *4 in CO a O' 05 o

cvj ro 4 in >0 N *4 *4 O' >0 oj o N ^4 00 O' ro <* a> 4 N o ro *0 CO o <M 4 'O <0

<\J (VJ <VJ CVJ CM CM cu O' in 4 4 4 m n r> CVJ cvj *4 *4 1 *4 CVJ CVJ ro ro ro lO <* 4 4 4 4

I I I ! I I I I I I ! I

ininmintninininininintninininininininininin ininininintninininininininininininininininininininininin

z

o o o o o o o o o o o © o o o o © © © o © © UJ

3: •4 o o o o o o o o o o o o o © © o o o o o o o o u

© o o o o o o o © o © o o o © © © © o o o o _J 4

O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' U. Uu.oo

oooooooooooooooooooooooooooo

oooooooooooooooooooooooooooo OOOOOOOOOOOOOOOOOOOOOOOOOOOO

'O'O'n'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'On'O'O'O'O'O'O'O

o

o>o-(\)rri<tir)iosooo'o — NP)<n/i'Oso»o

— — x-*(M(\jf\jrg(\jcgc\ir\IW(viro Q

UJ H

U

l— »— t— l— »— V— >— *— t— — t— I— »— h-»— I— I— »— I— »— I— I— UJ zzzzzzzzzzzzzzzzzzzzzz z

- z ooooouooooooooooaooooo o aaaaaaaaaa.aaao.aaQ.aaaaa u

— Pjr04in'0N000'O — CUI04in'0NC00'O - cu ro 4 in <0 s 4H4444HHH4(\| CVJ CVJ CVJ M cu cu cu cu

h h- K »- t- 1- - h- >- t- I- t- K H H h* 1- »- h- H- 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 o o o o o o u o u O U a o O o o o O o u O o o o o o o

a a a. a a a a a a Q. 0. 0- a a a a a a a a a a a. 0. Q. a a a

oooooooooooooooooooooo

o

UJ 3

X o « J

u. u. oooooooooooooooooooooooooooo

ao NOOI()tN®ffO-Nn — cvjcvicvjcg 4 in \0 f- ® o

csj c\j cm cm <\j c\j ro

<

>

UJ

-IUJ

a.o

o

_i

u. HNI<)4in«M00'O"NI,)ffl0OS10J'O-

~ — CM CM <U ® 4

CM CM CM

in io r- co

CM CM CM CM

C6
CM

— 0 in

O' s • •

ro cm * ® 0 a -

•0- CM

>->

n in

in in • • » •* O <\J 0 in 1 I

in ro 0 in

CMO • • 4 ro

>1 04 o in N 4 — O 4 «M • • 4 ro

in ro o o o o o o w

a.UJX

fl.

® <0

vo in

in ro • • 4 io f 4 <\j

C/) o • 0 to O' o <\j O' in CM® in

N 4 o Nm

H 4 ® CM MD O' ro

N ^4 in O'n

® o o CM a. IT) ro o ® 0 ro 44 O' 0 4 CM O' N in CM o CM 4 O' •4 m CO H m in ® o 10 r X a ro ro ro CM CM CVJ CM 44 *4 44 44 o O o o o o o o O ^4 ^4 CM CM CM CM ro ro < z 3 1 1 1 i 1 1 1 1 1 i i | | 1 o 4 O CM o O 44 o W a. i- 0

(M io z • • Q 4 w a ro Z o H o < J u. K a v0 O' ro 0 O' CM in ffl r4 4 4 o ro in ® 44 N m i*) *-4 O' N in W ^4 ® 0 CM ®

KH

M < < o o O' O' ro

N

N >0 in in 4 ro ro CM •4 o o o •4 CM ro ro 4 in

n

N

N

® O'

© CM 4 ro 44 z z UL 3 X > o r4 o o O O O o o o o o o O o o o o O O o o o o O O o © o ro 44

•-J 44

m in UJ o o o _! z 1 t i ! 1 1 1 i i 1 1 1 1 | a a UJ >- <n UJ 44 o £ V) O o 4% o H a O CM O' to o CM N 0 < • IT) CM O' in CM O' >0 ro o N <* CM O N in ro O' 44 ro 'D CD *4 <* N *4 •4 ® 44 in O' X a CD 0 o cc CM O' N in CM o CO 0 ro •<4 O' N <* CM o H <0 CO o ro in N o CM <t N O' 44 • • O O < to UJ Q ro ro CM CM CM CM CM •4 44 *4 o o O O o o o o o ^4 H CM CM CM CM CM ro z 0 • • a 1- 4 £ 3 3 1 1 1 1 1 1 i 1 1 1 i j | 1 w

ii. 1 10

O 4 CM o o l in

UJ -J 0 z a -I UJ 44 < a UJ a

»- * £ 0 O' <VI •0 O' CM >

O' ro

N o 4 CD ro

N 44 in CM CD in CM O' M3

m o 0D in roO 0 3

tn < _i

a

O O' CO C0N vO 0 in to to CM •4 o o o 44 *4 CM ro

m in <0 •a

N CD O' *4 if) * o < 44 O o o o o O o o c o o o o o o o o o o o o o O o o o o o o o ro if) CM o

O' o o M < O UJ *4

CM ro -1 cm - H z 1 1 1 i 1 i 1 i 1 1 1 i 1 1 a

a CM

®

u. Z in 44 a. • • UJ UJ 0 n J a a UJ a T u. o 44 O' < u. o o O' 00 N 0 in <r to CM *4 o O' ® N >0 in ro CM o O' <D K v0 tr> n CM •4 X 4 J ro CM cm CM CM CM CM CM CM CM CM UJ o IL

C7
2.600

18

2.450

19

2.300

20

2.149

21

1.999

22

1.849

23

1.699

24

1.549

1.399

26

1.250

27

1.100

28

.951

29

.801

30

.652

j © CO © N O 4

.. in o O ©

O' N © o s ©

a • a a • a © CM "* © CM r*

© <0 4 ffl O 4

© © © N ffl

*4 ro co CM S

vO o 4 in O ©

o a* N o O N

• • • • a a

CM CM ft N in ro K © ©

ro <• in O N oo

in © CM N in © O

o 4 O' o- a 4 o

CM o 00 N •4 O o

• • a • • a • © art 4 ©

'O 4 CM O © 4 CM

fO 4 © vD © S CO

o* r) o o

4 O' 14 O' IO - O ffl . . . • 4 © cm

in n - o

© 4 © ©

GO o* m —

«

© oo 4 o © -. o o>

• . . • 4 © CM

in ro — o* O MS O

© r>« o ©

O' ro © ro 4 to — O ....

4 © CM «.

4 CM © 05

ro 4 in in

w ro <o ro © ro co 4

4 ro — o

• . • • 4 t0 CVI «

4 CM O ffl © N ffl ffl

N N 4 in © ffl CM N

S 4 10 -

» a a a 4 10 CM ^

10 - 9 N

ro 4 4 in

© © O' <0 CM N CM ffl © 4 ro «. * » . • 4 ro CM «« © «a 0> N

© N N ©

— s © ©

ffl CM © -

N © 4 © ... a 4 to CM «*

O O CM O' N CM N CM

h- © 4 to • . a • 4 to CM ^

CM O ffl © CM O ffl © f044in © N N ffl J

-I a

ID o

3 t- a < ~ ffl CM © -i > 4 © © — <CMN-«©IU-<©^N HO'S©4_JO't'©4

I/I ... a IS a a . • 4 © CM — 4 © CM «

t" ©

— CM

-> O' N © — O' N © ©©4© © © N ffl

C8
APPENDIX D PROGRAM LISTING

APPENDIX D. PROGRAM LISTING

* MAIN PROGRAM

S>NBS*PL IBS.SHOW

C

A. MAI

C

C

c

c

c

c

cc

cc

c

ccc

c

cc

cc

cc

c

c

ccc

cc

cc

c

cc

c

cc

cc

c

cc

c

c

c

cc

c

c

COMPUTER PROGRAM FOR AIR FLOW ANALYSIS IN BUILDINGS

SPECIFICALLY FOR ANALYSIS OF PRESSURIZED STAIRWELLS

AND PRESSURIZED ELEVATORS

PROGRAM VARIABLES

C FLOW COEFFICIENT BETWEEN BUILDING POINTS

CO FLOW COEFFICIENT TO OUTSIDE

CS FLOW COEFFICIENT OF SHAFT

E LIMIT WITHIN WHICH CONVERGENCE IS ACCEPTABLE

F NET FLOW INTO POINT I FC FLOW BETWEEN INTERNAL POINTS

FF FIXED FLOW INTO POINT I FO FLOW TO OUTSIDE

FSS NET FLOW INTO SHAFT IS

H HEIGHT FROM GROUND TO MIDPOINT OF FLOOR

I BUG OUTPUT VARIABLE

ICONV INTEGER USED IN SUBROUTINES BLDGP AND SHAFTP

IF ICONV = 0 THEN THE PRESSURES WERE UNCHANGED

I FLOOR FLOOR LEVEL WHERE POINT IS LOCATED

IT POINTER TO TEMP PROFILE FOR POINT I ITS POINTER TO TEMPERATURE PROFILE OF SHAFT

JC POINT NO* CONNECTED TO POINT I JOC DIRECTION OF OUTSIDE CONNECTION

N NO* OF BUILDING COMPARTMENTS

NC NO. OF INTERNAL POINTS CONNECTED TO POINT I NCO NO. OF OUTSIDE CONNECTIONS

NFS1 BOTTOM FLOOR OF SHAFT

NFS2 TOP FLOOR OF SHAFT

NH NO. OF FLOORS

NPO NO. OF OUTSIDE PRESSURE PROFILES

NS NO. OF SHAFTS

NS1 I VALUE FOR START OF SHAFT

NS2 I VALUE FOR END OF SHAFT

NT TOTAL NO. OF POINTS (BLDG AND SHAFT)

NTP NO. OF TEMPERATURE PROFILES

P PRESSURE AT POINT I PFO OUTSIDE PRESSURE PROFILES

PO OUTSIDE PRESSURE

PS PRESSURE PROFILE OF SHAFT - WORKSPACE

PZ PRESSURE DUE TO ELEVATION DIFFERENCE

T TEMPERATURE PROFILE ARRAY

TITLE PROJECT TITLE

TITSH SHAFT TITLE

X FLOW EXPONENT TO INTERIOF POINTS

XO FLOW EXPONENT TO OUTSIDE

c PROGRA M PARAMETERS

c MB MAX NO. OF BUILDING COMPARTMENTS

c MM MAX NO. OF POINTS

c MS MAX NO. OF SHAFTS

c MC MAX NO. OF CONNECTIONS FOR ANY POI NT

c M PO MAX NO • OF OUTSIDE PRESSURE PROFIL ES

c M TP MAX NO. OF TEMPERATURE PROFILES

Dl
r>nnnn

ro»-

nonr>nr>nnnn

nnnnonnnnn

n

n

MAIN PROGRAM

C MFL MAX NO* OF FLOORS

PARAMETER MM =2 20* MS= 2 • MC=4 • MP0=2 . MTP=2 . MFL=1 05 , MB = 1 05

COMMON NT« P( MM) « C(MM,MC) ,NC( MM) , JC{ MM.MC) . ITS(MS )

1 FC ( MM« MC )cPZ!MM»MC)»POCMM, MPO ) ,CO( MM*MPO) , F ( MM ) , PFO! MF L, MPO) «, 2 FF(MM) «FO( MM, MPO) ,CS( MS) ,PS< MFL ) .NS1 ( MS) . NS2(MS)

,

3 FSS(MS) e N* NS » NPO » I CON V • E . I BUG , X ( MM , MC ) ,XO(MM,MPO),TITSH(MS.3)r

4 NH,H< MFL) . I FLOOR! MM) , T ( MTP • MFL ) . NF SI ( MS) * NFS2 (MS) . IT (MB) #NTP

5 9 NCO(MM) e JOC ( MM* MPO) e TOUT

DOUBLE PRECISION P»PO,PS

COMMON /RUN/I RUN

DATA MRUN/6HRUN /

DATA MR UN 1/6 HR UNI /

DATA NTEST /6HTEST /

DATA NFINDF/6HFIND F/

READ ITEST TO DETERMINE RUN TYPE

I TEST = RUN FOR NORMAL RUN

ITEST = RUN 1 FOR A NORMAL RUN OF 1ST DATA AND

A TEST OF ALL OTHER DATA

ITEST =TEST FOR A TEST RUN OF ALL DATA

ITEST = FIND F TO FIND FIXED FLOW » FF « FOR A SET PRESSURE

ACROSS A FLOOR OF A SHAFT

RE AD ( 5o 804 ) ITEST

NITER=500

I R U N = 1 CALL INPUT TO READ DATA

CALL INPUT

SET UP FOR FIND F RUN TYPE

IS = SHAFT NO*

IF = FLOOR NO.

DPS = SET PRESSURE ACROSS SHAFT IS AT FLOOR IF

EDP = ALLOWABLE CONVERGENCE LIMIT

IF( ITEST *NE. NF I NDF ) GO TO 10

READ! 5* 700) IS. IF. DPS. EDP

I DP = NS 1 ( I S ) + I F—

1

JDR = JC( I DP . 1

)

I 0=0

N 1 = NS 1 ( IS)

N 2 = NS2 ( IS)

WRITE!

6

e 807) IS. IF, DPS, EDP

0 E = 0 * 2

0 I CS=1

IF! ITEST • EQ* NTE ST ) GO TO 41

IF! I R U N • GT • 1 .AND* ITEST • EQ * MRUN11G0 TO 41

TEMPERATUPE CORRECTION

CALL CORR

CALL I NT T TO INITIALIZE PRESSURE ARRAY , P

D2
no

non

o

n

o

n

•£•

r>oo

r>

non

*

no

o

w

oon

non

* MAIN PROGRAM

C

22 CALL INIT

C

C

C DO LOOP TO 30 IS ITERATIVE SOLUTION TO PRESSURE ARRAY

C

24 DC 30 ITER=1, NITER

C

C CALL BLDGP TO SOLVE FOR BUILDING PRESSURES

C

CALL BLOGP

I CB= ICGNV

I F ( I CB • EQ • 0 .AND. ICS .EQ. 0 » GO TO 40

CALL SHAFTP TO SOLVE FOR SHAFT PRESSURES

CALL SHAFTP

I CS = I CONV

I F ( ICB * EQ » 0 .AND. ICS .EQ. 0)G0 TO 40

CALL PZAD TO CALCULATE PZ TERMS

CALL PZAD

CONTINUE

IF ROUTINE FAILS TO CONVERGE IN NITER

ITERATIONS PRINT ERROR MESSAGE

W R I TE (6c BOO

CONTINUE

WRITE(6« 80 1 ) I TER

GC TO 50 FOR FIND F RUN TYPE

IF ( ITEST .EQ. NFINDFSGO TO 50

CALL OUT TO OUTPUT SOLUTION

2 CALL OUT

I F ( I TE ST .EQ. NFINDF) WR ITE( 6*806) I S »FM

1 IF< ITEST .EQ® MRUN)WRITE( 6. 305

CALL INPUT2 TO CHECK IF ANOTHER RUN IS REQUIRED

WITH MODIFIED DATA

CALL INPUT2

IF( ITEST • NE« NF INDF ) GO TO 20

ZERO FF FOR SHAFT IS FOR 1ST ITERATION OF A NEW FIND F RUN TYPE

00 44 I=N1,N2

F F ( I ) = 0 * 4 CONTINUE

1 Q = 0

GO TO 20

D3
*

cc50

Cc

c52

51

53

54

56

C

C

c58

59

60

C

c

cccc700

800

801

804

805

806

807

808

809

900

MAIN PROGRAM

SELECTION OF NEW FIXED FLOW • FM . FOR FIND F RUN TYPE

DPM =P ( I DP )-P( JDP) — DPS

I Q- IQ+1

I F ( 1 Q .GT. 30 ) GO TO 42

I F ( I BUG • ECU 0 ) GO TO 52

WRITE!6»901 )D PM » DP 1 ,DP2 * F 1 ,F2 • I DP , J DP

CHECK CONVERGENCE OF DIFFERENTIAL PRESSURE , DPM

I F ( AB S ( DPM ) • LT • EDPJGO TO 42

I F ( I Q .NE. 1 ) GO TO 51

DPI =DPM

F1 = 0.

F M = 500

GO TO 58

IF( IQ .NE. 2 ) GO TO 53

I F ( DPM .LT. 0 • > GO TO 42

GO TO 54

IF(DPM*0P1 .GT. 0 • ) GO TO 56

DP2=DPM

F2 = FM

F M= SORT < (F1*F1 ) -DP 1 *( F 2 *F 2-F 1 *F 1 )/( DP2-DP1 ) ) GO TO 58

DPI =DPM

FI =FM

FM=SQRT! ! F1*F1 )-DP 1 *! F2*F 2-F 1 *F 1 )/( DP2-DP1 ) ) ASSIGN FM TO FF FOR SHAFT IS

DO 59 I=N I . N2

FF( I )=FM

I F ( I BUG .EQ. 0 ) GO TO 60

WRITE!6.900 )FM

I CS=1

IF( IQ .GT. 3) GO TO 24

GO TO 22

END OF SELECTION OF MEW FF FOR FIND F RUN TYPE

FORMAT STATEMENTS

FORMAT! ) FORMAT (/////5X, 35! 1H1 >//5X,

+ 3 5H F A I LURE OF MAIN PROGRAM TO CONVERGE //5X, 35 ! 1 HI )//

FORMAT! 10X.I5.5X.il HI TEP AT I O NS ) FORMAT ! A6

)

FORMAT! 1 H 1

)

FORMAT! 10X.20HFIXED FLOW IN SHAFT .I2.3H IS.F3.2/1H1)

FOFMAT! 10! /) .10X.28HRUN TO DETERMINE FIXED FLOW » 1 18HPEP FLOOP IN SHAFT .I2/10X.

2 41HSO THAT THE PRESSURE DIFFERENCE AT FLOOR .13,

3 3H IS.F8.2.7H PA SC AL/ 1 OX , 2 8H W I TH A CONVERGENCE LIMIT OF # 4 F8 • 5 o 7H PASCAL# 10! /) ) FORMAT! 10X# 10HNET FLOWS / /& ! 6X* 1HI .4X.3HNF )/)

F0PMAT!/10X»1 5HNET FLOW SHAFT .I3.3H IS.F6.1)

FORMAT! 10X ,4HFM =,F10.2)

D4
MAIN PROGRAM

901 F OR MAT ( // 1 0X» 5HDPM =* E 1 2. 4* 2X* 5H0P 1 =• E 1 2 • 4* 2X » 5HDP2 = *E12.4

+ / 1 OX® AHF 1 = 9 E 12. 4.2X* AHF 2 =. E l 2# 4

2X . 5HI DP =

I 4 • 2X , 5H JOP =.14)

END

®HDG 9 P * SUBROUTINE INPUT .L .

D5
u

uu

uu

u

u

* SUBROUTINE INPUT

EDNB

Cc

c

c

S*PLIB$.SHOW A. INPUT

SUBROUTINE INPUT

THIS ROUTINE READS AND PRINTS DATA

AND INITIALIZES PZ ARRAY

PARAMETER MM=220 , MS=2,MC=4 ,MP0=2 * MTP=2 , M FL=1 05 . M B= 1 0

COMMON /PZZ/ PGZ

COMMON /I 0/T I T LE ( 12 ) , I OUT, IUN I T , NCOMP ( MFL ) , S NCOMP ( MFL

)

COMMON NT o P ( MM ) , C ( M M , MC ) , NC ( MM ) , JC ( MM , MC ) , I T S( MS )

,

1 FC (MM* MC) e PZ ( MM, MC) ,PO(MM,MPO) ,CO( MM, MPO) , F ( MM ) , PFO( MFL

2 FF (MM) « FO( MM, MPO) ,CS( MS) ,PS( MFL ) ,NS1 ( MS) , NS2( MS)

3 FSS(MS) o N,NS ,NPO , I CONV , E , I BUG , X C MM , M C ) ,XQ(MM,MPO),TITSH«

4 NH,H( MFL ) , IFLOOR C MM) ,T( MTP , MFL ) , NF S 1 ( MS) , NFS2( MS). IT (MB]

5 c NCO( MM) c JOCC MM, MPO) . TOUT

DOUBLE PRECISION P,PO,PS

DIMENSION II (MFL) , TT ( MFL ) , P AR < 7 ) , C W ( MPO ) , P H( MFL ) ,NZZ(MM)

DATA PAR/3H MM , 3H MS,3H MC ,3HMP0 , 3HMT P, 3HMFL , 3 H MB/

DATA IBUG/0/

READ AND WRITE PROJECT TITLE

READ (

P 600 ) (T ITLE ( I ) , 1 = 1 , 12)

W R I TE (6.601) ( T I TLE ( I ) ,1=1 ,12)

READ GENERAL DATA

TOUT = OUTSIDE TEMPERATURE

I UN I T = 1 FOR SI UNITS

= 2 FOR ENG UNITS

I OUT = 0 FOR NO SUMMARY OUTPUT

OTHERWISE I OUT IS FILE NO. TO

WHICH SUMMARY OUTPUT IS WRITTEN

READ( 5, 700 )TOUT , I UN IT, I CUT

WR I TE(6 P 4 1 1 )TOUT, I UNI T, I OUT

I F ( IUN I T • GT • 2 .OR. IUNIT ,LT, 1 )G0 TO 105

READ HEIGHTS

NN=0 FOR INPUT OF ALL HEIGHTS

NN=1 FOR CALCULATION OF HEIGHTS

89

97

RE AD ( 5 9 700 ) NH, NN

WR I TE ( 6 9 4 1 2 ) NH , NN

I F ( NH • LE. MFL ) GO TO

I P A P = 6

GO TO 1 10

I F ( NN • EQ. 1 ) GO T 0 97

RE AD( 5 s 700 ) ( H ( I ) , I = 1. WR I T E ( 6 9 41 3 ) ( H( I ) e 1=1

GO TO 9 Q

RE AC ( 5 0 700 ) H ( 1 ) 9 D H

WR I TE ( 6 e 4 1 4 ) H ( 1 ) , D H

DO 98 I = 2. NH

I M = 1-1

i NH)

MPO) ,

MS, 3 )

„

« NTP

D6
* SUBROUTINE INPUT

98 H( I )=H( IM)+DH

C

C R E AO TEMPERATURE PROFILES

C

99 READ! 5* 700 )NTP

WRITE(6*415)NTP

I F ( NTP « LE o MTP )G0 TO 90

I PAR=5

GO TO 110

90 DO 3 IP=1 e NTP

READ!

C 700 ) NNN « ! I I ! J) ,TTI J) .J=l *NNN J WRITE! 6o 416) NNN, ! I I < J ) * TT ( J ) * J = 1 ,NNN)

IF( NNN • GT • 1 ) GO TO 2

DO 1 IFF=1,NH

1 T( IPc IFF )=TT I 1

)

GO TO 3

2 J = 1 JP1 =2

DO 4 I FF= 1 c NH

T ( IP* IFF) =TT( J ) + ! TT! JP1 ) — TT ( J ) ) * ( IFF— I I ( J ) ) / ( 1 1 ( J P 1 )- I 1 ( J ) ) IF( IFF .NE. 1 1 ( JP1 ) ) GO TO 4

I F ( JP 1 .EQ • NNN ) GO TO 4

J=JP1

JP1=J+1

4 CONTINUE

3 CONTINUE

C

c

C READ OUTSIDE PRESSURE PROFILES

C NN=0 FOR INPUT OF ALL PRESSURES

C NN=1 FOR CALCULATION BY POWER LAW

C

READ! 5o 700 ) NPO » NN

WRITE! 6c417)NPO»NN

IF! NPO *LE • MPO ) GO TO 91

IP4R=4

GO TO 110

91 IF!NN .EQ. 1 ) GO TO 81

C

C READ ALL OUTSIDE PRESSURES

C

DO t 1 = 1* NPO

6 R EAD (5* 700 ) PGZ , ! PFO ( J , I ) , J=1 ,NH)

WPITE!6o418)PGZ , ! PFO! J , I ) , J=1 ,NH)

GO TO 85

C

C CALCULATE OUTSIDE PRESSURES

C PATMOS IS ATMOSPHERIC PRESSURE !PA)

C

81 READ!5*700)VW»HW»XW,(CW!I ) » I =1 ,NPO)

WRITE!6*419)VW»HW»XW»!CW! I) *1=1 ,NPO

)

IFIIUNIT • EQ* 1 ) VW=VW*0.277e

IFdUNIT • EQ • 2 )VW=VW*0.4470

PATMOS=101 325.

T OQ=TOUT + 273.

IFIIUNIT .EQ. 2)T0C = !T0UT +4 60. ) / 1 •

PVA=176.4*VW*VW/T00

Z=— 0 .03417/T00

D7
SUBROUTINE INPUT

I F ( I UNI T .EQ. 2)Z=0.3048*Z

CWM=CW< 1

)

I F ( NPO • EQ* 1)60 TO 212

DO 211 1=1, NPO

I F ( CW < I ) .LT, CWM) CWM=CWI I

)

211 CONTINUE

212 PGZ=PATMOS*EXP(HI NH ) *Z ) +C WM*PV A* ( I H C NH ) /H W ) ** ( 2 .*XW ) )-100.

DO 210 1=1 »NH

PHI I ) = P ATMOS* EXP { HI I) *Z ) 210 CONTINUE

DO 82 I = 1 e NPO

DO 82 J=1 • NH

PFOIJe I )=PH(J)+CW( I)*PVA*(IHIJ) /HW ) ** ( 2 • * X W ) ) -PGZ

82 CONTINUE

C

c

C EUILDING DATA INPUT

C NFLS = NO. OF FLOORS IN BUILDING

C I F 1 = LOWER FLOOR IN SERIES OF SIMILAR FLOORS

C IF2 = UPPER FLOOR IN SERIES OF SIMILAR FLOORS

C NOC = NO. OF COMPARTMENTS PER FLOOR

C NZ = NO. OF CONNECTIONS TO COMPARTMENTS ON SAME FLOOR

C NA = NO. OF CONNECTIONS TO COMPARTMENTS ON FLOOR ABOVE

C

85 1=0

SNCOMPI 1

) =0

•

READI5.700) NFLS

WR ITEI 6» 420 )NFLS

IFINFLS .GT. NH ) GO TO 106

7 READ! 5e 700) IF1 , IF2, NOC

WRITE I6#400)IF1 , IF2.N0C

I F ( IF 1 .GT. I F2 ) GO TO 107

NCOMPI IF1 ) =NOC

I FP= I F 1 +1

SNCCMPI IFP)=SNCOMP( I FI )+NOC

DO 10 I Z= 1 c NOC

1 = 1+1

R E ADI 5c 700 ) N Z , N A , N NO , FF I I ) , IT! I

)

WRITE 16*401) NZ. NA , NNO , FF ( I) , IT I I

)

NZZ I I ) = NZ

NN=NZ+NA

I FL OOP ( I ) = IF

1

IFINN .LE. MC)GO TO 111

IPAF=3

GO TO 110

111 I F I NNO .LE. MPO ) GO TO 112

I P A R = 4

GO TO 110

112 IF(ITII) .GT. NTP .OR. ITII) .LT. 1 ) GO TO 102

NCI I )=NN

IFINZ .EQ. 0 ) GO TO 63

C

C INPUT CONNECTIONS TO COMPARTMENTS ON SAME FLOOR

C

READ(5t 700) I JCI I, J ) »C! I »J ) »X( I , J ) » J = l ,NZ)

WRITEI6 9 402)

WRITE(6®403) (JC(I*J),CI I, J) ,XII , J ) , J= 1 , NZ

)

DO 62 J=1 g NZ

D8
* SUBROUTINE INPUT

62 JC( It J)=JC( It J )+SNCOMP( IF 1

63 I F ( NA .EQ. 0)GO TO 8

C

C INPUT CONNECTIONS TO COMPARTMENTS ON FLOOR ABOVE

C

NP=NZ+1

READ! 5* 700 ) ( JC< I* J ) • C ( I , J ) . X ( I , J ) . J =NP . NA ) W R I TE ( 6 o 4 04

KRITE( 6. 403 ) ( JC( I , J ) .C < I , J) .X ( I , J ) , J=NP,NA ) 00 66 J = NP « NA

66 JC( Is J )=JC ( It J ) +NC CMP ( I FI )+SNCOMP( IF1

)

8 NCO < I ) =NNO

I F ( NNO .EQ, 0 ) GO TO 10

C

C INPUT CONNECTION TO OUTSIDE

C

RE AD (5c 700) (JOC(I,JJ),CO(I,JJ) ,X0< I , JJ ) . JJ = 1 ,NNO>

WRI TE (6,405)

WRITE(6,403)( JOC( ItJJ)tCO(I.JJ).XO(ItJJ)eJJ=lcNNO)

DO 9 J J = 1 o NNO

J=J0C(I 9 JJ)

9 P0( Ic JJ )=PF0( IF1 , J

)

10 CONTINUE

I F ( I F 1 ,NE. I F2 ) GO TO 11

I F { I F 1 .EQ. NFLS) GC TO 20

GO TO 19

C

C ASIGN DATA FOR FLOORS SIMILAR TO FLOOR IF1

C

11 I FP = I F 1 +

DO 17 I FF= I FP » I F2

NCOMP(IFF)=NOC

IFF F= I FF+

SNC CMP( IFFP) =SNCOMP( IFF) + NOC

DO 16 I Z= 1 1 NOC

1 = 1 + 1 11=1 Z+SNCOMP ( I F 1 ) I FLOOR ( I ) = I FF

FF ( I ) =FF ( I 1

I T ( I ) = I T( II ) NN=NC ( II

NNO=NCO (II)

NC ( I ) =NN

NCO ( I

) =NNO

I F ( IFF ,NE. NFL S ) GC TO 23

NN=NZZ( II ) NC ( I ) =NN

23 I F ( NN .EQ. 0 ) GO TO 14

DO 12 J = 1 » NN

C(It J)=C( lie J)

X ( I * J ) = X( II , J

)

J C ( I « J ) = JC (II * J) + SNC0MP(IFF)-SNC0MP(IF1

12 CONTINUE

14 I p( NNO .EQ. 0 ) G O TO 16

DO 15 J J= 1 » NNO

JOC ( I s JJ) =JOC( I 1 , J J

J=JCC( It JJ

C0( It J J )=CC( I 1 . JJ)

D9
SUBROUTINE INPUT

XO ( I o J

J

)=XO( I 1 » JJ

)

15 PO< Is JJ)=PFO( IFF, J)

16 CONTINUE

17 CONTINUE

18 I F( IF2 .EQ. NFLSJGC TO 20

19 CONTINUE

GO TO 7

20 N=I

N2 = N

I F ( N • LE. MB ) GO TO 114

I PAR=7

GO TO 110

C

c

c SHAFT DATA INPUT

114 R E AC ( 5s 700 )NS

I F ( NS .LE. MS ) GO TO 113

I PA R= 2

GO TO 110

113 DO 100 IS=1,NS

RE AD (5* 603) (TITSH( IS, I ) .1=1 .3

)

WRI TE(6 e 406 ) ( T ITSH ( I S . I ) • 1 = 1 * 3

)

READ(5*700)CS(IS) » NFS I ( IS) .NFS2< IS ) * ITS( IS ) WRITE(6*4Q7)CS( IS) » NF SI < IS) »NFS2( IS) • ITS< IS)

N 1 =N2+

1

N2=N1+NFS2 (IS )— NFS1 ( IS)

NS1 ( I S ) =N1

NS2 ( IS) =N2

I FF = NF S 1 ( IS)-1

READ (5s 700)NN0,FFF, JCP.CC »XX»NNN

WRI TE( 6. 408)NN0,FFF, JCP.CC.XX .NNN

I F ( NNO .EQ. 0 ) GO TO 21

RFAD( 5*700) ( JOC ( N 1 » J ) »C0( Nl.J) .XO(NltJ) • J = 1 »NNO)

WRITE(6*403)(JOC(N1*J)*CC(N1*J) ,X0( N1 » J ) . J=1 sNNO)

21 DO 24 I =N 1 « N2

NC< I )=1

NCO ( I )=NNO

I FF = I FF+ 1

I FLOOR ( I ) = I FF

J C ( I * 1 ) = JCP+SNCOMP ( IFF)

C ( I * 1 )=CC

X ( I * 1 )=XX

FF ( I )=FFF

I F ( NNO .EQ. 0 ) GO TO 24

DO 22 J = 1 » NNO

JJ=JOC(NlsJ)

P0< Is J ) =PFO ( IFF* JJ)

JOC ( I* J) = JJ

C0( Is J ) = C0(N1 , J ) 22 X0( Is J ) =X0 ( N1 » J

)

24 CONTINUE

Cc

cc

cc

C

EXCEPTIONS TO GENERAL SHAFT INPUT

NNN = NO. OF EXCEPTIONS

KE = 1 FOR FF EXCEPTION

KE = 2 FOR OUTSIDE CONNECTION

KE = 3 FOR INTERNAL CONNECTION

DIO
SUBROUTINE INPUT

I F ( NNN ,EQ. 0 ) GO TO 100

DO 69 I K = 1 » NNN

REAC(5«700)KE, IFF

WRITE(6c409)KE, IFF

I =N S 1 ( IS) +IFF-NFS1 ( IS)

I F ( KE • EQ* 1 ) GO TO 41

I F ( KE • EQ • 2 ) GO TO 42

I F ( KE .EQ. 3 ) GO TO 51

GO TO 104

41 READ(5»700) FF ( I

)

WRITE(6*410)FF( I

GO TO 69

42 READ(5o700) J.CCC.XXO

WR ITE(6,405)

WRITE (6, 403) J,CCO,XXO

NNC=NCO( I

)

I F ( NNC .EQ, 0 ) GO TO 44

DO 43 K= 1 « NNC

I F ( JOC ( It K) .EQ. J ) GO TO 46

43 CONTINUE

44 NJC=NNC+1

NCO ( I )=NJO

47 PO( I e NJO)=PFO( IFF, J)

JOC< I * N JO ) =

CO( I„NJO)=CCO

XO ( I« NJO) = XXO

GO TO 69

46 N JO =K

KK=K+

I F ( CCO .NE. O ) GO TO 47

N J0=NNC-1

NCO ( I )=NJO

I F ( N JO .EQ. 0 ) GO TO 69

DO 49 K=KK,NNC

KM=K-1

P0( Ic KM )=P0( I, K ) JOC ( I * KM ) = JOC ( I , K

C0( I c KM ) =CO ( I ,K

49 X0( I«KM)=XO( I ,K)

GO TO 69

51 READ(5o700)JCP»CC»XX

WR ITE(6»402)

WRITE<6« 403) JCPeCC«XX

J= JC°+SNCOMP( IFF)

NN=NC( I

I F ( NN .EQ. 0 ) GO TO 53

DO 52 K= 1 , NNC

I F ( JC ( T e K ) .EQ. J)GO TO 55

52 CONTINUE

I F( CC .NE. 0. ) G 0 TO 53

WRITE (6, 520) IS.KE, IFF

GO TO 69

53 N J= NN+

NC ( I )=NJ

54 JC( I e N J ) =

C( I t NJ )=CC

X ( I «N J ) = XX * GO TO 69

Dll
* SUBROUTINE INPUT

55 NJ=K

KK=K+1

I F ( C C .NE. O • ) GO TO 54

N J= NN—

1

NC ( I ) = N

J

I F ( N J • EQ. 0 ) GO TO 69

DO 61 K = KK * NN

KM=K-1

JCC Ic KM)=JC( It K)

C( I S KM)=C( I ,K>

61 X ( I « KM ) =X ( I , K

)

69 CONTINUE

100 CONTINUE

NT=N2

I F ( NT .LE. MM ) GO TO 160

I P A R=

1

GO TO 110

C

C PRINT OUTSIDE TEMPERATURE

C

160 WR ITE ( 6, 60 1 ) (TITLE! I ) * 1 = 1 *12)

IFCIUNIT • EQ • 1 )WR ITE( 6,800 )TOUT

IF( IUNIT .EQ, 2)WRITE( 6*5 00 ) TOUT

IFCIUNIT • EQ. 2 )TOUT=( TOUT-32. >/1.8

TOUT=TOUT+273t

C

C PRINT HEIGHT AND TEMPERATURE PROFILES

C

I F ( I UN I T .EQ. 1 )WRITE( 6,811 )( IP, IP=1 ,NTP>

I F ( IUN IT .EQ. 2 ) WR ITE(6,511 ) ( IP, IP=1 ,NTP)

WR I TE (6,813)

DO 30 I FF = 1 , NH

30 WRITE(6o312)H(I FF ) «(T( IP, IFF) , I P=1 , NTP

)

C

C CCNVERT TEMPERATURES TO DEG K

C

DO 33 I FF = 1 , NH

DO 33 I P= 1 , NTP

IFCIUNIT .EQ. 2)T( IP, IFF) =(T( IP, IFF)—32. )/l.8

33 T( IP, IFF)=T( IP, IFF ) + 273.

C

C PRINT OUTSIDE PRESSURE PROFILES

C

IF( IUNIT .EQ. 1 ) GO TO 79

WRITE(6,514)(IP,IP=1,NP0)

WRITEC6, 81 3)

DO 76 IFF=1,NH

DO 77 J = 1 » NPO

77 PFO(IFF, J)=PFO( IFF,J)/248.8

WRITE! 6, 51 5) H( IFF ) , (PFO( IFF , J ) , J=1 , NPO)

DO 78 J= 1 c NPO

78 PFOC IFF, J)=PFO( IFF, J)*248.8

76 CONTINUE

GO TO 83

79 WRITE(6#814) ( IP, IP=1, NPO)

WRITE(6 0 813)

DO 31 IFF=1,NH , WRITE(6#815)H( IFF) ,(PFO(IFF,J) ,J=1,NP0)

D12
* SUBROUTINE INPUT

31 CONTINUE

C

C CORRECT FOR CONNECTIONS ONLY INPUTED ONCE

C

83 00 60 1= 1 « NT

NN=NC ( I

IFCNN .EQ. 0 ) GO TO 60

DO 58 JJ=1.NN

J=JC< Ic JJ

IF( J .EQ. 0 ) GO TO 58

NNJ=NC ( J

)

IFCNNJ .EQ. 0 ) GO TO 57

DO 56 I A= 1 • NN

I F ( JC ( J o I A ) .EQ. I ) GO TO 58

56 CONTINUE

57 NNJ=NNJ+1

I F ( NNJ .LE. MCIGO TO 59

I P A R = 3

GO TO 110

59 NC(J)=NNJ

JC( Je NNJ )= I C( JcNNJ)=C( I* JJ)

X ( Jt NN J ) =X (

9 J J )

I F ( J .GT. N .OR. I .GT. N)GO TO 58

PZ( JsNNJ) =-PZ< I . J J ) 58 CONTINUE

60 CONTINUE

€

C CORRECT UNITS

C

I F ( I UN I T .EQ. 2 ) C A L L UNITS

C

C INITIALIZE PZ FOR BUILD COMPARTMENTS

C

87 DO 40 I =

9 N

NN=NC ( I

I F ( NN .EQ. 0 ) GO TO 40

I A= IT ( I)

I FI=IFLOOR( I

DO 38 J J=1 »NN

J = JC< Ie JJ)

I F J = I FLOOR ( JO

I F ( I F I .EQ. I F J ) GO TO 38

I 8= IT ( J ) TEMPA=0.5*(T(IA,IFI)+T(IB.TFJ))

PZ< I e JJ)=3A62.*<H( IFJ )-H< IF I ) ) /TEMP

38 CONTINUE

40 CONTINUE

C

C INITIALIZE PZ FOR SHAFTS

C

DO 50 IS=1,NS

N 1 = NS 1 { IS)

N2=NS2 ( IS )-l

I TT= I T S ( IS)

DO 45 I =N 1 . N2

I F I =1 FLOOR ( I

I F J= IF 1 + 1

D13
* SUBROUTINE INPUT

TEMPA=0.5*<T( ITTe I F I ) +T ( I TT . I F J )

)

PZ< I« 1 )=3462.*< H( IFJ )-H( IF I ) ) /TEMP

A

45 CONTINUE

50 CONTINUE

C

C CHECK SHAFT CONNECTIONS

C

00 240 I S= 1 » NS

N 1 =NS 1 ( IS)

N2=NS2< IS)

DO 239 I=N1»N2

NN=NC( I

)

I F ( NN .EQ. 0 ) GO TO 239

DO 236 J = 1 » NN

JJ = JC< It J

)

I F ( I FLOOR ( I ) • NE. I FLOOPC J J ) ) GO TO 103

236 CONTINUE

239 CONTINUE

240 CONTINUE

RETURN

C

C

C DIAGNOSTIC OUTPUT

C

102 WRITE (6.902) I , IT( I

)

GO TO 109

103 W R I TE ( 6c 90 3

)

GO TO 109

104 WP ITE ( 6« 904

)

GO TO 109

105 WR I TE ( 6e 905

)

GO TO 109

106 WR ITE( 6* 90 6

)

GO TO 109

107 WRI TE ( 6e 907

)

GO to 109

110 WPITE(6c910)PAR( IPAR)

C

C PRINT CORRECTED BUILDING DATA

C

109 WR ITE (6e 940

)

DO 70 1=1.

N

NN=NC ( I )

I F ( NN .GT. 0 ) GO TO 180

WRI TE ( 6o 94 1)1 , I FLOOR < I)eIT(I).FF(I)

GO TO 182

180 WRITE!

6

0 942)1 .IFLOOP(I).IT(I).FF(I).JC(I.l).C(I.l).X(I.l)

I F ( NN .EQ. 1 ) GO TO 182

WRITE(6 0 943)(JC(I.J)sC(I«J).X(I.J)» J = 2.NN)

182 NNC = NCO ( I ) I F ( NNO .EQ. 0 ) GO TO 70

WRI TEC 6,94 4) ( JOC< I . J) ,C0( I , J) ,X0( I . J) • J=1 .NNO)

70 CONTINUE

C

C PRINT CORRECTED SHAFT INPUT DATA

C

DO 80 I S= 1 » N

S

WRITEC6.816) ( TIT SH (IS«I)*I=1»3)

D14
*

72

74

75

80

Ccc400

401

402

403

404

405

406

40 7

408

409

410

411

412

413

414

415

416

417

418

4 19

420

500

511

514

515

520

600

601

603

700

800

SUBROUTINE INPUT

WRITE! 6* 806) IS* CS! IS)* ITS! IS)

N 1 = NS 1 ( IS)

N2=NS2 ( IS

WRITE!6*807)

DO 75 I =N 1 * N2

NN=NC< I

IF ( NN • GT • 0 ) GO TO 72

WPITE!6»801 )I FLOOR ( I ) * FF ( I

)

GO TO 74

WRITE! 6* 808) I FLOOR! I) • FF! I ) , JC ( I • 1 ) • C ( I • 1 ) , X ( I • 1

)

IF! NN .EQ. 1 )GO TO 74

WRITE!6»809)!JC!IeJ)*C! I * J ) X ( I * J ) * J=2 1 NN)

NNC=NCO! I

IF! NNO • EQ* 0 ) GO TO 75

WRITE!6»810) (JOC! I , J ) , C 0! I » J ) » XC ! I «J) . J=1 .NNO)

CGNTINUE

CONTINUE

STOP

FORMAT STATEMENTS

FORMAT! 5X* 5HIF1 =,I3.7H, IF2 =,I3,7H, NOC =.I3)

F OR MAT ! 5X» 4HN Z =,I3,6H NA = . I 3 , 7H , NNO = ,13, 6H. FF =,F8.1,

+ 7H, IT = • I 3

FORMAT !5Xe 25HCONNECTION ON SAME FLOOR ) FORM AT ! 5Xe 3HJ =*I3*5H, C =*F8.1*5H* X =,F5.2)

FOF MAT! 5X* 26HC0NNECT ION TO FLOOR ABOVE ) FORMAT !5X« 22HC0NNECT ICN TO OUTSIDE ) FOP MAT !5X*3A6)

FORMAT !5X. 4HCS =.F9.1,8H, NFS1 =,I3,8H, NFS2 =,I3,7H, ITS = ,I3>

FORMAT !5X» 5HNN0 =,I3,7H, FFF =,F8.1*5H* J = *I3*5H* C = * F8,l*

+ 5 He X =*F5.2*7H* NNN = *13)

FOPMAT (5X* 4HKE =,13. 7H, IFF =.13)

FORMAT! 5X« 4HFF =,F8.1)

FORMAT (5X* 6HT0UT =,F6.0.9H, IUNIT =,I3.8H, IOUT =.I3)

FOPMAT !5X« 4HNH =*I3*6H« NN =*I3)

FORMAT! 5X* 7HHEIGHTS /!10F8#2))

FOPMAT !5X* 6HH! 1 ) =,F8,2,6H, DH =,F8.2)

FOPMAT !6X* 5HNTP =.13)

FORMAT!5X* 20HTEMPERATURE PROFILE /I 5 , ! 1 0 ! I 4 , F7 . 1 ) ) ) FORMAT !5X,5HNP0 =,I3,6H, NN =,I3)

F0RMAT!5X*5HPGZ =*F12. 1/1 7HPRE SSURE PROFILE /!10F12«1))

FOPMAT (5X« 4HVW =*F6*1*6H* HW =.F6.1,6H. XW =.F4,2.6H. CW =. + ! 1 0F4 • 2 )

FOPMAT !/5X. 6HNFLS =,I3)

FORMAT I//10X.20HOUTSIDE TEMPERATURE .F6.1.2H F)

FORMAT! ///5X « 6HHE I GHT , 5X • 2 9H TEMPERATURE PROFILES !DEG F) /

+ 7X ,2HFT, 3X* 191 6)

F OP MAT ! ////5X * 6HHE I GHT * 5X* 26HOUTSI DE PRESSURE PROFILES

1 1 1 H ! IN H20 ) /7X , 2HFT .3X.8I10)

FORMAT! FI 1.2.3X.8F10.3)

FOF MAT!/// 5X . 1 5 HE PR OR IN SHAFT . I 2 . 1 5HEXCE PT I ON KE = ,12.

+ 2X.5HFL00R ,13//)

FORMAT! 1 2A6)

FORMAT! 1H1///10X, 12A6///)

FORMAT! 3A6* 12

FORMAT! ) F OR MAT ! //10X « 20H0UTSIDE TEMPERATURE .F6.1.2H C)

D15
SUBROUTINE INPUT

801 FORMAT! I13.F11.1

)

806 FORMAT! 1 OX * 1 2H SHAFT NUMBER . I 4/ 1 OX , 1 7H SHAFT COEFFICIENT .F10.1/

1

10X.20HTEMPERATURE PROFILE .14)

807 F0RMAT(/21X. 5HF IX ED* 25 X# 4HFL0W. 12X. 4HFL0W/ 1 OX . 5HFL00R . 6X

1 4HFLQW.5X.1 2H CONNECTED TO • 6X . 1 1 HCOEFF I C IENT . 6X . SHEXOONENT

2 /

)

8 08 FORMAT! I 13* FI 1. 1. 6X. 5HPOI NT . I 5* F 1 6. 1 • F 1 2. 2

)

809 FORMAT! 30 X. 5HPO INT .I5.F16.1 .F12.2)

810 FORMAT!30X. 7HOUTSI DE . 13 • FI 6* 1 • F 12 • 2

)

811 FORMAT! ///5X.6HHEIGHT. 5X * 29HTEMPER ATURE PROFILES !05G C) /

+ 7X.2HM « 3 X. 1916)

812 FOFMAT<F 1 1.2. 3X. 19F6. 1

)

813 FORMAT!/)

814 F OR MAT I////5X » 6HHE I GHT , 5X , 26HOUTS I DE PRESSURE PROFILES

1 1 1 H < PASCALS ) /7X « 2HM .3X.8I10)

815 FORMAT (FI 1 .2. 3X » 8F 10*1)

816 FORMAT! /// 10X * 3A6 ) 817 FOFMAT! 10X. 45HFL0W COEFFICIENTS CORRECTED FOR TEMPREATURE ) 902 FORMAT! 10! /) , 10X. 1 1 HCOMPA RTMENT .14/

1 10 X. 20 HTEMPERATURE PROFILE .I4.17H DOES NOT EXIST /

+ 10X* 16HPR0GRAM STOPPED ,10!/))

903 FORMAT! 10!/), 5X. 23HSHAFT CONNECTION ERROR .

1 /I OX. 16HPR0GRAM STOPPED .10!/))

904 FORMAT! 10!/). 10X, 40HINPUT ERROR IN EXCEPTIONS TO SHAFT DATA

1 /I OX. 16HPR0GRAM STOPPED .10!/))

905 FORMAT! 10(/)o 10X. 37HINPUT ERROR IN UNIT TYPE DESIGNATION /

1 1 OX. 16HPR0GRAM STOPPED .10!/))

906 FORMAT! 10!/). 10X. 37HINPUT ERROR NO. OF FLOORS EXCEEDS NH /

1 10X. 16HPR0GRAM STOPPED ,10!/))

907 FORM AT (10(/)»10X* 2 5H INPUT ERROR IF1 .GT. IF2 /

1 10X. 16HPR0GRAM STOPPED ,10!/))

910 FOF MAT! 10! /)» 10X.36HI NPUT EXCEEDS DIMENSION PARAMETER . A 3/

+ 10X. 16HPR0GRAM STOPPED .10!/))

930 FORMAT! 10X. 3A6)

935 FORMAT!// 10X.26HFL0W COEFFICIENTS AS READ ) 940 FOF MAT! 10X. 15HBUILDING DATA //34X * 1 1 HTEMPER ATURE . 4X . 5HFIXED.

1 1 2X, 2! 1 1 X,4HFLOW)/10X, 1 1 HCCMPA RTMENT . 4X® 5HFL00R. 6X. 7HPR0F ILE.

2 6X.4HFLOW. 5X. 1 3HCONNEC T I ON TO . 4X. 1 1 HCOEFF IC IE NT .4X,

3 8HEXP0NENT )

» I3.F1 1. 1.F13. 1

)

941 FORMAT

942 FORMAT

943 FORMAT

9 44 FORMAT

END

7DHDG.P * SUBROUTINE CORR L, 1

D16
n

n

\o

oo

>i

o>ui

4

>oj

ro

-

on

* SUBROUTINE CORR

3NBS*°L IBJ.SHOW A. CORR

SUBROUTINE CORR

C

C THIS ROUTINE CALCULATES FLOW COEFFICIENTS

C WHICH ARE CORRECTED FOR TEMPERATURE

C

°AR AMETEP MM =2 20. MS= 2* MC=4 , MP0=2 * MTP = 2 • MFL= I 05 . M B= 1 05

COMMON /CORR/C1 ( MM. MC ) , C2 (MM ,MC ) .CO 1 < MM.MPO) .C02( MM.MPO)

COMMON NT. PC MM) ,C( MM.MC) ,NC( MM) , JC( MM.MC) . IT3( MS ).

1 FCCMM, MC) .PZCMM, MC) , POCMM. MPO) .CO( MM.MPO) , F ( MM ) , PFOC MF L, MPO) . 2 FF ( MM) .FO( MM, MPO) ,CS( MS) ,PS < MFL ) • NS1 ( MS) , NS2 ( MS)

,

3 FSS( MS ) . N. NS . NPO , I CON V. E , I BUG , X ( MM , M C ) , X O ( MM , MPO ) • TI T SH< MS • 3 ) • 4 NH.HC MFL ) , IFLOORC MM) ,T< MTP.MFL) ,NFS1 ( MS) . NFS2C MS). IT C MB) ,NTP

5 . NCOC MM) . JOC( MM. MPO) .TOUT

DOUBLE PRECISION P.PO.PS

DO 12 1=1 . NT

C

CORRECT C

NN=NC ( I ) I F ( I .GT. N ) GO TO 1 IP= IT ( I ) GO TO 4

DO 2 I S=1 . NS

I F ( I ,LE. NS2 (IS) .AND. I • GE . NSl(IS))GO TO 3

CONTINUE

WP ITE (6. 700

STOP

IP= IT S ( IS)

I FF = I FLOOR ( I ) T 1 =T ( IP, IFF)

I F ( NN .EQ. 0 ) GO TO 10

DO 9 J=1 , NN

J J = JC ( I. J ) Cl(I.J)=C(l9J)*((294./Tl)**<2.5*X(I.J)-0.75>)

I F ( J J .GT. N ) GO TO 5

I P= IT ( J J

GO TO 8

DO 6 I S = 1 9 NS

I F ( J J .LE. NS2 (IS) .AND. JJ .GE. NS1(IS))G0 TO 7

CONTINUE

WR I TE ( 69 700

STOP

IP= IT S ( IS)

IFF=IFLOOR( JJ)

T2=T( IP, IFF)

C2(I.J)=C(I.J)*((294./T2)**(2.5*X(I,J)-0.75))

CONTINUE

CORRECT CO

C

10 NNC=N CO ( I

I F ( NNC .EQ. 0 ) GO TO 12

DO 11 J = 1 , NNC

COl (I.J)=CO(I.J)*( (294./T1) **(2.5*XO(I.J)-0.75))

C02 (I 9 J ) = CO( I , J )*( ( 294 ./TOUT) ** (2.5*X0( I , J )— 0.75)

)

11 CONTINUE

D17
* SUBROUTINE CORR

12 CONTINUE

RETURN

700 F0RMATC///10X, 36HPR0GRAM STOPPED IN SUBROUTINE CORR

END

3) HOG 9 P * SUBROUTINE INIT .L.

//)

1

Dl 8
* SUBROUTINE INIT

®NBS*PL IBS. SHOW A. INIT

SUBROUTINE INIT

C

c

C THIS ROUTINE INITIALIZES THE PRESSURE ARRAY

C

PARAMETER M M =2 20 . MS= 2 * MC =4 * MP0=2 . MTP=2 • MFL=1 05.MB=105

PARAMETER M8P=MB+1

COMMON NT. P(MM).C(MM.MC).NC(MM). JC(MM.MC). ITS(MS).

1 FC(MM.MC).PZ(MM.MC ) oPO(MM. MPO) *C0( MM. MPO) . F ( MM ) , PF 0{ MFL . MPO)

2 FF (MM) «F0( MM, MPO) ,CS< MS) . PS ( MFL) .NS l ( MS ) . NS2(MS)

3 FSS(MS).N.NS , NPO . I CON V. E , I BUG . X ( MM , MC ) , X 0 ( MM , MPO) . TI T SH ( MS . 3 ) . 4 NH » H ( MFL ) . I FLOOR ( MM) . T( MTP.MFL) .NFS1 ( MS) . NFS2( MS) . IT (MB) . NTP

5 . NCO< MM) . JOC ( MM. MPC) . TOUT

DOUBLE PRECISION P.PO.PS

DIMENSION SC( MS) .SCO(MS)

COMMON /MAT/A(MB*M BP ) . XX ( MB ) » NNN

DOUBLE PRECISION A, XX

NNN= N e

C CALCULATE AVERAGE OUTSIDE PRESSURE

e

SUM =0

DO 10 J= 1 . NPO

DO 10 1=1. NH

10 SUM=SUM+PFO( I , J

PA=SUM/(NPO*NH)

C

cc

C THE DO LOOP TO STATEMENT 30 ESTIMATES

C SHAFT PRESSURES

C

DO 30 I S= 1 » NS

C

C CALCULATE SHAFT PRESSURE DIFFERENCE . DP

C

SUM =0

SUMN=0

SUM X=0

NX = 0

N 1 = NS 1 ( IS)

N2=NS2 (IS)

DO 18 I=N 1 . N2

SUM=SUM+FF< I ) NN=NC< I )

I F ( NN .EQ. 0. ) GO TO 16

DO 15 J= 1 . NN

SUMX=SUMX+X ( I . J

)

NX=NX+1

SUMN=SUMN+C( I. J)

15 CONTINUE

SC( IS ) =SUMN

16 NNC=NCO(I)

I F ( NNO .EQ. 0 ) GO TO 18

DO 17 J= 1 . NNO

S UM X = SUMX+ XO ( I . J)

D19
n

n

oj

ro

non

* SUBROUTINE INIT

N X= NX+1

SUMN=SUMN+CO< I « J

)

17 CONTINUE

SCO ( IS) =SUMN— SC ( IS ) 18 CONTINUE

A X = NX/SUMX

DP2=SUM/SUMN

S I GN= 1

IP( DP 2 .LT. 0 • ) S I GN=— 1

DP= S IGN*(DP2**AX)

C

C CALCULATE AVERAGE TEMP OF SHAFT

C

SUM =0

IP= ITS ( IS)

DO 20 I=N1,N2

IFF=IFL00R< I

)

20 SUM=SUM+T ( IP, IFF)

TA=SUM/(N2— Nl+1

)

C

C ESTIMATE PRESSURE AT BOTTOM OF SHAFT • C

HH=0.5*(H(NH)-H<1 ) )+H( 1

)

NF 1 =NF S 1 < IS)

P80T=PA+DP+3462.* ( HH-HINFl ) )/TA

ESTIMATE OTHER SHAFT PRESSURES

P ( N 1 ) =PBOT

NM=N2—

1

DO 24 I=N1»NM

IP1=I+1

4 P ( I PI )=»< I )-PZ( 1,1)

0 CONTINUE

END OF SHAFT PRESSURE ESTIMATES

C

C SET UP MATRIX FOR BUILDING COMPARTMENTS

C

NP 1 =N+

1

DO 50 1 = 1 o N

NN = NC( I ) SUM I 1 = 0.

SUM NP = 0

I F ( NN .EQ. 0 • ) GO TO 42

DO 40 JJ = 1 , N N

J = JC( I, JJ

)

I F (

J

.GT. N ) GO TO 34

A ( I

«

J )=C( I • J

J

) SUM I I = SUM I I-C( I , JJ ) SUMNP=SUMNP-C( I , J J )*PZ ( I , J J ) GO TO 40

34 SUM II=SUMI I-C< I , JJ ) S UMNP = SUMNP— C ( I » J J ) *P< J)

40 CONTINUE

42 NNC=NCO( I

)

I F ( NNO .EQ. 0 ) GO TO 46

DO 45 K=1 , NNO

PBOT

D20
* SUBROUTINE INIT

SUM 1 1 = SUM I I — CO ( I.K)

45 SUMNP = SUM NP—CO (

t K)*PO( I «K)

46 A ( I « I ) = SUM I I A ( I . NP 1 )=SUMNP-FF( I)

50 CONTINUE

C

C WRITE MATRIX

C

I F ( I BU G • EQ • 0 ) GO TO 84

WR I TE ( 6 « 802

DO 52 1=1.

52 WRI TEC 6. 803) ( A< I . J ) . J=1 ,NP1

)

C

C

C CALL ROUTINE TO SOLVE FOR INITIAL BUILDING PRESSURES

C

84 CALL SIMEQ

C

C OUTPUT INITIAL PRESSURES

C

I F ( IBUG .EQ. 0 ) GO TO 89

WRITE(6» 800

WRITE (6*801)(I. XX (I)*I=1.N)

NN=NS 1 ( 1

)

WPITE(6.801 ) ( I»P( I ). I=NN» NT)

C

C

C ASSIGN BUILDING PRESSURES

C

89 DO 90 1=1.

90 P< I ) =XX { I ) RETURN

800 FORMAT! ///8<6X • 1HI , 4X.3HP )/)

801 FORMAT (8(17. F7.1 )

802 FORMAT ( ///l OX . 20HMATR IX COEFFICIENTS /)

803 FOF MAT ( 10X. 1 1 FI 1 . 1 ) END

SHOGs P * SUBROUTINE BLDGP

D21
nnnnronn

nn

n

r>

r>n

* SUBROUTINE BLDGP

$NBS*PL I B$. SHOW A. BLDGP

SUBROUTINE BLDGP

C

c

C THIS ROUTINE CALCULATES STEADY STATE PRESSURES

C FOR BUILDING COMPARTMENTS

C

cc PARAMETER MM=220

MS=2 * MC=4 , MPO =2

«

MTP= 2

•

MFL=1 05

.

M B= 1 0

5

COMMON NT p P( MM) #C(KM« MC) .NC(MM) . JC( MM.MC) , ITSCMS ),

1 FC <MM« MC) « PZ (MM, MC) ,PO( MM, MPO) ,CO< MM, MPO) , FC MM ) , PF 0( MFL, MPO) , 2 FF ( MM) c FO( MM, MPO ) ,CS< MS) ,PS < MFL ) , NS 1 < MS ) ,NS2(MS),

3 FSS ( MS ) * N, NS ,NPO, ICONV.E, I BUG , X ( MM , MC ) , X 0( MM , MPO ) , T I T SH ( MS , 3 ) , 4 NH , H ( MFL ) ,IFLOOR( VM) , T< MTP , MFL ) ,NF SI (MS) ,NFS2(MS>* IT(MB).NTP

5 c NCO( MM), JOC( MM, MPO) » TOUT

DOUELE PRECISION P,PO,PS,PI

I F ( I BUG • GT• 0 )WRITE( 6,806)

ITM=20

I CONV=0

DO 15 1=1,

N

CALCULATE NET FLOW ,FI, INTO POINT I F I = PFLOW ( I , P ( I ) ) CHECK MAGNITUDE OF FI

IF(AeS(FI) .LT. E ) GO TO 15

I C0NV=IC0NV+1

SET UP PARAMETERS FOR ITERATION

3D=1.0

I PH A SE = 1 DP 1=0.

EE=0.2*ABS< FI

)

I F ( EE .LT. E ) EE=E

S IGN=1

I F ( F I .LT. 0.)SIGN=-1

I K =0

I F ( I BUG .GT. 0) WRITE(6. 802)

ITERATION TO RE DU I CE MAGNITUDE OF FN

I K= IK+1

NEW ESTIMATE OF PRESSURE ,PI, AT POINT I P I=P< I )+SIGN*DP

CALCULATE NET FLOW , FN , INTO POINT I USING PI

FN = FFLOW ( I , P I

)

I F ( I BUG. GT.O) WRITE (6, 80 4) I , I K , F I , FN ,FP , DPI , DP , DPP , P I , I PH A SE

C

C CHECK MAGNITUDE OF FN

I F ( ABS ( FN ) .LT. EE ) GO TO 10

C

C CHECK NUMBER OF ITERATIONS

I F ( I K .GT. I TM ) GO TO 25

C

C CHFCK PHASE

D22
* SUBROUTINE BLDGP

IF( IPHASE .EQ. 2 ) GO TO 6

C

C CHECK FOR TRANSITION FROM PHASE 1 TO PHASE 2

I F( F I *FN • LT• 0 • ) GO TO 4

ec

C PHASE 1 DPI=DP

DP=5.0*DP

F I = FN

GO TO 2

C

C °HASE 2

4 IPHASE=2

GO TO 9

6 IF(FI*FN #GT. O.IGO TO 8

C

C NEW DP BETWEEN DPI AND DP

9 DPP=DP

FP=FN

DP= DPI 4- (DPP—DP I >*FI/IFI-FN)

GO TO 2

Cr

C NEW DP BETWEEN DP AND DPP

8 FI=FN

DPI =DP

DP=DPI+ (DPP-DPI )*FN/(FN—F°)

GO TO 2

10 P(I)=PI

15 CONTINUE

C

RETURN

25 WPITE<6*800>

STOP

C

C FORMAT STATEMENTS

C

800 F OF M AT ( /// 1 OX * 20( 1 H* ) //! OX • 22HEXCESSI VE ITERATIONS /

+ 1 0 X* 8H IN BLDGP /// 10X * 20 < 1 H* ) //////)

002 c OP MAT ( //I IX* 1HI*2X»2HIT*12X »2HF I * 1 3X « 2HFN * 1 3X * 2HFP « 1 2X * 3 HOP I

*•1 3X • 2HDP • 1 2X * 3HDPP « I3X.2HPI «3X«5HPHASE /)

604 FOFMAT( 9X* 214, 3E1 5,4,4F15,6* 15)

806 FORMAT! /// 1 0 X, 6HBLDGP ) END

5)HDGo R * SUBROUTINE SHAFTP ,L,1

D23
nnnnnn

n

n

o

o

i\> n

o

on

no

on

* SUBROUTINE SHAFTP

<DNBS*°L IB$, SHOW A. SHAFTP

SUBROUTINE SHAFTP

C

C

C THIS ROUTINE CALCULATES STEADY STATE PRESSURES

C FOR SHAFTS

C

c PARAMETER MM =220 , MS=2 , MC =4 , MPO=2 , M TP=2 , MFL =1 05 ,M8=105

COMMON NT, P< MM) ,C< MM ,MC) ,NC(MM), JC< MM, MC)# ITS(MS )•

1 FC(MM»MC) *PZ(MM*MC) • PO ( MM* MPO ) , CO ( MM • MPO ),F( MM ) , PFO< MFL, MPO) , 2 FP (MM),FO(MM,MPO) ,CS( MS) . PS C MFL ) , N SI ( MS) . NS2(MS)

.

3 FSS( MS) • N* NS ,NPO, ICONV,E.I8UG,X<MM,MC),XO(MM,MPO),TITSH<MS.3),

4 NH,H< MFL) * I FLO OR (MM) »T(MTP , MFL ) ,NF SI ( MS) , NFS2C MS ) , IT < M B ) , N TP

5 ,NCO(MM), JOC(MM,MPO)*TOUT

DOUBLE PRECISION P,PO,PS,PI

I F ( I BUG • GT 9 0 )WRITE(6,806>

I TM =2 0

I C0NV=0

DO 15 1=1, NS

CALCULATE NET FLOW ,FI, INTO POINT I N 1 = NS 1 ( I

)

FT = SFLOW( I » P ( N 1 )

)

CHECK MAGNITUDE OF FI

IF( ABS(FI ) ,LT. E ) GO TO 15

ICCNV=ICONV+l

SET UP PARAMETERS FOR ITERATION

DP= 1 •

0

I PH ASE = 1 D P I =0

•

EE=0.2*ABS(FI

)

I F( EE .LT, E)EE=E

S I GN=

1

I F ( F I , LT • 0.)SIGN=—

1

IK=0

I F ( I FUG • GT • 0) WRITE(6,802)

ITERATION TO REDUICE MAGNITUDE OF FN

I K= IK+1

NEW FST 1 MATE OF PRESSURE ,PI, AT BOTTOM OF SHAFT I P I=P(N1 )+SIGN*DP

CALCULATE NET FLOW , FN , INTO SHAFT I USING PI

FN=SFLOW< I , P I

)

IF< IBUG. GT.O ) WRITE < 6,804) I, I K , F I , FN , FP , DPI , DP , DPP , P I , I PHASE

CHECK MAGNITUDE OF FN

I F ( ABS ( FN ) *LT • EE ) GO TO 10

CHECK NUMBER OF ITERATIONS

IFC1K , GT, I TM ) GO TO 25

CHECK PHASE

D24
no

o

o

* SUBROUTINE SHAFTP

I F ( IPHASE .EQ. 2 ) GO TO 6

CHECK FOR TRANSITION FROM PHASE 1 TO PHASE 2

I F( F I *FN ,LT* 0 • ) GO TO A

PHASE 1 DPI=DP

0P=5.0*DP

F I =FN

GO TO 2

C

C PHASE 2

4 I PH A SE = 2

GO TO 9

6 IF<FI*FN .GT. 0 . ) GC TO 8

C

C NEW OP BETWEEN DPI AND DP

9 DPP=DP

FP= FN

DP=DPI+(DPP-DPI )*FI/(FI-FN)

GO TO 2

C

C NEW DP BETWEEN DP AND DPP

8 F I = FN

PI =DP

D P= DPI + ( DPP—DP I )*FN/(FN-FP)

GO TO 2

10 N2=NS2 ( I

)

DO 11 I F =N 1 1> N2

I 1= IF+1-N1

11 P ( I F ) =PS ( II)

15 CONTINUE

C

RETURN

25 WPITEC 6s 800

STOP

C

C FORMAT STATEMENTS

C

800 FORMAT (///10X»20( 1 H* )///10X* 2 2 HE XCE SSI VE ITERATIONS /

+ 1 OXo 9HIN SHAFTP /// 1 OX . 20 ( 1 H* )//////

)

8 02 FOFMAT< //I IXo lHIo2X„2HIT# 1 2 X » 2HF I « I 3X.2HFN, 13X, 2HFP « 1 2X * 3 HD PI

+ 1 3X o 2HDP* 1 2X « 3HOPP « 13X, 2HPI ,3X.5HPHASE /)

804 FOP MAT <8X. 2I4 # 3E15.4,4F15*6tI5)

806 FORMAT ( ///10X* 6 H SHAFTP

END

tDHDGsP * SUBROUTINE PZAD .L«l

D25
SUBROUTINE PZAD l *

aNBS*PLIBS.SHOW

SUBROUTINE

A. PZAD

PZAD

C

C

c

2

1

8

10

1 5

20

800

801

802

803

THIS POUTINE CORRECTS PZ TERMS FOR PRESSURE

PARAMETER MM=220 . MS=2 , MC=4 . MP0=2 . MTP=2 * MFL=1 05 . MB= 1 05

COMMON NT, P( MM ) «C( MM»MC) »NC( MM) • JC( MM* MC) * ITS(MS1»

1 FCtMMg MC) « PZ (MMC MC)*PO(MM, MPO) tCOC MM, MPO) *F( MM ) # PFO( MFL# MPO)

«

2 FF ( MM ) e FO ( MM , MPO ) ,CS( MS) ,PS< MFL ) ,NS1 ( MS) , NS2 ( MS )

3 FSS(MS),N*NS »NPO , I CONV , E , I BUG, X ( MM , MC ) » XO ( MM , MPO) , TIT SHI MS ,3 )

,

4 NH,H< MFL) ,IFLOOR( MM) , T ( MTP , MFL ) , NFS1 C MS) , NFS2<MS)t IT (MB) ®NTP

5 pNCO(MM)# JOCC MM, MPO) # TOUT

COMMON /PZZ/ PGZ

DOUBLE PRECISION PgPOpPS

I F ( IBUG .GT. — 2 ) GO TO 1 WRITE( 6, 800)

DO 2 1 = 1 p N

NN=NC .( I ) I F ( NN .EQ. 0 ) GO TO 2

W R I TE ( 6 , 80 1 )( I , J , P Z ( I • J ) • J=1 ,NN

)

CONTINUE

NP 1 =N+

1

WRITE <6© 80 2) ( IL*PZ( IL* 1 ), IL=NP1 , NT

)

DO 10 1=1 c N

NN=NC< I

)

I F ( NN • EQ* 0 ) GO TO 10

I A= I T ( I ) I F 1= IFLOOR ( I

)

DO 8 J J=1 » NN

J=JC( Ip JJ)

I FJ = IFLOOR ( J

)

I F( I F I • EQ • I F J ) GO TO 8

I B= IT ( J

)

TEMPA=0.5*(T( IA»IFI)+T(IB»IFJ) ) PAVE=0.5*(P(I )+P( J) )+PGZ

PZ( I. JJ)=(0.03416*PAVE/TEMPA)*(H< IFJ)-H( IF I )

)

CONTINUE

CONTINUE

DO 20 I S=1 , NS

N 1=NS 1 ( IS)

N2=NS2( I S ) —

1

I TT= ITS ( IS)

DO 15 I =N 1 , N2

I FI= IFLOOR ( I

)

IFJ=IF 1+1

TEMPA=0.5*(T( ITT, I FI )+T< ITT.IFJ )

)

J= I + 1 PA = 0. 5* <P< I >+P< J) ) +PGZ

PZ( I, 1 ) = <0.03416*PA/TEMPA)*(H(IFJ)-H< IF I) ) CONTINUE

RETURN

p OP MAT (/10X, 10HINITIAL PZ /)

FORMAT ( lOXp 3HPZ (, 12, 1H, 12, 4H) = ,F12.4)

p OR MAT( lOXp 3HPZ ( , I 2 , 6H . 1 )• ,F12.4)

FORMAT ( / 10X , 1 1 HADJUSTED PZ /)

END

D26
* SUBROUTINE OUT

5>NBS*PLIB$.SH0W A. OUT

SUBROUTINE OUT

Cc

C THIS ROUTINE OUTPUTS FLOWSAND DIFFERENTIAL PRESSURES

C FOR ALL SHAFTS AND BUILDING COMPARTMENTS

C

c PARAMETER MM=220 * MS=2 . MC =4 ,MP0=2 « MTP=2 , MFL = 1 05 . MB= 1 05

COMMON /CORR/C1 (MM , MC I , C2 < MM , MC ) • CO 1 ( MM, M PO) . C02 ( MM. MPO

COMMON /IO/T ITLE( 12 ). I OUT. IUN I T. NC OMP ( MFL ) . SNC OMP ( MFL

COMMON NT* P(MM) . C(MM« MC) sNC(MM). JC( MM. MC) . ITS( MS).

1 FC < MM. MC ) .PZ(MM. MC) . PO( MM. MPO) .CO(MM.MPO) ,F( MM ) • PFO( MFL. MPO)

2 FF (MM ) , FO ( MM , MPO ) *CS( MS) .PS(MFL)»NS1 ( MS) » NS2 ( MS ) . 3 FSS(MS).N.NS »NPO .ICONV.E.I BUG* X(MM»MC) »XO( MM « M PO )»TITSH(MS»3)

4 NH,H( MFL ) , IFLOOR ( MM) ,T( M TP, MFL ) .NFS1 ( MS) . NFS2< MS ) . IT ( M 8) . NTP

5 cNCO(MM). JOC(MM, MPO). TOUT

DOUBLE PRECISION P.PO.PS

INTEGER COM

C

C I UN IT = 1 FOR SI UNITS

C I UN I T = 2 FOR ENG UNITS

C WHEN IUN I T = 2 GO TO 100

I F ( I UNIT .EQ. 2 ) GO TO 100

C

C BUILDING COMPARTMENT OUTPUT

C

1

1=0

I L= 0

WRITE (6.800) ( T I TLE ( I) .1=1 .12)

DO 30 IFF=1,NH

NNN=NCOMP( IFF

I F ( NNN .EQ. 0 ) GO TO 30

DO 29 I C=1 . NNN

1 = 1 + 1 NN=NC ( I

NNC=NCO( I

)

I L= IL+NN+NNO+2

I F ( IL *LT. 51 ) GO TO 2

WRITE( 6. 800 ) (TITLE! I ) .1 = 1 .12)

I L= NN + NNO+ 2

2 IF(NN .GT. 0)GO TO 3

WR ITE(6.801 ) IFF ,IC,P( I ) , IT( E ) ,FF(I

GO TO 21

3 DO 20 J=1 . NN

J J = JC ( I. J ) DP = P ( J J ) — P ( I ) + PZ< I , J)

CC=C2 ( I. J ) I F ( DP .LT. 0. ) CC=C 1 ( I . J

)

I F ( J J .LE. N ) GO TO 10

DO 5 I S= 1 . NS

I F ( J J .GE. NSl(IS) .AND. JJ .LE. NS2(IS))G0 TO 6

5 CONTINUE

6 IF(J .GT. l)GO TO 7

WR ITE(6,802) IFF,IC.P( I ) , I T( I ) ,FF( I ) . (TlTSH( IS.K ) . K= 1 .3)

+ ,DP.CC«X( 1,1 ) * FC (1,1)

GO TO 20

D27
* SUBROUTINE OUT

7 WRITE (6 9 80 3) ( TI TSH ( I S* K ) • K= 1 1 3 ) • DP. CC • X ( I . J ) . FC ( l • J

)

GO TO 20

10 IFJ=IFLOOP( JJ)

COM=J J-SNCOMP ( IFJ

)

I F ( J • GT • 1 ) GO TO 12

WRITE(6,804)IFF»IC»P(I )• IT( I) »FF ( I ) . I F J ,C OM«DP , CC , X ( I . 1 ) ,FC( I .1

)

GO TO 20

12 WRITE (6. 805) I F J .COM , DP , CC tX ( I , J) ,FC( I. J)

20 CONTINUE

21 I F ( NNO • EQ • 0 ) GO TC 29

DO 23 J=1 « NNO

JJ=JCC< It J)

DP=PO( Ic J )-P( I

)

CC= C02< It J)

I F ( OP .LT. O. )CC=C01 ( I • J

)

23 WRITE(6»806)JJ» DP ,CC,XO(I «J)«F0(I,J)

29 WRITE(6 # 807)F( I

)

30 CONTINUE

C

C SHAFT OUTPUT

C

IL= IL+6

DO 60 IS=1,NS

N1=NS1 ( IS)

N2=NS2( IS)

IL= IL+N2+ 1 8— N

1

I F ( IL .LT • 56 ) GO TO 32

WP I TE (6*814)

IL=N2+18-N1

32 WRITE (6, 308) (TITSH( IS.K). K=1 .3)«ITS(IS).CS(IS)

DO 50 I=N1,N2

NN=NC( I

)

IF ( NN .GT, O ) GO TO 35

WRI TE( 6,809) I FLOOR ( I ) * P ( I ) , FF ( I

)

GO TO 41

35 DO 40 J = 1 9 NN

JJ=JC( Ie J)

DP=P( J J )—P ( I

)

CC=C2(

1

9 J

)

I F ( DP .LT. 0. )CC=C1(I,J)

I F J = I FLOOR ( J J

)

COM= J J— SNCOMP ( IFJ

)

I F ( J .GT. 1 ) GO TO 36

WRITE(6»810) I FLOOR ( I ) , P ( I ) , FF ( I ) , I F J , COM, D P ,CC , X ( I , 1 ) , FC ( I , 1 ) GO TO 40

36 WRI TE (6o 81 1 ) IF J,COM,DP,CC , X ( I , J ) , F C ( I , J

)

40 CONTINUE

41 NNO =NCO ( I ) I F ( NNO .EQ. 0 ) GO TO 50

DO 46 J = 1 , NNO

JJ= JCC( I. J ) DP=PO ( I c J >-P( I ) CC=C02 ( I . J ) I F ( CP .LT. 0. )CC=C01 ( I , J

)

4 6 WRITF(6,812) J J , DP , CC , XO ( I ,J) ,FO(I, J)

50 CONTINUE

WRITF(6c813)FSS( IS)

60 CONTINUE

D28
SUBROUTINE OUT

GO TO 165

c

C BUILDING DATA OUTPUT FOR IUNIT = 2

C

100 1=0

I L = 0

WRITE! 6,800) (TITLE! I) ,1=1 ,12)

DO 130 IFF = 1 « NH

NNN=NCOMP< IFF

IF! NNN • EQ • 0 ) GO TO 130

DO 129 I C= 1 » NNN

1=1 + 1 FF I =F ( I )/0.47l9

PI I I = P( I ) /248 •

FFF =FF ( I J/0.4719

NN=NC ( I

NNC=NCO ( I

)

I L= IL+NN+NNO+ 2

I F ( IL ,LT • 51 ) GO TO 102

WRITE (6,800) (TITLE! I) ,1=1 ,12)

I L = NN+NNO+ 2

102 IF! NN • GT • 0 ) GO TO 103

WP ITE(6, 60 1 ) I FF , I C , P 1 1 I , I T( I) ,FFF

GO TO 121

103 DO 120 J=1 ,NN

F CC C=FC ( I e J ) /O • 47 1

JJ=JC( I, J)

0P= ( P ( J J ) — P( I )+PZ( I, J))/248.8

CC=C2 ( I* J

IF! DP ,LT. 0, ) CC=C 1 ( I , J

CC=CC*2, 119/(0«004019**X( I , J )

)

I F ( J J ,LE. N ) GO TO 110

DO 105 IS=1 ,NS

IF! J J .GE, NSl(IS) .AND. JJ • LE • NS2(IS))G0 TO 106

105 CONTINUE

106 IF! J • GT • 1 ) GO TO 107

WRITE(6,602) I FF , IC,PIII ,IT( I ) , FFF ,(TITSH( IS,K) ,K=1 ,3)

+ »DP,CC*X(I ,1 ) ,FCCC

GO TO 120

107 WRI TE! 6*603) (TIT SH ( I S® K ) , K= 1 , 3 ) ® DP , CC • X ( I , J ) , FC CC

GO TO 120

110 IFJ=TFLOOR( JJ)

COM=J J—SNCOMP ( IFJ

)

IF! J • GT • 1 ) GO TO 112

W R ITE! 6, 60 4)IFF,IC»PIII,IT( I), FFF ,IFJ,C0M,DP,CC,X(I,1),FCCC

GO TO 120

112 WPITE(6» 605 ) IFJ, COM, DP, CC »X(I,J),FCCC

120 CONTINUE

121 IF! NNO • EQ • 0 ) GO TO 129

DO 123 J=1 ,NNO

FCO=FO( I, J ) /O ,471

J J= JCC! I • J

0°= (PO( I« J)-P( I ) )/248,8

CC = C02( I® J

I F ( CP ,LT, 0 • ) CC=C01 ( I , J

)

CC=CC*2. 119/! 0, 0 04019*^X0(1 ,J) ) 123 WRITE!6,606)JJ,DP,CC.X0(I,J) ,F00

129 WRITE(6»807)FFI

D29
* SUBROUTINE OUT

130 CONTINUE

C

C SHAFT OUTPUT FOR IUNIT = 2

C

IL = IL + 6

00 160 I S= 1 , NS

CSS=CS( IS) /O • 02992

FF I=FSS( IS) /O *4719

N 1 = N S 1 C IS)

N2=NS2< IS)

IL= IL+N2+18-N1

I F ( I L .LT. 56 ) GO TO 132

WP ITE<6,81 4)

I L = N2+ 1 8— N

1

132 WRITE (6s>808)< TITSHI IS B K)#K=1*3)«ITS(IS),CSS

DO 150 I =N 1 « N2

FFF =FF ( I ) /O • 4719

PII I = P( I)/248.8

NN= NC ( I

)

I F ( NN » GT • 0 ) GC TO 135

WRITE ( 6*609) I FLOOR! I ) *°I 1 I * FFF

GO TO 141

135 DO 140 J=1 iNN

FCCC=FC( Is J)/0.4719

J J = JC ( I e J ) DP=(P( JJ )-P( I ) )/248.3

CC=C2 C I« J)

I F ( DP .LT, 0 • ) CC=C 1 ( I * J

)

CC=CC*2. 1 19/< 0. 00401 9**X( I, J)

)

I FJ= I FLOOR ( J J

)

COM= J J— SNCGMP ( I F J

)

I F ( J .GT. 1 ) GO TO 136

W D ITE(6»61 0) IFL 00 R ( I ) » P 1 1 I » FFF , I F J , COM, DP , CC , X ( I , 1 ) , FCCC

GO TO 140

136 WRI TE(6c61 1 ) IFJcCOMo DP c CC « X(I*J),FCCC

140 CONTINUE

1 41 NNC = NCO( I ) I F ( NNO .EQ. 0 ) GO TO 150

DO 146 J= 1 , NNO

F 00 =F0 ( I, J )/0.4719

J J= JOC ( Is J ) DP=(PO(I*J )—P ( I ) ) /248.

8

CC = C02 ( Io J

)

I F ( DP .LT. 0. ) CC=C01 ( I , J)

CC = CC*2. 1 19/( 0. 004019**XO( I • J) ) 146 WPITE(6»612)JJ»DP»CC»X0(I »J) ,F00

150 CONTINUE

WRITE(6« 813) FF

I

160 CONTINUE < C

C SUMMARY OUTPUT

C USER INSERTS WRITE STATEMENTS TO FILE I OUT

C

165 CONTINUE

RETURN

C

c

C FORMAT STATEMENTS

D30
* S7URHI 3yk§tt§d?‘*filE OUT

C

60 1 602

603

604

60 5

606

609

6 10

611

612

800

FORMAT !/4X, I3,I10.F13.3,I8.F12,0) STURHI .A W0H8 .$ P I JR*8 GUG

FOR MAT ( /4X * I3,110,F13.3,I8,F12.0,3X ,3 A6 , i F ?8iI)VP f 8^ £ »F 1 1 • 1 ) FORM AT { 53X, 3 A6 , FI 6. 3, FI 5.0, FI 0.2, FI 1. 1

FC5GH»?3/3*«f 3#f fe&M^ARTMENT* 13

i fi i.3 a fT6tio9?i0J€9Pi¥Iiy obriudbr si hur rbhtoha 31 ro

FOR MAT ( 53X» 5HFL00R « 13* 12H COMPARTMENTS 3 ?^**!' 1 ?§ .Pfd.tyjP 10 ,2 »F1 ? .1 ) F OR MAT ( 53 X » 1 7 H OUTS I DE D IRECT ION • 13 , FI 4 . 3# F ?3 .80^1 6 • 2# 1 • 1 ) FORMAT < 4X® I 3» F 1 0. 3 *F 1 1 • 0 ) FORMAT! 4X, I 3 . F 1 0 .3 , F 1 1 • 0 , 3X . 5HFL00R , 1 3 • 1 aH^Otf^RfMEN?'? 13 ,F11 .3,

1 F15.0.F10.2.FU.1 ) YJMQ 2THI0R S H33WT38 3 ROB 4 = D»

F0RMAT(31Xt 5HFLOOR *13*1 2H COMPARTMENT, I3,F11.3.F15.0,F10.2,Fll2l)

fop mt r^5MPl?Xg§*2 8x«

* IGR* •43HMi*RH*eWt¥A49SH*«Mr*QawtemiWe^4R¥9MNL§0Rl2H»¥tfi>«gaM*4RTMENT

2 ,2x.8HP*eBfcdfte«2Ritf?APMcMettua«4fti*&e9f s, To ,

. <c*23*eRT«BPAt9!Swms;9*#itHcHB^iWei^T?aXV§fl£xP0NEfl¥#‘SxUMi>&0W S') 8 0pTM^dR»i*tt^4Ma$3?¥l<jSFi3Sf^lfeVP¥29^>T, CMM >fl00J3I. ( J3M >H,HH 4

80 2 FORMAT !/4X» 13, 110, F 13.1, 1 8 » F 1 2, <6 V§* 9 i ?<?£ i 2 SF 1 1 .1 ) F0RMAT(53X,3A6,F16.1 ,F15.1 , FI 0 .2 * I Of

q H0ISI03RR 3J0UOO

FOi*Ma¥IS»¥«?6*3# ^9M^»6tment # I3,

1 FI 1. Is F15.1* F10.2«F11 . 1 ) HURI\HUR\ HDMMOD

FORMAT! 53X C 5HFLOOR, I 3, 12H COMPARTMENT, I 3, FI 1 . 1 tPVS i Pf b .2 ,F1 1 .1 ) FORMAT (53X* 17H0UTSI0E DIPECTI0N,I3,F14.1,F15.1, FI 0 . ¥S i 1

FORMAT! 1 15X,F8,1 .4H NET) MHW( OOY «5 ) 0 A3

FOFMAT(////20X, 3A6//20X ,2 OHTEMPERATURE !£? • T v> i 3^2 d H i

1 238-SHAFT FLOW COEFF I Cl 0N?I • 1 £ f 9 *?8W AB.il&SfgOi* £3*^ SNF I XED , 2 2 8 X ® 1 2 HD I FFERENT I AL » 5 X* 4 HFLOW, 8Xe 4HFLQW/4 WSHpEdSR&g*, 9GPRESSURE,

3 5X« 4 HFLOW «3X« 1 6HCCNNECT I ON TO , 1 2X , 88Pfc£s3tf fcEQSx S i ?At0 EFF I C I EN T 4, 2X ,8HEXP0NENT,5X,4HPQ§w0?)0D! 4 .TO, OX .RO . I .TJ. 0X)3I

FORMAT !4X, 13, F10, 1

,

FI 1 .0) 3=88

FORMAT! 4X S I 3, FI 0, 1 , F 1 1 • 0 • 3X , 5HFL00R , f> 3

,

91H°§<*>tfp a6¥MElQ¥ ) ?3 , FI 1.1,

1 FIE. 1 .F10.2.F1 1 .1 ) °S OT 00 ! S .03. OX)BI

FORMAT !31 Xs 5HFLC0R* 13 s 12H C OMPA RTME Ntf I SjrF??^ fc F?§ • 1 id J 2, F 1 1.1)

F OP MAT !31X® 1 7H0LT SIDE DIRECTION s I 3 , P? 4S?f , P? 4$ 1 tPl © • 2$ I* I i • 1 ) FOR MAT! 93X„F8. 1 ,4H NET) ROTS

FORMAT! 1 HI ) END 21 T3AHS HI 33 HI 30HAHD

803

804

805

8 0 6

807

808

909

810

81 1 812

813

814

SHQGsP * SUBROUTINE IN°UT £2 I ) Igfclf!

(SI )£3H=SH

P IY4. 0*88=08 ( S .03. TIHUI)3I

SH , I H= I ei 00

88= ( I ) 33

3UHITH00

E ,3 I ( 008 ,8 )3T IRW

001 OT 00

D

D

0

0 I

8 I

D SI T3AH3 R03 0 HI 30HAH0 D

D (SI)ISH=IH OS

D31
n

r>

*-*

*-• no

o

* SUBROUTINE INPUT2

£>N3S*PL IB*. SHOW A.INPUT2 « 0 • S I B , 8 I * F . 0 I -I , Z> i I » "

< I . I : SUBROUTINE Xygyi^.dAE. XE» 0 .£ 13 , B I • C *E 13 *0 1 I* £ I *XAN ) TAM

c (X .1 13, S. 013,0.6 I3,0.dI3.dAE,X06> TAMPO^ EOd

C THIS ROUTINE DETERMINES IF THE PROGRAM • IS TO EE STOOGED

C OR IF ANOTHER RUN IS REQUIRED WITH MODIFIED DATA

C ilWW. £9*1* if ifN'MJf e ?th 3MTPAPM00 HS f *EI . POO JBH6 *XE6 ) TAM < c r * I «E. A 13 ,01 .MOITDBPIC BO I 3TU0HT I * XE8) TAMPQ3

C KC = 2 FOR C ( 0 • 1 1 3 » 0 • 0 I 3 o F I »

• '. ,€.1X3, 0I$fH5M?p£9Bo§oHSI .0 I* P00J3HB . XO, 0.113,0.0X3, 01 »XA ) T A M

C KC = 4 FOR C BETWEEN 2 POINTS ONLY ( I . I IB ,S* 0!

<lClI3,S.0I3,0.6I3,0.1I3.0I .TMBMTBAPMOD HS 1 ,0 I , BOOJ3H6 * X I E > T AM 303 X Id

PARAMETER MM=220 , MS=3f».*IC=* ;.MPE)=a» MTReSsMfrfeR* ©S«M8b 106

COMMON NT* P(MM| x §fMMfM^,NE<;MMi#UG(MM*MC>f fTS(MS)*

1 FC(MMoMC>* PZ<MM# MC*»P9(;MM»MPOJf§9ftMM«MRQiAf < MM >* PFO< MFL* MPO) • 2 FF (MM)oFO<MM,MPOI f CS(iMS) ,R§ ( MRMa NSli t MS* t MSB <MS» * 3 FSS<MS)»NoNS *NPO* ICON V* E * I BUG > X ( MM «MC ) *X 0( MM,MP0> .TITSHCMS.3)

.

A NH * H( MFL ) . IFLOOR ( MM ) * T < f f MSk I W?8? t MS ) IT ( M B > . NTR 0 8

5 * NCG< MM), JOCCMM* MPP.£t*eM¥.S 13,81 , I „EIB>OX I :

’ I v DOUBLE PRECISION P « ?0 f»?$ , S . 0 13 , I. 3X3, I

d 13 • dAE# XEB> TAMPO

B

COMMON /IO/THiyE(J ta^»XeW¥»l9NI¥#NC0MR€MFW) *SNCOMP(MFL)

COMMON /RUN/I RUN (IeII3»S»OI3*I.6I3oX.I IB

X. X XB, s. 89M9?169Ntl . I ib,oi , tmbmtpapmod HS f ,0 I ,RQOJ3HS «X05 itampo

I lf¥^T§*OIB ,I.5I3,i.AI3,0I,M0ITD3BIG HCI6TU0HYI «XE6)TAM 303

5 READC 5*700 )NNN ( T3H H A » X • 8B » X 6 X X ) T A lRtHN«iE«9T *B9!850R 3PUTAP3qM3TM0S, XOSWdAO , XOSWW ) TAM BOB

WRI TE<6, 600) (TITLE! »lBl *12) • I RUN D IBB BOD WCJB / -i

“

IU8 2BPBR8 **§Qb$&3%WB9Wa\W0J3HA »X8 « WO JBH A»X3#JAI TH3B333 I CHS I B X 8<

T MB

I

D IBBBPF AC (5c7OO)KC»IS*08,XSX,OT HOI TDBMHOOHd X *XE«W0J3HA » X? 0

I F ( KC »LT • 1 • OR* KC • GT • 4 ) GO ( ^OwgftOlHA , XE , TM3H0PX3H8 , XS » * BB=e (0* X XB» X .0X3 „0I ,XA ) Tt

, X . X X^tfcft ^ft3MfflARfcd§0Hl9,l;Q. ROOJ3HE,XO.O • 113,1.013,01 ,XA ) TAMPO B

I F{ KC .EQ. 2 ) GO TO 20 ( X . i X B , S • 0 XB , X . 3 X B X

X • X X3,sl6f*tr«§9R.n§93't@l39MBMTPAqMOO HS X *01 »ROO JBHE 8 X IE) TAMPO

• If I 9 HOI T 03 3 1 Q BO I 2 T JOHY I a X IE) TAM 3

STOP (T3H HA, X

83 * XO 9 ) TA M 3Q3

( XHI )

T

AMPOB AIR

CHANGE IN FF IN SHAFT IS

0 N1=NS1(IS) 3HITU0P3US

N2=NS2( IS)

I F ( I UN I T .EQ. 2 ) BB=BB* 0.4719

DO 15 I=N1,N2

FF ( I )=BB

5 CONTINUE

WRITE! 6c 800) I S,

B

GO TO 100

CHANGE IN C FOR SHAFT IS

C

20 N 1 =NS 1 ( IS

)

N2=NS2< IS)

DO 28 I=N1,N2

NN= NC ( I ) I F ( NN .EQ. 0 ) GO TO 28

DO 24 J=1,NN

IFCIUNIT .EQ. 2 ) BB = B*0.47 19/( 2A8.8**X( I • 1 ) ) C ( I c J )=BB

D32
* SUBROUTINE INPUT2

JJ=JC«I*J)

NM=NC( JJ)

30 22 K = 1 e NM

IF< JC< JJ.K ) .EQ. I ) GO TO 23

22 CONTINUE

23 CIJJ*K)=8B

24 CONTINUE

26 CONTINUE

WRITE<6* 80i)IS,B

GO TO 100

c

C CHANGE CO IN BUILDING

C

30 IF I IS • GT * 0 ) GO TO 40

30 38 1=1.

NNO=NCO ( I

I F < NNG • EQ. 0 ) GO TO 38

00 34 J= 1 . NNO

IF«IUNIT ,EQ. 2)BB=B*0 .4719/(248.8* *XO<I»J))

COd.J ) = BB

34 CONTINUE

38 CONTINUE

WRITE(6 e 802)B

GO TO 100

C

C CHANGE CO IN SHAFT IS

C

40 N 1 = N S 1 ( IS)

N2=NS2( IS)

DO 48 I=N 1 • N2

NNO =NCO ( I )

1 F ( NNO .EQ. O ) GO TO 48

00 44 J = 1 » NNO

IF* IUNIT .EQ. 2 )BB=B*0.4719/(248.8**X0< I, J))

C0{ Ic J)=BB

44 CONTINUE

48 CONTINUE

WP I TE ( 6 ® 8 03 ) IS.B

GO TO 100

C

C NEW C BETWEEN POINTS IS AND JS

C

50 PFAD(5»700)JS

NN=NC* IS)

DO 51 J = 1 e NN

IF* JC ( IS. J ) .EQ. J S ) GO TO 52

51 CONTINUE

GO TO 200

52 IF( IUNIT .EQ. 2 ) BB=B*0 • 47 19/< 248 • 8* *X( I S . J )

)

C( I Sc J )=BB

NN=NC( JS)

DO 53 J= 1 c NN

I F ( JC ( JSc J ) .EQ. IS) GO TO 54

53 CONTINUE

GO TO 200

54 C(JS«J)=BB

WRITE < 6c 805) IS, JS,

100 CONTINUE

D33
* SUBROUTINE INPUT2

200

600

700

800

801

802

803

804

805

RETURN

WR I TE (6# 804

)

STOP

FORMAT! ///I OX. 12A6//10X. 3 HR UN. I 4///

)

F OF MAT ( ) FORMAT! 10X.25HNEW

FORMAT! 1 OX, 16HNEW

FORMAT! 10X.23HNEW

F OF MAT! 10X. 17HNEW

FIXED FLOW FOR SHAFT ,I2,4H

C FOR SHAFT *12, 4H IS.F8.1//)

CO FOR BUILDING IS »F8*1//)

CO FOR SHAFT .I2.4H IS.F8.1//

F0RMAT!//10X. 2 1 HERR OR IN INPUT 2 DATA /)

FORMAT! 10X. 14HNEW C BTTWEEN ,I3.4H AND.I3.3H IS

END

3) HD Go ° * SUBROUTINE UNITS

IS.FS. 1/)

),F8.1//>

• L , 1

D34
* SUBROUTINE UNITS

S)NBS*PL IBS. SHOW A. UNITS

SUBROUTINE UNITS

C

C THIS ROUTINE CONVERTS VARIABLES H,FF,C,CO,CS TO SI UNITS

C

c PARAMETER MM = 220 . MS=2 « MC=4 , MP0=2 • MTP = 2 • M FL=1 0 5 . MB= 1 0

COMMON NT, P( MM) ,C( MM,MC) ,NC( MM) . JC( MM. MC) , ITS( MS ),

1 FC< MM e MC) e PZ(MM» MC ) , PO ( MM» MPO ) * CO ( MM. MPO) *F( MM). PFO( MFL. MPO)

2 FF (MM) ,FO< MM, MPO) ,CS( MS) , PS ( MFL ) ,NS1 ( MS) ,NS2(MS)

,

3 FSS(MS). N. NS ,NPO, I CONV , E , I BUG , X ( MM , MC ) , X O ( MM , M PO ) , T I T SH ( MS , 3 ) . A NH *H ( MFL ) , I FLOOR (MM) . T ( MTP , MFL ) . NF SI (MS) ,NFS2(MS), IT(MB).NTP

5 . NCO(MM ). JOC( MM, MPG) • TOUT

DOUELE PRECISION P,PO,PS

DIMENSION B ( 4

DATA B/0* 30 48 ,248. 8,0.4719, 0. 02992/

DO 10 1 = 1. NH

10 H( I )=H( I )*B( 1

DO 20 1=1, NT

FF ( I )=FF( I )*B(3)

DO 16 J= 1 , MC

IF(X(I e J) .EQ. 0 ) GO TO 16

BB=E (3 )/(B(2 )**X( I, J) ) C( I * J ) =C ( I, J)*BB

16 CONTINUE

DO 18 J= 1 « MPO

IF(XO(I®J) .EQ. 0 ) GO TO 18

BB=B(3)/(3<2)**X0( I , J)

CC( I g J ) =C0 ( I , J )*BB

18 CONTINUE

20 CCNTINUE

DO 22 I S= 1 . NS

22 CS( IS) = CS(IS)*B(4)

RETURN

END

£HDGe =» * SUBROUTINE SIMEQ L, 1
* SUBROUTINE SIMEQ

©NB S*PL I B $ • SHOW A* SI MEQ

SUBROUTINE SIMEQ

C

C CHOLESKY* S METHOD OF SOLUTION OF

C SIMULTANEOUS LINEAR ALGEBRIC EQUATIONS

C

PARAMETER MM=220, MS=2* MC=4* MPO=2 • MTP = 2 ,MFL=1 05 • M8=l 05

PARAMETER MBP=MB+

1

DOUBLE PRECISION A*X

COMMON /MAT/ A< MB • MBP ) , X< MB ) «

N

NP1=N+1

ZEP.G= 1 • 0E—35

K = Q

C

c

C SEE IF A ( 1 • 1

)

IS ZERO

C IF SO ADD ANOTHER ROW TO ROW 1 IF( ABS( A( 1 1 )

)

.GT. ZERO) GO TO 40

DO 31 1=1,

N

1 F( A( I«

1

) • NE • 0* ) GO TO 32

31 CONTINUE

12 WRITE( 6, 804)K

STOP

32 DO 33 J= 1 « NP

1

33 A ( 1 o J )=A< 1, J ) +A ( I, J

)

C

C CALCULATE UPPER AND LOWER

C TRIANGULAR MATRICES OVER ORIG

C MATRIX A 40 A A= A ( 1 c 1

)

DO 2 J=2, NP1

2 M1«J)=A(1,J)/AA

DO 10 1 = 2* N

K = 0

C

C STORE A ( I o

1

) ••• A (

I

, I

)

IN X ARRAY

C IN CASE NEW A(I,I) IS ZERO

C ROW I CAN BE RECALCULATED

4 DO F J = 1 * I 5 X ( J ) = A ( I * J ) K = K+ 1 DO 10 J=2*NP1

I F( J • GT • I ) GO TO 8

J Ml = J1

A 4=0.

DO 3 IR= 1* JM1

3 AA=AA + A( I* IR) * A ( I R , J ) A ( I * J ) = A ( I, J )-AA

C

C CHECK IF A ( I , I

)

IS ZERO

C IF SO MULTIPLY OLD ROW I BY 2*

C

I F < I .NE, J ) GO TO 1

0

IF( ABS(A(II) ) .GT. ZERO ) GO TO 10

DO 6 JJ=1,I

6 A ( I * J J ) =X ( J J

)

DO 7 J J= 1 * NP

1

D36
* SUBROUTINE SIMEQ

Cc

c

cr

c

7 M I o J )=2.*A( I, J

I F ( K • GT • 3 ) GO TO 12

GO TO 4

8 I M 1 = I— 1 A A — 0 * DO 9 I R= 1 * 1 M

9 AA = AA + AU« XR)*AC IR» J)

A (I# J)=( A(I ,J)-AA)/A(1 , I)

10 CONTINUE

END OF CALCULATION OF TRIANGULAR MATRICES

BACKWARD SUBSTITUTION

X(N ) = A(N*NP1

DO 20 I 1=2 •

A A=0«

I =NP1-I

IP1=I+1

DO 15 J=IPl,N

15 AA = AA + A< I® J)*X( -J) 20 X ( I ) = A( I • NP1 ) — A A 804 ^ OF. MAT( //////l OX » 1 6HPR0GR AM FAILURE ,13/////)

END

fO HD Go P * FUNCTION FLOW «L,1

D37
uu

* FUNCTION FLOW

5>NB S*°L I B $• S HOW A. FLOW

FUNCTION FLOWCPI ,PJ,PZ.C,X)

OOUELE PRECISION PI.PJ

THIS FUNCTION CALCULATES FLOWS BETWEEN TWO POINTS

C

I F ( ABS(C) »LT • OtOODGO TO 10

C>o=p j_pi +pz

S I GN= 1*0

I F ( OP .LT. .0)SIGN=-1.

I F (

X

.NE. *5 ) GO TO 11

FLOW=SIGN*C*SQRT< SIGN* DP)

RETURN

10 FLOU=0.0

RETURN

11 FLOW=SIGN*C*< < SIGN*DP)**X>

RETURN

END

SHDG.P * FUNCTION PFLOW • Lei

D38
ro

FUNCTION PFLOW

3)NBS*PL I 8$ • SHOW A. PFLOW

FUNCTION PFLOW ( I * P I

cc

C THIS FUNCTION CALCULATES NET FLOWS INTO POINT I C

PARAMETER MM=220 . MS=2. MC=4 ,MP0=2,MTP=2,MFL=105 ,MB=105

COMMON /CORR/C1 (MM.MC) ,C2 (MM * MC ) • COl ( MM* MPO) , C02 ( MM, MPO

)

COMMON NT o P ( MM ) C ( MM • MC) • NC( MM ) • JC( MM* MC ) ITS ( MS )

1 FC(MMoMC) * PZ ( MM * MC) »PO(MM»MPO) *CO( MM, MPO) ,F ( MM ) , PFO( MF L * MPO) , 2 FFCMM) ,FO( MM.MPO) ,CS( MS) .PS(MFL) ,NS1 ( MS) ,NS2(MS)

,

3 FSS ( MS ) * N» NS ,NPO, ICONV, E » I BUG X ( MM , MC ) ,X 0< MM , M PO) , TI TSH< MS , 3 )

,

A NH , H ( MFL ) , IFLOOR(MM) , T ( MTP , MFL ) , NFS1 (MS) , NFS2(MS), IT(MB),NTP

5 e NCO< MM) , JOCC MM, MPO) ,TOUT

DOUBLE PRECISION P,PO,PS.PI

NN=NC( I

SUM=Q

IF(NN «EQ. 0 ) GO TO 3

DO 1 JJ=1,NN

J=JC< I, JJ

CC= C 1 ( I , J J

IF(PI ,LT, P ( J ) ) CC=C2( I » J J

PZZ=PZ( I , J J)

I F ( I • GT • N ) PZZ=0

F C ( I , J J ) =FLO W (P I « P ( J ) »PZZ»CC»X( I,JJ )

1 SUM=SUM+FC( I, J J

NNC=NCO( I

I F ( NNO • EQ o 0 ) GO TO 4

DO 2 K= 1 , NNO

CC=C01 ( I,K

)

I F ( P I .LT. PO (I , K ) ) CC=C 02 ( I * K

F0( IcK ) = FLOW(P I ,P0 U ,K ) ,0 ,CC, X0( I ,K)

SUM = SUM+ FO( I * K)

PFLOW=SUM+FF( I

)

I F ( I ,LE, N)F( I )=SUM+FF( I ) RETURN

END

BHDG,P * FUNCTION PSFLOW ,L,1

D39
* FUNCTION PSFLOW

a)N9S4'PL I B $• SHOW A.PSFLOW

FUNCTION P SFLO W ( I » P I ) C

c

C THIS FUNCTION CALCULATES NET HORIZONTAL FLOWS

C INTO A FLOOR OF A SHAFT

C

PARAMETER MM=220 , MS=2 . MC=4 ,MP0=2, MTP=2 . MFL=1 05 .MB=105

COMMON /C0RR/C1 (MM , MC) * C2 (MM »MC ) • C01 C MM* MPO) • C02( MM, MPO

)

COMMON NT, P(MM)*C( MM* MC) * NC (MM), JC( MM* MC), ITS(MS),

1 FC(MM,MC) *PZ ( MM,MC) ,PO( MM, MPO) ,CO( MM, MPO) , F ( MM ) , PFO( MFL, MPO )

,

2 FF(MM) e FO ( MM * MPO ) ,CS( MS) *PS(MFL)*NS1 (MS) ,NS2(MS)

3 FSS(MS) *N*NS,NPO. I CON V , E , I BUG , X ( MM • MC ) • X 0 ( MM , M PO ) • TI T SH( MS ,3 )

,

A NH ,H ( MFL ) , IFLOOR ( MM ) ,T ( MTP , MFL ) . NFS1 (MS) ,NFS2(MS), IT(MB)*NTP

5 cNCO(MM)* JOC ( MM* MPO

)

DOUELE PRECISION P,PO,PS,PI

NN=NC( I

)

SUM=0*

I F ( NN *EQ* 0 ) GO TO 3

DO 1 J J= 1 * NN

J = JC( Ie JJ)

CC = C1 ( I* JJ

)

I F ( P I ,LT, P ( J ) ) CC =C2( I »JJ)

1 SUM=SUM+FLOW(PI ,P ( J ) ,0 , CC ,X ( I • J J )

)

3 NNC =NCC ( I

)

I F ( NNO • EQ* 0 ) GO TO A DO 2 K= 1 , NNO

CC=C01 ( I,K

)

I F ( P I ,LT* P0( I ,K) )CC=C02(I ,K)

2 SUM=SUM+FLOW(PI ,P0 ( I , K ) , O , CC , XO ( I • K )

)

A P SFLO W=SUM+FF ( I

)

RETURN

END

5>PRKPT PRINTS

D40
FEDERAL INFORMATION PROCESSING STANDARD SOFTWARE SUMMARY

01. Summary date

Yr. Mo. Day

04. Software date

Yr. Mo. Day

3 ) 0 5 2 7

02. Summary prepared by (Name and Phone)

John Klote 921-3387

05. Software title A Computer Program for Analysis of Pressurized

Stairwells and Pressurized Elevator Shafts

03. Summary action New Replacement

o

Previous Internal Software ID

07. Internal Software ID

Deletion

06. Short title 08. Software type

Automated Data System

[3] Computer Program

Subroutine/Module

09. Processing mode

Interactive

:jT] Batch

| | Combination

10. General Computer Systems

Support/Utility [~

Scientific/Engineering Q

Bibliographic/Textual Q

Application are

Management/

Business

Process Control Other

a

Specific Smoke control in

buildings

11. Submitting organization and address Fire Protection Systems Research

Center for Fire Research

National Bureau of Standards

Washington, D.C. 20234

12. Technical contact(s) and phone

John Klote

921-3387

13. Narrative Pressurized stairwells and pressurized elevators can be used as a means of providing a smoke free exit route during fire situations. This computer program analyzes systems

intended to pressurize stairwells and elevator shafts.

k

4ir

Ke$ovement ; computer programs; egress; elevator shafts; escape means; pressurization;

smoke control; stairwells

15. Computer manuf'r and model 16. Computer operating system 17. Programing language(s) 18. Number of source program state- ments Univac 1108 Fortran V 1900

19. Computer memory requirements 20. Tape drives 21. Disk/Drum units 22. Terminals

23. Other operational requirements

24. Software availability 25. Documentation availability

Available Limited In-house only Available Inadequate In-house only t] ©

26. FOR SUBMITTING ORGANIZATION USE

185-101 Standard Form 185

1974 July

U.S. Dept, of Commerce—NBS

(FIPS. Pub. 30)
NBS-114A (REV. 9-76)

4. TITLE AND SUBTITLE

1 COMPUTER PROGRAM FOR ANALYSIS OF PRESSURIZED

! TAIRWELLS AND PRESSURIZED ELEVATOR SHAFTS

U.S. DEPT. OF COMM. BIBLIOGRAPHIC DATA

SHEET

1. PUBLICATION OR REPORT NO. NBSIR 80-2157

7. AUTHOR(S)

John H. Klote

5. Publication Date January 1981

8. Performing Organ. Report No.

9. PERFORMING ORGANIZATION NAME AND ADDRESS

NATIONAL BUREAU OF STANDARDS

DEPARTMENT OF COMMERCE

WASHINGTON, DC 20234

11. Contract/Grant No.

12. SPONSORING ORGANIZATION NAME AND COMPLETE ADDRESS (Street, City, state, ZIP) 13. Type of Report & Period Covered

15. SUPPLEMENTARY NOTES

I | Document describes a computer program; SF-185, FIPS Software Summary, is attached.

16. ABSTRACT (A 200-word or tees factual summary of moat significant information. If document includes a significant bibliography or

literature survey, mention it here,) Pressurized stairwells and pressurized elevators can be used as a means of

providing a smoke free exit route during fire situations. This paper describes

a computer program which analyzes systems intended to pressurize stairwells or elevator shafts. The basic assumptions and limitations of the program are also

discussed. The appendices contain a program listing and examples.

17. KEY WORDS (six to twelve entries; alphabetical order; capitalize only the first tetter of the first key word unless a proper name;

separated by semicolons)

Air movement; computer programs; egress; elevator shafts; escape means;

pressurization; smoke control; stairwells.

18. AVAILABILITY §Unlimited

| | For Official Distribution. Do Not Release to NTIS

| | Order From Sup. of Doc., U.S. Government Printing Office, Washington, DC

20402, SD Stock No. SN003-003-

[X] Order From National Technical Information Service (NTIS), Springfield, VA. 22161

19. SECURITY CLASS

(THIS REPORT)

UNCLASSIFIED

20. SECURITY CLASS

(THIS PAGE)

UNCLASSIFIED

21. NO. OF

PRINTED PAGES

85

22. Price

$ 8.00

USCOMM-DC