A11100 ^5631

NBS

PUBLICATIONS

NAT'L INST. OF STAND & TECH

I III III II II II Hill III II HI A 111 Ob M 0 7 4 T

NBSIR 82-2512

A Computer Program for Analysis

of Smoke Control Systems

U.S. DEPARTMENT OF COMMERCE

National Bureau of Standards

National Engineering Laboratory

Center for Fire Research

Washington, DC 20234

June 1 982

Final Report

Sponsored in part by:

•QC—

100

, U56

02-2512

2 partment ot Health and Human Services

ashington, DC 20201

NATIONAL MVKHAl

OW 0TAKDAAU8

UBXAIT

JUl 7 1982

NBSIR 82-2512 Ato-L A COMPUTER PROGRAM FOR ANALYSIS

OF SMOKE CONTROL SYSTEMS

John H. Klote

U.S. DEPARTMENT OF COMMERCE

National Bureau of Standards

National Engineering Laboratory

Center for Fire Research

Washington, DC 20234

June 1982

Final Report

Sponsored in part by:

Department of Health and Human Services

Washington, DC 20201

U.S. DEPARTMENT OF COMMERCE, Malcolm Baldrige, Secretary

NATIONAL BUREAU OF STANDARDS. Ernest Ambler, Director
' ' • - A*' ; . -A " ' ' * •

* v:.

•• S' ' •<

•' < i

f

:U

•:c ^

.«

: .

1 ) •

•' 1 :

’

> „

• : 5

... ,aO

• : ’ i

•

*

‘ Wh V j'sO’ '

"• f'«v .n. •

: •• \

•• n,.- V /
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 con- sists 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

LIST OF FIGURE V

Abstract 1

1. INTRODUCTION 1

2. PROGRAM CONCEPT ...... 2

3. ASSUMPTIONS AND LIMITATIONS 3

4. EQUATIONS 3

5. PROGRAM DESCRIPTION 6

5 . 1 Main Program 7 5 . 2 INPUT Subroutine 7 5.3 CORR Subroutine 7 5.4 INIT Subroutine 7 5.5 BLDGP Subroutine 9 5.6 SHAFTP Subroutine 9 5.7 PZAD Subroutine 9

5 . 8 OUT Subroutine . 9

6. FUTURE DIRECTION 9

7. REFERENCES 10

APPENDIX A. DATA INPUT DESCRIPTION FOR COMPUTER PROGRAM 11

APPENDIX B. INPUT EXAMPLES . 18

APPENDIX C. EXAMPLE OUTPUT 21

APPENDIX D. PROGRAM LISTING 26

IV
LIST OF FIGURE

Figure 1. Flow chart for main program

A COMPUTER PROGRAM FOR ANALYSIS OF SMOKE CONTROL SYSTEMS

John H. Klote

Abstract

This paper describes a computer program developed to

analyze systems intended to control smoke in building

fires. These systems include pressurized stairwells,

pressurized elevator shafts, zone smoke control systems,

and pressurized corridors. This program calculates air

flows and differential pressures throughout a building

in which a smoke control system is operating. 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

The majority of fire fatalities result from smoke inhalation. As a result of

this, a number of systems have been designed and built to control smoke movement in

building fires. The most common smoke control systems are pressurized stairwells

and zone smoke control systems

1

. These systems are intended to control smoke movement in a building by use of air flows and by differential pressures. The computer

program described in this paper provides a means to calculate the air flows and

differential pressures throughout a building (either real or conceptual) in which a smoke control system is operating.

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

2 throughout a building [1,2] . Other programs go beyond this to calculate smoke

concentrations throughout a building that would be produced in the event of a fire

[3-7]. In general, most of these programs are capable of analyzing smoke control

1 The concept of extending the use of smoke control to protect elevators is currently

being investigated at NBS.

2Numbers in brackets refer to the literature references listed at the end of this

paper.

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

for analysis of smoke control systems and is an extension of a program specifically

written for analysis of pressurized stairwells and elevators [8]

.

While the basic

theory of this program is the same as that of the stairwell program it has been

extended to include analysis of (1) stairwells with vestibules, (2) elevators with

elevator lobbies, (3) zone smoke control systems, and (4) pressurized corridors.

The data input has 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

pressure differences across all of the building shafts, as well as the flows and

pressures throughout the building.

This program was originally intended primarily as a research tool to investigate

the feasibility of specific smoke control systems and to determine the interaction

between these systems and the rest of the building. The predecessor [8] of this

program has already been used to analyze pressurized stairwells without vestibules

and to evaluate factors which affect the performance of these systems [9]

.

And,

this program has been used to generate data for an National Bureau of Standards

(NBS) Handbook on Smoke Control Design which is being developed. This paper is not

intended to be a design guide for smoke control systems. The state-of-the-art of

these systems is still under development 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 pressure difference 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 pressurization of any other building space.

In addition, any building space can be exhausted. This allows analysis of zone smoke control systems where the fire zone is exhausted and other zones are pressurized.

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 to 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

3 pressure (101325 P ) . 3

4

.

EQUATIONS

A. Flow equation

F = CA /2pAP (3.1)

where

F = mass flow rate

C = flow coefficient

A = flow area

p = density of air in flow path

AP = pressure difference across flow path

The flow coefficient is dimensionless and for smoke control analysis it is

generally taken to be in the range of 0.6 to 0.7. Because of the large number of

flow calculations performed during the computer analysis the flow equation is

rewritten in the program as F = C' /AP. Using the ideal gas law, the adjusted flow

coefficient, C', can be expressed as

C' CA J

2 P

atm~

1 RT

where:

P atm

= ab s °lute barometric pressure at ground level

R = gas constant of air

T = absolute temperature of air in flow path

(3.2)

3 The results of the program are not 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

and one in the subroutine CORR.

3
B . Mass Balance Equations

For building compartment i

S F

(i,j)

+ 2

j=l k=l

F + F o(i,k) *f(i) 0 ( 3 . 3 )

and for shafts

where

[i, j) + Fo(i, k) + Ff(i)

k=l

0 ( 3 . 4 )

F

(i i) = mass flow ra te from space j to space i. For building

'1 compartments this flow can be either horizontal or vertical,

however for shafts this flow can only be horizontal.

F

. i k) = mass flow rate from direction k outside of the building

' to space i.

F

f

,., = net mass flow rate of air due to the air handling system

1 or due to a pressurization system.

N

c = number of building spaces connected to space i.

N

q

= number of connections to the outside from space i.

is the space number at the bottom level of the 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 relationship is used to calculate the gage pressure, P^, at

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

P... = P ,. . > -P

(i) (l-l) z P

f

where

P

z = hydrostatic pressure difference

P^ = pressure loss due to friction

( 3 . 5 )

The following equation is used to calculate the hydrostatic pressure difference.

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

than in a shaft.

4
p

RT

h

(i)

h

( i-1) 3.6

where

= height of point i

h^_jj = height of point i-1

g = gravitational acceleration

R = gas constant

T =

T

(i) + T(i-1)

P , . . + P , . , . P =

(l)

- X ^ --- + P. 2. D

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,

f. \ 2 Im a P

f = s lc

H

l (3.7)

where

= upward flow from i-1 to i in shaft

C = shaft flow coefficient

sS = sign of m

u

D . Outside Pressures

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

following method.

o (i)

P, , . . + C P ,. . h ( l) w v ( l)

(3.8)

where

Po (i)

h ( i)

'v(i)

Cw

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

at ground level

hydrostatic pressure difference between h(i) and ground level

velocity pressure due to the wind at height h(i)

pressure coefficient

Because the outside temperature is constant

P h(i) * p atm exp I - qhU) \

l RTout )

- p. (3.9)

5
where

P atm

= absolute barometric pressure at ground level

TQut = outside absolute temperature

When the outside pressures are calculated by the computer the wind velocities

are assumed to be described by the power law.

where:

V

Q

= wind velocity at height h

Q n = wind exponent

This relationship has been extensively used to describe the boundary-layer

velocity profile of the wind near the surface of the earth. It 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 appropriate 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 = |pvM.

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 z pressure coefficient depends upon building geometry and varies locally over the wall

surface. Numerical values for C and n as well as practical engineering information w

are available from a number of sources [10-13].

5. PROGRAM DESCRIPTION

This program is written in ANSI FORTRAN on the UNIVAC 1100/82 and a program

listing is provided in appendix D. The following is a detailed description of the

main program and the major subroutines.

6
5.1

Main Program

The main program calls the subroutines which read the data, calculate the

adjusted flow coefficients, calculates the initial values of pressures and interatively solves for the pressures according to the logic illustrated in the flow

chart of figure 1.

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 areas

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.

Appendix A contains a detailed description of the data input method.

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 adjusted flow coefficients for all flow paths using eq.

(3.2). Two sets of these 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 [1] . 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

.

7
Figure 1. Flow chart for main program

8
5.5

BLDGP Subroutine

The iterative solution for the building pressures and flows consists of t e 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 com- partment 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 [14] . 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 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 require modification 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 appropriate variables are converted to the engineering system before

output.

6. FUTURE DIRECTION

It is planned to use this computer program in a project at NBS to study the

feasibility of protected elevators as a means of fire escape for handicapped

individuals. Consideration is being given to further development of the program for use as a design tool. Also, a program may be developed for microcomputers which can

be used interactively.

9
7 . REFERENCES

[1] Sander, D. M. , FORTRAN IV Program to Calculate Air Infiltration in Buildings,

National Research Council Canada, DBR Computer Program No. 37, May 1974.

[2] 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.

[3] 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.

[4] 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, pp. 70-75, 1969.

[5] 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, pp. 78-87,

1969.

[6] Evers, E. and Waterhouse, A., A Computer Model for Analyzing Smoke Movement in

Buildings, BRE Fire Research Station, Borehamwood, Hertsfordshire, England,

November 1978.

[7] Wakamatsu, T., Calculation Methods for Predicting Smoke Movement in Building Fires

and Designing Smoke Control Systems, Fire Standards and Safety, ASTM STP 614,

A. F. Robertson, Ed., American Society for Testing and Materials, pp. 168-193,

1977.

[8] Klote, J. H., A Computer Program for Analysis of Pressurized Stairwells and

Pressurized Elevator Shafts, Nat. Bur. Stand. (U.S.), NBSIR 80-2157, January

1981.

[9] Klote, J. H. , Smoke Control by Stairwell Pressurization, Engineering Applications

of Fire Technology, SFPE and Nat. Bur. Stand., April 1980.

[10] Sachs, P., Wind Forces in Engineering, Pergamon Press, New York, 1972.

[11] Houghton, E. L. and Carruther, N. B., Wind Forces on Buildings and Structures,

John Wiley & Sons, New York, 1976.

[12] Simiu, E. and Scanlan, R. H., Wind Effects on Structures: An Introduction to

Wind Engineering, John Wiley & Sons, New York, 1978.

[13] MacDonald, A. J., Wind Loading on Buildings, John Wiley & Sons, New York, 1975.

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

John Wiley & Sons, Inc., New York, 1969.

10
APPENDIX A. DATA INPUT DESCRIPTION FOR

COMPUTER PROGRAM

Data input consists of the following elements:

1. Initial data

2. Building heights

3. Temperature profiles

4. Outside pressure profiles

5. Building data

6. Shaft data

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

Elements 1 through 6 are always required. In the following sections the input

required for each of the six 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

and integers can be written with or without decimal points. 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

project title (col. 1-72)

outside unit indication summary output

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

1The user must assign this file before program execution.

11
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

entered as follows:

h

(l)

h

( 2

)

where h

(i)

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

For each temperature profile the following data must be supplied.

no. of temp. level temperature level temperature

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

level temperature

no. (°C,°F)

12
4 . Outside pressure profiles

no. of outside

pressure profiles

input parameter

(either 0 or 1)

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

Po(l) P

o (2)

Po(3) P o(i) P° ( V

where P ,.

.

o (l)

is the outside pressure at the center of level i.

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

data are required.

V

o wind velocity

(mph)

h

o height at

velocity is

which

measured

n wind

exponent

pressure coefficients for each pressure profile

CW(1) C

W(2) CW(Npo)

o o o o

5 . Building data

N

f no. of levels

(or floors)

13
All of the following data in this input element are supplied for each level, or consecutive groups of similar levels.

I 1 Starting floor

I

2 Ending floor

Ncom

No. of compartments per floor

(Floor data is entered in ascending order of levels or floors. When data

is for only one level then 1^ = and the same number is supplied for

both.

)

For each compartment on a level the following data are supplied:

NCS

No. of connections

to other compartments on the same level

NCA

No. of connections

to compartments on

the level above

NC0

No. of connections

to the outside

Net flow^

(1/s, cfm)

Temperature

profile

number

For each connection between this compartment and another on the same floor the

following data are required.

Other compartment

number on the same level

C

flow coefficient flow

(m 2

,

A

area

ft 2

)

For each connection between this compartment and one on the level above the following

data are required.

2All net flows are at standard conditions of 21°C (70°F) and one atmosphere.

14
Other compartment

number on floor

above

C A

flow coefficient flow area

(m 2

, ft 2

)

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

outside pressure

profile number

C

flow coefficient

A

flow area

(m 2

, ft 2

)

6

.

Shaft data

no. of shafts

All of the following data in this input element are required for each shaft,

shaft title (col 1-20)

shaft flow

coefficient

bottom

level of shaft

top level

of shaft

temperature

profile number

Enter the following typical data which applies to each level of the shaft.

Exceptions can be entered later.

15
no. of connections

between typical

level of shaft and

outside

net flow into

typical level

of shaft

(1/s, cfm)

The connection data to the building for a typical level are required.

compartment no. to which shaft is

connected

C A

flow coefficient flow area

(m 2

, ft 2

)

For each connection to the outside, the connection data for a typical floor are required

.

outside pressure

profile

C

flow coefficient

A

flow

(m 2

, area

ft 2

)

The number of exceptions to the typical data is required.

no. of exceptions

All of the following data in this input element are required for each exception.

exception type level of shaft

(1, 2 or 3)

16
The next card depends on the exception type. For exception type = '

.

to the net flow into the floor of the shaft is defined.

net flow

(1/s, cfm)

For exception type = 2 , an exception to an outside connection for this shaft is

defined

.

outside pressure

profile number

C

flow coefficient

A

flow area

(m 2

, ft 2

)

For exception type = 3, an exception to the connection between the shaft and the

building is defined.

compartment no. to which shaft

is connected

C

flow coefficient

A

flow area

(m2

, ft 2

)

17
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 550 cfm1 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 wind pressure coefficients are 0.7 for the

windward wall and -0.7 for the leaward wall. The flow areas are the same vertically

and are listed in Table Bl. The flow coefficient is taken to be 0.65 for all

connections

.

Table Bl. Flow areas for example 1

Connection location Area (ft 2

)

Between stairwell & building 0.42

Between building & outside into

the wind 0.75

Between building & outside away

from the wind 0.75

1.1 Data for Computer Input

initial ( TEN STORY BUILDING WITHOUT VERTICAL LEAKAGE

data ( -20 20

bui Iding i 10 1 heights \ 5 10

temperature

f 2 profi les 1 70

( 2 1 50

outside

( 2 1 pressure < 30 30 .15

profi les (0.7 -0.7

1At standard conditions of 21°C ( 70 °F) and one atmosphere.

18
1

2 0 1 0.75

0.75

10 2

.42

2. Example 2

This is a 10 story building which is 70°F inside. Outside the air temperature

is -5°F and there is no wind. This building has a stairwell and an elevator. The

flow areas which are generally the same vertically are listed in table B2 and the

flow exponents are taken to be 0.5. The stairwell is pressurized by a net 550 cfm

per floor. The elevator shaft has a 4 ft 2 vent to the outside at the top. On

floors 2 through 10 the elevator lobby separated from the building by doors that

automatically close in the event of a fire. The flow coefficient is taken as 0.65

in all connections.

bui Iding

data

10

1

0

1

2

10

0

.65

.65

shaft

data

STAIRWELL

80000 1

0 550

1 .65

0

Table B2. Flow areas for example 2

Connection location Area (ft 2

)

Between stairwell & building 0.42

Between building & outside 1.5

Between elevator & elevator lobby 0.65

Between elevator lobby & building 0.55

2.1 Data for Computer Input

initial (TEN STORY BUILDING WITH ELEVATOR

data (“5 2 0

bui Iding fiO 1 heights 1 5 10

temperature i 1 profile 1 1 1 70

outside

( 1 1

pressure < o 1 1

data 1 1

19
bui Iding

data {10

shaft

data

1st { 1 1

floor < 0 0 1

( 1 0.65 0.75

2 10 2 2nd 1 0 1 through 2 0.65 0.55

10th floors 1 0.65 0.75

0 0 0

{ 2

0

0

0

1

1

1

shaft l

shaft 2

STAIRWELL

80000 1 10

0 550

1 .65 .42

0

ELEVATOR

2.7E6 1 10

0 0

2 .65 .65

2

2 10

1 .65 4.0

3 1

1 .65 .65

1

1

2 . 2 Example 2 Output

The output

for Cases 2 and

for example

3) is given

2 case 1 (the data above not including modifications

in appendix C.

20
APPENDIX C. EXAMPLE OUTPUT

21
OUTSIDE

TEMPERATURE

00

CNO

o

uJ<

13

Q.

ZDO

UJ

Q-

<r>cUJ

O

UJo

l/D UJO

a

a.

UJa3

6-<a

UJ

CLs oooooooooo

oooooooooo

h- oooooooooo X oooooooooo OH-

•-« u. in in in m in to id in m in

uj *-oj(r)*^incot^a)(j) x

a»^^iooi»-^ioo>riQ: O^NOCONOPHOOUJ QMD’J-WOOJMn^H

oooooooooo

22

THE

FOLLOWING

UNITS

ARE

USED

FOR

OUTPUT
TEN

STORY

BUILDING

WITH

ELEVATOR

ADJUSTED

TEMP

FIXED

DIFFERENTIAL

FLOW

FLOW

FLOOR

COMPARTMENT

PRESSURE

PROFILE

FLOW

CONNECTION

TO

PRESSURE

COEFFICIENT

AREA

FLOW

h~

UJZ

UJz

in n oj t- co co t t n in r- o cn cm co ooo cn co cm co cn o — cj) in co co cn — cn — t*** o m n cv

r- t- « 0(0(0* oo co cm in oo oo o coco T (O — CO T — COCO

00 cn CM 0-0(0

co in —

co oo t o co • T in t t — co — — —

co co co in cm co co o cn co - o —

id t id co oo CM in CM CM CM in'oo T (O CM

cm in co

cr> o o — o o

co m in m cm m

o co r-~ o o- c* - o - o o

in m m cm in oco o t- - o - o oin in uo cm in 0(0 o t Is- — o - o oin m in cm in 0(0 o T t" - o - o cin in in cm in oco o t c"-

•n <3- o

cn cn — O -D - — — CM

in cl t

oo cn in 'no (TOD T (0

in in

co c* in TOO

T T

CO O

T -3 t in lo t <? oo cn in oo cn TOO ^ vO t m in t «r

oo cn m oo cn n O D ST O

t- in in

co cn in T O CD

CM (0 (O

00 CO O

- - O N CO o O

t" cn o in o-o o o

0)Nh <n CM

in — o in t o cm o o o

oo in cm oo in T co — T oo o cm o o oCD CM cn CO 00

CO If; — CO CM o CM o o o

<n rCM (D CM O CM O

CM -I

UJ UJ 3 O

CL CL *-* O *-* (0 < »- -M- D

U- L/D O

cr < o > o UJ

o *-« -J o

CO -jUJ UJ 3 O

cc ct •—

<

O CO O < H

-I H D

u. m o

co a:o»-

cc < o >

O UJ

o a

u -*

-i a

T -*

UJ UJ 3 O

cc o: m

o *-* in o < »- -H-D

u- in o

C < o > O UJ

3 O

o:

*-* in < th- ZD m o

cn < o > o ui

s i O i

< -J l CD —1 UJ i 3 l o: ck * O -* c < »

COo

cn T

(O

CM

CMin

oo <n oocno

CM CM CM CM

CM CM CO CO T t in in co (0

23

.069

1

0.

FLOOR

6

COMPARTMENT

1

.029

143-1.

.051

?4b

ELEVATOR

-.021

169-1.

.650

-24 3:
TEN

STORY

BUILDING

WITH

ELEVATOR

ADJUSTED

TEMP

FIXED

DIFFERENTIAL

FLOW

FLOW

FLOOR

COMPARTMENT

PRESSURE

PROFILE

FLOW

CONNECTION

TO

PRESSURE

COEFFICIENT

AREA

FLOW

ID2

ID2UJ2UJ2

cm cm cm <t cm co 'T to cm co «t cm to o to cm to co cm co o o co in cm

CO O CM • CO 00 O Q) t" O O

CM If) CO CM CM

CO CO if) . CO 00

to O) CM to U)

lO <T

'T’TO)

CM O P»

tO <T

<r « in cm • min • CM CM (''**- co t>*

• to in •

- o o

Lf) CM a") O'TN

— o - o o o

Lf) lo if) cm in m in o to o <r o to - o oin cm m

o r- - o o o

in in m cm in o to O^TN

- o

in in o to

*r if) if) rr *r CO CD If) co O) 'TOO to

«r in in

•*o o f) %* o J)

<r ’T ld if) co o co o', in N to ^ O CD

*r ^ in in o- CO O CO CD if) ro CD

or to or o cd to

«- - to «- in

CM CO CO CM *- O CM o o o

<T 'TO

»- CD <T *- -r- O CM o o o

CO T O

o o to O CO o

co in o o

o o

co co in O *“ f** CO oCO CM o o

o o

CM CM CM CM

2

IDS 2

H- O

CL •— < hCL U

S ID O cl (J —

-J Q

-iID ID 3 Q

CL CL o *-• tn O < •- JHD u_ lO O

z

ID 2ID zID 2 2 z 2

H- f- o )-

CL cc *—

«

CL < < h- <

CL CL (J CL 2 2 UJ 2 O O CL O

(J CJ U

r- cl -J oCD -J co tr o

hCL < a

ID ID 3 Q

cr —

«

O

hcr < O > o -« l/> o > o UJ o < »- o UJ

-J _J _J —J _J

li_ UJ lx. L/) o U. UJ

z

ID5 z

- o

o: *-• < HQ- U

5 ID O a u

-i co O) -JUJ UJ 3£ Q

a: a: *-• o ^ to o < h- DH 3

u. t/> O

Z

UJ2

CL<

CL2

O

u

O) a

o

a: < O >

o UJ

2 2

»- O

cr

O a

o -•

-l Q

O -J

c- UJ ID 3 a

a a «

o if) o < »- -Jh D

u. D) O

O

CJo a

tr- o>-

o: < O >

uj

I

s* (D CD CD CD

24
TEMPERATURE

PROFILE

1

SHAFT

FLOW

COEFFICIENT

80000.

incncr, com — r^<roo*~ c cn o ^ o O n cn • o r c t 1 C 3) t o —

,?<Tininminininioio

Mnoo-f'^^rD(niiV7 c o - in co cd m cn • <% ro r - 3) ^ o ^ fN > s

lonnncNCiCN’-- r-*

3 < O ul 00)(DCT>(X>00)0)<T>0) nnncnnnnnnn oooooooooo2 < O UJ ooooooooooo lOlOlOlOUDiDCOCOlOiOO OOOOOOOOOOO

O UJ UJ —

*- 2 O

in O ~

O -J u_

"D u. U- Q UJ < Oo

inminminininininin 000000000)0 oooooooooo »- 2 o

o o

DJU.

“3 U. U. O UJ < o

u

^r*r<T<3

, «T<r<T^*T*rco ooooooooooo* iD.OlDlOlDlOiOlOlDlO^ — «- *- o

Z D

uj cn a cn UJ UJ

u. cc u_ a

cNr^r-mcNC"- CDO— OLOCOCDOOr-

* CNCNCNCNCNCNCOCO z o

uj cn

cc cn

UJ Ul

u. cc

U_ CL

(00«Nina3»-inoinojtcnin^rmcNCN — — oor* -oooooooooo

— CNCNCNCNCNCNCNCNCN

z z

UJ UJ S 2

I- »- a cc <

cl a S 2 o o o o

z z z z z z

1 1

1

III III 111 III III ssssss »— I— !— •— >—— a CL oc cc cc cc <<<<<< a o. o. o. a a ssssss o o o o o o

o o o o o o

z z

Ul Ul s s

I- »-

cc cc < <

CL CL

z z : UJ UJ I 2 2 : H- — I cc cc < < a cl 2 2 o o o o

cc cc < <

CL CL 2 2 o o o o

*-• — CN CO ^ in. CO MDOO •-• *- c* O

z z z z z z III 111 111 III III III 2 2 2 2 2 2 Z

r— h- h- — I— O

CC CC CC OC CC CC *-* <<<<<< bcl a a a o. n u

2 2 2 2 2 2 Uj O O O O O O cc O O O O O u

o

If) ID N CD O O

o

Ul 2

X o

crccccccoccccccro'cc oooooooooo oooooooooo

oooooooooo inmininincninininin

inininininininminin

O u

CC UJ a oo

UJ

OC 2

=> o »- —

J

< a Q

UJ 2

X O

ccaccccccccacccccc-* OOOOOOOOOOO) OOOOOOOOOOi- _J_J_J_J-J_l_»-J-J-JO UUU.U.UUU.U.UUO

OOOOOOOOOO

0 CN CD n cun »- co ’j o

co^cninocDCNt^na)

OO)C^C0 incOCNOO>C^

CN — — — — — — —

r^coom — co^rococo (OCNMDCMTO^^N

f'-m<Tco — oo>r^cn<y

-cNco^rincoc^cocno — oi<o<rtn(o^©a>o

25

NET
. .. ~ - - -

.... . : r .• ,j ‘ • • • • * . • *5 -,

--S ' '
APPENDIX D. PROGRAM LISTING

26
/
(*iAnnnftnnnnr\nnAnrir>rinnnnnAononAriAr\ftnnAAnAnonnAnAn

MAIN PROGRAM

£NBS*PLI BS.SHOW A. MAIN

C

C COMPUTER PROGRAM FOR AIR FLOW

C SPECIFICALLY FOR ANALYSIS OF

CC

c

cc

c

ANALYSIS IN BUILDINGS

SMOKE CONTROL SYSTEMS

FROGRAM VARIABLES

IS ACCEPTABLE

OF FLOOR

AI LEAKAGE AREA OF INTERNAL CONNECTION

AO LEAKAGE AERA OF CONNECTION TO OUTSIDE

C FLOW COEFFICIENT BETWEEN BUILDING POINTS

CO FLOW COEFFICIENT TO OUTSIDE

CS FLOW COEFFICIENT OF SHAFT

E LIMIT WITHIN WHICH CONVERGENCE

F NET FLOW INTO POINT I FC FLOW BETWEEN INTERNAL FOINTS

FF FIXED FLOW INTO POINT I FO FLOW TO OUTSIDE

FSS NET FLOW INTO SHAFT IS

F HEIGHT FROM GROUND TO MIDPOINT

IB U G OUTPUT VARIABLE

ICON V INTEGER USED IN SUBROUTINES BLDGP AND SHAFTP

IF ICONV = 0 THEN THE PRESSURES WERE UNCHANGED

IFLCCR 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 NC. OF OUTSIDE CONNECTIONS

NFS 1 BOTTOM FLOOR OF SHAFT

NFS 2 TOP FLOOR OF SHAFT

NH NO. OF FLOORS

NPC NO. OF OUTSIDE PRESSURE PROFILES

NS NO. OF SHAFTS

NS 1 I VALUE FOR START OF SHAFT

NS2 I VALUE FOR END OF SHAFT

NT TOTAL NO. OF POINTS (BLDG ANO SHAFT)

NT P NC. OF TEMPERATURE PROFILES

P PRESSURE AT POINT I PFO OUTSIDE PRESSURE FROFILES

FO OUTSIDE PRESSURE

PS PRESSURE PROFILE OF SHAFT - WORKSPACE

PZ PRESSURE DUE TO ELEVATION DIFFERENCE

T TEMPERATURE FRGFILE ARRAY

TITLE PROJECT TITLE

TITSH SHAFT TITLE

FROGRAM PARAMETERS

MB MAX NO. OF BUILDING COMPARTMENTS

MM MAX NO. OF POINTS

MS MAX NO. CF SHAFTS

MC MAX NO. OF CONNECTIONS FOR ANY POINT

MPO MAX NO. OF OUTSIDE PRESSURE PROFILES

MTP MAX NO. CF TEMPERATURE PROFILES

MFL MAX NO . OF FLOORS

27
n

n

a

»

h

inn

n

in

m

or.r>r»ro

r>

n

n

n

MAIN PROGRAM

PARAMETER (MM=140 > MS = 8 , MC = 9 , MFC=2 . MT P=2 » MFL= 25 * MB=50

)

CCMMCN NT * P( MM ) ,C(MM, MC) »NC (MM) , JC ( MM, MC ) , ITS(MS)

,

1 FC(MM.MC) » P Z ( M M , MC ) , PO ( MM , 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) «N » NS * NPO • ICONV, E, I BUG, AI ( MM, MC ) , AO( MM, MPO) , T I T SH (M S , 5 ) , 4 NF, KMFL), IFLOOR (MM) ,T( MTP , MFL ) ,NFS1 ( MS) ,NFS2( MS) , IT ( MB ) .NTP

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

DOUBLE PRECISION P.PO.PS

COMMCN /RUN/I RUN

Cl MEN S ION Bl( MM,MC) ,B2(MM,MF0)

N I TER=5000

IR UN = l CALL INPUT TO READ DATA

CALL INPUT

E= 0,2

1CS= 1 SAVE A I ( I , J ) IN BKI.JI AND FIND

MAX VALUE OF A I ( I , J )

AZ Z = C AMAX = 0 DO 10 1=1, NT

CO 8 J = 1 , MC

El ( I « J > = A I ( I , J)

IF(AI(I,J) .GT. A MAX ) AMA X=AI ( I • J

)

CONTINUE

DO S J = 1 » MPO

E2( I ,J ) = AO( I, J)

IF (AC ( I, J ) .GT. AMAX) AMAX=AC( I , J

)

CONT INUE

0 CONTINUE

ADJUST FOR LARGE VALUES OF FLOW AREA

IF (AMAX .LT. 0.2) GO TO 25

AZ Z = 1 AM= 0 .2/ ( AMAX— 0 • 1

)

EB = 0 • 1 * ( 1.0-AM)

CO 15 1=1, NT

CO 12 J= 1 » MC

IF(A1(I«J) .LT. 0.1) GO TO 12

AI ( I . J ) = AM*AI ( I « J ) +BB

2 CONTINUE

DO 14 J= 1 , MPO

IF(AC(I,J) .LT. 0.1) GO TO 14

AO < I , J ) = AM*AO( I . J ) +BB

4 CONTINUE

5 CONTINUE

TEMPERATURE CORRECTION

25 CALL CCRR

28
MAIN PROGRAM

C

C CALL INIT TO INITIALIZE PRESSURE ARRAY , p

C22 CALL INIT

C

c

C DO LOOP TO 30 IS ITERATIVE SOLUTION TO PRESSURE ARRAY

C24 DO 3C ITER=1. NITER

CC CALL BLDGP TO SOLVE FOR BUILDING PRESSURES

C CALL BLDGP

IC E = ICON

V

IF ( I CB .EQ. 0 .AND. ICS .EQ. 0 ) GO TO 40

CC CALL SHAFTP TO SOLVE FOR SHAFT PRESSURES

C CALL SHAFTP

ICS= 1C0NV

IF ( I C B .EQ. 0 .AND. ICS .EQ. 0 ) GO TO 40

CC CALL PZAD TO CALCULATE PZ TERMS

C CALL PZAD

30 CONTINUE

C

C IF ROUTINE FAILS TO CONVERGE IN NITER

C ITERATIONS PRINT ERROR MESSAGE

C HR I TE ( 6.800)

AO CONTINUE

HRITE(6»801)ITER

IF ( A 2 Z .EQ. 0 . ) GO TO 42

AZZ = 0.

DO 6C 1=1. NT

DO 5 C J=1.MC

50 AI ( I , J )= B 1 ( I, J ) DO 55 J=1 , MPO

55 AO ( I , J ) = B2( I, J

)

60 CONTINUE

CALL CCRR

GO TC 24

C

cc

C CALL OUT TO OUTPUT SOLUTION

C

A2 CALL OUT

C wRiTE(e.eos)

STOF

c

c

C FORMAT STATEMENTS

C£00 FORMAT (/////5X. 35 < 1H1 )//5X.

+35HF A I LURE OF MAIN PROGRAM TO CONVERGE / /5X , 35 ( 1 Hi )//

)

29
MAIN PROGRAM

801 FORMAT! I OX , 1 5. 5X , 1 1

H

I 7ERAT I C NS ) £05 FORMAT (11-1)

END

aFDG.P SUBROUTINE INPUT. L.l

30
c\

r>

rf\ri(^nrtr>nr\r\r>nr\

r*r>r>

r

n

n

n

SUBROUTINE INPUT

$NBS*PLI B$ .ShGW A. INPUT

SUBROUTINE INPUT

THIS ROUTINE READS AND PRINTS DATA

ANC INITIALIZES PZ ARRAY

PARAMETER ( MM = 14 0 . MS = 8 . MC=9 . MFO=2 , MT P=2 , M FL= 25 , M B=50

)

COMMCN /PZZ/ PGZ

COMMON / IO/T ITLE ( 18) , I OUT, IUNIT,NC0MP< MFL) ,SNCOMP( MFL ) COMMCN NT, P ( MM ) ,C(MM , MC ) ,NC (MM) , JC ( MM, MC) , ITS( MS) , 1 F C ( M M *MC) , PZ (MM, MC),PC( MM, MFC), 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.IBUG,AI(MM,MC),AO(MM,MPO),TITSH(MS,5),

4 NH, h( MFL) , IFLOQR ( MM) , T ( MT P , MFL ) , NFS 1 ( MS) , NFS 2 (MS), IT(MB)»NTP

5 , NCC (MM ) , JOC ( MM, MPO) , TOUT

COLELE PRECISION P.PO.PS

CHARACTER PAR *6

DIMENSION II( MFL) ,TT( MFL) , PAR < 7 ) , CW ( MPO ) , PH( MFL) ,NZZ(MM

)

CATA FAR/' MM',' MS*,' MC • ,

• MFO * ,• MTP • , • MFL* ,

• MB'/

IBUG=0

READ AND WRITE PROJECT TITLE

READ (5,600) (T ITLE ( I ) . 1 = 1 . 18)

WRITE(C,e01)( TITLE ( I) ,1 = 1 ,18)

READ GENERAL DATA

TOUT = OUTSIDE TEMPERATURE

IUN IT = 1 FOR SI UNITS

= 2 FOR ENG UNITS

ICLT = 0 FOR NO SUMMARY OUTPUT

OThERWlSE IOUT IS FILE NO. TO

WhICH SUMMARY OUTPUT IS WRITTEN

RE AC <5,7 00) TOUT, I UNIT, I OUT

Wrt ITE(€,41 1 )TOUT» IUNIT.IOUT

IF ( ILNIT .GT. 2 .OR. IUNIT .LT. 1 ) GO TO 105

READ hEIGHTS

N N = 0 FOR INPUT OF ALL hEIGHTS

C NN = 1 FOR CALCULATION OF hEIGHTS

C READ (5 ,700) NH , NN

WR ITE (C,412)NH,NN

I F ( N h »LE. MFL ) GO TO 89

IP AR =6

GO TO 110

69 IF (NN .EG. 1 ) GO TO 97

READ (5 ,700) (H( I ), I =1 ,NH)

WR ITE(6«413)(H(I) , 1 = 1 , NH

)

GO TO 99

97 READ (5 ,700 )H( 1 ) ,DH

wR I T E ( e,414)H( 1 ),DH

CO 9 6 1= 2 , NH

31
r>

r>

r>

SUBROUTINE INPUT

I M = I— 1 5 8 MI) = MIM)+DH

READ TEMPERATURE PROFILES

59 REAC (5»700)NTP

WRITEC 6.415JNTP

IF ( N T P .LE. MTP ) GO TO SO

IPAR=5

GO TC 110

50 DO 2 IF=1,NTP

REAC (5 ,700)NNN, (I I ( J ) , TT ( J ) , J = 1 ,NNN)

WR ITE(6»416)NNN»( I I(J ) »TT(J) , J = 1 ,NNN)

IF ( NHN .GT. 1 ) GO TC 2 CO 1 I FF = 1 » NH

1 T( IP .IFF )=TT( 1 ) GO TC 3

2 J= 1 JP 1 = 2 CO 4 I FF = 1 9 NH

T( IP , IFF ) =TT ( J ) + ( TT ( JP 1

) — TT ( J ) )*( IFF- I I( J) )/( II ( JPI )-I I ( J))

IF ( I FF ,NE. I I ( JP I ) ) GC TC 4

IF ( JF 1 .EQ. NNN )G 0 TO 4

J=JP 1 JPl = J+l

4 CONT 1NUE

2 CONTINUE

CC

C READ OUTSIDE PRESSURE PROFILES

C NN=0 FOR INPUT OF ALL PRESSURES

C NN=1 FOR CALCULATION BY POWER LAW

C REAC (5.700 JNPO.NN

WR XTE( 6»417)NPC »NN

IF ( N FO .LE. MPO )GO TO SI IP AR = 4

GO TC 110

51 IF ( N h .EQ. 1 ) GO TO 81

CC READ ALL OUTSIDE PRESSURES

C CO 6 1=1. NPO

€ READ (5,7 00 ) PGZ . ( PFO ( J . I ) , J= 1 , NH

)

WR ITE(€.413)PGZ,( PFO ( J . I ) » J= 1 » NH

)

GC TC es

CC CALCULATE OUTSIDE PRESSURES

C PATMCS IS ATMOSPHERIC PRESSURE (PA)

C

8 1 REAC(5.7CO)VW,hW,XW,(CW< I ) .1 = 1 ,NPO)

WRITE( e.419)VW,HW,XW,(CW(I) *1=1, NPO

)

IF ( IUN IT .EQ. 1 )VW = VW*0 .2778

IF ( I UN IT .EQ. 2)VW = VW*0. 4470

PATMCS=10 1325.

TO 0=1CLT + 273.

IF ( IUN IT .EQ. 2 )TOO=( TOUT +460 . ) /I .6

PVA= 17 8 . 4*VW*VW/T0C

32
<-» V

SUBROUTINE INPUT

Z=-0 .034 17/TOO

IF ( I UN IT .EQ. 2)Z=0.3048*Z

CWM=CW ( 1 ) IF ( N FQ .EQ. 1 ) GO TO 212

DO 211 1 = 1. NP O

IF(CMI) .LT. CWM)CWM=CW( I)

211 CONTINUE

212 PG Z= F ATMOS* EXP < H( NH ) *Z ) +CW M*PV A* C ( H (NH ) /HW )** ( 2 . *X W ) )-100.

CO 210 1=1. NH

PH ( I )=PATMOS*EXP( H ( I ) 4Z

)

210 CONTINUE

DO 82 1=1. NPO

CO 8 2 J= l . NH

PFO ( J * I ) =PH( J) +CW( I ) *PVA*( ( H ( J )/HW)** (2 .*XW ) )-PGZ

82 CONTINUE

C

cC BUILDING DATA INPUT

C NFLS = NO. OF FLOORS IN BUILDING

C IF 1 = LOWER FLOOR IN SERIES OF SIMILAR FLOORS

C I F 2 = UPPER FLOOR IN SERIES OF SIMILAR FLOORS

C NCC = 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

C85 1 = 0 SNCCMP ( 1 ) =0 * READ (5 .700) NF L

S

WRITE(6.420)NFLS

IF ( NFL S .GT. NH )GO TO 106

7 REAC (5.700) IF1 , IF2.N0C

WR I T E ( 6 « 40 0 ) I F 1 .IF2.N0C

IF ( I f 1 .GT. IF 2 )GO TO 107

NCOM F ( IF 1 ) =NOC

IFP= IF 1+1

SNCO M P ( I FP ) =SNC0MP ( I F 1 >+NOC

CO 1C IZ=1 . NOC

1= 1 + 1 READ (5 .700) NZ . NA. NN0.FF( I) .IT < I)

WR IT E( 8.401 )NZ,NA .NNO.FF (I ) .IT (I

)

NZ Z ( I ) =NZ

NN=N 2 + NA

IFLOCR ( I ) =1 F

1

IF (NN .LE. MC)G0 TO 111

IP AR = 3 GO TC 110

111 IFCNNC .LE. MPO )G 0 TO 112

IP AR =4

GO TC 110

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

NC ( I ) = NN

IF (NZ .EQ. 0) GO TO 83

INFUT CONNECTIONS TO COMPARTMENTS ON SAME FLOOR

REAC (5.7C0) ( JC( I. J ) .C< I * J) .AI< I . J) . J=1 .NZ)

WRITE (6 .402)

WRITE(6.403)(JC(I.J).C(I.J).AI(I,J).J=1.NZ)

33
SUBROUTINE INPUT

CO 6 2 J=I.NZ

62 JC< I . J )= JC( I. J)+SNCCMP< IF1

)

63 IF <NA .EC. 0) GO TO 6

CC INFUT CONNECTIONS TO COMPARTMENTS ON FLOOR ABOVE

C NP = N 2+1

READ (5 ,700 H JC ( I , J ) ,C ( I . J ) . AI ( I . J ) , J= NP, NN)

WR ITE( 6,404)

WRITE (6,403) ( JC<I,J),C(I.J),AI(I.J) . J=NP, NN)

DC 66 J=NP,NN

6 6 JC ( I , J ) = JC( I , J ) +NCCMP ( I F 1 ) + SNCO MP ( IF1 ) e NCC ( I

) =NNO

I F ( N NO .EQ. 0 ) GO TO 10

C

C INFUT CONNECTION TO OUTSIDE

C REAC <5,700) (JOC(I * J J ) »CQ(I»JJ) «AO( I * J J ) , J J = 1 , NNO)

WR I T E ( 6 , 40 5 ) WR ITE< 6,40 3 ) ( JCC( I , J J ) „ C C( I , J J ) , AO ( I , J J ) , JJ = 1 ,NNO)

CO 9 J J= 1 * NNO

J= JOC( I . JJ)

9 FO ( I , J J ) =PFO( IF 1, J ) 10 CONTINUE

I F ( I F 1 .NE. IF 2 )GO TO 11 IF ( I F 1 .EQ. NFLS) GO TO 20

GO TC 19

CC ASIGN CATA FOR FLOORS SIMILAR TO FLOOR IF1

C11 IFP=IF1 + 1 CO 17 IFF=IFP,IF2

NCOMF ( IFF )=NOC

IFFP= IFF+1

SNCOMP! I FFP ) =SNCO MP ( IFF ) +NOC

CO 16 IZ=1.NQC

1= 1+ 1 I1 = I2 + SN CO M P ( I F 1

)

IFLCCR ( I )= IFF

FF ( I )=FF( II)

IT< I )= I T < I 1 ) NN=NC< 1 1 ) NN0=NCC< ID

NC ( I ) = NN

NCC ( I ) =NNO

IF ( IFF .NE. NFLS ) GO TO 23

NN = N 2 Z ( I 1)

NC ( I ) = NN

23 IFINN .EG. 0 ) GC TO 14

DO 12 J= 1 , NN

C< I.D=C<I1,J)

A I (I , J ) = A I ( 1 1 , J ) JC ( I , J )= JC( II , J )+ SNCOMP { IFF ) — SNCOMP ( I F 1 ) 12 CONTINUE

14 IF (NNO .EQ. 0 ) GO TO 16

DO 16 JJ=l,NNO

JOC ( 1 , J J ) = JOC< 1 1. J J)

J= JCC( I , J J)

34
SUBROUTINE INPUT

CO ( I ,JJ )=CO( II , JJ ) AO ( I « J J ) =AQ( I 1 . JJ ) 15 FO ( I ' J J )=PFO( IFF, J ) 16 CONTINUE

17 CONTINUE

18 IF ( I F 2 .EQ. NFLSIGC TO 20

19 CONTINUE

GO TC 7 20 N= I N2 = N IF (N .LE . MB) GO TO 114

IP AR=7

GO TC 110

C

C SF AFT CATA INPUT

C114 READ (5 .700 )NS

IF ( N S .LE. MS ) GC TC 113

IP AR =2

GO TC 110

113 DO ICO I S= 1 » NS

READ <5.603 ) <T I T SH < I S . I ) • 1= 1 . 5

)

WRITE(6,406)( TITSH{ IS. I ) .1=1 .5)

READ (5.7 0 0>CS< IS) .NFS1 < I S ) .NF S2( I S ) , I T S< I S

)

WR I TE ( 6 • 407 ) CS( IS) » NFS1 ( IS ) ,NF S2 ( IS ) . ITS( IS)

N 1 =N 2 + 1 N2 = N 1 + NF S2 ( IS) -NFS 1C IS)

NS 1 ( IS ) = N1

NS2 ( IS ) =N2

IFF = NFS1 C IS)-1

READ (5,700) NNO . FF F.JCP.CC.AA

WR ITE(6.408)NNC.FFF, JCP.CC »AA

IF ( NNO .EQ. 0 ) GC TO 21

RE AC (5.7 CO) (JOCCN 1 .J) ,C0(N1 ,J) .AO(Nl.J) ,J = 1 ,NNO)

WRITE (6 ,403) ( J0C( N 1 , J ) . CO ( N1 . J ) , AO ( N1 , J ) , J= 1 , NNO)

21 DO 24 I=N1.N2

NC( I )=1

NCO ( I ) =NNO

IF F = 1FF+ 1 IFLOCRC I )=IFF

I F ( IFF .GT. NFLS) GO TO 25

FF ( I )=FFF

IF ( J CP .GT. NCOMP ( IFF) ) GO TO 25

JC ( I . 1 ) = JCP+SNCOM P ( IFF)

C( I, 1 )=CC

AI (I .1 ) =AA

26 IF ( NNO .EQ. 0 ) GO TC 24

CO 22 J= 1 , NNO

JJ= JCC(N1 , J)

FO ( I , J ) = PFO( IFF , J J ) JOC< I , J) = J J CO ( I .J ) = C0(N1 , J)

22 AO (I ,J )=AO(Nl , J)

GO TC 24

25 NC(I)=0

GO TC 26

24 CONTINUE

C

35
na

n

SUBROUTINE INPUT

EXCEPTIONS TO GENERAL SHAFT INPUT

NNN = NO. OF EXCEPTIONS

KE = 1 FOR FF EXCEPTION

C KE = 2 FOR OUTSIDE CONNECTION

C KE = 3 FOR INTERNAL CONNECTION

C READ {5 .7 00) NNN

IF (NNN .EQ. 0 ) GO TO 100

DO 69 IK=1,NNN

READ <5 .700 )KE. IFF MR I T E ( 6 • AO 9 ) KE t IFF

I=NS 1 ( IS )+ IFF-NFS 1 ( IS)

IF (KE .EQ. 1) GO TO 41

IF ( K E .EQ. 2 ) GO TO 42

I F ( K E .EG. 3) GO TO 51

C-0 TC 10 4 41 READ (5 .700)FF ( I ) WR ITE ( 6. 410 )FF ( I)

GO TC 69

42 READ (5»700)J . CCG* A AO

MR IT E ( 6. 405)

WRITE1 6 ,403) J .CCO.AAO

NN C= N CO ( I) IF ( NNC .EQ. 0 ) GO TQ 44

CO 43 K=1,NNC

IF(JCC(I.K) .EQ. J ) GO TO 46

43 CONTINUE

44 NJC=NNC+

1

NCO ( I ) =N JG

47 PO( I ,N JO )=PFO( IFF , J)

JOC ( I .NJO)=J

CO( I .NJO)=CCO

AO ( I «N JO )=AAO

GO TC 69

46 NJO -K

KK = K4 1 IF (CCC .NE. 0 ) GO TC 47

NJO =N NC— 1 NCC ( 1 )=N JO

IF ( NJO .EQ. 0 ) GO TO 69

DO 49 K= KK . NNC

KM=K—

1

FO ( I ,KM)=PO ( I . K ) JOC ( I « KM ) = JOC ( I , K

)

CO (I .KM)=CO(I.K)

49 AO ( I iKM =AO( I . K ) GO TC 69

51 READ (5 .700 ) JCP .CC . AA

WRITE (6. 402)

WRITE( 6 ,403) JCP.CC .AA

J= JCF+SNCOMP< IFF)

NN = N C ( I ) IFCNN .EG. 0 ) GO TO 53

DO 52 K=1,NN

IF(JC( l.K) .EQ. J ) GO TO 55

52 CONTINUE

IF (CC .NE. 0. >GC TC 53

MR ITE ( 6. 520 ) IS .KE . IFF

36
SUBROUTINE INPUT

S3

GO TC 69

NJ=NN+1

£4

NC ( I 1 =N

J

JC ( I ,N J )=J

55

C< I .NJ )=CC

A I ( I »N J ) =A

A

GO TC 69

NJ =K

KK=K4

1

IF ( A A .NE. O.iGO TO 54

N J = N N— 1 NC ( I )=NJ

IF ( N J .EQ. 0) GO TO 69

61

CO 61 K=KK,NN

KM =K— 1 JC( I . K M ) = JC ( I . K ) C ( I .KM ) = C< I . K

1

AI ( I ,KM ) =A I ( I . K ) 69 CONTINUE

100 CONT 1NUE

NT =N 2 IF (NT .LE. MM ) GO TO 160

c

c

IP AR = 1 GO TC 110

PRINT OUTSIDE TEMPERATURE

C160 WR ITE<6. 601 M TITLE! I) .1 = 1 . 12)

IF ( I ON IT .EQ. I )WRITE(6,800)TOUT

IF ( ION IT .EQ. 21WRITE! 6.5C01TOUT

IF ( I UN IT .EQ. 2 )T0UT = { TOUT-32 . )/l .8

TOUT =TCUT +273.

CC PRINT FEIGHT AND TEMFERATURE PROFILES

C IF < I UN IT .EQ. 1 )WR ITEC 6, 81 1 ) ( IP. IP=1 .NTP)

IF ( I UN IT .EQ. 2)WRITE(6.511 )(IP.IP=1»NT°)

HR I T E ( 6 , 8 1 3 1 CO 3 C 1FF= 1 »NH

20 HR ITE( 6.8121H! I FF ) «(T( IF»I FF ) * I P= 1 « NT P

1

C

C CONVERT TEMPERATURES TO DEG K

C DO 32 IFF=1.NH

CO 32 IF = 1 » NT P

IF ( I UN IT .EQ. 2)T( IP, IFF)=(T( I P . IFF )- 32 . ) /I .8

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

.

C

C PRINT OUTSIDE PRESSURE PROFILES

C IF I I UN IT .EQ. 1 ) GO TO 79

HRITE<6,514)( IP » I P= 1 * NPC ) HR I T E ( 6 • 8 1 3 ) CO 7 € IFF= 1 »NH

DO 77 J= 1 » NPO

77 PFO ( IFF. J )=PFO( IFF , J )/248.

8

HR ITE(6.515)H( IFF ) »(PF0C IFF » J ) »J = 1 . NP 0

)

CO 76 J= 1 . NPO

37
SUBROUTINE INPUT

7e PFC< IFF, J)=PFO( IFF,J)*24E.8

76 CONT IN UE

C-0 7C 6 3

79 WR ITE(6.814)( IP.IP=l,NPO)

WR I T E ( 6 , 6 1 3 ) DO 31 IFF=I ,NH

WR ITE(C» 615)H( IFF ) , (PFO< IFF ,J ) . J=1 ,NPO)

31 CONTINUE

CC CORRECT FOR CONNECTIONS ONLY INPUTED ONCE

C63 DO 6C 1=1, NT

NN=NC< I) IF { N K .EC. 0 ) GO TO 60

CO 56 J J = 1 , NN

J = JC ( I , J J) IF ( J .EQ. 0 ) GO TO 58

NN J = NC ( J ) IF < N N J .EQ. 0 ) GO TO 57

DO 56 IA =1 , NN

J

IF(JC(J, IA ) .EQ. 1)60 TO 58

56 CONTINUE

57 NNJsl'NJ+l

IF ( NN J .LE. MC)GO TO 5S

IPAR = 3

GO TO 110

£9 NC(J) = NNJ

JCC J ,NNJ )=I

C ( J , NN J ) =C ( I , J J ) AI ( J ,NNJ ) = AI ( I , JJ

)

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

PZ ( J *NN J ) =— PZ ( I , J J ) 58 CONTINUE

60 CONTINUE

CC CORRECT UNITS

C IF ( I UN I T .EQ. 2)CALL UNITS

CC INITIALIZE PZ FOR BUILD COMPARTMENTS

C67 CO 4 C 1=1, N NN = NC ( I) IF (NN .EG. 0 } GO TO 40

IA = I 1(1)

IF 1= IFLOOR ( I)

CO 36 J J = 1 , NN

J=JC ( I ,J J)

IF J= I FLOOR! J)

IF ( 1 F I .EQ. IF J )GQ TO 36

IB=IT( J)

TEMPA=0.5*(T( IA,IFI)+T< IB.IFJ) ) PZ ( I , J J) =3462. *(H( IF J)-H(IFI ) ) /TEMP

A

38 CONTINUE

40 CONTINUE

C

C INITIALIZE PZ FOR SHAFTS

C

38
SUBROUTINE INPUT

DC SC IS = 1 »NS

N1=N£1 (IS)

N2=NS2( IS)-I

IT T= IT S ( IS)

DO 4 S I=N1,N2

IF 1= IFLOOR ( I ) IF J= IF 1+ 1 TEMPA=0.5* (T( ITT, IFI)+T( ITT, IF J) ) PZ ( I . 1 )=2462.*( H< I FJ)-H( IFI ) ) /TEMP

A

45 CONTINUE

50 CONTINUE

cC CHECK SHAFT CONNECTIONS

C CO 240 IS=1.NS

N1 = N S 1 (IS)

N2=NS2( I S)

CO 229 I =N 1 . N2

NN=NC( 1 ) IF ( N F .EQ. 0 ) GO TO 239

CO 226 J=1.NN

JJ = JC( I, J) IF ( I FLOOR ( I ) .NE. IFLOOR ( JJ )) GO TO 103

236 CONTINUE

239 CONTINUE

240 CONTINUE

RETURN

C

C

C DIAGNOSTIC OUTPUT

C102 WRI 1E( 6.902 ) I . IT( I ) GO TO 109

103 KR ITE( 6 .903)

GO TC 109

104 taRITE(6.904)

GO TC 109

105 \»R ITE(6*905)

GO TC 109

106 1*R ITE(6.906)

GO TC 109

1C7 1*R ITE( 6.907)

GO TC 109

110 \#RI TE( 6.910 )PAR( I PAR)

C

C PRINT CORRECTED BUILDING DATA

C109 *RI 7E( 6.940)

DO 7 C 1=1. N NN =NC ( I) IF ( N N .GT. 0 ) GO TO 180

*RITE( 6. 941 ) I . IFLOOR ( I ) , IT ( I ) , FF ( I

)

GO TC 182

180 *R ITE( 6.942) I , I FLOOR! I ) . I T( I ) , F F ( I ) , JC( I . 1 > ,C ( I . 1 > . A I < I , 1 ) IF ( N F .EG. 1 ) GO TO 182

V*R ITE (6 »942)( JC( I , J ),C (I . J ) .AI ( I . J) .J = 2 ,NN)

182 NNC =FCC ( I ) IF ( NFO .EQ. 0 ) GO TC 70

39
70

CC

c

72

74

75

£0

C

CC400

401

402

403

404

405

406

407

408

409

410

411

412

413

4 14

415

4 16

417

418

419

420

500

£11

£14

£15

SUBROUTINE INPUT

WRITE(6,944){ JOC( I , J ) , CO ( I , J ) , AO ( I , J ) , J=1 ,NN0)

CCNT INUE

PRINT CORRECTED SI-AFT INPUT DATA

CO 8 C IS = 1 , NS

kR ITE ( 6,816 ) ( TITSH (IS, I ) • I L ,5 )

Vi/R ITE ( 6.806 )IS.CS( IS), IT c

( I

:

S)

N1 =N£ 1 ( I S)

N2 = NS 2 (IS)

HR ITE ( 6,807 ) CO 7 £ M

II z « N2

NN = NC ( I) IF (NP .GT. 0 ) GO TO 72

HR ITE ( 6 ,801 ) IFLCGR ( I) , FF ( I ) GO TC 74

HR ITE ( 6,808 ) IFLCOR ( I ) , FF ( I ) . JC (I , 1) ,C( I .1 ) , A I (1,1 ) IF (NN .EG. 1 ) GO TO 74

HR ITE ( 6 . 809 )( JC(I , J) ,C ( I * J ) i A I (I , J) , J=2 .NN)

NNC = NCC( I ) IF (NNO .EQ. 0 ) GO TO 75

HRITE(6,810)( JOC( I . J ) . CC ( I , J ) . AO ( I . J ) , J=1 * NNO

)

CONT 1NLE

CONT INUE

STOP

FORMAT STATEMENTS

FORM AT (5X, 5HIF1 =. 13,71-1, IF2 =,I3,7H, NOC =,I3)

FORMAT (5X » 4HNZ =,I3,6H NA = ,I3,7H, NNO = .I3.6H, FF =,F8«1,

+ 7 P , IT = , 13)

FORM AT (5X, 25HCONNECTION ON SAME FLOOR ) FORMAT (5X.3HJ =,I3.5H, C =.F10.3,5H, A =,F9.4)

FORM AT (5X, 26HCGNNECT ION TO FLOOR ABOVE ) FORM AT <5X, 22HC0NN ECT ION TO OUTSIDE ) FORMAT (5X.5A4)

FORMAT (5X.4HCS =.F9.l,8H, NFS1 =,I3,8H, NFS2 =,I3,7H, ITS=,I3)

FORMAT <5X, 5HNN0 =,I3,7H, FFF =,F8.1,5H, J =,I3,5H, C =, F10.3,

+ £ H , A = , F9 • 4

)

FORMAT (5X.4PK.E =,I3. 7Pi, IFF = ,I3)

FORMAT (5X ,4HFF - , F 8 • 1

)

FORMAT (5X.6PT0UT =,F6.C,9H, IUNIT =,I3.8H, IOUT =,I3)

FORMAT (5X,4 HN H =,I3,6H, NN =,I3)

FORMAT (5X »7HHE IGHT S /(10F8.2))

FORMAT (5X,6HH( 1 ) =,F8.2,6H, DH =,F8.2)

FORM AT (6X, 5HNTP =,I3)

FORMAT (5X.20HTEMPERATURE PROFILE / 15 , ( 1 0 ( 14 , F7 . 1 ) )

)

FORMAT (5X.5PNPO =,I3,6h, NN =,I3)

FORMAT (5X,5HPGZ = , F12 . 1 /17HPRESSURE PROFILE /(10F12.1))

FORMAT (5X ,4HVW =,F6.1,6H, HW = ,F6.1,6H, XH =,F4.2,6H, CW =, + ( 10 F4 .2 ) ) FORMAT </5X, 6HNFLS =,I3)

FORMAT (//10X,20HOUTSIDE TEMPERATURE .F6.1.2H FI

FORMAT! ///5X ,6liFEICHT ,5X» 29 H TEMPERATURE PROFILES (DEG F) / + 7X, 2FFT ,3X, 1 916)

FORMAT (////5X .6HHEIGHT , 5X , 26H0UTSI DE PRESSURE PROFILES

1 IIP (IN H20 ) / 7X , 2P.F T«3X«8I10)

FORMAT (Fll ,2,3X,8F10«3

)

40
SUBROUTINE INPUT

520 FORMAT !///5X» 15FERRGR IN SHAFT , I 2 . I 5HEXCEPT I ON KE - ,12,

2 X » EH FLOOR ,13//)

£00 FORMAT ( 1 8A4

)

£01 FOR M AT ( l F 1/// 1 OX » 18A4///)

£03 FORMAT (5A4)

700 FORMAT ( ) £00 FORM AT (//I OX » 20HOU T SIDE TEMPERATURE ,F6.1,2H C)

£01 F0RMAT(I13iFll.l)

£06 FORMAT < 10X * 1 2HSFAFT NUMBER , I 4/ 1 OX , 1 7H SHAFT COEFFICIENT .F10.2/

1 10X .20HTEMPERATURE PROFILE .14)

007 F0RMATC/21X.5HFIXED.25X. 4HFL0W * 12X . 4HFL0W / 1 OX , 5HFL00R , 6 X

,

1 4FFL0Vy,5X,12HCGNNECTED TO *6X . 1 1 HCOEFFIC IE NT .6X.6H AREA

2 / ) e08 FORMAT ( I 13, FI 1.1, 6X ,5HPG INT , I 5 ,F 16 . 1 , F 1 5. 4)

ECS FORMAT ( 3 OX , 5HP0 IN T . I 5 . F 1 6 . 1 .FI 5. 4)

eiO FORMAT (30X .7H0UTS IDE , I 3 , FI 6 . 1 . F 1 5 . 4

)

£11 FORM AT ( ///5X.6FFEIGHT . £ X , 29HTEMPERA TURE PROFILES ( DEG C) /

+ 7 X , 2HM ,3X, 1 9 16)

£12 F0RMATCFH.2.3X.19F6. 1 ) £13 FORMAT!/)

£14 FORMAT (////5X .6HHEIGHT . 5X . 26H0UTSI DE PRESSURE PROFILES

1 1 1 F (PASCALS) /7X.2HM .3X.8I10)

£15 FORM AT ( F 1 1 • 2 . 3X ,8 F 1 0 . 1

)

£16 FORMAT <///10X . 5A4

)

817 FORMAT! IOX.45HFLOW COEFFICIENTS CORRECTED FOR TEMPREATURE ) 502 FORMAT (10(/)«10X, 1 1HCC MPARTMENT ,14/

1 1 OX ,2 OHTEMPERATURE PROFILE .I4.17H DOES NOT EXIST /

+ 1 OX , 1 6 H PRO GRAM STOPPED .10!/))

503 FORMAT (10!/) ,5X .23HSHAFT CONNECTION ERROR * 1 / 1 0 > » 1 6 F PRO GRAM STOPPED ,10!/))

504 FORMAT (10!/) .10X.40HI NPUT ERROR IN EXCEPTIONS TO SHAFT DATA

1 / 1 0 X , 16 FPROGRA M STOPPED ,10!/))

505 FORMAT ( 1 0!/) , 10X, 37HINPUT ERROR IN UNIT TYPE DESIGNATION /

1 1 OX , 16FPR0GRAM STOPPED ,10!/))

5C6 FORMAT (10!/), 10X, 37HINPUT ERROR NO. OF FLOORS EXCEEDS NH /

1 10X « 1 6H FRO GRAM STOPPED .10!/))

5 C 7 FORMAT! 10!/) . 10X .25HINPUT ERROR IF1 .GT. IF2 /

1 1 OX , 1 6HFR0GRAM STOPPED ,10!/))

510 FORMAT (10!/) « 10X, 36HI NPUT EXCEEDS DIMENSION PARAMETER , A3/

+ 1 0 X , 1 6F PROGRA M STOPPED ,10!/))

530 FORM AT ! 1 OX ,3A6)

535 FORMAT!// 10X.26HFL0W COEFFICIENTS AS READ ) 540 FORM AT ( 1 OX, 15HBUILDING DATA //34X , 1 1 HTE MPERATURE .4X.5HFIXED,

1 12X ,2! 1 IX , 4HFL0 W ) / 1 OX , 1 1H CO MPARTMENT , 4X , 5HFL00R , 6X, 7HPR0F ILE

.

2 6X.4FFLQW, 5X « 1 3HCCNNECT ION TO » 4X , 11 FiCOEFFIC I ENT ,4X,

3 8 F AREA ) 541 F0RMAT(/4X.3I12.F14.i

)

5 42 FORMAT !/4X. 3112. F 14.1 .4X.5HP0INT ,I7»F 11 .2 ,F1 5.4)

5 43 FORMAT (58X ,5HPC INT , I7.F11.2.F15.4)

544 FORMAT (58X,9HOUTS IDE • I 3 . F 1 1 . 2 ,F 1 5. 4

)

END

aHDG.P SUEROUTINE CORR.L.l
nrir>ift

no

-*i

rh

<n

x>

(u

w

-

r\

r\

r\

c\

r>

r>n

SUEROUTINE CORR

£NBS*PLI E$ .SEC IN A. CORR

SUERCUTINE CORR

THIS ROUTINE CALCULATES ADJUSTED FLOW COEFFICIENTS

(C1.C2.CC1 . CC2

)

PARAMETER (MM = 14 0 ,MS=e . MC=9 , MFG=2 , MT P=2 , MFL=25 . MB=50

)

CCMMCN /CORR/C1 (MM , MC ) . C2 (MM , MC ) .CO 1 ( MM , MPO ) . C02(MM,MP0

)

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),FO( MM , MFO > ,CS( MS) . PS ( MFL ) .NS1 (MS) ,NS2(MS)

.

3 FSS (MS ) .N.NS ,NPO, ICONV ,E, I BUG .AI (MM, MC) . AO( MM, MPO) , TITSH(MS .5) , 4 NH, M MFL) , IFLOOR ( MM ) , T ( MTP , MFL ) .NFSl ( MS) ,NFS2( MS) , IT( MB) ,NTP

5 ,NCC( MM ) . JOC( MM, MFC) , TOUT

COUELE PRECISION P.FO.PS

CO 12 1=1. NT

CORRECT C

PATMCS= 101325.

EB= 1 CO 0.*SQRT(2.*PATMO S/287. )/ 1.2

NN=NC( I) IF ( I . CT. N ) GO TO 1 IP=IK I)

CO TC 4 CO 2 I S= 1 , NS

IF ( I .LE. NS2 (IS) .ANC. I . GE . NS1(IS))G0 TO 3

COM INUE

WRITE(6»700)

STOP

ip=nsas>

IFF= IF LOGR ( I)

T1 =T ( IF, IFF)

IF (Nh .EQ. 0 ) GO TO 10

CO 9 J = 1 ,NN

JJ = JC ( 1, J)

Cl (I , J )=BB*C( I, J) *AI ( I , J ) /SORT ( T1

)

IF ( J J .GT. N) GO TO 5

IP=IT ( JJ

)

GO TC 8 DO 6 I S= 1 , NS

I F ( J J .LE. NS2 ( IS ) .ANC. JJ .GE. NS 1 ( I S ) ) GO TO 7

CCNT INUE

WRITE (6,700)

STOF

IP=ITS( IS)

IFF = IFLOOR( JJ)

T2=T ( IP, IFF)

C2 ( I , J )=BB*C( I . J) *AI ( I . J)/SQRT(T2)

CONT INUE

CORRECT CO

10 NN C = M CO ( I) IF ( NhC .EO. 0 ) GC TC 12

CO 11 J= 1 , NNC

C0 1 ( I, J)=BB*CC( I, J )*AG( I , J)/SGRT(T1

)

42
SUBROUTINE CORR

CO 2 ( 1. J)=6B*C0< It J )*AC< ItJ)/SQRT(TOUT)

11 CONTINUE

12 CONTINUE

RETURN

700 FORM AT (/// 1 OX » 36HPR0GRAM STOPPED IN SUBROUTINE CORR

END

£ F DG P SUBROUTINE INIT.L.l

//)

43
A

A

*

AA

SUERCUTINE INIT

SNBS*PLIES .SHCfc A. INIT

SUERCUTINE INIT

C

CC THIS ROUTINE INITIALIZES Tl-E PRESSURE ARRAY

C PARAMETER <MM= 140 , MS = 8 , MC=9, MPG=2 , MTP=2 ,MFL=25 , MB=50)

PARAMETER (MEP=MB+1)

COMMON /C0RR/C1 (MM ,MC) . C2 ( MM , M C ) * CO 1 ( MM , M PO ) , C02( MM.MPO

)

COMMCN NT. P(MM ) ,C(MM,MC) ,NC(MM) , JC ( M M , MC ) . ITS(MS) . 1 FC( KM . MC) .PZ ( MV, MC) .PC( KM, MPC ) ,CO( MM , MPO ) ,F( MM ) , PFO( M^L, MPO )

,

2 FF ( MM ) ,FC( MM.MPO ) ,CS( MS) »PS(MFL) ,NS1 ( MS) , NS2 ( MS) , 3 FSS (MS) ,N,NS , NFO , ICCNV.E, I EU G . AI ( MM , MC ) .AO(MM.MPO) ,TITSH(MS,5)»

4 NH , K MFL ) . IFLOCR ( MM) ,T< MTP, MFL) .NFS1 ( MS) ,NFS2( MS ) . IT(MB ),NTP

5 ,NCC ( MM ) , JOC( MM, MPC) , TOUT

COUELE PRECISION P.PO.PS

CIMENS1CN SC ( MS ) , SCO (MS)

COMMCN /MAT/A( ME. MEP) ,XX ( MB) ,NNN

COUELE PRECISION A, XX

NNN = N

C CALCULATE AVERAGE OUTSIDE PRESSURE

SUM=C . CO 1C J=l,NPO

CC 1C 1=1, NH

SUM = SU M-f PFQ( I , J)

FA=SLM/( NPO*NH)

C

C Tt-E CO LOOP TO STATEMENT 30 ESTIMATES

C SHAFT PRESSURES

C DO 30 IS = 1 , NS

CC CALCULATE SHAFT PRESSURE C I F FERENCE , DP

C SUM = C . SUMN = 0 . N 1 =N S 1 ( I S

)

N2=NS2 ( IS)

CO 16 I=N1 , N2

SUM=SUM+FF( I)

NN =NC ( I) IF ( N N .EC. 0 . ) GO TO 16

DO IE J= 1 , NN

SUMN = SUMN+CK I , J)

15 CONTINUE

SC ( IS )=SUMN

16 NN C = NCC( I) IF (NNO .EQ. 0 ) GO TO 18

DO 17 J= 1 ,NNO

SUM,N = SUMN + C01 ( I , J ) 17 CONTINUE

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

44
SUBROUTINE INIT

CP 2= SUM/SUMN

SI GN = 1 . IF C CF2 .LT. 0.)SIGN=-1.

CP = S 1GN*CSIGN*DP2)**2

CC CALCULATE AVERAGE TEMP CF SHAFT

C

SUM=C.

ip=itscis)

DO 20 I = N 1 • N2

IFF = IF LOCR( I

)

20 SUM=SUM+TC IP, IFF)

TA = SLM/C N2-NI + 1

)

C

C ESTIMATE PRESSURE AT BCTTCM OF SHAFT . C HH=0.5*CHCNH)-HC1 ) )+HCl ) NF 1 = N FS 1 ( IS)

PB0T=PA+DP+3462.*C HH-F(NF1 ) )/TA

C

C ESTIMATE OTHER SHAFT PRESSURES

C

PCM ) = PBCT

NM =N 2 - 1 CO 24 I=N1,NM

IP1=1+1

24 PC IP 1 ) =P( I )-PZ( I. 1

)

3 0 CONTINUE

C

C ENC CF SHAFT PRESSURE ESTIMATES

C

C SET UP MATRIX FOR BUILDING COMPARTMENTS

C NP 1 = l CO 50 1=1, N NN = NCC I)

SUMI 1 = 0.

SUMNF=0.

IFCNN .EC. 0 . ) GC TO 42

CO 4 C JJ=1,NN

J= JCC I ,J J)

IF C J .CT. N ) GO TO 34

AC I ,.)=C1C I.JJ)

SUMI 1=SUMII-C1C I.JJ)

SUMNF=SUMNP—Cl C I, JJ)*PZC I.JJ)

GO TC 40

34 SUMI 1=SUM1I-C1C I. JJ>

SU MNF= SUMNP—C 1 C I* J J ) *P C J ) 40 CONTINUE

42 NNC=NCCC I

)

IF CNNO .EQ. 0 ) GO TO 46

CO 45 K= 1 , NNO

SUMI I = SUM I I—CO 1 C I , K)

45 SU MN F = SUMNP—CO 1CI »K)*FOCI»K)

46 AC I, I ) = SUM 1

1

AC I * NP 1

) = SUMNP—FF C I)

50 CONTINUE

C

PBOT

45
c*

r.n

SUBROUTINE INIT

C WF ITE MATRIX

C IF < I EU 0 .EQ. 0 ) GO TO e4

V»RITE(€.802)

CO £2 1=1. N

£2 WRITE(£,803)(A( I, J) .J=l ,NP1>

Cc

C CALL ROUTINE TO SOLVE FOR INITIAL BUILDING PRESSURES

C84 CALL SJMEQ

C

C OUTPUT INITIAL PRESSURES

C IF ( I EU G .EQ. 0 ) GO TO 85

WRITE( 8,800)

WRITE (6, 801 )( I,XX< I). 1 = 1 ,N>

NN=NS1 ( 1

)

WR ITE (6. 801 )( I.P( I ) « I=NN . NT

)

ASSIGN EUILDING PRESSURES

C89

SO

DO SC 1=1,

N

pc i >=xxc i)

RETURN

800

801

FORMAT C///8C6X , 1HI .4X.3HP

FORM AT (8(17. F7.1) ) )/)

802

6 C3

F OR MAT( /// 10X, 20 H MATRIX

FC RM AT ( 1 OX . 1 1 F 1 1 . 1 ) ENC

COEFFICIENTS / )

a " DG • P SUBROUTINE BLDGP.L.l

46
AA

A

A

A

A

A

A

A

A

M

A

A

A

A

A

A

SUBROUTINE BLDGP

ENBS*PLI E$ .SHOW A. BLDGP

SUERCUTINE BLDGP

Cc

C THIS ROUTINE CALCULATES STEADY STATE PRESSURES

C FOR BUILDING COMPARTMENTS

C

c

c PARAMETER (MM=140 , MS= e , MC=9 , MPC = 2 . MTP=2 , MFL=25 , MB=50

)

COMMON NT, P(MM ) ,C( MM ,MC ) ,NC (MM) . JCC 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 • NFO , I CON V , E , I BUG , AI ( MM , MC ) , AO( M M, MPO ) , TI T SH (M S ,5 )

,

4 Nh.KMFL), IFLOOR ( MM ) , T ( MTP , M FL ) .NFS! ( MS ) ,NFS2( MS ) . IT < MB ) , NT

P

5 , NCC( MM ) , JOC( MM, MFC) .TOUT

CCUELE PRECISION P.FO.PS.PI

IF ( I E U G .GT. 0 ) WR I TE( 6,806

)

IT M = 100

ICON V=0

CO IS 1=1, N

C

C CALCULATE NET FLOW .FI. INTO POINT I FI=PFLCW ( I ,P( I )

)

CHECK MAGNITUDE OF FI

IF (AES(F I) .LT. E ) GO TO IS

ICCNV=ICCNV+1

SET UP PARAMETERS FOR ITERATION

CP = 1 .0

IPHA SE = 1 CP I = C.

EE = 0 .2*ABS(FI

)

IF ( EE .LT. E ) EE=E

SI GN = 1 IF ( F 1 .LT. 0. >SIGN=-1

IK = 0 IF ( I EU G .GT. 0 ) MR I TE (6.602)

ITERATION TO RECUICE MAGNITUDE OF FN

IK = I K+ 1 NEW ESTIMATE OF PRESSURE .PI. AT POINT I PI = F ( I ) + S I GN*QP

CALCULATE NET FLOW , F N , INTO POINT I USING PI

FN = F FLCW ( I , P I

)

IF ( IEUG.GT.O) WRITE! 6.8 04)1, IK , F I , FN , FP , DP I , DP , DPP , P I , I PHASE

CHECK MAGNITUDE OF FN

IF(AES(FN) .LT. EE ) GO TO 10

CHECK NUMBER OF ITERATIONS

IF( IK .GT. I T M ) GO TO 25

CHECK PHASE

47
<u

r>

r>

SUBROUTINE ELDGP

IF ( IFHASE .EQ. 2) GO TO 6

C

C CHECK FOR TRANSITION FROM PHASE 1 TO PHASE 2

IF ( F I *FN .LT. 0.)GC TO A

CC PHASE 1 CP I = CP

CP = 5 .0 *DP

F I =F H

C-0 TC 2

CC PHASE 2

4 IPHASE=2

GO TC S

t IF ( F I* FN .GT. 0.)GC TO 8

CC NEW CP EETWEEN DPI AND DP

5 CPP=CP

FP = FH

CP = DF !< CPP-DP I )*F I/( FI-FN

)

GO TC 2

CC NEW CP EETWEEN CP AND CFF

€ FI = FN

CP I = CP

CF=DFI4< CPP-DP I )*FN/(FN-FP)

GO TC 2 10 P( I ) =P I 15 CONTINUE

C RETURN

25 WR I T E ( 6 t 80 0 ) STOF

C FORMAT STATEMENTS

00 FORMAT (///10X, 20( 1 H* >/// 1 OX, 22HEXCESSI VE ITERATIONS /

+ 10X.8HIN ELDGP ///10X»20<1H*) //////)

£02 FORMAT (//l IX, 1HI , 2X , 2H I T . 12 X, 2HF I , 1 3X « 2HFN* 13X, 2HFP, 12X , 3HDPI,

+ 13X, 2H CP. 1 2X. 3H DPP* 13 X, 2 HP I ,3X,5HPHASE /)

£04 FORMAT (8X,2I4,3E1 5. 4, 4F 15. 6. 15)

£06 FORM AT ( /// 1 OX, 6HBLDGF ) END

aFDG.P SUBROUTINE SHAFTP.L.l
n

n

n

n

n

n

n

n

m

n

n

n

a

n

a

a

a

a

aa

a

a

SUBROUTINE SHAFT®

aNES*PLI ES . SFO I* A.SHAFTF

SUERCUTINE SFAFTP

C THIS F 0 U T I NE CALCULATES STEADY STATE PRESSURES

FCR SHAFTS

PARAMETER ( M M= 14 0 . MS=8 • MC = 9 » MFC = 2 , MTP=2 . MFL=25, MB=50

)

COMMON NT. F< MM ) ,C( MM ,MC ) ,NC(MM) JC<MM, MC > . ITS(MS) .

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

,

2 FF(MM)»FO( MM . MPO ) .CS( MS) »PS(MFL) .NS1 ( MS) .NS2(MS)

.

3 FSSCMS) .N.NS . NPO . ICONV.E, IBUG.AI (MM.MC) . AO( MM* MPO )»TITSH(MS,5) . 4 Nt-.KNFL), IFLCCR ( MM ) ,T( MTF. MFL) .NFS1 ( MS) .NFS2( MS ) . I T ( M B ) .N TP

5 ,NCC(MM), JOCCMM, MPO). TOUT

OOUELE PRECISION P.PO.PS.PI

IF ( I EU C .GT. 0 )WR ITE( 6,806)

ITM= 100

ICGN V = 0

DO IE 1=1, NS

CALCULATE NET FLOW ,FI, INTO POINT I N1=NS 1 ( I ) FI=SFLCW ( I ,P( N1 )

)

CHECK MAGNITUDE OF FI

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

ICON V = ICCNV + 1

SET UP PARAMETERS FOR ITERATION

CP = 1 .0

IPHA SE = 1 DP I=C.

EE =0 «2*ABS( FI

)

IF ( EE .LT. E) EE=E

S I GN = 1 IF (F I .LT. 0. )SIGN=—

1

IK =0

IF ( I E U G .GT. 0 )WR ITE( 6 ,802

)

ITERATION TO REDUICE MAGNITUDE OF FN

IK = IK+ 1 NEW ESTIMATE OF PRESSURE .PI, AT BOTTOM OF SHAFT I PI = F ( N 1 )+SIGN*DP

CALCULATE NET FLOW , F N , INTO SHAFT I USING PI

FN= S FLOW ( I ,P I ) IF(IEUG.GT.0)WRITE(6,804)I.IK,FI,FN.FP,DPI , DP , DPP, PI , I PHASE

CHECK MAGNITUDE OF FN

IF (AES (FN ) .LT. EE ) GO TO 10

CHECK NUMBER OF ITERATIONS

I F ( I K .GT. IT M ) GO TO 25

C

C CHECK FHASE

49
r\

n

in

r>

n

ft

SUBROUTINE SHAFTP

IF ( I F F AS E .EQ. 2)G0 TO 6

C

C CFECK FOR TRANSITION FROM PHASE I TO PHASE 2

IF(FI*FN .LT. 0 • ) G 0 TO 4

C

C FHASE 1 CP I=CP

CP = 5 .0*0P

FI =Fh

GO TC 2 C

C FHASE 2

4 IPHA SE =2

GO TC 9 IF ( F 1* FN .GT. 0 • ) GO TO 8

NEW CP BETWEEN CPI AND DP

CP P=CP

FP = FN

DP=DFI+( CPP-DPI )*FI/(FI-FN)

GO TC 2 NEW CP E ET WEEN DP AND DPP

E F I =F N CP I=CP

DP=CFI+( DPP—DPI )*FN/(FN-FP>

GO TC 2 10 N2 = N £2(1)

CO 11 IF=N1,N2

II = IF+ 1-N1

11 F( IF 1 = PS (II)

15 CONTINUE

C RETURN

25 WRITE(C.800 ) STOP

C

C FORMAT STATEMENTS

C€00 FORMAT (///10X »20( 1 H* > /// 10X , 22HEXCESS I VE ITERATIONS /

+ 1 OX .9 H IN SHAFTP /// 1 OX , 20 ( 1 H* )//////

)

€02 FORMAT {//l IX. I H I , 2 X . 2H I T • 12X . 2HF I . 1 3X • 2HFN . 1 3X . 2HFP . 1 2X . 3HD PI.

+13X.2HCP . 1 2X . 3FDPP .13X.2HPI .3X.5HPHASE /»

€04 F0RMAT(8X,2I4,3E1 5.4.4F15.C .15)

€06 FORMAT ( ///10X, CHSHAFTP)

END

£ H DG .

P

SUBROUTINE PZAD L.l
SUBROUTINE PZ AD

£NBS*PLI E$ .SHOW A.PZAD

SUERCUTINE PZAD

C

C THIS ROUTINE CORRECTS PZ TERMS FOR PRESSURE

C PARAMETER ( M M= 1 40 . MS=8 . MC = 9 , MFC=2 , MTP=2 , MFL=25 , MB=50

)

CQMMCN NT, P ( MM ) , C { MM , MC ) , NC ( MM ) , JC ( MM , MC ) , ITS( MS ) ,

1 FC( MM ,MC> ,PZ ( MM, MC) ,P0( MM ,MPC) ,CO( MM , MPO ) , F( MM ) , PFO ( MFL » MPO )

,

2 F F ( M M ) , FO( MM , MPO ) ,CS( MS), PS (MFL) ,NSl ( MS) ,NS2(MS) , 3 FSS (MS ) ,N,NS ,NFO , ICONV, E. I EUG , AI ( MM , MC ) , AO( MM, MPO) ,TITSH(MS,5).

4 NH, F( MFL) , IF LCCR ( MM) , T ( MTP , MFL) ,NFSI ( MS) ,NFS2(MS), IT(MB)»NTP

5 ,NCC (MM ) , JOC( MM, MPO) , TOUT

COMMCN /PZ Z/ PGZ

COUELE PRECISION P,F0,PS

IF ( I EUG .GT. — 2 ) G O TO 1 WRITE (€.800)

DO 2 1 = 1 ,N

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

WRITE (€.801 ) ( I.J.PZd, J),J=1.NN)

2 CONTINUE

NP 1 = N + 1 WR ITE ( €, 802 )( IL,PZ( IL, 1 ) ,IL = NP 1 , NT)

1 CO 1C 1=1,

N

NN = NC( I)

I F ( N N .EC. 0 ) GO TO 10

IA= IT ( 1) IF I = IF LOCR( I

)

DO 8 JJ=1,NN

J= JC ( I , J J)

IF J= IFLOOR ( J ) I F ( IFI .EQ. I F J ) GO TO 8

IB= IT ( J

)

TE MP A = 0 .5*(T(IA,IFI) + T( I B , I F J ) ) FAVE = 0 .5* (P (I )+P ( J ) ) +PGZ

PZ ( I «JJ)=(0. 03416 *PA VE/TEMPA )* (H( IF J ) — H ( IFI )

)

€ CONTINUE

10 CONTINUE

CO 20 IS= 1 , NS

N l =N £ 1 ( I S

)

N2=NS2(IS)-1

ITT= 1TS( IS)

CO IS I=N 1 ,N2

IF I = IFLOQR( I

)

IF J = IF 1+ 1 TEMPA=O.S*(T(ITT,IFI)+T( ITT. IFJ)

)

J= 1+ 1 PA= 0 • 5 * ( P( I ) +P ( J ) ) +PGZ

15 PZ ( I , 1 )=(0.034l€*PA/TEMPA)*(H( IFJ)- H( IFI)

)

20 CONTINUE

RETURN

E00 F0RMAT(/10X ,1 OHINI TIAL FZ /)

601 FORMAT ( 10X ,3HPZ( , 12, 1H, 12. 4H ) = ,F12.4)

£02 FORM A T ( 1 OX , 3HPZ ( I2.6H.1 )

• ,F12.4)

£03 F0RMAT(/10X ,HHAD JUSTED PZ /)

END

51
rr>

(n

uro

^an

^

an

<*>

SUBROUT INE OUT

aNES*PLIE$ .SHOW A. OUT

SUERCLTINE OUT

C

cC THIS ROUTINE OUTPUTS FLOWS AND DIFFERENTIAL PRESSURES

C FOR ALL SHAFTS ANO BUILDING COMPARTMENTS

C

C PARAMETER ( MM= 14 0 , MS=8 . MC=9 . MPC=2. MTP=2 . MFL=25 . MB=50

)

CO M VC N /CORR/C1 (MM ,MC ) ,C2 (MM , MC ) ,COi ( MM, MPO) » C02( MM, MPO

)

COMMCN / IO/T ITLE ( 18) , IOUT , I UN I T , NCOMP ( MFL ) ,SNCOMP(MFL)

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) ,FO( MM.MFO) ,CS< MS) , FS( MFL ) *NS1 ( MS) ,NS2 (MS )

,

3 FSS (MS ) ,N,NS ,NFO , ICONV ,E , I BUG . AI (MM.MC), AO (MM, MPO) , TITSH(MS .5 ) , 4 NH, H( MFL ) , I FLOOR ( MM) , T( MTP .MFL ) ,NFS1 ( MS) , NFS2( MS) , IT (ME) ,NTP

5 , NCC ( MM ), JOC( MM, MFC) .TOUT

CCUELE PRECISION P.FC.FS

INTECER COM

C IUN IT = 1 FOR SI UNITS

I U N I 7 = 2 FOR ENG UNITS

WHEN I UN IT = 2 GO TO 100

IF ( I UN IT .EQ. 2 ) GO TO 100

BUILDING COMPARTMENT OUTPUT

1 = 0 IL = 0 WRITE(e,800)(TITLE( I) .1 = 1.18)

CO 2 C 1FF= 1 . NH

NNN=NCCMF( IFF)

IF (NNN .EQ. 0 ) GO TO 30

DO 29 IC =1 .NNN

1= 14 1 NN =NC ( I)

NNC=NCC ( I) IL = I L 4 HN 4NN04 2

IF ( IL .LT. £ 1 ) GO TO 2

taR I T E ( €.80 0 )( TI TLE ( I ) ,1=1 , 18

)

IL = NN4NN042

IF ( N N .GT. 0 ) GO TO 3

WR ITE ( 8, 801 ) IFF, IC,P(I ) , I T ( I ) ,FF( I ) GO TC 21

CO 2 C J= 1 , NN

JJ= JC ( I . J)

CP= F (JJ)-P(I)4PZ( I.J)

CC = C 2 ( I.J)

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

IF ( J J .LE. N ) GO TO 10

CO 5 I S= 1 . NS

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

CONTINUE

IF ( J .C-T. 1 ) GO TO 7

WR ITE(C, e02) IFF, IC .P ( I ) . IT ( I ) . FF ( I ) , ( TI TSH( IS .K ) ,K=1 ,5)

4 , DF ,CC. AI ( 1,1) ,FC ( 1,1)

GO TC 20

52
SUBROUTINE CUT

7 WR IT E ( £ . 603)(TITSH( I S . K ) « K= 1 S ) .DP . CC • A I ( I.J) .FC(I.J)

C-0 TC 20

10 IF J = IF LOOK ( J J ) CO*=JJ-SNCOMP( I F J

)

IF ( J .GT. 1 )G0 TO 12

WRITE! £.804) IFF,IC.P( I).IT(I),FF(I) , I F J , COM. DP. CC * A I ( I . 1 ) .FC( I • l ) GO TC 20

12 WR ITE(6.805)IFJ,CGM»DP. CC * A I ( I ,J),FC(I. J)

20 CONTINUE

21 IF ( NNG .EQ. 0 ) GO TO 29

CO 2 2 J-l.NNO

JJ = JCC (I.J)

CP=PC( I * J )-P( I

)

CC=CC2 (I.J)

IF ( 0 F .LT. 0 . ) CC=C01 (I.J)

23 WRITE! £.806) JJ.DP .CC • AO (I.J) *F0( I.J)

29 WR ITE (£. 807 )F( I ) 30 CONTINUE

WR IT E ( £ • 900

)

C

C SFAFT OUTPUT

C CO 6 C IS = 1 . NS

N1=NS1 ( IS)

N2=NS2 (IS)

22 WRITE(£.ei4)(TITLE(I) .1 = 1.18)

WR IT E ( £ . 80 8)(TITSH(IS*K).K = 1,E) . ITS( I S ) *CS( IS

)

CC SC I=N1 • N2

NN = NC( 1)

IF( NN .GT. 0 ) GO TO 35

WR ITE(£. 609)1 FLOOR (I).P(I) .FFC I)

GO TC 41

25 CO 4 C J= 1 . NN

JJ=JC( I.J)

CP=P ( J J )—P ( I

)

CC =C 2 ( I.J)

IF ( OF .LT. 0. )CC=C1(I.J)

IF J=IFLOCR( JJ)

COI»= JJ-SNCOMP( IFJ ) IF ( J .C-T. 1 ) GO TO 36

WR ITE(6. 810) I FLOOR (I).P(I).FF(I) , I F J, COM , DP ,CC * AKI, l ) .FC(I.l)

GO TC 40

2£ WRITE (£.81 1 )IFJ«CCM»DP»CC» AI (I.J). FC (I.J)

40 CONTINUE

41 NNC = NCC ( I ) IF (NNO .EQ. 0 ) GO TO 50

DO 4£ J= 1 . NNO

JJ = JCC (I.J)

CP = FC ( I.

J

)-P( I

)

CC =CC2 ( I .J)

IF ( CF .LT. 0 • ) CC=CO 1 ( I . J ) 46 WRITE(6.812)JJ.DP,CC.A0(I,J),FG(I,J)

£0 CONTINUE

WRITE(£ ,813)FSS( IS)

WR ITE (£.900

)

£0 CONTINUE

GO TC 165

C

53
SUBROUTINE OUT

C EUILCING CATA OUTPUT FOR IUNIT = 2

C100 1=0

IL = 0 WR ITEtC.eOO )( TITLE! I) .1=1 . 18)

CO 120 IFF=1,NH

NNN=NCCMF { IFF ) IF ( NNN .EQ. 0 ) GO TO 130

CO 129 IC=1 »NNN

1=1+ 1 FF 1 = F ( 1 )/ 0 . 47 19

PI II=P(I )/248.8

FF F = FF ( I )/ 0.4719

NN=NC ( I ) NNC= N CO C I) IL = I L + NN+NNO+2

IF ( IL .LT . 51 ) GO TO 102

WR ITE( €.80 0 )( TI TLE( I ) .1 = 1 , 18 ) IL = N N +NN0+2

102 I F ( N N .GT. 0 ) GO TO 103

WR ITE( €, 601 ) IFF, I C. PI 1 1 , IT< I > »FFF

GO TC 121

103 CO 120 J = 1 » NN

FCCC=FC( I, J )/0 .47 1

9

JJ = JC ( I , J)

CP= ( F ( J J )-P(I )+PZ ( I,J» )/24 8.8

AA 1=41 ( I » J ) /O • 092 9

CC = C 2 ( 1. J)

IF ( D F .LT. 0. )CC=C1( I.J)

CC=CC*33 .43

IF ( J „ .LE. N ) GO TO 110

CO 1C5 I£=l .NS

IF ( J J .GE. NS 1 ( IS ) .AND. JJ .LE. NS2(IS))GO TO 106

105 CONTINUE

106 IF ( J .GT. 1 ) GO TO 107

WR ITE< 6.60 2) IFF.IC.PIII . ITC I I .FFF »(TITSH( IS «K ) , K= 1,5)

+ .CP.CC.AAI, FCCC

GO TC 120

1 C 7 WR IT E ( €

,

602 )(TITSH( IS.K ) ,K=1 .5) .DP ,CC ,AAI .FCCC

GO TC 120

110 IF J= 1 FLO OR ( JJ ) CO ^ = J J - S NCC MP ( IF J ) IF C J . GT . 1 ) GO TO 112

WR ITE(£,604)IFF .IC.PIII , IT (I ) » FFF , I F J , COM, DP , CC . AAI . FCCC

GO TC 120

112 WRITE! 6, €05) I F J ,CC * » DP ,CC ,AAI .FCCC

120 CONTINUE

121 IF ( NNO .EQ. 0 ) GC TO 129

CO 123 J = 1 » N NO

FOC = F 0 ( I , J )/0 • 471

9

JJ = JCC ( I • J

)

DP= (FC (I.J )-P( I ))/248.8

AAC= AC ( I , JJ/0.0929

CC = CC 2 ( I.J)

IF ( D F .LT. 0 • ) CC = C0 1 (I.J)

CC=CC*22 .43

123 WRITE(€.€06) JJ.CP »CC , A AC , FOO

129 WR ITE ( €, 807)FF

I

54
SUBROUTINE OUT

130 CONTINUE

*R 1TE(£.901)

C

C SI-AFT OUTPUT FOR IUNIT = 2

C DO ISO IS=1,NS

CSS=CS< IS)/0. 02992

FF 1= F S S ( IS)/0.4719

N1=NS1 (IS)

N2=N£2 ( IS)

132 IaR I T E < £,814)( TITLE! I ) «I=1 . 18)

IaR IT E ( £ , 80 8 ) ( T ITSH < I S . K ) , K=1 , 5 ) * IT S< I S) .C SS

DO ISO I=N1,N2

FFF=FF ( I )/0.4719

FI I I=P(I )/248 .

8

NN=NC( 1) IF ( N N .GT. 0 ) GO TO 135

IaR ITE(£.€09) I FLOOR ( I ) * PI II »FFF

GO TC 141

135 DO 1 A 0 J=1,NN

FCCC=FC( I • J )/0 .4719

JJ=JC ( It J)

DP = ( F ( J J ) — P C I ) )/24 8.8

A A 1= A I ( I , J )/0 *0929

CC=C 2 ( I» J) IF ( C F .LT. 0 • )CC=C 1(1 t J

)

CC=CC*23.43

IF J= IF L0CR( JJ)

CCM= JJ-SNCOMP! IF J ) IF ( J .GT. 1 ) GO TO 136

WRITE(€,610) I FLOOR (I).PIIItFFF , IF J , COM . OP ,CC , AAI .FCCC

GO TC 140

136 <aRITE(£,€11)IFJ,C0M,DP,CC, AAI, FCCC

140 CONTINUE

141 NNC=NCO( I)

IF ( NNO .EQ. 0 ) GO TO 150

CO 146 J= 1 , NNO

FOC=FO(I , J ) /0 . 471

9

JJ=JCC ( I , J

)

CP = (FO(I ,J)-P( I }) / 24 8.8

AAG = AO ( I , J )/0.0929

CC=CC2 ( I ,J)

IF ( CP .LT. 0. )CC=C01( I . J

)

CC =C C* 23 .4 3

146 \aR ITE(£,£12)JJ, DP , CC , A AO , FOO

150 CONTINUE

IaR ITE ( £. 813 )FF I IaR I T E ( € , 90 1 ) ICO CONTINUE

C

C SUMMARY OUTPUT

C USER INSERTS \aR I TE STATEMENTS TO FILE IOUT

C

1 £ 5 CONTINUE

RETURN

CC

C FCFMAT STATEMENTS

55
SUBROUTINE OUT

c

d 0

1

FORM AT </4X

eo2 FORMAT (/4X

€03 FORM AT (53X

d 04 FORMAT </4X

1 F 1 1 .3 ,F 15

€05 FORMAT (53X

€06 FORM AT (53X

€09 FORM AT(4X,

€ 1 0 FORMAT (4 X,

1 F15.0.F10

€1 1 FORMAT ( 3 IX €12 FORMAT (3 IX eoo FORM AT (1 HI

FOR M AT (/4X, 1 3, I 10 ,F13.3,I8,F12.0,3X,5A4,F14.3,F15.0,F10.3,F11.1)

4,F14.3,F15.C,F10.3,Fll.l)

FORM AT </4X, 13, 1 10 ,F13.3 , I 8 . F I 2 . 0 , 3X , 5HFLC OR , I 3, 12H COMPARTMENT, 13,

0 , F 1 0 ,3 ,F11 .1 ) 5HFLC0R • 13 , 12H COMPARTMENT, 13 ,F1 1 • 3 ,F1 5 . 0 , F 1 0 .3 , FI 1 .1 ) 17HCUTSIDE DIRECTION, I3,F14.3,F15.0„F10.3,F11.1>

3, F 1 0. 3 »F 1 1 • 0 ) FORMAT (4X,I3,F10.3,F11 . 0 , 3X , 5HFLOOR , I 3 , 1 2 H COMPARTMENT, I3.FH.3,

FI 1 . 1 ) IFLOOR, 13 , 12H COMPARTMENT ,I3,F11.3,F15.0,F10.3,FI1.1)

HOUTSIDE DIRECTION , I 3 ,F 1 4 . 3 ,F 1 5 . 0 ,FI 0 . 3 , F 1 1 . 1 ) FORMAT (1 HI , 20X , 18A4 »/S4X » 8H A D J LSTED/35X , 4 HTEMP » 7X , 5HF IXED,28X«

1 I 2H DIFFERENT I AL, 5X.4HFL0W ,8X , 4HFL0W/4X, 5HFL00R.2X, 1 1 HCQMPARTMENT

2 » 2X » 8 F PRESSURE , 2X , 7H PROFILE » 5 X » 4HFL0W , 3X , 1 fcH CONNECT ION TO , 3 12 X * E,4X , 1 1HC0EFFICIENT ,2X ,8H AREA ,5X,4HFL0W /)

£01 FORM,

£02 FORM

£03 FORM

£04 FORM

1 F 1 1 £05 FORM,

£06 FORM

£07 FORM,

£C8 FORM

l 23F

2 2 ex

3 5 X , 4 ,2 X , £09 FORM

£ 1 0 FORM

1 F 15

811 FORM

£12 FORM

£13 FORM

814 FORM

9 0 0 FORM

1//5X

2/5X,

SOI FORM

I .F15.1.F10.4

FLOW

13.1 , I8,F12.0,3X ,5 HFLOOR ,13, 12H COMPARTMENT, 13,

Fll.l ) 00 R , 1 3 • 1 2H COMPARTMENT, 1 3 , F 1 1 . 1 ,F1 5 . 1 , F 1 0 .4 , FI 1 . 1 ) UTS IDE DIRECTION. I3,F14.1,F15.1,F10.4,F11.1)

1.4H NET)

£A 4//20X , 2 OH TEMPERATURE PROFILE ,13/ 20X,

COEFFICIENT , F 1 0 .0//7 2X , 8 HADJ USTED/24X , 5HFI XED

,

€h AREA , 5X , 4HFL0W /)

TO, 12X .8HPRESSURE.4X ,1 1 H COEFFICIENT

COMPARTMENT, 13. FI 1.1

,

1 //5X , 'FLOW

2 /5X ,* PRESSURE

END

5HFLC0R , 13 , 12H COMPARTMENT ,I3.F11.1,F15.1,F10.4,F11.1)

1 7H0UT S IDE DIRECTION , 1 3 ,F 1 4. 1 ,F1 5 . 1 ,Fl 0 .4 , FI 1 .1 ) F8.1 ,4 H NET ) 20X , 18A4)

iX.'TFE FOLLOWING UNITS ARE LSED FOR OUTPUT*

IN LITERS PER SECOND AT 21 DEG C AND 1 ATM*

IRE IN PASCALS */EX, ‘AREA IN METERS SQUARED')

5X ,

' THE FOLLOWING UNITS ARE USED FOR OUTPUT*

IN CFM AT 70 DEG F AND 1 ATM*

IN INCHS H20 ' / 5X ,

• ARE A IN FEET SQUARED*)

SHDG ,

P

SUERCUTINE UNITS. L.l

56
SUBROUTINE UNITS

5NBS*PLI ES .SHOW A. UNITS

SUERCUTINE UNITS

CC THIS ROUTINE CONVERTS VARIABLES H ,FF , Al * AO. CS TO SI UNITS

Cc PARAMETER (MM= 1 40 . MS= 8 , MC=9 , MFC=2 , MTP=2 , MFL=25 , M8=50

)

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

,

1 FC(MM.MC), PZ ( MM, MC) »PG( MM.MPO) . CO ( MM , MPQ ) , F ( MM ) , PFO ( MFL , MPO )

,

2 FF ( MM ) , FO ( MM , MPO ) ,CS ( N S ) , PS ( MFL ) ,NS1 ( MS) ,NS2 < MS) , 3 FSS(MS) ,N ,NS ,NPO , ICONV, E, I BUG, AI (MM, MC) ,A0( MM, MPO) ,TITSH(MS,5)

4 NH,M MFL) . IFLOCR ( MM) , T ( MTP , MFL ) ,NFS1 (MS) ,NFS2( MS ) , I T ( M B ) ,N TP

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

COLELE PRECISION P.PO.PS

C IMEFS ION B( 5

)

CATA B/0.3048 .248 .8,0.4719,0.02992,0.0929/

CO 10 1=1, NH

10 F( I)=H(I >*B(1 ) CO 2 C 1=1, NT

FF ( I } = FF ( I )*B( 2)

DO 1C J= 1 , MC

AI ( I ,J ) = A I ( I , J )*B ( E)

16 CONTINUE

DO 18 1 , MPO

AO (I ,J ) = AO ( I, J)*B( E)

18 CONTINUE

20 CONT 1NUE

DO 22 IS=1 ,NS

22 CS(IS)=CS(IS)*E(4)

RETURN

END

£FDG,P SUEROUTINE SIMEQ.L.l
a

r.A

n

o

r>

r>

r>

a

a

an

r» A

r>

aA

a

SUBROUTINE SIMEG

$NBS*PLIB$ .SHOW A.SIMEQ

SUBROUTINE SIMEG

C CFCLESKY'S METHOD CF SOLUTION OF

SIMULTANEOUS LINEAR ALGEBRIC EQUATIONS

PARAMETER (MM= 1 4 0 ,MS= e . MC=9 . MFC=2 . MTP=2 . MFL=25 * MB=50

)

PARAMETER (MBF=ME+1)

COUELE PRECISION A.X

COMMON /MAT/ A(MB,MBPJ.X(MB) ,N

NP 1 = N + 1 2ERO= 1 .OE-35

K = 0

SEE IF A ( I • 1 ) IS ZERO

IF SC ADC ANOTHER ROW TO ROW I IF (AES (A (1 , 1 ) ) .GT. ZERO ) GO TO 40

DO 31 I=1,N

IF<A(I.l) a NE • 0 o ) GO TO 22

31 CONTINUE

12 WRITE(e,804)K

ST C F 32 CC 32 J= 1 , NP

1

33 A ( 1 , j )=A (1 . J)+A ( I , J ) CALCULATE UPPER AND LOWER

TRIANGULAR MATRICES OVER ORIG

MATR IX A 40 AA = A (1.1)

CO 2 J=2,NP1

2 A( 1 *U)=A(1»J) /A

A

DO 1C 1=2. N K=0

STORE A(I.l) ... A(I.I) IN X ARRAY

IN CASE NEW A(I. I) IS ZERO

ROW I CAN BE RECALCULATED

4 CO 5 J = 1 .

1

5 X( J ) = A ( I ,J)

K=K+ 1 CO 10 J=2. NP1

IF ( J .GT. I) GO TO 8 JM 1 = 1 A A = C . DO 3 IR= 1. JM1

3 AA = AA+A( I. IR)*A(IR. J)

A ( I ) =A ( I . J ) — AA

CHECK IF A(I, I) IS ZERO

IF SC MULTIPLY OLD ROW I BY 2.

I F ( 1 .NE. J ) GO TO 10

IF ( AE S ( A ( I » I

)

) .GT. ZERO) GO TO 10

CO 6 JJ=1,I

6 A( I, JJ ) = X(JJ)

DO 7 J J= 1 » NP1

58
r\

r\

SUBROUTINE SIMEQ

7 A( I. w)*2.*A( I. J ) IF ( K . CT . 3) GO TO 12

GO TC 4

8 IM 1=1-1

AA = C . CO 9 IR= 1 « I Ml

9 AA=AA+A( I,IR)*A(IR. J» A< I, J )=< A< I, J )-AA ) /A (1*1)

10

CONTINUE

C ENC CF CALCULATION OF TRIANGULAR MATRICES

CC EACKkARD SUBSTITUTION

X(N)=A<N iNP 1

)

00 2 C 1 1 =2 * N

AA =0 . I=NF 1- II IP 1 = 1+1

00 IE J= IP 1 »N 15 AA=AA+A ( 1 1 J)*X( J)

20 X( I ) = A ( I iNPl ) -A

A

C 804 FORMAT (////// 10X, 1CHPR0GRAM FAILURE ,13/////)

END

ihDGiP FUNCTION F LO W • L 1
r>

r\

FUNCTIGN FLOW

aNES*PLI ES .SHOW A.FLCW

FUNCTION FLOW (PI.PJ.PZ ,C

COUELE PRECISION PI.FJ

THIS FUNCTION CALCULATES FLOWS EETWEEN TWO PO INTS

C IF C C .LT. 0«001 J G 0 TO 10

CP=PU-FI+PZ

SI CN= 1 .0 IF ( D F .LT. • 0 ) S I GN=— 1 . FLCW=SIGN*C*SQKT( SIGN*OP)

RETURN

10 FLOW = 0.0

RETUFN

END

a HOG »P FUNCTION PFLOW.L.l
r\

r\

FUNCTION PFLOW

SNBS4PLI ES • SFC W A. PFLOW

FUNCTION PFLOW ( I* P I )

C THIS FUNCTION CALCULATES NET FLOWS INTO POINT I C PARAMETER ( MM=140 , M S=e , MC=9 . MFC=2 . MTP=2 * MFL= 25 , MB=50

)

COMMON / CORR/C 1 (M M , MC ) , C2 (MM , MC ) , CO 1 < MM , M PO ) , C02<MM,MP0

)

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

,

1 FC< MM ,MC) ,PZ< MM, MC) , PC ( MM, MPC ) , CO < MM.MPO ) ,F< MM ) , PFO< MFL * MPO )

2 FF ( MM )»FO(MM»M PO ) * C S ( M S > , PS ( MFL ) ,NSi < MSI ,NS2(MS) * 3 F SS (MS ) »N»NS , NFO , I CON V , E , I EUG , AI ( MM, MC ) , AO( MM , MPO ) , TI T SH (MS ,5)

4 Nh, M MFL) , IFLOOR ( MM) , T( MTP.MFL) ,NFS1 < MS) ,NFS2(MS ),IT(MB),NTP

5 , NCC< MM ), JOC( MM, MPC) ,TOUT

COUELE PRECISION P,PO,PS,PI

NN=NC ( I)

SUM=C.

IF ( N M .EG. 0 ) GO TO 3

CO 1 J J= 1 , NN

J= JC ( I , J J)

CC = C1( 1. JJ)

IF ( F I .LT. P( J) )CC=C2< I.JJ)

PZZ=FZ (I ,JJ)

IF (I . C-T. N)PZZ=0.

FC ( I , J J )=FLOW(PI,P( J ) .PZZ.CC)

1 SUM=SUM+FC( I , J J

)

2 NNC=NCC ( I

)

IF ( NNO .EQ. 0 ) GO TO 4

CO 2 K=l.NNO

CC = CO 1(1 , K

)

IF ( F I .LT. PO ( I «K ) )CC=CG2 ( I , K ) FO ( I ,K )=FLOW( PI ,PO(I,K) ,0,CC)

2 SUM= SUM + FO ( I , K

)

4 FFLOH=SUM+FF( I

)

IF ( I .LE. N ) F ( I )= SUM+FF ( I ) RETURN

ENC

aFDG.P FUNCTION SFLOW.L.l

61
ft

n

o

-

<n

ft

ft

ft

ft

ft

FUNCTION SFLCW

SNBS + PLI E$ .SNOW A.SFLUW

FUNCTION SFLOW ( IS * PI

)

< THIS RCUTINE CALCULATES NET FLOW INTO A SHAFT AND

SI-AFT PRESSURE PROFILE

PA R AN E TER ( M M = 1 4 0 , NS-=8 * MC=9 * MPC=2 , MTP=2 * MFL=25 , MB=50 ) COMMCN NT, P( MM ) ,C( MM ,MC) ,NC(MM) , JC< MM, MC> . ITS(MS),

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

2 FF < NM ) ,FO(MM,MFO) ,CS( NS) , PS (MFL) ,NS1 (MS) ,NS2(MS)

,

3 FSS (MS > ,N,NS ,NPO. ICCNV , E , I BUG . AI ( MM , MC) , AO(MM.MPO) , T ITSH(MS ,5 ) , A NF , F< NFL) , IFLCOR ( NM) , T ( MTP ,MFL ) ,NF SI ( MS) ,NFS2( MS) , IT(MB) ,NTP

5 ,NCC ( NM ), JGC< NN, MFC) .TOUT

COUELE PRECISION P.PO.PS.PI

IF ( I £ UG . GT. 1 )WRITE(6.800)IS

SUM=C.

N1 =NS1 ( IS)

N2 = N S2 ( I S

)

FS( I ) = PI FUF = C . CSS=CS( IS)

CO 1C I=N1.N2

1 1 = I 4 1 —N 1 FLG=FFLOW< I,PS( II ) ) FUP = FLC4-FUP

SUM=«UM+FLO

IF ( I .EQ. N2)GC TO 5

IIP1=I I+l

SI GN = 1 IF ( FLP .GT. 0 . ) SI GN =— i

.

FS ( I IP 1 )=PS( 1 1 )-PZ ( I, 1 ) +SIGN*FUP*FUP/( CSS*CSS

)

IFCIEUG .GT. 1 ) WR I TE( 6,801)1,11 ,PS( II) ,FLO, FUP* SUM

0 CONTINUE

FS S ( IS) = SUM

SFLO\»=SUM

RETURN

FCRNAT STATEMENTS

600 FORMAT(///5X,17HFLCW - SHAFT NO ,15/)

601 FORM AT (5X.3HI =,I3.5X,4HII =,I3,5X.4HPS =.

4- E15.7.5X, 5HFL0 = , E 1 0 .

4

, 5X , 5HFUP = , El 0 .4 , 5X, 5HSUM =,E10.4/)

END

aERKPT PRINTS

62
FEDERAL INFORMATION PROCESSING STANDARD SOFTWARE SUMMARY

01. Summary date

Yr. Mo. Day

04. Software date

Yr. Mo. Day

2 0 5 2 8

08. Software type 09. Processing mode 10. Application area General Specific Automated Data Computer Systems Management/

| | System | | Interactive n Support/Utility Business Smoke control in

PH Computer Program jx] Batch pq Scientific/Engineering [~ Process Control buildings

| | Subroutine/Module | | Combination | | Bibliographic/Textual Q

Other

11. Submitting organization and address 12. Technical contact(s) and phone

02. Summary prepared by (Name and Phone)

John Klote 921-3387

05. Software title A Computer Program for Analysis of Smoke Control

Systems

06. Short title

03. Summary action New Replacement Deletion Q

Previous Internal Software ID

07. Internal Software ID

Suppression & Extinguishment Group

Center for Fire Research

National Bureau of Standards

Washington, D.C. 2023*+

John Klote

921-3387

13. Narrative Pressurized stairwells and pressurized elevators can he 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.

14.. Keywords Air movement; computer programs; egress; elevator shafts; escape means; modeling;

pressurization; simulating; smoke control; stairwells

15. Computer manuf’r and model Univac 1100/82

16. Computer operating system 17. Programing language(s)

ANSI FORTRAN

18. Number of source program state- ments 1756

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 G

26. FOR SUBMITTING ORGANIZATION USE

In-house only Available Inadequate In-house only

Standard Form 185

1974 July

U.S. Dept, of Commerce—NBS

(FIPS. Pub. 30)

185-101
INSTRUCTIONS

01. Summary Date. Enter date summary prepared. Use Year, Month, Day format: YYMMDD.

02. Summary Prepared By. Enter name and phone number (including area code) of individual who prepared this summary.

03. Summary Action. Mark the appropriate box for new summary, replacement summary or deletion of summary. If this software summary is a replace- ment, enter under “Previous Internal Software ID” the internal software identification as reported in item 07 of the original summary, and enter the new internal software identification in item 07 of this form; complete all other items as for a new summary. If a software summary is to be deleted, enter under “Previous Internal Software ID" the internal software identification as reported in item 07 of the original summary; complete only items

01, 02, 03 and 11 on this form.

04. Software Date. Enter date software was completed or last updated. Use Year, Month, Day format: YYMMDD.

05. Software Title. Make title as descriptive as possible.

06. Short Title. (Optional) Enter commonly used abbreviation or acronym which identifies the software.

07. Internal Software ID. Enter a unique identification number or code.

08. Software Type. Mark the appropriate box for an Automated Data System (set of computer programs), Computer Program, or Subroutine/Module,

whichever best describes the software.

09. Processing Mode. Mark the appropriate box for an Interactive, Batch, or Combination mode, whichever best describes the software.

10. Application Area. General: Mark the appropriate box which best describes the general area of application from among: Computer Systems Support/Utility Process Control Management/Business Bibliographic/Textual

Scientific/Engineering Other

Specific: Specify the sub-area of application; e.g.: “COBOL optimizer” if the general area is “Computer Systems Support/Utility”; “Payroll” if the general area is "Management/Business”; etc. Elaborate here if the general area is “Other."

11. Submitting Organization and Address. Identify the organization responsible for the software as completely as possible, to the Branch or Division level, but including Agency, Department (Bureau/Administration), Service, Corporation, Commission, or Council. Fill in complete mailing address, including mail code, street address, city, state, and ZIP code.

12. Technical Contact(s) and Phone: Enter person(s) or office(s) to be contacted for technical information on subject matter and/or operational aspects

of software. Include telephone area code. Provide organization name and mailing address, if different from that in item 11.

13. Narrative. Describe concisely the problem addressed and methods of solution. Include significant factors such as special operating system modifications, security concerns, relationships to other software, input and output media, virtual memory requirements, and unique hardware features.

Cite references, if appropriate.

14. Keywords. List significant words or phrases which reflect the functions, applications and features of the software. Separate entries with semicolons.

15. Computer Manufacturer and Model. Identify mainframe computer(s) on which software is operational.

16. Computer Operating System. Enter name, number, and release under which software is operating. Identify enhancements in the Narrative (item 13).

17. Programing Language(s). Identify the language(s) in which the software is written, including version; e.g., ANSI COBOL, FORTRAN V, SIMSCRIPT 11.5, SLEUTH II.

18. Number of Source Program Statements. Include statements in this software, separate macros, called subroutines, etc.

19. Computer Memory Requirements. Enter minimum internal memory necessary to execute software, exclusive of memory required for the operating system. Specify words, bytes, characters, etc., and number of bits per unit. Identify virtual memory requirements in the Narrative (item 13).

20. Tape Drives. Identify number needed to operate software. Specify, if critical, manufacturer, model, tracks, recording density, etc.

21. Disk/Drum Units. Identify number and size (in same units as "Memory”—item 19) needed to operate software. Specify, if critical, manufacturer, model, etc.

22. Terminals. Identify number of terminals required. Specify, if critical, type, speed, character set, screen/line size, etc.

23 Other Operational Requirements. Identify peripheral devices, support software, or related equipment not indicated above, e.g., optical character devices,

facsimile, computer-output microfilm, graphic plotters.

24. Software Availability. Mark the appropriate box which best describes the software availability from among: Available to the Public, Limited Availability

(e.g.: for government use only), and For-ln-house Use Only. If the software is “Available", include a mail or phone contact point, as well as the price and form in which the software is available, if possible.

25. Documentation Availability. Mark the appropriate box which best describes the documentation availability from among: Available to the Public, In- adequate for Distribution, and For In-house Use Only. If documentation is “Available", include a mail or phone contact point, as well as the price and

form in which the documentation is available, if possible. If documentation is presently "Inadequate", show the expected availability date.

26. For Submitting Organization Use. This area is provided for the use of the organization submitting this summary. It may contain any information deemed useful for internal operation.

GPO : 1974 O - fih-l-G55
NBS-114A (rev. 2-ec) U.S. DEPT. OF COMM. 1. PUBLICATION OR

REPORT NO.

2. Performing Organ. Report No. 3. Publication Date BIBLIOGRAPHIC DATA

SHEET (See instructions) NBSIR 82-2512 June 1982

4. TITLE AND SUBTITLE

A COMPUTER PROGRAM FOR ANALYSIS OF SMOKE CONTROL SYSTEMS

5. AUTHOR(S)

John H. Klote

6. PERFORMING ORGANIZATION (If joint or other than NBS, see instructions) 7. Contract/Grant No. NATIONAL BUREAU OF STANDARDS

DEPARTMENT OF COMMERCE 8. Type of Report & Period Covered

WASHINGTON, D.C. 20234 Final Report

9. SPONSORING ORGANIZATION NAME AND COMPLETE ADDRESS (Street. City. State, ZIP)

Department of Health and Human Services

Washington, D.C. 20201

10. SUPPLEMENTARY NOTES

J Document describes a computer program; SF-185, FIPS Software Summary, is attached.

11. ABSTRACT (A 200-word or less factual summary of most significant information. If document includes a significant

bibliography or literature survey, mention it here)

This paper describes a computer program developed to analyze systems

intended to control smoke in building fires. These systems include pressurized

stairwells, pressurized elevator shafts, zone smoke control systems, and

pressurized corridors. This program calculates air flows and differential

pressures throughout a building in which a smoke control system is operating.

The basic assumptions and limitations of the program are also discussed.

The appendices contain a program listing and examples.

12. KEY WORDS (Six to twelve entries; alphabetical order; capitalize only proper names; and separate key words by semicolon s) Air movement; computer programs; egress; elevator shafts; escape means; modeling;

pressurization; simulation; smoke control; stairwells

13. AVAILABILITY

[jj^Unl imited

| | For Official Distribution. Do Not Release to NTIS

f | Order From Superintendent of Documents, U.S. Government Printing Office, Washington, D.C.

20402.

[ST] Order From National Technical Information Service (NTIS), Springfield, VA. 22161

14. NO. OF

PRINTED PAGES

69

15. Price

USCOMM-DC 6043-P80
t