^-) rl)

.^,y

■J

t I

STATION BULLETIN 5U>%

p:^March, 1979

^'S'f

W^

> ^'*^7^^**xr^tSr7^F^^ ■-

NVENT: FOREST INVENTORY SYSTEM
FOR NEW HAMPSHIRE LANDOWNERS

>fEW HAMPSHIRE

AGRICULTURAL EXPERIMENT STATION

UNIVERSITY OF NEW HAMPSHIRE

DURHAM, NEW HAMPSHIRE

ACKNOWLEDGEMENTS

This publication is a result of the research program of the Institute
of Natural and Environmental Resources. The Institute is a multi-disciplinary
group of scientists involved in a coordinated program of research, teaching
and extension. The research effort encompasses investigations of: problems
affecting the quality of the environment, economics of agriculture, forest
and wildlife resources, the efficient use and conservation of water and soil,
and the regional and community planning and development .

Development of the forest inventory program INVENT was supported by the
Pinchot Consortium for Environmental Studies, the Mclntyre-Stennis 19 project
entitled "Social and Economic Opportunities of Town Forests in New Hampshire",
and the Hatch 149 project entitled "Effects of Site and Stocking on the Growth
of Eastern White Pine".

The authors would like to express their appreciation to the following
people whose help was instrumental in making INVENT a viable system: Jean Brien,
Henry Corrow, Gibb Dodge, Jeff Gove, Roger Langevin, Bill Leak, Kevin McCarthy,
Ken Norcott, Peter Pohl, Don Quigley and Brad Wyman. Special thanks goes to
Phyllis Groves for laying out and typing the final manuscript.

The program was written in FORTRAN-10 on the DECsystem-10 Model 1090
Computer at the University of New Hampshire. A card deck of the program may
be obtained by contacting:

Dr. James P. Barrett

Institute of Natural and Environmental Resources

University of New Hampshire

Durham, New Hampshire 03824

PKogfuxmi o{i thz Mcu> Hamp6kln.e. AgfilcuZtuAxti Expexlrmnt Station oAz ope.n to
alt peA6on6 iMithout izgafid to fuict, colon, nationcit oKlgtyi ok izx. Thz
UnlveAilty 0^ Mew) Hamp^kOiz -Li an A^^-iAmattve. Actyion/EquaZ OppohtixyUty
Employed.

ABSTRACT

A computerized forest inventory system has been developed to meet the
needs of New Hampshire landovmers and foresters. The system will analyse
all sampling designs currently in wide use throughout the State producing
tables easily understood by landowners. The system includes a rigorous
statistical analysis of volume standard errors as well as a newly developed
height double sampling technique. The system is available through the
Institute of Natural and Environmental Resources.

KEY WORDS: Forest inventory. Computer inventory. Sampling. Double sampling.

TABLE OF CONTENTS

Page

Introduction

Section 1. Field Forester Information 3

Section 2. Data Coding and Program Execution 17

Section 3. Program Documentation 25

Section 4. Statistical Analysis 37

Appendix A Data Coding Information ^3

Appendix B Statistical Formulas and Derivations 51

Appendix C Volume Formulas 55

Appendix D Inventory Design and Sample Size Determination . . 59

Appendix E Program Listing 65

Literature Cited

86

ii

INVENT:
FOREST INVENTORY SYSTEM FOR NEW HAMPSHIRE LANDOWNERS

by
D. S. Linden and J. P. Barrett*

INTRODUCTION
INVENT Is a computerized forest inventory processing system designed
for New Hampshire landowners. INVENT was designed with the following five
goals in mind:

1. Produce a system that will analyse all sampling designs and methods
presently in wide use throughout the State.

2. Produce output tables easily understood by foresters and landowners —
output tables that look as if they were designed by a forester, not

a computer.

3. Find and implement volume equations that are well adapted to the
New Hampshire species — calculate both board-foot and cubic-foot
volumes .

A. Indicate the reliability of volume estimates.

5. Produce a system that any forester can use and understand — a system
available to all foresters at low cost.

The realization of these goals is the result of two years work in close coopera-
tion with consulting foresters, county foresters, and forest industry. As a
result. New Hampshire foresters have "been emancipated from the drudgery and
limitations of computation and can concentrate on forest data needs." (Davis
1966, p. 282).

The use of INVENT is made available by the Institute of Natural and
Environmental Resources at the University of New Hampshire. Institute personnel
will process data from tally sheets submitted by subscribing foresters. Persons
interested in using the system should contact Professor James Barrett at I.N.E.R.,
University of New Hampshire.

This publication is written in four separate sections. Section 1 contains
all the information a field forester needs to know to utilize the service.
Section 2 describes how the data is computer coded and how the program is run.
Section 3 documents the computer program itself with subroutine descriptions

*Graduate Research Assistant and Professor of Forest Resources, Institute of
Natural and Environmental Resources, respectively.

-2-

and variable dictionaries. Section 4 discusses the statistical techniques
used within the program. The appendices include tally forms, species code
list, volume equations, a guide to sample size determination, a program
listing and other supportive material.

-3-

SECTION 1. FIELD FORESTER INFORMATION

The first step in a forest inventory is to systematically arrange the
forest into groups or categories based upon some definite scheme. We shall
refer to these groups as either compartments or strata, using these terms
interchangeably. The boundaries of these compartments may be determined by
forest type, stocking level, merchantability class or administrative conveni-
ence. Whatever the criteria actually used, INVENT will assume that each
compartment should be processed individually and will produce a complete set
of summary tables for each compartment. Furthermore, INVENT will produce a
set of summary tables for all compartments combined assuming stratified
random sampling (Cf. Freese 1962). INVENT will handle up to 999 compartments
although in practice one would seldom have more than 20 compartments within
any forest.

For each forest to be processed, the forester must supply the informa-
tion requested on the INVENT Processing Request Form (Appendix A, page 48),
This information includes :

1. Name and address of forester making request.

2. Number of compartments in the forest.

3. Alternate form class specification.

4. Diameter class size.

5. Confidence level.

6. Tract description-
Alternate^ Form Class Specification. If heights are measured to the

merchantable top, form class equations are used to determine board foot
volume. INVENT uses the following form classes as a default:

Hemlock 70 All other softwoods 76 All hardwoods 74
Alternate form classes may be specified using the space provided on the INVENT
Processing Request Form. These form classes will be used for all compartments
in the forest.

-4-

Dlameter Class Size. INVENT assumes that trees will be grouped Into
diameter classes. The size of the class Is up to the user within the follow-
ing limitations. The class size must be a whole number. The smallest diameter
class allowable is the class size; i.e., if the class size is 2 in. then the
first class must be 2 inches. All other diameter classes must be integer
DHjltiples of the class size. The maximum diameter class must be less than or
equal to 40 inches.

Confidence Level. INVENT uses the 95 percent confidence level as a
default. Using this level, there is only a one in twenty chance that the
true volume will be outside the confidence interval computed by INVENT (Cf.
Freese 1967, p. 11). This is the most common confidence level used in
forestry in the Northeast. If the user wishes to use a different level, he
may specify that level as long as it is in the range from 65 percent to
99 percent.

Tract Description. A tract description of up to 70 characters should be
specified for the forest. This description will be printed at the top of
each page of the combined compartment summaries.

For each compartment within the forest, the forester must supply the
information requested on the INVENT Compartment Information Form (Appendix A,
page 49). This information includes:

1. Compartment name and/or description.

2. Sampling method and sample unit size.

3. Tree height measurement units.

4. Tree top specification.

5. Whether or not the heights of all volume trees were
measured.

6. Size of compartment in acres.

Compartment Name. A compartment name and/or description should be
specified for each compartment. The description may be up to 70 characters
in length and is printed at the top of each page of the compartment output
tables.

-5-

Sampling Method. INVENT will process point cruises (variable plot) ,
plot cruises (fixed radius plot), strip cruises, and 100 percent tallies.
The sampling method need not be the same for each compartment. When using
point or plot sampling, a complete statistical analysis of the reliability
of the volume estimates will be calculated only if the data for each point
or plot is tallied separately. If only summary data is available, no such
statistical analysis can be made. When point sampling, the BAF of the angle
used must be specified. When plot sampling, the plot size in acres must be
specified. Strip cruises are analysed as a plot cruise where only one large
plot is sampled. In this case, the plot size is the total acreage of all
strips. There is no statistical analysis of reliability on strip cruises.
One hundred percent tallies require no statistical analysis of reliability
since all trees in the population have been measured. One hundred percent
tallies are also treated as a plot cruise where only one large plot is sampled.
Here, however, the plot size is equal to the total acreage of the compartment.

Tree Height Measurement. Tree height may be measured in feet or sixteen
foot logs to the nearest half -log.

Tree Top Specification. Trees may be measured to a merchantable top or
measured by total height. There is also a mixed option where softwoods are
measured by total height and hardwoods by merchantable height. As mentioned
earlier, form class volume equations for board foot volume are used when
merchantable height is specified.

Were Heights of All Volume Trees Measured? INVENT allows double sampling
for height (Cf. Freese 1962, p. 43). Using this technique, the forester only
measures height on a subsample of the trees. INVENT develops a height-dbh
relationship from the height trees sampled and uses this relationship to
estimate height for the remaining trees. This method must be used with some
caution. On at least one-third of the trees in every species, the height
should be measured. At least three height trees must be measured for every
species within each compartment. An analysis table of the height-dbh relation-
ship is printed for each compartment using the option. A warning message is

-6-

printed if an insufficient niimber of trees within any species was sampled.
(See table 17).

Size of Compartment . The size of each compartment in acres must be
specified. If not known exactly it must at least be estimated.

Tally Forms and Product Specification

Tally forms are available for use with INVENT (see Appendix A, page 46) .
Trees are tallied using a two letter species code (Appendix A, page 45),
diameter class, height and product. The product column is used as follows:
INVENT allows trees to be graded as sawlog or pulpwood. If you do not wish
to grade trees simply ignore the product column and all trees will be assumed
sawlog. If you do wish to grade trees there are two methods available. You
can grade the entire tree as sawlog or pulpwood by entering an S in the product
column for sawlog, or entering a P for pulpwood. Alternatively, you may spec-
ify the percent sawlog in the tree. If you put 60 in the product column, 60
percent of the volume will go into sawlog and 40 percent will go into pulpwood.

The column marked # is used if two or more trees of the same species, size
and product are tallied on a given plot.

Output Tables
Seventeen tables are presented as an example of INVENT 's output. Tables
1-16 are part of a 59 page inventory report based on a 3 compartment cruise
of Whitaker Woods in North Conway, N.H. Tables 1-9 are for compartment 1
while tables 10-16 are for all 3 compartments combined. Table 17 is taken
from another cruise where the height double sampling option was used. This
table is a summary of that relationship. The output tables are further dis-
cussed in section 3 on page 34. However, the tables should be self-explanatory
to most foresters.

-7-

Table 1. Job Control Information (Individual Stratum)

INVENT VER.

DiTE OF EON:

11/1/78 I.N.E.R.

7-Dec-78

O.N.H.

B.A.F. = 10.00

DIAHETER CLASS SIZE = 2

HEIGBTS OF ALL VOLOHE TREES HEASOBED

INPUT FILE NAME: CONilT.CI

CONFIDENCE LEVEL = .95

NOBBEE OF POINTS SAHPLED = 20

HEIGHTS HEASORED BY HERCH HEIGHT IN LOGS

STSATOB I 1

ACREAGE OF TRACT = a7.0

Table 2. Species Stand and Stock Table

WHITE PINE FC = 76 *

*
**************************************************************************************************

BASAL AREA

TREES

SAHLOG

SAHLOG

POLPHOOD

DIAHETER

PER

PER

COBIC FOOT VOLDME

BOARD FOOT VOLOHE

CUBIC FOOT VOLUHE

CLASS

ACRE

ACRE

PER ACRE

PER ACRE

PER ACRE

8

1.5

4.3

0.0

0.0

30.7

10

1.0

1.8

16.8

82.8

0.0

12

«.5

5.7

91.1

535.4

0.0

14

5.0

4.7

109.9

705.8

0.0

16

7.5

5.4

168.7

: 1115.0

16.4

18

9.0

5.1

221.1

1545.0

9.6

20

7.0

3.2

196.5

1405.3

0.0

22

3.5

1.3

101.6

747.0

0.0

2H

2.5

0.8

70.3

: 538.0

0.0

26

a. 5

1.2

134.1

1035.6

0.0

28

1.0

0.2

33.7

258.3

0.0

30

1.0

0.2

28.4

233. 1

0.0

32

2.0

0.4

58.4

483.9

0.0

TOTAL

50.0

34.3

1230.6 ♦

35)1

8685.2 i 35X

56.7 1 83X

HEAN STAND DIAHETER = 16.3

BOARD FOOT COEFFICIENT OF VARIATION = 78X

MERCHANTABLE H.S.D. = 16.3

POLP-HOOD COEFFICIENT OF VARIATION = 183%

PERCENT CRUISE = 6.2H

-8-

Table 3, Species Stand and Stock Table

*****•••*•**•***•••••••*•*•••***••****••**••***•••••••••••**•«•••*••***•••*•••**••*********••*•****

•

BED

OAK

FC = 74 »

•

•

*•*********•*•**•****•*****************•«••*•••*•**•**••••*«•*•••*••**•**»********•**•*•******•**•*

: EASiL AfiEA

TBEES

SAWLOG

SAHLOG

PULPBOOD

DliHETEB

PES

PEB

CUBIC FOOT yOLDHE

BOABD FOOT VOLOHE

CUBIC FOOT VOLOHE

CL&SS

ACRE

ACSE

PEB ACBE

PEB ACBE

PEB ACBE

6

1.0

5.1

0-0

0.0

19.4

8

1.5

14.3

6.1

26.3

16.7

10

1.5

8.3

5U. 1

: 273.4

13.2

12

5.5

7.0

72.11

417.1

0.0

14

2.0

1.9

30.9

183.9

0.0

16

2.0

l.t

28.2

179.2

0.0

18

3.0

1.7

U1.8

: 277.1

0.0

20

1.5

0.7

25.1

165.4

0.0

22

0.5

0.2

U.6

: 35.6

0.0

TOTAL

21.5

30.5

263.3 *

58X

1557.9 1 56JI : 49.3 ♦151X

HEAN STAND DIADETEB = 11.4

BCAED FOOT COEFPICIEMT OF VABIATIOH = 122*

HEBCHAHTABLE H.S.D. = 12.2

PDLP-BOOD COEFFICIENT OF TABIATION = 327%

PEBCENT CEOISE = 3. OX

Table ^. Softwood Summary Stand and Stock Table

« AIL SCFTBOODS

*

BASAL ABEA

TBEES

SA8L0G

SAMLOG

POLPHOOD

DIABETEB

PEE

PEB

COBIC FOOT VOLUHE

BOABD FOOT VOLOHE

COBIC FOOT VOLOBE

CLASS

ACBE

ACBE

PEB ACBE

PEB ACRE

PER ACRE

6

0.5

2.5

0.0

0.0

8.1

8

3.5

10.0

0.0

0.0

57.3

10

7.0

12.8

86.1

392.6

4.0

12

10.5

13.4

169.7

9 56.7

0.0

14

9.0

8.4

183.7

1147.4

0.0

16

11.5

8.2

238.4

1574.8

16.4

18

13.5

7.6

307.3

2079.7

9.6

20

8.5

3.9

226.6

1593.2

0.0

22

4.5

: 1.7

125.5

915.2

0.0

24

2.5

0.8

70.3

538.0

0.0

26

: 4.5

1.2

134.1

1035.6

0.0

28

: 1.0

0.2

33.7

258.3

0.0

30

1.0

: 0.2

28.4

233.1

0.0

32

2.0

0.4

58.4

483.9

0.0

TOTAL

: 79.5

: 71.5

1662.2 1

29X

11208. 3 ♦ 30X

95.4 ♦ 63%

SEAN STAND DIAHETEB = 14.3

BOARD FOOT COEFFICIENT OF VARIATION = 65X

BEBCBANTABLE B.S.D. = 14.5

PDLP-BOOD COEFFICIENT OP VARIATION = 138X

PEBCENT CROISE = 4.7X

-9-

Table 5, Hardwood Summary Stand and Stock Table

• AIL BIBDROODS •

: BASAL AREA

TSEES

SAflLOG

SAHLOG

POLPiOOD

OI&HETEB

PER

PES

COBIC FOOT VOLOHE

BOARD FOOT VOLOHE

COBIC FOOT TOIDBE

CLASS

ACRE

ACRE

PER ACRE

PEB ACRE

PER ACBE

6

3.0

15.3

0.0

0.0

48.6

8

3.5

10.0

6.1

26.3

47.8

10

6.0

11.0

60.0

302.8

25.2

12

7.5

9.5

77.3

448. 5

14.4

^n

3.0

2.8

1*3.8

260.4

0.0

16

2.5

1.6

35.8

224.4

0.0

18

3.0

1.7

41.8

277. 1

0.0

20

1.5

0.7

25.1

165.4

0.0

22

0.5

0.2

4.6

35.6

0.0

TOTAL

30.5

53.0

29U.7 1 52%

1740.5 + 5U

136.0 1 72*

HEAH STAHD DTAHETEB = 10.3

BOARD FOOT COEFFICIEBT OF TABIATIOB = 109X

HEBCBAHTABLE H.S.D. = 11.6

PDLP-BOOD COEFFICIEHT OF VABIATIOM = 157X

PEBCEMT CBOISE = 2. 4 J

Table 6. All Species Summary Stand and Stock Table

ALL SPECIES

BASAL AREA

DIAHETEB

PEE

CLASS

ACRE

6

3.5

8

7.0

10

13.0

12

18.0

14

12.0

16

14.0

18

16.5

20

10.0

22

5.0

24

: 2.5

26

: 4.5

28

1.0

30

: 1.0

32

: 2.0

TBEES

FEB
ACRE

SARLOG

COBIC FOOT VOLOHE

PER ACRE

SANLOG

BOARD FOOT VOLOHE

PEB ACRE

POLPROOD

COBIC FOOT VOLOHE

PEB ACBE

17.8

0.0

20.1

6.1

23.8

146.1

22.9

246.9

11.2

227.6

10.0

274.3

9.3

349.1

4.6

251.7

1.9

130.1

0.8

70.3

1.2

134.1

0.2

33.7

0.2

28.4

0.4

58.4

0.0

26.3

695.4

1405.2

1407.8

1799.2

2356.7

1758.6

950.8

538.0

1035.6

258.3

233.1

483.9

56.7
105.1

29.2

14.4
0.0

16.4
9.6
0.0
0.0
0.0
0.0
0.0
0.0
0.0

TOTAL

110.0

124.5

1956.8 + 21)1

12948.9 * 22X

231.4 * 50%

HEAN STAHD DIAHETER = 12.7

BOABD FOOT COEFFICIEHT OF VABIATIOH = 48X

HEBCBAHTABLE H.S.D. = 13.5

POLP-ROOD COEFFICIEHT OF VASIATIOH = 109X

PERCEHT CRDISE = 3. 8X

-10-

Table 7, Species Composition Table

• *

* SPECIES COHPOSITIOH BT PEBCEHT *

S&WLOG

SAWLOG

PULPWOOD

COBIC FOOT

BOARD FOOT

: COBIC FOOT

SPECIES

BASAL ABEA

as. 45

TREES

VOL ORE

VOLOHE
67.07

VOLOHE

UHITE PINE

27.59

62.89

24.50

RED PINE

5.91

5.15

6.90

6.76

0.00

HEMLOCK

20.91

211.67

15.15

12.72

16.75

SDGAB HiPLE

0.91

3.20

0.00

0.00

6.82

RED BAPLE

0.91

1.53

0.26

0.25

3.68

WHITE BIBCH

0.91

2,56

0.25

: 0.24

3.92

BEECH

5.U5

10.80

1.09

0.91

23.03

RED O&K

19.55

2«.51

13.46

12.03

21.31

ilL SOFTBOODS

72.27

57. U1

61.94

86.56

41.25

ALL HARDWOODS

27.73

U2.59

: 15.06

13.44

58.75

_„_—_«

_________«___«.

Table 8. Volume Summary for All Species

***************************************************************************************************

* *

* VOLOHE TOTALS FOR ALL SPECIES ♦

* *

SAWLOG

COBIC FOOT VOLOHE

PER ACRE

SAWLOG :

BOARD FOOT VOLOHE :

PEE ACRE :

8685.

2 :

876

0 :

1647.

2 :

0

0 :

32.

8 :

31.

4 :

118.

4 :

1557.

9 ;

PDLPWOOD

COBIC FOOT VOLOHE

PEE ACRE

WHITE PINE
RED FINE
H EH LOCK
SOGAE HAPLE
RED HAPLE
WHITE BIRCH
BEECH
RED OAK

1230.6

135.0

296.5

0.0

5.1

4.8

21.4

263.3

56.7

0.0

38.8

15.8

8.5

9.1

53.3

49.3

ALL SOFTWOODS
ALL EABDWOODS

1662.2
294.7

11208.3
1740.5

95.4
136.0

ALL SPECIES

1956.8

12948.9

-11-

Table 9. Volume Summary for All Species Expanded by Acreage

t^t ********************** *************************************************************************

*

VOLUHE TOTALS EXPSNDED BY iCBEiGE *

*

ttttt ******************************** *************************************************************

S&RLOG
CUBIC FOOT VOLDBE

SKHLOG
BOARD FOOT VOLOBE

POLPHOOD
CUBIC FOOT TOLOSE

WHITE PIHE
BED PINE
HEMLOCK
SDGAB SAPLE
BED HAPLE
WHITE BIBCH
BEECH
BED OAK

57838

63a7

13936

0

2«0

228

looa

12377

408202

U1170

77al9

0

1542

1478

5563

73222

2664
0

1822
742
400
426

2505

2317

ALL SOFTWOODS
ALL HARDWOODS

78122
13849

526792
81805

4486
6390

ALL SPECIES

91971

608597

10876

Table 10. Job Control Information (Stratified Total)

IHVENT VER. 3 11/1/78 I.H.E.E. O.N.H.

DATE OF BON: 7-Dec-78

DIAHETEE CLASS SIZE = 2

TOTAL ACEEAGE SAMPLED = 105.0

C08FIDEBCE LEVEL = .95

NOHBEE OF STRATA SAMPLED = 3
TOTAL POINTS SAMPLED = 56

-12-

Table 11, Species Volume Table (Stratified Total)

* »

* WHITE PINE PC = 76 »

* *

DIAHETES
CL&SS

SiHLOG
CDBIC FOOT VOLOBE

SAHLOG
BOARD FOOT VOLOHE

POLPWOOD
COBIC FOOT VOLBHE

6

8
10
12

ia

16
18
20
22
2"t
26
28
30
32
3U
36

0

823

1«99

5111

5550

9600

1«5ai

141)68

10177

7902

92117

2506

2074

4228

568

455

0

2718

7465

30058

35604

63484

100826

103635

75100

60263

71309

19544

16148

34874

4555

3920

914
1716
0
0
0
769
453
0
0
0
0
0
0
0
0
0

TOTAL

88750 + 26X

629503 ♦ 26X

3853 ♦ 73*

Table 12. Species Volume Table (Stratified Total)

BED OAK FC = 74

DIAHETER
CLASS

SAHLOG
CDBIC FOOT VOLOHE

SAHLOG
BOARD FOOT VOLOHE

POLPHOOD
COBIC FOOT VOLOHE

6
8
10
12
14
16
18
20
22
24
28

0

1234

23270

42632

16995

29345

19598

9308

1674

0

2035

1635

2981

622

310

0

0
664
294

0
291

0

24968 ♦ 37X

146091 ♦ 36X

6799 ♦ 61X

-13-

Table 13. Softwood Summary Volume Table (Stratified Total)

• *

* AIL SCFTHOODS •

DI4METEE
CLASS

SAHLOG
CDBIC FOOT VOLDIJE

SAILOG
BOAED FOOT VOLOHE

PULPHOOD
CDBIC FOOT VOLOBE

6
8
10
12
14
16
18
20
22
24
26
28
30
32
34
36
40

0

1053

6105

10944

12957

22644

26139

18292

13103

9151

10723

2799

3291

4407

915

45S

459

0

3451

28274

61345

79177

141792

172099

127282

95876

68560

81395

21620

24339

36297

7039

39 20

3204

1874

4161

320

107

0

769

591

0

0

0

0

0

0

0

0

0

0

TOTAL

143437 + 19X

955670 ♦ 19X

7822 + 50X

Table 14. Hardwood Summary Volume Table (Stratified Total)

AIL EAROHOODS
*************************************************************************************************

DIASETEE
CLASS

SAHLOG
COBIC FOOT VOLOHE

SAWLOG
BOARD FOOT VOLOHE

PDLPWOOD
COBIC FOOT VOLOHE

6
8
10
12
14
16
18
20
22
24
26
28

607

1234

48605

56781

24463

34023

19598

9308

2825

0

2009

2035

7887
11231

2785

2805
826
488

1212
«58
610
607
237
0

35115 + 26X

201490 + 26%

29146 + 27X

-14-

Table 15, All Species Summary Volume Table (Stratified Total)

* ilL SPECIES »

DIAHETEB
CL&SS

SAiLOG
COBIC FOOT VOLDHE

SAWLOS
BOAED FOOT VOLOHE

POLPUOOD
COBIC FOOT VOLOHE

6
8
10
12
14
16
18
20
22
24
26
28
30
32
34
36
40

366

1341

15522

20955

17124

28116

29144

19695

13484

9151

11038

3089

3291

II407

915

455

459

607

4685

76879

118127

103641

175815

191697

136589

98701

68560

83404

23655

24339

36297

7039

3920

3204

9760

15393

3104

2912

826

1257

1803

458

610

607

237

0

0

0

0

0

0

178552 + 14)1

1157160 ♦ 14J

36968 ♦ 22X

Table 16. Volume Summary for All Species (Stratified Total)

» •

» VOLOBE TOTALS FOE ALL SPECIES *

« *

SPECIES

SAHLOG
COBIC FOOT VOLOME

SARLOG
BOARD FOOT VOLOHE

PDLPHOOD
COBIC FOOT VOLOHE

OBITE PINE
BED PINE
HEMLOCK

ED-BL-WH SPEOCE
SOGAB MAPLE
BED HAPLE
OBITE BIBCH
VELLOS BIRCH
BEECH
BED OAK

88750

10094

44593

0

0

1072

4474

89

4513

24968

629503

65651

260516

0

0

6271

25213

649

23266

146091

3853

0

3746

224

742

6150

3255

0

12200

6799

AIL SOFTWOODS
ALL BABDHOODS

143437
35115

955670
201490

7822
29146

ALL SPECIES

178552

1157160

36968

-15-

Table 17. Regression Summary (If Option Used)

JIHILISIS OF THE HEIGHT - DBH RELiTIOHSHIP
aODEL: HEIGHT = BO + B1/0BH

_

BO

.

B1

.

B

STAHOABD

ATEEAGE

VOLOBE

HEIGHT

SPECIES

IHTEECEPT :

SLOPE

:

SQOABED

EBBOB %

HEIGHT

TBEES

TBEES

8HITE PIHE

3.296

-21.321

':"

54.0

5.1

1.5

92

40

HEBLOCK

1.799

-8.528

;

16.2

: 9.6

1.1

48

23

BED BIPLE

3-106

-29.436

:

62.1

13.5

0.9

39

11

BLICK CBERBI

• ««

HABHIliG

***

0.0

0.0

3

0

SHEET BIBCB

• *•

RABHIBG

• *•

: 0.0

0.0

1

0

B»SS»OOD

0.0

2.0

2

2

BED OilK

1.1«0

:

-5.U71

;

30.1

: 19.7

0.8

: 5

4

WHITE OAK

*•*

KABNIBG

•••

0.0

0.0

3

0

Bl&CK OXK

••*

RABHIBG

• ••

0.0

0.0

2

0

ELH

1.157

z

-«.992

;

18.7

: 15.0

0.8

9

6

BICKOBT

2.865

•

-23.076

•

73.0

9,2

1.3

6

3

ilL SOPTBOODS

: a.6

1.4

140

63

ALL HABOROODS

s 8.1

1.0

70

26

ALL SPECIES

4.1

1.3

210

89

••* RABHIBG •••

IHSDFFICIEHT BDHBEB OF HEIGHT TBEES BEASOBED FOB THIS SPECIES!

-17-

SECTION 2, DATA CODING AND PROGRAM EXECUTION

Program Operating Environment

INVENT was developed on the DEC-10 system at the University of New Hamp-
shire. INVENT is an interactive system acquiring all of its job control
related information from the user through a dialogue. Unit 5 (TTY on DEC-10)
is used as the input and output device for all of INVENT 's interactive
dialogue.

For each stratum the plot and tree data must reside on a disk file with
a unique name. The name of the file is specified by the user at execution
time as INVENT processes each stratum. Unit 21 is used to input each disk
file. The file uses a fixed format as described below.

If the user wishes to use alternate form classes, these must reside on
a disk file also. The user specifies the name of this file at execution
time and unit 1 is used to input this file. The file uses a fixed format as
described on page 18.

All output is written to unit 3 (LPT on DEC-10) to a file named INVENT. LPT.

Data File Format and Coding
An individual data file must exist for each stratum. The file consists
of card images coded as follows:

Card Type Columns Format Variable Description

Title 1-70 12A5 TITLE Stratum title

Plot 1-2 12 POINT Plot number (1)

Number of trees on plot (1)

Tree 1-2 A2 SPEC Species code

Diameter at breast height

Product code (2)

Tree height

Number of trees of this

description (3)
Percent sawlog (4)

(1) When using the summary option where trees are not tallied by Individual
plots, only one plot card is used and it is coded as follows: POINT =-9,
TREES = number of plots sampled .

(2) Product codes: Sawlog (S or blank); Pulpwood (P) ; Percent Sawlog (%) .

(3) If IREP is coded as 0 or blank, the program assumes IREP = 1.

(4) PCTSAW is ignored unless PROD = %. PCTSAW is coded as a percent (0-100)
and the F3.2 format reads it as a decimal (0-1.00).

1-2

12

POINT

3-5

13

TREES

1-2

A2

SPEC

4-5

12

DBH

6

Al

PROD

7-9

F3.

0

HEIGHT

10-11

12

IREP

12-14

F3.

2

PCTSAW

-18-

Examples of various tree card codings may be found in Appendix A,
page ^^. The species code list is found in Appendix A, page 45.

The data file is constructed by combining the title, plot, and tree
cards as follows. The first card in the file is the title card. This is
followed by the plot card for plot number 1. The plot card is then followed
by all the tree cards for that plot. Next, the plot card for plot number 2
is included followed by its tree cards; and so on for the remaining plots.
The actual value of POINT (plot number) need not be consecutive.

An example of an INVENT tally sheet and the resulting data file listing may
be found in Appendix A, pages 46 and 47.

Alternate Form Class File Coding
The user has the option to specify a set of alternate form classes.
Form class only affects volume when tree heights are measured to a merchant-
able top (see page 3) . A file containing these alternate form classes must
reside on a disk file. This file consists of four rows and eight columns

organized as follows:

Row 1
Row 2
Row 3
Row 4

Form classes for species 1-8 (see species list in
" " " " 9-16 Appendix A, page 45)

17_24
It II II II 2S— '^2

Each row is formatted 8(12, X).

-19-

As an example, assume the following form classes were desired:

All pines 80

Hemlock, Spruces, Ash, 78

Black Cherry

Balsam Fir, Tamarack, Cedar

Other Softwoods , Sugar Maple 76

White Birch, Yellow Birch

All other species 75

The data file that would implement these form classes is:

80 80 80 78 76 76 78 78

80 76 76 76 75 78 75 75

75 78 76 76 75 75 75 75

75 75 75 75 75 75 75 75

Running the Program

Running INVENT consists of interactively supplying the program all job
control information. This information includes: alternate form class
specification (if desired), diameter class size, number of strata, confidence
level, input file name for each strata, sampling method (point or plot),
basal area factor or plot size, method of height measurement, acreage of
each stratum, and title of combined strata output (if more than one stratum
sampled) .

This interactive portion of the program is self-explanatory. The program
asks the user certain questions while listing the possible answers. The
number of questions asked depends on the options used and the number of strata
sampled. The questions are asked as each data file is read and processed.
Under conditions of heavy computer usage there may be a time lapse of as much
as two minutes between certain questions. This process is best explained by
example.

Tables 1-16 are part of a 59 page inventory report based on a 3 compart-
ment cruise of Whitaker Woods in North Conway, N.H. Listed below is the
entire user-program dialogue that produced that output. Each of the questions
asked has been numbered for reference and user responses are underlined.
Most questions are totally self-explanatory and therefore only a few will be
discussed.

-20-

♦RUN INUENT

1, DO YOU WANT TO READ IN ALTERNATE FORM CLASSES FROM DISK?
NO

2.
3.

6.

7.
8.

9.

10.
11.
12.

ENTER DIAMETER CLASS SIZE
2

ENTER NUMBER OF STRATA
3

INVENT USES THE 95% CONFIDENCE LEVEL AS A DEFAULT.
DO YOU WISH TO USE A DIFFERENT LEVEL?

5 ENTER INPUT FILE NAME FOR STRATUM * 1
CONWIT.Cl

WHITAKER WOODS 1978 COMP. 1

SPECIFY SAMPLING METHOD

TYPE POINT OR PLOT
POINT

ENTER B.A.F.
10

ARE HEIGHTS IN FEET OR SIXTEEN FT. LOGS?

TYPE FEET OR LOGS
LOGS

ARE HEIGHTS TOTAL OR MERCHANTABLE?

TYPE TOTAL f MERCH» OR MIXED
MERCH

WERE THE HEIGHTS OF ALL VOLUME TREES MEASURED?

ENTER ACREAGE OF TRACT.
42

ENTER INPUT FILE NAME FOR STRATUM ♦ 2
C0NWIT.C2

WHITAKER WOODS 1978 COMP. 2

13.

14.
15.

16.

-21-

SPECIFY SAMPLING METHOD

TYPE POINT OR PLOT OR SAME
SAME

ENTER ACREAGE OF TRACT.
42

ENTER INPUT FILE NAME FOR STRATUM * 3
C0NUIT.C3

UHITAKER WOODS 1978 COMP. 3

SPECIFY SAMPLING METHOD

TYPE POINT OR PLOT OR SAME
SAME

17 ENTER ACREAGE OF TRACT,

18 WHAT DO YOU WANT TO TITLE THE COMBINED STRATA OUTPUT?
WHITAKER WOODS 1978 STRATIFIED TOTAL CQMPS. 1 - 3

END OF EXECUTION

CPU time: 4.32 ELAPSED TIME.* 5:31.23

EXIT

Answering NO to questions 1 and 4 indicates the user will use the default
form classes and the default confidence level of 95 percent.

Question 9 asks the user if heights are total or merchantable and offers
three responses: TOTAL, MERCH, or MIXED. The response MIXED would have
indicated that softwoods were measured by total height while hardwoods were
measured by merchantable height. In this case, however, MERCH was answered
indicating merchantable height was used on all trees.

There would be a pause of about thirty seconds between questions 11 and
12 as the program processes stratum 1.

Question 13 asks the user to specify the sampling method as did question
6. However, question 13 offers an additional response: SAME. Answering SAME
indicates that in stratum 2 the sampling method, basal area factor, and method
of height measurement are identical to those in the previous stratum. There-
fore questions 7-10 are not repeated.

-22-

Once again there would be a pause of about thirty seconds between
questions lA and 15 and between questions 16 and 17.

Question 18 asks for a title for the combined strata output. This
title may be up to 70 characters in length and is printed at the top of each
table in the combined strata output section.

After an additional pause of about thirty seconds the end of execution
message appears and the output is sent to the line printer queue.

The following example is Included to Illustrate the alternate form
class option and the alternate confidence level option.

♦RUN INVENT

1.

lA.
iB.

2.

3.

4.

4a.

DO YOU WANT TO READ IN ALTERNATE FORM CLASSES FROM DISK?
YES

ENTER FORM CLASS FILE NAME
UHIT.FRM

DO YOU WANT A LIST OF THE FORM CLASSES JUST READ?
YES

UP =80 RP = 80 PP =80 HM = 78

BF = 76 TA = 76 SP = 78 NS = 78

SC =80 CE = 76 OS =76 SM = 76

RM =75 UA = 78 WI =75 AS = 75

YP = 75 BC = 78 WB = 76 YB = 76

SB = 75 GB = 75 BE = 75 BW = 75

RO = 75 WO = 75 BO =75 EL = 75

(3U = 75 HI = 75 HH = 75 OH = 75

ENTER DIAMETER CLASS SIZE

ENTER NUMBER OF STRATA
3

INVENT USES THE 95% CONFIDENCE LEVEL AS A DEFAULT.

DO YOU WISH TO USE A DIFFERENT LEVEL?

YES

ENTER CONFIDENCE LEVEL AS A DECIMAL IN THE RANGE .65 TO .99

.90

-23-

Answering YES to question 1 envokes the alternate form class option.
Question la asks for the file name while question lb gives the user the
option to list the new form classes. These form classes are the ones de-
scribed on page 19. Answering YES to question 4 allows the user to specify
an alternate confidence level.

-25-

SECTION 3. PROGRAM DOCUMENTATION

The FORTRAN program consists of the main program: INVENT; nine sub-
programs: FRMCLS, LINREG, GUTS, VOLMER, VOLTOT, STRAT , OUTPl, MDSTI, and
MDNRIS; and five named common storage areas: COMUNE, FORMCL, REG, OUTl,
and STOUT. A description of each along with a variable dictionary are
presented below. The program listings are found in Appendix E.

Common Area COMUNE

COMUNE transfers job control information from the main program to the

subroutines.

Description

Acreage of compartment

Prism factor or plot size

Diameter class size

Sampling method

Input file name

Height measurement length technique

Conversion factor to transform logs

to feet
Two-tailed area under t-curve
Regression option indicator
Output title indicating sampling

method used
Compartment title
Height measurement top technique

Variabl(

2 Type D
Real

imension

ACRES

BAF

Real

-

CLASIZ

Int

-

FRAME

Real(C*5)l

-

IFILE

D.P.(C*10)

-

LOGANS

Int (C*4)

-

MULT

Int

-

PROBLV

Real

_

REGANS

Int (c*3)

-

SAMMTD

D.P.(C*8)

-

TITLE

Int (C*5)

14

TOPANS

Int (C*5

-

Common Area FORMCL
FORMCL transfers the form class array.
Variable Type Dimension Description

FC Int 32 Form class array

Common Area REG
REG transfers information to and from the regression subroutine.

Description

Average height

Minimum diameter sampled for height
Number of height trees sampled
Regression slope
Regression Intercept
Coefficient of determination
Standard error of the regression
Species code

Number of non-cull trees above
minimum d.b.h. sampled

■(C*5) indicates the variable is used to store 5 characters.

Variable

I Type

Dimension

MEANY

Real

35

MINDBH

Int

32

NTR

Int

35

REGCOF

Real

33

REGINT

Real

33

RSQR

Real

33

SEREG

Real

35

SPECIE

Int (C*2)

32

TOTNTR

Int

35

-26-

Common Area OUTl

OUTl transfers processed data to the output subroutine.

Variable

. Type

Dimension

BASPAC

Real

40,35

BASTOT

Real

35

BDTOT

Real

35

BDVLSM

Real

40,35

CBTOT

Real

35

CBVLSM

Real

40,35

FINPOP

Real

—

NPNTS

Int

—

PLTOT

Real

35

PLVLSM

Real

35

SEBRD

Real

35

SECUB

Real

35

SEPUL

Real

35

TREPAC

Real

40,35

TRETOT

Real

35

VARBRD

Real

35

VARPUL

Real

35

Description

Basal area/acre by d.b.h.

Basal area/acre

Sawlog bd.ft. volume

Sawlog bd.ft. volume by d.b.h.

Sawlog cub.ft. volume

Sawlog cub.ft. volume by d.b.h.

Finite population correction

Number of sample units sampled

Pulpwood cub . ft . volume

Pulpwood cub.ft. volume by d.b.h.

Sawlog bd.ft. volume standard error

Sawlog cub. ft. volume standard error

Pulpwood cub.ft. volume standard error

Trees/acre by d.b.h.

Trees/acre

Sawlog bd.ft. volume variance

Pulpwood cub . ft . volume variance

Note; All volumes are per acre

Common Area STOUT
STOUT transfers combined strata information to the output subroutine.

Description

Sawlog bd.ft. voltime

Sawlog bd.ft. volume by d.b.h.

Sawlog cub.ft. volume

Sawlog cub.ft. volume by d.b.h.

Sawlog bd.ft. volume effective d.f.

Sawlog cub.ft. volume effective d.f.

Pulpwood cub.ft. volume effective d.f.

Pulpwood cub.ft. volume

Pulpwood cub.ft. volume by d.b.h.

Sawlog bd.ft. volume standard error

Sawlog cub.ft. volume standard error

Pulpwood cub.ft. volume standard error

Total acres sampled

Note : All volumes are total

Variable

Type

Dimension

BDTST

Real

35

BDVLST

Real

40,35

CBTST

Real

35

CBVLST

Real

40,35

EFDFBD

Real

35

EFDFCB

Real

35

EFDFPL

Real

35

PLTST

Real

35

PLVLST

Real

40,35

STSEBD

Real

35

STSECB

Real

35

STSEPL

Real

35

TACRES

Real

—

Main Program INVENT
INVENT provides interactive user communication and subroutine control.
Input and output are to unit 5 which is the user's TTY on the DEC-10. The
following job control information is acquired interactively from the user:

(1) Alternate form class specification.

(2) Diameter class size.

(3) Nvnnber of strata.

-27-

(4) Confidence level.

(5) Input file name.

(6) Sampling method and sampling frame.

(7) Method of height measurement.

(8) Acreage of strata.

(9) Title of combined strata output (if more than one strata) .
INVENT opens the output file on unit 3 (LPT on DEC-10) naming it

INVENT. LPT. INVENT also sets the default confidence level to 95 percent.

The following subroutine calls are made in INVENT: FRMCLS, LINREG,
GUTS, OUTPl, 0UTP2, STRAT, and STRAT2.

Common Area Used: COMUNE.

Other

Variables Type Dimension Description

ALPHA Real — Confidence level as a decimal

ALPHCH Real(C*3) — Confidence level change flag

FRMANS Real(C*3) — Form class change flag

III Int — Loop counter for strata number

ISTRAT Int — Number of strata

REPEAT D.P.(C*8) ~ User message

SVFRAM Real(C*5) — Saves the contents of FRAME

Subroutine FRMCLS

FRMCLS is called by INVENT if the user specifies that he wishes to use
alternate form classes. FRMCLS interactively asks the user to specify the
name of the file where the alternate form classes are stored. This file is
read from unit 1 (disk on DEC-10). As in INVENT, unit 5 is used to communicate
with the user. The user is given the option to list the new form classes after
they are read from the file.

Common Areas Used: FORMCL and REG.

Other

Variables Type Dimension Description

FCFILE D.P.(C*10) — Input file name

I Int — Implied do loop index

LSTANS Int — Form class list option flag

Subroutine LINREG
LINREG is called by INVENT if the user specifies that only the heights
of a subsample of the volume trees were measured. LINREG then reads through

-28-

the data file specified by the variable IFILE and located on unit 21 (disk
on DEC-10) . For every species with three or more height trees measured,
LINREG fits the following regression model: HEIGHT = B + B * 1/DBH. The
following statistics are calculated and stored in common area REG:

(1) Regression intercept.

(2) Regression slope.

(3) Standard error of the regression.

(4) Number of height trees measured.

(5) Total number of merchantable trees sampled.

(6) Average height.

2

(7) Coefficient of determination (R ) .

(8) Mlnimtim diameter of all height trees sampled if less than 8 inches.
LINREG ignores cull trees (HEIGHT = 999) . It also checks that input parameters
are within the following ranges :

1 <_ DBH <_ 40 0 <_ HEIGHT £ 7 if height measured by logs.

0 <_ HEIGHT <^ 200 if height measured in feet.
If data is found outside these ranges, the value is set to the minimum or
maximum of the range. If an illegal species code appears, the tree is ignored
and an error message is written to unit 5.

LINREG also calculates certain statistics for the Softwood, Hardwood and
All Species accumulations. These statistics include (3), (4), (5) and (6)

above .

Common

Area Used:

REG

,

Other

Variables

Type

Dimension

CORMNY

Real

_.

CORNTR

Real

—

COVXY

Real

33

DBH

Int

~

HEIGHT

Real

— -

IDUMMY

Int

— -

IREP

Int

—

ISPEC

Int

—

ITYPE

Int

—

JTYPE

Int

—

LOOP

Int

—

MEANX

Real

35

MEANY

Real

35

MNYCOR

Real

3

Description

Corrected mean height for species

Corrected number of trees sampled

Covariance of X and Y

Diameter at breast height

Tree height

DO loop index for reps on input

Number of reps on input

Species loop counter

Forest type code

Forest type code

Loop counter

Mean of X

Mean of Y

Corrected mean height for forest types

-29-

Other

Variable

'.s Type

Dimension

Description

MSE

Real

—

Mean square error

NTRCOR

Int

3

Corrected number of trees for forest type

NUMREG

Int

3

Number of regressions run

POINT

Int

—

Point number

SPEC

Int(C*2)

—

Species input code

SSE

Real

—

Sum of the squares of error

TREES

Int

—

Number of trees on a sample unit

TYPSSE

Real

3

Forest type error sum of squares

VARX

Real

33

Variance of X

VARY

Real

33

Variance of Y

X

Real

33

1/D.B.H.

XSUM

Real

33

Sum of X

XSUMSQ

Real

33

Sum of X^

Y

Real

33

Height

YSUM

Real

35

Sum of Y

YSUMSQ

Real

33

Sum of y2

Subroutine GUTS

GUTS is called by INVENT to read the input data and calculate the basal
area per acre, trees per acre, and all of the volumes along with their stan-
dard errors. GUTS also calculates the number of units sampled and the percent
cruise. Either subroutine VOLMER or VOLTOT is called by GUTS repeatedly to
calculate tree volume.

GUTS reads the input file specified by the variable IFILE and located on
unit 21 (disk on DEC-10) . It checks the range of input data in the same way
LINREG does except that the minimum allowed DBH is set at the value of CLASIZ
Instead of 1.

GUTS checks the value of POINT on the first plot to see if POINT = -9.
Such a value indicates the data is presented in summary mode (not tallied by
individual plot) and therefore no statistical analysis of the sampling error
can be made. It also indicates that the value of TREES is the number of plots
sampled .

GUTS reads all the data accumulating sums of basal areas, number of trees,
sawlog board foot volume, sawlog cubic foot volume, and pulpwood cubic foot
volume. After all the data is read, the accumulated sums are expanded into
per acre estimates based on the sampling technique used (point or plot) . The
variances and standard errors of these volumes are also calculated. Before

-30-

attempting to calculate variances, GUTS checks to see if more than one plot
has been sampled. If not, GUTS sets itself into summary mode and does not
attempt to calculate any variances or standard errors.

Common Areas Used: COMUNE, REG, and OUTl.
Other

Variables

L Type

Dimension

DBVOL

Real

^_

BRDPT

Real

35

BRDSQ

Real

35

CUBPT

Real

35

CUBSQ

Real

35

CUBVOL

Real

—

DBH

Int

—

FACTOR

Real

—

HEIGHT

Real

—

I

Int

—

IREP

Int

—

I SPEC

Int

—

I START

Int

2

ISTOP

Int

2

LOOP

Int

—

NPLOTS

Int

—

PCTCRZ

Real

—

PCTSAW

Real

—

PLTSIZ

Real

—

POINT

Int

—

PROD

Real(C*l)

—

PULPT

Real

35

PULSQ

Real

35

PULVOL

Real

—

REGWGT

Real

—

SPEC

Int(C*2)

—

TRECNT

Int

40,32

TREES

Int

—

VARCUB

Real

—

Description

Board foot volume of a tree

Sawlog bd.ft. volume of a sample unit

Sum of BRDPT^

Sawlog cub.ft. volume of a sample unit

Sum of CUBPt2

Cubic foot volume of a tree

Diameter at breast height

Sample method expansion factor

Tree height

DO loop index

Number of reps on input

Species loop counter

First of a forest type

Last of a forest type

Loop counter

Number of plots sampled

Percent cruise

Percent sawlog

Fixed radius plot size

Point number

Product code

Pulpwood cub.ft, volume of a sample unit

Sum of PULPt2

Pulpwood cub.ft. volume of a tree

Weight for regression standard error

Species input code

Number of trees tallied by d.b.h.

Number of trees on a sample unit

Sawlog cub.ft. volume variance

Subroutine VOLMER

VOLMER is called by GUTS to calculate the merchantable cubic and board
foot volume of individual trees when tree height is specified to a merchant-
able top. The input parameters are d.b.h., merchantable height, and species.

The cubic foot volume equations are from Barnard et al. (1969). To
use these equations all species must be placed in one of seventeen groups.
The species grouping is as follows:

-31-

Spec

ies

in

Group

WP,

RP,

OH

SP

BF

HM

PP,

TA,

NS, SC

CE

SM

RM,

YP

WA,

AS

BC

WB,

YB,

SB, GB

BE

B

RO,

BO,

GU

WO,

EL

HI

WI,

HH,

OH

Species Group
Number

1

2

3

4

5

6

7

8

9
10
11
12
13
14
15
16
17

_ The board foot volume equation is from Wiant and Castenada (1977),

This equation is a formulation of the Mesavage and Girard (1946) form class

tables. A form class must be specified for each of the 32 species. VOLMER

sets the form classes to the following defaults :

Hemlock 70

All other softwoods 76
All hardwoods 74

The user can specify alternate form classes at execution time.

The cubic foot equations are based on a 4" merchantable top while the
board foot equations are based on a 10" merchantable top. The equations
are presented in full in Appendix C, pages 56 and 57.

Common Area Used: FORMCL.

Other

Variabl

es Type

D

imension

BDHT

Real

BDVOL

Real

—

CBHT

Real

—

CUBINT

Real

17

CUBSLP

Real

17

CUBVOL

Real

—

DBH

Int

—

GRP

Int

32

HEIGHT

Real

—

I

Int

—

ISPEC

Int

—

WIANT

Real :

Functii

on

Description

Board ft. merchantable height in logs

Board ft. volume of tree

Cubic ft. merchantable height in feet

Cubic ft. volume equation intercept

Cubic ft. volume equation slope

Cubic ft. volume of tree

Diameter at breast height

Species group for cub.ft. volume equation

Tree height

Implied DO loop index for data statements

Species loop counter

Computes board foot volume

-32-

Sub routine VOLTOT

VOLTOT is called by GUTS to calculate the merchantable cubic and board
foot volume of individual trees when total tree height is measured. The
input parameters are d.b.h., total height, and species.

Both the cubic and board foot equations are from Honer (1967). In

applying these formulas, the gross cubic foot volume of the tree is first

calculated. The merchantable cubic and board foot volumes are then derived

from the gross volume. To use these equations it is necessary to place all

species into one of 16 groups. The species grouping is as follows:

Species Group Species in

Number Group

1 WP

2 RP, PP. TA, SC, OS

3 SP, NS

4 BF

5 CE

6 HM

7 AS

8 WB, SB, GB

9 YB

10 SM, RM, YP

11 BW

12 BE

13 BC

14 EL

15 WI, HH, OH

16 RO, WO, BO GU, HI, WA

Stump height and merchantable top are variables in the equations.

VOLTOT sets these to the following defaults:

Stump Height 0.5 ft.

Cubic Foot Top 4 in.

Softwood Board Foot Top 6 in.

Hardwood Board Foot Top 8 in.

The equations are presented in full in Appendix C, page 58.

No common area used.

-33-

Variables

. Type

Dimension

ASDRBF

Real

3

ASDRCF

Real

3

BDVOL

Real

—

BRDTOP

Real

2

CUB INT

Real

16

CUBSLP

Real

16

CUBTOP

Real

2

CUBVOL

Real

—

DBH

Int

—

GRP

Int

32

HEIGHT

Real

— -

I

Int

—

ISPEC

Int

—

ITYPE

Int

—

STUMP

Real

—

TCBVOL

Real

—

XBRD

Real

_

XCUB

Real

Description

Board foot volume equation coefficients

Cubic foot volume equation coefficients

Sawlog bd.ft. volume of tree

Bd.ft. merchantable top diameter

Cub.ft. volume equation intercept

Cub.ft. volume equation slope

Cub.ft. merchantable top diameter

Cubic ft. volume of tree

Diameter at breast height

Species group for volume equations

Tree Height

Implied DO loop index for data statements

Species loop counter

Forest type indicator

Stump height

Total cubic foot volume

Intermediate step in bd.ft. volume

calculation
Intermediate step in cub.ft. volume

calculation

Subroutine STRAT

STRAT is called by INVENT if more than one stratum are being processed.
STRAT calculates combined strata statistics for sawlog cubic foot volumes,
sawlog board foot volumes and pulpwood cubic foot volumes. These statistics
include total volume over all strata along with their standard error and
effective degrees of freedom. STRAT also calculates the total acreage of
all the strata combined.

STRAT has two entry points, STRAT and STRAT2. STRAT is used as an
entry point after each individual stratum is processed. Sums are accumulated
for volumes, variances, and squared variances. STRAT2 is called after all
of the strata have been processed. The previously accumulated sums are used
to calculate total volumes, standard errors, and effective degrees of freedom.

Common Areas Used: OUTl and STOUT.

Other

Variables

Type

Dimension

ACRES

Real

—

CLASIZ

Int

—

DBH

Int

—

DFBD2

Real

32

DFCB2

Real

32

DFPL2

Real

32

ISPEC

Int

—

STVMBD

Real

35

STVMCB

Real

35

STVMPL

Real

35

I

Int

—

J

Int

_

Description

Acreage of compartment
Diameter class size
Diameter at breast height
Intermediate calculation in EFDFBD
Intermediate calculation in EFDFCB
Intermediate calculation in EFDFPL
Species loop counter
Sawlog bd.ft. variance
Sawlog cub.ft. variance
Pulpwood cub . ft . variance
Implied DO loop index
Implied DO loop index

-34-

Subroutine OUTPl

OUTPl is called by INVENT to produce all of the program's output. The
subroutine has two entry points, OUTPl and 0UTP2. OUTPl is used to produce
the individual stratum output while 0UTP2 is used to produce the stratified
total output. All output is written to unit 3 (LPT on DEC-10) creating a
file named INVENT. LPT.

OUTPl produces the following output tables for each stratum:

1. Summary of job control related information (Table 1).

2. Analysis of the height-dbh relationship if regression option
used (Table 17) .

3. Individual stand and stock table for each species sampled
(Tables 2, 3).

4. Stand and stock tables for Softwood summary. Hardwood summary,
and All Species summary (Tables 4, 5, 6).

5. Table of species composition by percent (Table 7).

6. Table of volume totals for all species (Table 8).

7. Table of volume totals for all species expanded by acreage
(Table 9).

The stand and stock tables include basal area/acre, trees/acre, sawlog cubic

and board foot volumes per acre, and pulpwood cubic foot volume per acre,

all presented by diameter class as well as a total over all diameters. A

confidence interval expressed as a percent is given for each volume total.

Also listed are mean stand diameter, merchantable m.s.d. (mean stand diameter

of trees with dbh greater than or equal to 8 inches), percent cruise, sawlog

board foot coefficient of variation, and pulpwood cubic foot coefficient of

variation. Furthermore, the title of each compartment is printed on the top

of each page and the pages are consecutively numbered.

0UTP2 is used as an entry point to produce the following stratified total

output if more than one stratum were sampled:

1. Summary of job control related information (Table 10).

2. Volume tables for each species sampled (Table 11, 12).

3. Volume tables for Softwood summary. Hardwood summary, and All Species
summary (Tables 13-15).

4. Table of volume totals for all species (Table 16).

-35-

All volume figures from 0UTP2 are total volumes expanded over the entire
acreage of all the strata combined. The individual species tables and the
three summary tables list volumes by diameter class as well as giving a
total over all diameters. A confidence interval expressed as a percent is
given for each volume total. Each page of the output is consecutively
numbered and a specified title is printed at the top of each page.

Both OUTPl and 0UTP2, when producing individual species tables, will
print the form class used to calculate board foot volume next to the species
title if the trees were measured to a merchantable top. Otherwise, no form
class information is printed.

Common Areas Used: COMUNE, FORMCL, OUTl, REG, and STOUT.

Other

Variable

es Type

Dimension

AFCNUM

ReaKC

*2)

AFCTIT

ReaKC

*5)

—

ALPHA

Real

—

BASCOM

Real

—

BDCOM

Real

—

CUBCOM

Real

—

DBH

Int

—

DP

Real

—

I

Int

—

IBDCV

Int

—

IBDTST

Int

—

IBDVST

Int

—

ICBTST

Int

ICBVST

Int

—

IDATE

int(c*:

5)

2

lERBRD

Int

—

lERCUB

Int

—

lERPUL

Int

—

IFLAG

Int

—

IPAGE

Int

—

IPLCV

Int

—

IPLTST

Int

—

IPLVST

Int

—

ISPEC

Int

ISTRA

Int

—

J

Int

—

MSD

Real

—

OFRAME

ReaKC^

*5)

—

PCTCRZ

Real

—

PULCOM

Real

—

RBSTOT

Real

—

RECUSE

ReaKC^

*A)

—

RMSD

Real

—

RTRTOT

Real

—

SPETIT

Int(C*

5)

35,3

Description

Form class number as a literal

Output title

Confidence level as a decimal

Percentage of basal area

Percentage of sawlog bd.ft. volume

Percentage of sawlog cub.ft. volume

Diameter at breast height

Degrees of freedom

Implied DO loop index for data statements

Board ft. coefficient of variation

Total sawlog bd.ft. volume all strata

Total sawlog bd.ft. volume all strata by

dbh
Total sawlog cb.ft. volume all strata
Total sawlog bd.ft. volume all strata by

dbh
Date

Percent error in sawlog bd.ft. volume
Percent error in sawlog cb.ft. volume
Percent error in pulpwood cb.ft. volume
Warning message flag
Page number

Pulpwood coefficient of variation
Total pulpwood cb.ft. volume all strata
Total pulpwood cb.ft. volume all strata

by dbh
Species loop counter
Number of strata sampled
Implied DO loop index
Mean stand diameter
Saves the contents of FRAME
Estimated percent cruise for point sampling
Percentage of pulpwood volume
Basal area of trees with dbh <^ 8"
Output message

Merchantable mean stand diameter
Trees/acre of trees with dbh <^ 8"
Species titles for output

-36-

Other

Variabl

.es

Type

Dimension

TOTABD

Int

__

TOTACB

Int

—

TOTAPL

Int

—

TOTPTS

Int

—

TRECOM

Real

—

TVAL

Real

—

Description

Total sawlog bd.ft. volume
Total sawlog cb.ft. volume
Total pulpwood volume
Total number of sample units
Percentage of trees/acre
t -value

Subroutines MDSTI and MDNRIS

MDSTI and MDNRIS are adapted from the International Mathematical and
Statistical Libraries, Inc. (IMSL) . Both subroutines have been edited
removing capabilities not needed in this application. The remaining code
is reproduced in Appendix E in accordance with the IMSL policy on research
work as stated on page INTRO-22 of the July, 1977 IMSL library manual.

MDSTI is called by OUTPl to calculate the percentage point of the
Student's t distribution. MDNRIS is called by MDSTI and is used to evaluate
the inverse normal probability distribution function. MDSTI gives 5 signif-
icant digit accuracy (IMSL 1977). Furthermore, table 26.10, p. 999 of the
Handbook of Mathematical Functions edited by Abramowitz and Stegun was
duplicated by MDSTI (IMSL 1977) .

No variable dictionary is presented. The interested reader may refer
to pages MDSTI-1, MERFI/MERFCI/MDNRIS-1, and MERFI/MERFCI/MDNRIS-2 of the
1977 IMSL Library manual.

-37-

SECTION 4. STATISTICAL ANALYSIS

Overview

The purpose of INVENT 's statistical analysis is to indicate the reli-
ability of the volume estimates. This is accomplished by computing a
confidence interval for all volume estimates that are summed across the
diameter classes. This confidence interval is specified as a percent of
the respective volume estimate. The user specifies the confidence level
as a decimal in the range .65 to .99. All statistical analysis assumes a
stratified random sample where each sample unit is a cluster. There is an
option to double sample for height within the clusters.

INVENT calculates confidence intervals for sawlog cubic foot volume,
sawlog board foot volume, and pulpwood cubic foot volume. The analysis
for each volume is carried out independently of the other volumes. Since
the method is exactly the same for each volume, the calculation of error
for an arbitrary volume V is presented. Throughout the discussion we will
assume the reader is familiar with the techniques for calculating the mean
volume per acre when the heights of all volume trees are measured. (Cf.
Barrett and Nutt, 1975; Avery, 1975; Husch, Miller, and Beers, 1972;
Dillworth and Bell, 1976.)

The statistical analysis is presented in the same sequence as the
program performs it. For each individual stratum the standard error of the
volume must be calculated. If the height subsampling option is used, the
program first calculates the standard error due to height estimation. Next
the percent cruise of the stratum is calculated. This is used to refine
the standard error due to the variation between clusters which is calculated
next. This standard error is then combined with the height estimation
standard error to give the volume standard error of the stratum. Finally,
the volume standard errors of all the strata are combined to yield the
standard error of the stratified volume estimate along with the effective
degrees of freedom of that estimate.

-38-

Standard Error of Height Estimation
INVENT allows the user to double sample within each cluster for tree
height. For each species where three or more height trees are measured
within the stratum, a regression relationship between height and the inverse
of d.b.h. is analysed. When calculating volume for a tree whose height was
not measured, the regression equation is used to generate a height for that
tree. If the user were required to measure the heights of all trees on
some of the clusters while measuring no heights on the remaining clusters,
the volume standard error could be calculated using techniques developed by
Donald Bruce (1961) and Floyd A. Johnson (1958). Bruce 's method is outlined
in Appendix B, page 53. Johnson's method is an application of double
sampling with ratio of means estimation and this method is presented in
Appendix B, page 52. Neither method Is applicable if the user is free to
subsample a portion of the trees within each cluster. A new statistical
technique was developed to approximate the standard error of the resulting
volume estimate.

The authors believe that this new technique is far more efficient in
field application. It allows the user to sample heights more intensely on
the more valuable species. It provides far more flexibility than the
traditional approach.

H = height D = d.b.h. T = // trees sampled t = // height trees sampled
Regression Model: H = B^ + B^ ^"n"^ "*" ^

„ ^ /-u^ c2 MSE
Variance (H) = S„ = — — -
u t

where MSE = mean square error of the regression (for

computational formula see Neter and Wasserman,
1974, p. 45).

^H
Standard Error (H) = SE„ = — 2_ (iqO)

" H

SE is calculated for every species where three or more height trees
H

are measured. It is also necessary to calculate SE^ for the Softwood, Hard-
wood and All Species summary volumes. The method used is from Neter and

-39-

Wasserman pp. 163-4 (1974). It is equivalent to fitting a regression for

each of the three groups using pooled data from each species in the group.

However the method used is far more efficient. We simply define T , t , and

MSE for each summary group as a function of T, t, and MSE for every species

in that group. For instance, suppose we want SE for hardwoods where the

n

hardwood group in this stratum consists of three species. Then,

*

T = T + T + T
1 2 ^3

. MSE, (t,-2) +MSE_(t,-2) + MSE-(t_-2)
MSE = ^^ 2 2 3 3

(t^-2) + (t2-2) + (t3-2)

Using these new values, SE is calculated as before.

H

Percent Cruise
Percent cruise is simply the percent of the population actually sampled.
When sampling with fixed-radius plots, percent cruise is easily calculated
as follows:

Let P = percent cruise

PLT = plot size in acres

n = number of plots sampled
AC = tract size in acres
Then. P = " '^J^'^ (100)

When sampling with variable-radius plots (point sampling) the true value
of P is unknown. However, based on our sample estimate of MSD we can estimate
the value of P (Cunia, 1959).

Let P = percent cruise

n = number of sample points
AC = tract size in acres
BAF - basal area factor of angle used

B = total basal area of all trees on tract
T = total number of trees on tract

MSD = mean stand diameter (quadratic mean)

100 • n • B _ 0.5454 • n • MSD^
' T • AC • BAF ~ AC • BAF

-40-

It should be noted that in a variable-radius cruise, the percent cruise may
be different for each species. See tables 2 and 3.

In forest sampling, we deal with finite populations. The formula we
use to estimate the variance of the mean volume per acre for each stratum
includes a term known as the finite population correction (FPC) . The FPC
reduces the variance of the mean when a significant (usually > 5 percent)
portion of the population has been sampled.

Let N = number of sample units in population
n = number of units actually sampled

N - n
Then, FPC = " "^

N

AC
In plot sampling N = and n = number of plots sampled. However,

if Lt L

in point sampling there is no easy was to calculate N. Therefore, mensura-
tionists have traditionally ignored the FPC when analysing point cruises.
This can be quite inefficient when tract sizes are small as is often the
case in New Hampshire. This problem can be easily overcome using the
previous result permitting the calculation of percent cruise when point
sampling. It can easily be shown that the FPC is simply a function of P
(percent cruise) .

FPC = N - n ^ _N D_ = 1 _ ^ = 1 P

N N N N 100

INVENT uses this relationship to refine the variance estimates when point

p
sampling as well as plot sampling. INVENT calculates FPC as 1 - y^ .

The standard error due to the variation in volume between clusters, SE^,

is calculated as follows:

Y = mean volume /acre of stratum

2
S = sample variance of volume

n = number of sample units

1

then, SE^ =

(100) (FPC)

If the height of all volume trees was measured, then SE— is the standard
error of the volume estimate for the stratum, SE . If, however, the

-41-

height double sampling option was used, SEL must reflect both SE— and SE„,
the standard error due to height estimation. In this case, SE^ is
calculated as follows :

^^voL = ^/ 4 ^ -H^ ^'4 )

Combined Strata Statistics
We have seen how INVENT calculates the mean volume/acre with its standard
error for each stratum in the sample. The final step in the statistical
analysis is to use the individual stratum statistics to calculate the total
volume across all strata along with the standard error.
Let L = number of strata sampled

AC. = size of stratum i in acres

N = number of sample units in stratum i

/

n = number of units sampled in stratum i

Y = mean volume per acre of stratum i

VOL = total volume of stratum i = AC " Y

^%

OL = volume standard error of stratum i as a percent

2

2 2 2
VOL = (SE^Q, • VOL ) / 100 = variance of total volume

1
then we calculate the stratified total statistics as follows:

AC = Z AC = total acreage sampled
i=l ^

L

VOL = Z VOL = stratified total volume
i=l ^

J I 'K

»^0L "°«

Having calculated VOL and SE^ , all that remains in calculating the
confidence interval is determining the appropriate percentage point of the

Student's t distribution. In order to do this we must calculate the

2
degrees of freedom associated with S . Satterthwaite (1946) states that

2
the exact distribution of the complex estimate S is too involved for

VULi

-42-

everyday use. He suggests approximating the true distribution with a Chi-
square distribution of equal variance. Cochran (1977 p. 96) uses this
approach and defines the degrees of freedom of the approximating Chi-square
as^

DF

eff

i=l ^

i=i ^i"^

where g .

"i ^^i " ^i^

S. = sample variance

A more appropriate form for our purposes is given by
L

DF

eff

ifl '^°^i

L VOL .
i=l ^i- ^

The equivalence of the two formulas is shown in Appendix B, page 54.

Using this result the appropriate value t is derived from the Student's
t distribution and the confidence interval is given by:

VOL ± fSE.

VOL

-43-

APPENDIX A

DATA CODING INFORMATION

-44-

EXAMPLES OF TREE CODING OPTIONS

■3-

o

o

>— 1

hO

LO

lO

■—1

CNJ

I— 1

I— 1

(N

rH

.-H

ro

H

i-H

CD

■—1

CT>

o

00

o

o

d-

o

OO

•

d-

CN

«

ID

•

r^

CN

CM

H

UD

CO

w

Pi

*<>

dO

LO

UD

00

o

ID

CM

CV4

.=r

rH

I-H

H

.H

CM

H

N-\

CNI

pL-

a,

W

O

Oi

S

1— 1

s

a;

m

a;

CO

X

Z

(0

•p

■M

W

p

1 w

CO Uh (X

bO

0)

>i 0)

■P bO

•p (1)

bO •,-) <u

z

0

0)

M (I)

O 0

O 01

O O fj •

3

>,M

>,4h

C M-1

3 -i

D iP

>i H 0) P

-I

o

O

T3

T3

.H O, Mh O

c c

c c

C

0 c

O C

c c t/1 o ;3

O -H

O -H

T3 -H

fl -rH

fn -H

O -rH T3

O

cu

a,

p n 0

bO -M

W) -P

0 +->

■p

p

bO P rH O f<

O JZ

0 ^

3 J=

XI JH

T3 JZ

o j2 3 -H a.

■H bO

f-{ M

a bo

dJ bO

0) bO

rH bO (0 P

S -H

S -H

H -H

X -H

X -H

s -H 4h m xi

Ifl 01

(0 0)

„^ I*

•H QJ

•H QJ

to 0) 0) o c

w a:

w ac

CL, X

s as

S X

00 a: Q -H m

-A5-

INVENT Species List

White Pine

WP

Red Pine

RP

Pitch Pine

PP

Hemlock

HM

Balsam Fir

BF

Tamarack

TA

Spruce

SP

Norway Spruce

NS

Scotch Pine

SC

Cedar

CE

Other Softwoods

OS

Sugar Maple

SM

Red Maple

RM

Ash

WA

Willow

WI

Aspen

AS

Yellow Poplar

YP

Black Cherry

BC

White Birch

WB

Yellow Birch

YB

Sweet Birch

SB

Grey Birch

GB

Beech

BE

Basswood

BW

Red Oak

RO

White Oak

WO

Black Oak

BO

Elm

EL

Gum

GU

Hickory

HI

Hophornbeam

HH

Other Hardwoods

OH

-46-

INVENT TALLY SHEET

COMR

^

POINT

±

ff

TREES

1 1

SPECIES

DBH

HEIGHT

PRODUCT

i=t:

f^P

tW

l.O

wP

ao

H.O

iRP

/^

z.^

WjP

/9r

H.O

wP

o^

s.o

wP

I"?

l.iT

wp

l(o

z.o

wp

0.0

c/.O

wp

XX

v.r

U^\P

-x^

V'^T

RP

l(o

3.0

COMR

/

L

POINT

^

^ /^

TREES

^

SPECIES

DBH

HEIGHT

FRODUC T

i=tr

luP

207

-i.O

^uP

^2.

^. s

RP

2^

3.r

RP

l^

2>.S"

({o

N

LO

(lo

:i^

\.o

Ro

1?

1. o

wp

Q.L,

?'C

^

COMR
POINT

a

^ /^

TREES

/^

SPECIES

DBH

HEIGHT

PRODUCT

:«r

vue

/z

a r

LV-'P

9.0

3.6~

^'P

3a

v.o

(ip

/6

3.0

a^

/ o

/.o

Ce

OG

• r

P

A

ujp

/(:?

3>-S

UjP

/G

^.o

wP

O

ao

Rp

/<^

a.o

wp

^(o

'V.O

COMR
POINT

i

V ^

■TREES

IV

SPECIES

DBH

HEIGHT

PRODUCT

:fct

u^P

2-:i

"i.O

wP

2 o

H-O

u^f

3 0

"i-o

(

wP

:^o

Z.'T

ucP

ic

H^O

3.

U-P

a.(^

9.S-

Cc/P

/t-

"i.O

U.P

^?

H. r

u/P

/§■

l.o

U/ P

:ikc

w. s~

wP

?.?

^. s~

wP

5:i,

w.o

^p

3LC

^.s~

-47-

INVENT DATA FILE

WHITAKER WOODS

01 11

RP 14

3.0

WP 20

4,0

RP 14

3.5

WP 18

4.0

WP 12

3.0

WP 18

3.5

WP 16

3.0

WP 20

4.0

WP 22

4.5

WP 22

4.5

RP 16

3.0

02 12

WP 12

2.5

WP 20

3.5

WP 32

4.0

RP 16

3.0

BE 10

1.0

BE 06P0.5 2

WP 16

3.5

WP 16

4.0

WP 12

2.0

RP 16

3.0

WP 26

4.0

03 09

WP 26

4.0

WP 32

4.5

RP 22

3.5

RP 18

3.5

RO 14

1.0

RO 22

1.0

RO 18

1.0

WP 26

5.0 2

04 14

WP 22

4,0

WP 30

4.0

WP 30

4.0

WP 20

3.5

WP 20

4.0 2

WP 26

4.5

WP 16

1 . 0

UP 18

4.5

WP 18

3.0

WP 20

4.5

WP 28

4.5

WP 32

4.0

WP 20

4.5

05 13

WP 24

3.0

RP 16

2.0

WP 32

4.0

WP 26

3.5

RO 12

1.0

HM 18

1.5

HM .1.6

1.0

RP 12

1.0

RP 10

1.0

RO 16

1.0 2

UP 14

2.5

HM 18

2.0

1978 COMP. 1

-48-

INVENT Processing Request Form

Forester Information:

Name:

Firm Name;_
Address:

Tract Information:

Number of compartments:

Do you wish to use alternate form classes? ( yes or no )

If yes enter fona classes below.

WP

RP

PP

HM

BF

TA

SP

NS

sc

CE

OS

SM

RM

WA

WI

AS

YP

BC

WB

YB

SB

GB

BE

BW

RO

WO

BO

EL

GU

HI

HH

OH

Diameter class size

INVENT uses the 9^% confidence level as a default. If you wish
to use a different level, please specify (65 - 99)

Tract Description (up to 70 characters)

-49-

IKVENT Compartment Information Form

Compartment Name and/or Description (up to TO characters)

Sampling Method ( Point, Plot, Strip, or 10055)

If Point: BAT =

If Plot : Plot size = acres

If Strip: Total acreage of all strips = acres

Tree Height Measurement Units ( l6 ft. logs or feet )

Tree top specification ( merchantable, total, or mixed )

Were the heights of all volume trees measured? ( yes or no )

Acreage of compartment : acres

-51-

APPENDIX B

STATISTICAL FORMULAS AND DERIVATIONS

-52-

DOUBLE SAMPLING WITH RATIO OF MEANS EXPANSION

X. = Basal area per acre at cluster 1

Y = Volume per acre at cluster i,
N = Number of possible sample units
n = Total number of units sampled

n = Number of volume units sampled
o

X = Average basal area of all units sampled

X = Average basal area of the volume units sampled

Y = Average volume of the volume units sampled
o

Y = Ratio estimate of the mean volume per acre

Y

B - -T-

O

Y = B • X

r

t - C-^X^l (-ii^^^ 4(-^

-53-

STATISTICAL ANALYSIS OF A VARIABLE PLOT CRUISE USING BRUCE 'S METHOD

(Bruce 1961)

Notation
y - volume points

a - all points
n - number of points

TC

tree count

VBAR

volume to basal area ratio

VOL

volume per acre

SE

standard error expressed as a percent

CF

correction factor

Example of Notation;

SE - standard error of the tree count on all points expressed as a
a percent

Estimated volume = VOL

TC /n

g g

TC /n
Y Y J

= (VOL • TC • n )/(TC • n )
Y a Y Y a

/ 2 2

Uncorrected standard error = / SE,„.„ + SE

. / VBAR TC
V Y a

Corrected standard error = CF- (uncorrected standard error)

/OL

^hc

where CF =

<AR +<
Y Y

-54-

FORMULA FOR EFFECTIVE DEGREES OF FREEDOM

DF

eff

i=l "i"^

where g .

^i ^^i-\^

L N. (N.-n^) ^^y
LV n, / ^1

1=1

n^-1

L S^ _ (N,-n,) n2 Y
1=1 "l ^1 J

1=1

n^-1

a \- »o '

a ^vo.)^

4 4

J. 1

i-l "l"^

1=1 ^r^

-55-

APPENDIX C

VOLUME FORMULAS

-56-

CUBIC FOOT VOLUME EQUATIONS

BASED ON MERCHANTABLE HEIGHT

(Barnard et al. 1973)

2

General form of equation: VOLUME = a + b (DBH -HEIGHT)

where height is measured to a A" top in ft.

Species Group

a

b

1

3.5142

0.00236

2

2.1998

.00257

3

1.4793

.00272

4

2.4784

.00242

5

-0.0496

.00303

6

1.5817

.00259

7

1.6823

.00310

8

1.1763

.00310

9

1.0067

.00293

10

1.1809

.00292

11

1.1339

.00281

12

1.2851

.00326

13

0.8976

.00349

14

1.2027

.00301

15

1.0009

.00300

16

1.4438

.00316

17

0.8591

.00309

For the species composition of each group see page 31.

-57-

BOARD FOOT VOLUME EQUATION

BASED ON MERCHANTABLE HEIGHT

(Wlant and Castenada 1978)
(INT. 1/4-lnch Rule)

D = diameter at breast height

H = merchantable height to a 10" top in 16 foot logs

FC = Girard form class

VOLUME = [ (aQ+a^H+a2H^) + (bQ+b^H+b2H^)D + (Cq+c^H+C2H^)D^] [ (FC-78) (.03)+!]

^0

ss

-13.35212

^

=

9.58615

^2

=

1.52968

^0

=

1.79620

^1

=

-2.59995

^2

=

-0.27465

•=0

=

0.04482

'^l

=

0.45997

<=2

=

-0.00961

Equations:

-58-

CUBIC AND BOARD FOOT VOLUME EQUATIONS

BASED ON TOTAL HEIGHT

(Honer 1967)
(INT. 1/4-inch Rule)

D = diameter at breast height

H = total height in feet

GCVL = gross cubic foot volume

MCVL = merchantable cubic foot volume

MBVL = merchantable board foot volume

STH = stump height

CTOP = cubic foot merchantable top diameter

BTOP = board foot merchantable top diameter

GCVL = °

H

X =

,_CTOPJ_,

H + STH

H

MCVL =

[b + b, X + b. X J • GCVL
o i i.

Y =

P BTOP , p H + STH ,

MBVL

[c + c- Y + c„ Y^] • GCVL
o i /

where

STH

=

0.

,5 ft.

CTOP

=

4,

.0 in.

for all species

BTOP

=

6,

,0 in.

for softwoods

8,

.0 in.

for hardwoods

t'l

=

0,

.9604

=

-0,

.1660

b2

=

-0,

.7868

<^0
'^l
^2

=

5,

.4332

=

-1,

.6281

=

-4,

.4710

Species

Group

^0

*1

1

0.691

363.676

2

0.710

355.623

3

1.226

315.832

4

2.139

301.634

5

4.167

244.906

6

1.112

350.092

7

-0.312

436.683

8

2.222

300.373

9

1.449

344.754

10

1.046

383.972

11

0.948

401.456

12

0.959

334.829

13

0.033

393.336

14

0.634

440.496

15

1.877

332.585

16

1.512

336.509

For the species composition of each group see page 32.

-59-

APPENDIX D

INVENTORY DESIGN AND SAMPLE SIZE ESTIMATION

-60-

Inventory Design

INVENT analyzes data based on two sampling designs:

(1) simple random sampling,

(2) stratified random sampling with simple random sampling within
each stratum.

Sampling Is assumed without replacement — that Is, sample units (plots
or points) are far enough apart so that a tree is sampled on only one sample
unit.

INVENT allows for the possibility of double sampling regression estima-
tion of volume. The error associated with regression estimation is small
relative to the error among sample units. Thus, it is ignored in the
discussion given below on sample size approximation.

Foresters commonly use systematic sampling rather than simple random
sampling. Generally, estimates from systematic sampling are somewhat more
precise than from simple random sampling provided the systematic sample
units are spread well over the tract sampled. Therefore, confidence
intervals based on simple random equations provide somewhat conservative
estimates for well-representative systematically collected data. For the
same reason, sample size is usually overestimated.

Sample Size
Sample size is treated briefly. More detailed explanations are given
in Barrett and Nutt (1975) , Cochran (1977) , Freese (1962) , and Dllworth and
Bell (1976).

If a tree is sampled on more than one sample unit, the estimated means are
still unbiased but the variance is underestimated.

-61-

Slmple Random Sample

Sample size from simple random sampling is estimated by the following

equation: ^

r t CV ,^
n = I — ] Equation 1

where t is a t-value, CV is the coefficient of variation, and E is the
allowable error in percent.

The t-value depends on the confidence level desired. The approximate
t-values for the most widely used confidence levels are:

Confidence level

(percent) t-value

80 1.3

90 1.7

95 2.0

99 2.7

The coefficient of variation is a measure of the variability among
volume estimates in a forest based on a particular technique of sampling.

Assume that from past experience a forester estimates that the
coefficient of variation is 50 percent when he samples with a 10 B.A.F.
prism. He wants to estimate the volume on a tract with an allowable error
of 10 percent at the 95 percent confidence level. That is, he wants the

chances to be 95 in 100 that his estimate after completing the inventory

3
is within 10 percent of the actual volume.

The sample size by Equation 1 is estimated as follows:

, = [_E_C^L]' = [^iP^]' = 100 points

Thus, the forester must select 100 points.

2

The estimate of sample size using simple random and the subsequent estimates
using stratified random sampling can be adjusted by the equation

n =

a 1 + P/100

where P is the percent cruise. This additional refinement in estimating
sample size is generally not necessary except on small tracts.

3
To be within 10 percent of the actual volume implies that errors such as tree

measurement errors, errors in estimating form class, and errors in volume

equations are negligible.

-62-

Stratified Random Sampling

Often a forester divides the area into compartments or stands before
beginning the inventory. In other words he stratifies the area.

If the forester's primary interest is in the estimates within each
stratum, then he can specify the allowable error for each stratum and apply
Equation 1 to estimate sample size.

If his primary interest is in the mean volume per acre or total volume
over all strata, then he can estimate sample size by one of two techniques.
Each technique depends on how sample units are allocated to strata. The
two techniques will be illustrated with an example.

First consider proportional allocation. By proportional allocation,

the forester determines sample size and allocates sample units to strata in

proportion to the acreage of each stratum. The equation for sample size is

t 2 L
n = [-^] I W.CV7 Equation 2

1=1

where t is the t-value, E is the allowable error, W. is the proportion of
acreage in stratum i, CV. is the coefficient of variation in stratum i,
and L is the total number of strata.

Assume a forester has the following information on an 80-acre tract he
wishes to sample with a 10-B.A.F. prism.

Stratum

Acreage

\

CV^

W^CV2

1.

2.
3.

small pines
large pines
mixed pine-hardwood

36
18
26

.450
.225
.325

30
50
80

405.00

562,50

2080.00

Totals

80

1.000

—

3047.50

Suppose he wants to estimate the mean volume per acre within 10 percent

with 95 percent confidence. Then, by applying equation 2, he estimates the

number of points required.

2 L 2

n = [-^] Z W^CV^ = [j|-] (3047.50) = 122 points

He would allocate points to strata as follows:

-63-

Stratum W n'W

1

.450

55

2

.225

27

3

.325

40

Totals 1.000 122

By the second technique of allocation, called Neyman allocation, the forester
takes into account the coefficient of variation as well as the relative
acreage in each stratum in assigning sample units to strata. Neyman allocation
minimized the number of sample units required to achieve a specified precision.
The equation for sample size is

2

n = i-^]'

T. W CV
1=1

-i2

Equation 3

To apply this equation to the data given above we require the following
computations :

Stratum

w.

1

CV^

W^CV^

1
2
3

.450
.225
.325

30
50
80

13.50
11.25
26.00

Totals

1.000

—

50.75

Assuming an allowable error of 10 percent with 95 percent confidence, the

forester estimates the number of points required.

-i2

n = [^]'

1 W CV
i=l ^

= [j|-] (50.75)^ = 103 points

The points are allocated to strata as follows ;

L

Stratum

W.CV^

n-W CV / l W.CV.
i 1 i=l ^

1

13

.50

27

2

11

.25

23

3

26

.00

53

Totals

50

75

103

Proportional and Neyman allocation are two of the most common allocations.
Other techniques of allocation are given in the references cited.

-65-

APPENDIX E

PROGRAM LISTING

-66-

c
c*

c

c
c
c
c
c
c
c*

PROGBAH INVERT
INVENT VEB. 3 1/1/79 DAVID S. LINDEN I. N. E. R. O.N. H.

THIS IS THE HAIN PBOGRAH OF INVENT

IT PROVIDES THO HAT OSER COHflDNICATION AND SOBBOOTINE CONTROL

DOOBLE PBECISION IFILE, SAHHTD, REPEAT
INTEGER LOG ANS, TOP ANS.CLASIZ, REGANS, TITLE (14)

COHHON/COHONE/CLASI Z, TITLE, HO LT,PBOBLV,BAF, TOP ANS.BEGANS, ACRE
1, FRAllE,SAnHTD, LOGANS, IFILE
DATA PHOBLV/.05/
OPEN OOTPOT FILE

OPEN (0NIT=3 , ACCESS= • SEQOOT* , FILE=« INVENT. LPT' )

OSER COHBONICATION PORTION
****** ********* •#««****«•******

ASK OSER IF HE WANTS TO OSE ALTERNATE F.C. •$
10 »RITE(5,20)

20 FOBHAT(*0DO ¥00 WANT TO READ IN ALTERNATE FORM CLASSES •,
1'FROa DISK?*)
READ(5,270) FBHANS
IF(FRHANS.BQ. 'YES') CALL FRllCLS
IF<FBHANS.NE. 'NO*. AND.FRBANS.NE.'IES') GOTO 10
HAVE OSER SPECIFY DBH CLASS SIZE
WRITE(5,30)
30 FORHAT(«0 ENTER DIAHETER CLASS SIZE')
READ(5,*) CLASIZ
HAVE OSER SPECIFY NOBBEB OF STRATA
NRITE(5,i(0)
40 FORHATCO ENTER NOHBER OF STRATA')
READ(5,*) ISTRAT
ASK OSER FOR NEW VALDE OF ALPHA IF DESIRED
50 HBITE(5,60)

60 FOHMATCOINVENT OSES THE 95X CONFIDENCE LEVEL AS A DEPAOLT. ' /
1' DO YOO WISH TO OSE A DIFFERENT LEVEL?')
READ(5,270) ALPHCH

IF(ALPHCH.NE. ' YES*. AND. ALPHCH. NE. 'NO') GOTO 50
IF (ALPHCH. Eft. 'NO') GOTO 90
70 HRITE(5,80)

eo FORHATC ENTER CONFIDENCE LEVEL AS A DECIHAL IN THE RANGE ',
1'.65 TO .99')
READ(5,*) ALPHA

IF (ALPHA.GT. 0.99. OB.ALPHA.lt. 0.65) GOTO 70
PROBLV=1.0- ALPHA
90 CONTINOE

BEGIN PROGBAB CONTROL LOOP

DO 310 111=1, ISTRAT

HAVE OSER SPECIFY INPOT FILE

WRITE(5,100) III
100 FORBATCO ENTER INPOT FILE NABE FOB STRATOB • *,I3)

READ(5,110) IFILE
110 FORBAT(AIO)
OPEN INPOT FILE

0PEN(0NIT=21,ACCESS='SEQIN',FILE=IFILE)
BEAD INVENTORY TITLE FROB FILE

HEAD(21,120) TITLE
120 FORBAT(14A5)
WRITE TITLE ON TEBBINAL

WRITE (5, 130) TITLE
130 FORBATCO', 14A5/)

SVFRAH=FRAaE
HAVE OSER SPECIFY SAHPLING HETHOD

IF(III.GT. 1) BEPEAT=' OR SABE'
140 WRITE(5,150) REPEAT
150 FOBBATCOSPECIFY SABPLING HETHOD'/' TYPE POINT OR PL0T',A8)

READ(5,240) FBABE

IF(FRAHE. EQ.' SABE'. AND. III. GT.1) 160 , 170

BAINOOIO
BAINa02a
IIAIN0030

♦♦♦♦•BAINOOHO
HAIN0050
HAIN0060
BAIN0070
HAIH0080
HAIN0090
HAINOIOO
HAIH0110

♦•***BAIN0120
BAIN0130
HAIN0149
HAlei0150
BAIN0160
HAIH0170
HAIH0180
BAIN0190
HAIN0230
HAIN0210
BAIN9220
HAIN0230
BAIH0240
HAIN0250
HAIN0260
NAIN0270
HAIN0280
BAIN0230
aAIN0300
BAIN0310
BAtN032Q
HAIN0330
HAIN03II0
HAIN0350
BAIN0360
HAIN0370
aAIN0380
aAIN0390
HAIN0400
BAIN3419
BAIN0420
BAIN0430
BAIN044Q
HAIN0450
BAIN04&0
BAIN0470
BAINa43Q
aAIN0490
BAIN0500
BAINQ510
HAIN0520
BAIN0530
HAIN0540
BAIN0550
aAIN0560
HAIN0570
aAIN0580
aAIN05)Q
BAIN0600
BAIN0610
HAIN0620
NAIN0630
NAINOenO
aAIH0650
HAIN0660
BAIN0670
BAIN0680
BAIN063O
HAIN0700
BAIN0710
HAIH0720
HAIN0730
aAIN0740
BAIN0750
BAIN0760
BAINa770
HAIN0780

-67-

169 FEAHE=SVFEAH

GOTO 260
173 IF (FRiHE.NE-« POINT* -AND. FBAflE. HE. 'PLOT') GOTO 1U0

IF{FBAHE. EQ. 'POINT') S AHHTD= 'B. A. F. •

IF(FEAaE-EQ. 'PLOT') SABJITD=' PLOT SIZE'
HAVE USER SPECIFY BAF OB PLOT SIZE

HBITE(5,180) SAHHTD
180 FOBHATCOENTEB ' , A9)

READ(5,*) BAF
ASK USEE IF HEIGHTS ABE IN LOGS
190 HBITE (5,200)

200 FOaaAT('0ARE HEIGHTS IN FEET OB SIXTEEN FT. LOGS?'/
1' TYPE FEET OE LOGS')

BEAD(5,210) LOSANS
210 FOBHAT(A(t)

IF(LOGANS.NE. 'FEET'. AND.LOGANS.NE. 'LOGS') GOTO 190
ASK OSES IF HEIGHTS ABE HEBCHANTABLE OE TOTAL
220 WBITE(5,230)

230 FOBHAT('0ARE HEIGHTS TOTAL OE MEfiCH ANTABLE?'/
1' TYPE TOTAL, HEBCH, OB flIXED')

EEAD(5,240) TOPANS
2U0 F0BBAT(A5)

IF (LOGANS. EQ. 'FEET') H0LT=1

IF(LOGAHS.EQ. 'LOGS') JJ0LT=16

IF (TOPANS- NE. ' TOT AL' . AND. TOPANS- NE. • HEECH' . AND. TOPANS. NE. 'MIXED')
1 GOTO 220
ASK OSES IF BEGBESSION BOOTIBE IS NEEDED
250 HEITE(5,260)
260 FOBBATCOHEBE THE HEIGHTS OF ALL TOLOBE TBEES BEASOEED?')

BEAD(5,270) BEGANS
270 F0EHAT(A3)

IF (BEGANS. NE. 'YES' .AND. EEGANS- NE. • NO') GOTO 250
ASK USEE FOB ACBEAGE
280 HBITE(5,290)
290 FOBHATCOENTEB ACBEAGE OF TEACT. ' )

BEAD(5,*) ACEES

IF NEEDED BON THE HEIGHT BEGBESSION BOOTINE

*** **4>4< ****«**•«««**** ******««4i**4i 41 4<*** ****** *4>*4<4i****4>*4i*4i«

IF(BEGANS. E0.'HO') CALL LIHEEG (IFILE, TOPANS)

******************************************

BEAD IN TBEE DATA AND PBOCESS
******************************************

300 CALL GUTS

**********************************
CALL THE OOTPOT EOOTINE
**********************************

CALL 00TP1

****************************************

ACCOHOLATE TOTALS FBOB EACH STBATA
****************************************

IFIISTBAT.GT. 1) CALL STBAT (ACEES, CLASIZ)

310 CONTINUE

END PEOGRAB CONTROL LOOP

*************************************
CALCULATE COHBINED STEATA STATISTICS
*************************************

IF(ISTEAT. GT. 1) CALL STBAT2

********************************
OOTPOT STBATA COBBINED TOTALS
********************************

IF(ISTEAr.GT. 1) CALL O0TP2

EHD

BAIN0790
HAIN0800
BAINOBIO
HAIN0820
flAIN0830
BAINOSltO
HAIN0850
BAINOSSO
aAIN0870
HAIN0880
flAIN0890
HAIN0900
BAIN0910
BAIN0920
BAIN0930
BAIN09ltO
HAIN0950
BAIN09Sa
HAIN0970
HAIN0980
BAIN099a
HAIN1000
BAIHIOIO
HAIN1020
HAIN1030
BAIN10II0
BAINIOSO
HAINIOSO
HAIN1070
BAIN1080
HAIN1090
HAIN1100
BAIN1110
HAIN1120
BAIN1130
HAINIIUO
HAIN1150
BAIN1160
BAIN1170
HAIN11S0
flAINl190
aAIHl200
flAIN1210
HAIN1220
BAIN123Q
BAIN1240
BAIN1250
BAIN125Q
HAIN1270
HAIN1280
aAIN1290
BAIN1300
BAIN1310
aAINl320
aAIN1330
aAIN13l»0
BAINISSO
BAIN1360
aAIN1370
aRIH1380
aAIN1390
BAIN1430
BAINKtlO
HAIN1420
BAXNIItSO
BAINItlO
BAIN1450
flAIN1460
aAIN1<«70
BAIN11130
BAIN1490
HAIN1500
BAINISIO
aAIN1520
UAIN1530
UAIN1540
BAIN1550
HAIN1560

-68-

C FBHC0013

SOBROOTIBE FRUCLS FBBC0020

C FFMC0030

C*********************************************************** ******* *****FB.tlCQOHO

C PEBC0050

C THIS SOBPEOGRAH READS IN AN AREA? OF FOEH CLASSES FEOH A FEnC0050

C OSER SPECIFIED DISK FILE. THIS DISK FILE IS FREE FOBHATTED FRHC0070

C OSING LIST DIRECT ED INPUT. THE USEE IS GIVEN THE OPTION FEHCOOSO

C TO LIST THE FOBB CLASSES AFTER THEY AEE READ. FEMC0090

C FEMC0 190

C** ***«******«>******•***********«********«****«************«***** *******FRHCO 110

C FRHC0120

INTEGER FC(32) ,SPECIE(32) PBMC3130

DOOBLE PRECISION FCFILE FRBC0140

REAL fiEGCOF(33) ,BEGINT(33) ,SEBEG(35) .BSQB (33) ,HERNY(35) FBHCOISO

INTEGEB aiNDBH(32) ,NTR(35) ,TOT»TR(35) FRHC0150

COMaON/SEG/REGCOF,BEGINT,SPECIE,HIHDBH,SEREG,NTR, RSQR.HEANY FRI1C0170

I.TOTNTB FBHC0130

COHBON/FORMCL/FC PRaC0190

C HAVE OSER SPECIFY FOBH CLASS FILE NABE FBHC0230

«EITE(5,10) FBHC0210

10 FOEBAT('0 ENTER FOEH CLASS FILE NAME') FEI1C0220

Z BEAD FOBM CLASSES FBOH FILE FEMC0230

BEAD(5,20) FCFILE FRMC02aO

20 FORHAT(AIO) FRHC0250

OPEN(FILE=FCFILE,ACCESS=«SEQIN',0HIT=1) PRaC0 260

READ(1,30) FC FRWC0270

30 FOBHAT(8{I2,X)) FBHC0230

C ASK USEE IF HE WANTS TO LIST FOEH CLASSES JUST BEAD FBBC0290

UO BEITE{5,50) FEHC0300

50 FORHAT('0DO YOO WANT A LIST OF THE FOEH CLASSES JOST READ?') FBaC0310

READ(5,60) LSTANS FBaC0320

60 F0RHAT(A3) FBBC0330

IF(LSTANS.NE. ■ YES' .AND.LSTANS. NE. 'NO') GOTO HO FBHC0340

IF (LSTANS. EQ. 'NO') EETOEN FEHC0350

WRITE(5,70) (SPECI2(I) ,FC(I) ,1=1,32) FRaC0350

70 FOBBAT('0',8(/' •,4(A2,' = ',I2,5X))) FRaC0370

EETOBN FEHC0330

END FBHC0390

-69-

SOBROOTIHE LIHREG (IPILE.TOPfiNS)
C

C

THIS SOBPBOGRAH READS IB DATA FE3H ALL THE SAHPLE TREES

IT FINDS THE TREES HITH HEIGHT HEASOREMENTS

IT THEN FITS A REGRESSION HODEL FOR EACH SPECIE SAHPLED

HODEL: HEIGHT = BO + B 1 * 1/DBH

C

c
c
c

c
c

DOUBLE PRECISION IFILE

REAL SEREG(35) ,RSQR(33) ,MSE,TYPSSE (3)

REAL VARX(33) , V ARY (33) , COVXY (33) , REGCOF (33) ,REGINT{33)
REAL X(33),Y(33) ,XS0a(33) ,YS0B(35) ,XS0HSQ(33) ,YS0HSQ(33)
REAL XYS0H(33) , BEANX (33) , tlEANY (35) , HEIGHT, HNYCOR (3)
INTEGER SPECIE(32),NTR(35) , POINT, TREES, DBH, SPEC .HINDBH (32)
INTEGER TOTNTR(35) ,N0IiREG(3) ,NTRC0E(3) ,TOPANS
COHHON/REG/REG CO F, REG INT, SPECIE, BIN DBH,SE8EG,NTR,BSQB,HEA NY
1,T0TNTR
Z SET HINiaOB DBH'S TO 8
DO 10 ISPEC=1,32
10 aiNDBH(ISPEC) =8

ZERO ALL ACCOHOLATOES

DO 20 JTyPE=1,3

NTRCOR(JrYPE)=0

HNYCOR (JTYPE)=0

NOHREG(jrYPE)=0
20 TYPSSE(JTYPE)=0
DO 30 ISPEC=1,35

NTR (ISPEC) =0 ; YSOM (ISPEC) =0

SEREG(ISPEC) =0; BEANY (ISPEC) =0

IF(ISPEC.GT.32) GOTO 30

XSOa(ISPEC) =0;XSOHSQ(ISPEC)=0

¥SOBSQ(ISPEC)=0:XYSUB(ISPEC) =0

R EGCOF( ISPEC) =0 ;EEGINT (ISPEC) =0 ; RSQR (ISPEC) =0
30 CONTINUE

INPUT POINT INFOEBATION

«0 READ(21,50,END=130) POINT, TREES
50 FORBAT(I2,X,I2)
CHECK TO SEE IF DATA IS IN SOHHAEY FORB
IF(POINT. NE.-9) GOTO 60
TREES=1000000
INPUT INDIVIDUAL TREE DATA
60 LOOP=0
70 IF (LOOP. GE. TREES) GOTO «0

READ{21,80, END=«0) SPEC, DBH, HEIGHT, lEEP
80 FORMAT(A2,X,I2,1X,F3.0,I2)
IF(IREP.EQ.O) IEEP=1
LOOP=LOOP+IREP
CHECK TO SEE IP TREE IS A CULL
IF (HEIGHT. EQ. 999) GOTO 70
CHECK DBH

IF(DBH. LT. 1) DBH=1
IF(DBH.GT.«0) DBH=<*0
CHECK HEIGHT

IP(HEIGHT-LT.O.) HEIGHT=0.

IF (TOPAHS„ EQ. • FEET' . AND. HEIGHT. GT. 200. ) HEIGHT=200.
IF (TOPANS. EQ. ' LOGS •• AND. HEIGHT. GT. 7.) HEIGHT=7.
CHECK TO SEE IF HEIGHT OF THIS TREE WAS BEASURED

IF(HEIGHr. EQ. OJ GOTO 70
DECODE SPECIE CODES
DO 90 ISPEC=1,32
90 IF(SPEC.EQ. SPECIE (ISPEC)) GOTO 110
PRINT EEEOE BESSAGE

WBITE(5,100) SPEC, POINT
100 FOEBAT(«0***ERROR*** ILLEGAL SPECIE CODE: ',A2,' AT POINT:
l/'ODATA FOR THIS TREE IGNORED! EXECUTION CONTINUES.')
GOTO 70
DETERBINE HINIBOB BEECHANTABLE DBH
110 IF(DBH,LT.BINDBH(ISPEC) ) flIHDBH (IS PEC) =DBH

LNBG0010
LNBG0020
LNRG0030

♦»»»»LBBGOO<»0
LNRG0050
LNRSOOSO
LNBG0070
LNBG0030
LNSG0090
LNRGOIOO
LNBG0110

♦•♦♦♦LB8G0120
LKRG0130
LNBG01U0
LNRG0150
LNBG01S0
LNBG0170
LNRG0130
LNBG0190
LNRG0230
LNBG0210
LNSG0220
LNRG023a
LNBG02'tO
LKBG0259
LNBG0260
LNBG0270
LNBG0280
LNBG02J0
LNBG0300
LNBG0319
LNBG0320
LNBG0330
LNBG0 3ltO
LNEG0350
LNBG0360
LNBG0370
LNBG0330
LNBG0390
LNBGOtOO
LNBG0410
LNBG0(»20
LNBSOU30
LNBG0440
LNBG01I50
LNRG0t60
LNBG0470
LNRGOaSO
LNRG0490
LNB60530
LNBG0510
LNBG0523
LNRG0530
LBBGOSUO
LNBGOSSO
LNBG0560
LHBG0570
LNRG0580
LNBG0590
LNBG0600
LNBG0610
LNBG0620
LHBG0630
LNRG06t»0
LNRG0650
LNRG0660
LNRG0670
LNBG06aO
LNRG0690
LNBG3700
LNBG0710
,12, LNBG0720
LNRG0730
LBBG07ltO
LNEG0750
LNBG0760

-70-

C COnPOTE ACCOHOLATED SOHS FOB INDIVIDOAL SPECIE

ITYPE=33

IF(ISPEC. GE. 12) ITYPE=3IJ

DO 120 ID0HMY=1 ,IREP
X(ISPEC) = 1.0/DBH
T(ISPEC) =HEIGHT
NTB(ISPEC)=NTB(ISPEq +1
BTB(ITYPE) =NTR(ITyPE) +1
NTB{ 35 )=NTE( 35 ) +1
XSOJ!(ISPEC) =XSOH(ISPEC) *t (ISPEC)
TSO H (IS PEC) = Y SO H ( IS PEC) ♦Y ( ISP EC)
YSDH(ITyPE) =YSOH{ITYPE)*Y (ISPEC)
¥SOa( 35 )=YSOH( 35 )+Y(ISPEC)
XSOMSQ (ISPEC) =XSOMSQ (ISPEC) +X (ISPEC) **2
YSOHSQ(ISPEC) =YSOaSQ(ISPEC) ^Y (ISPEC)**2
XYSOM(ISPEC)=XYSOa(ISPEC) +X(ISPEC) *Y(ISPEC)

120

130

COHTINOE
GOTO 70
CONTINUE

COUPOTE REGRESSION PABAHETEBS FOR ALL SPECIES

DO 170 ISPEC=1,32
JTYPE=1

IF(ISPEC-GE.12) JTYPE=2
IF(NTR(ISPEC) .LT. 1) GOTO 170
flEaNY(ISPEC) =YSOH (ISPEC) /HTR(ISPEC)
C CHECK TO SEE IF ENOOGH TREES EXIST FDR REGRESSION

IF(NTR(ISPEC) .GE.3) GOTO 140
C CALOLATE CORRECTION SOHS FOR SPECIES WHERE NO REGRESSION HAS BUN
HHYCOR (jrYPE) =BNYCOB (JTYPE) +YSOM (ISPEC)
NTHCOR(JTYPE) =NTRCOR (JTYPE) +NTR (ISPEC)
GOTO 170
1«0 NOaREG (JTYPE) =N0aBEG (JTYPE) +1

aEANX (ISPEC) =XS0a (ISPEC) /NTB(ISPEC)
VARX (ISPEC) = ( XSOaSQ (ISPEC) -

(XSOa{ISPEC) **2/NTB (ISPEC) )) /(NTB (ISPEC) -1)
VARY (ISPEC)=(YSOBSQ(ISPEC) -

(YSOa(ISPEC) ♦♦2/HTB (ISPEC) ))/(NTB (ISPEC) -1)
COVXY (ISPEC) = (XYSOB (ISPEC) -

XS OH (ISPEC) *YS0a (ISPEC) /NTH (ISPEC) ) /(NT R (IS PEC) -1)
TO SEE IF VARIANCE IS ZERO

IF (VARX (ISPEC) -GT. 0.00000 LAND. VARY (ISPEC). GT.O.OOOOOI)
GOTO 150

REG INT (ISPEC) =HEANY (ISPEC)
GOTO 160
150 REGCOF(ISPEC)=COVXY (ISPEC) /VARX (ISPEC)

REGINT(ISPEC) =MEANY (ISPEC) -REGCOF (ISPEC) *H£ANX (ISPEC)
160 IF(VABY (ISPEC). LT. 0.0000001) GOTO 170

SSE = YSOHSe (ISPEC) -REGINT (ISPEC) ♦YSOa(ISPEq -
1 BEGCOF(ISPEC) *XYS0H (ISPEC)

TYPSSE(jrYPE) =TYPSSE(JTYPE) ♦SSE
aSE=SSE/(NTB(ISPEC)-2)

SEBEG(ISPEC)=SQBT(HSE/NTB(ISPEC) ) * 100/HEANY (I SPEC)
RSQR (ISPEC) =100* (COVXY (ISPEC) *»2) / (VA RX (ISPEC) ♦VABY (ISPEC))
CONTINUE

1

1

CHECK

1

170

COHPOTE STANDARD ERROR FOB SOFTWOOD, HAR DHOOD, ALL SPECIES

NTRC0B(3) =NTRC0B(1) +NTBC0B(2)
HNYC0E{3) =aNYC0B(1) +aHYC0E(2)
TYPSSE(3) =TYPSSE(1) ■••TYPSSE(2)
NUHREG(3)=N0HREG(1) +N0HBEG(2)
DO 180 JTYPE=1,3

ISPEC=JTYPE+32

IF (NTR (ISPEC) .LE.O) GOTO 180

HEANY (ISPEC) =YSOa(ISPEC) /NTR(ISPEC)

CORNTB=NTR(ISPEC) -NTBCOE (JTYPE)

IF(COBNTB.LE. 0) GOTO 180

COBHN¥= (YSOH(ISPEC) -HNYCOB (JTYPE) ) /CORNTB

HSE=rYPSSE(JrYPE) /(CORNTB- 2*N0aEEG(JTYPE) )

SEREG (ISPEC) =SQRT(aSE/CORNTB) ♦lOO/CORBNY
180 CONTINUE

CLOSE (0BIT=21)

0PEN(UNIT=21,ACCESS=« SEQIH* ,FILE=IFILE)

SKIP RECORD 21

RETURN

END

LNRG0770
LNBG0780
LNBG0790
LNBG0800
LKRGOSIO
LNRG0820
LNBG0B3a
LNBG08U0
LNBG0850
LNRG0860
LNRG0870
LNRG0880
LNRG0890
LNRG0900
LNRG0910
LNRG0920
LNR60930
LKBG09l«0
LNRG0950
LNRG0960
LNRG0970
LNRG0990
LNRG0990
LNRG1030
LNRG1010
LNBG1020
LNRG1030
LNRGIOHO
LNRG10S0
LNRG1060
LNRG1070
LHRG1080
LHRG1090
LNRG1100
LNRG1110
LNRG1120
LKRG1130
LNRGIIitO
LNRG1150
LNRG1160
LNRG1170
LNRG1130
LtlRG1190
LNRG1200
LKRG1210
LUeG1220
LNBG1230
LNRG12U0
LNEG1250
LNBG1260
LNRG1270
LNRG1230
LNRG1290
LNRG1300
LNRG1310
LNBG1320
LNBG1330
LNRG1340
LNBG1350
LNRG1360
LNfiG1370
LNRG1380
LNBG1390
LNRGItOO
LNRGIltIO
LNBG1420
LNRG1(t30
LNRGIUaO
LNRG1450
LNBG1460
LNRG1470
LNRGIdSO
LNRG1<(90
LNRGISOO
LNRG1510
LNRal520
LNRG1530
LNRG1540

-71-

c

c*

c

c

c

c

c

c

c

c

c

c*

SOBBOOTIHE GOTS
THIS SOBPBOGBAH IS THE 'GOTS' OF INVEST.

GOTS PEBFOBaS THE FOLLOWING:

1. INPUTS PLOT AND TREE DATi

2. CALCOLATES BASAL AREA, TREES PER ACRE AND ALL VOLOHES

3. CALCULATES VABIAHCES AND STAHDABD ERRORS

**«««***««**«4< **•*«*«*********«*«** 4i**4i***4i*********4<«* **********

DOUBLE PRECISION SAatlTD.IPILE

INTEGER REGANS, ISTART (2), IST0P(2) , DBH.SPEC, TOTNTB (35)

INTEGER SPECIE (32) , POINT, TREES, TBECNT (UO, 32)

INTEGER CLASIZ,TITLE(14) .MOLT.TOPANS.HINDBH (32) .NrR(35)

REAL VABBRD(35) ,VABP0L(35) ,SEREG(35) ,BSQR (33) ,HEAN¥(35)

REAL BEGINT(33) ,REGCOF(33) ,SECDB(35) ,SEBRD(35)

REAL C0BPT(35) ,BRDPT(35) ,P0LPT(35) ,C0BSQ(35) , BRDSQ(35)

REAL BASPAC(U0,35) , BASTOT(35) , rREPAC(«»0,35) ,TBETOT(35)

REAL CBVLSa(l»0,35) ,CBTOT(35) ,BDVLSM (40 , 35) , BDTOT (35) , HEIGHT

REAL SEP0L(35) , P0LSQ(35) , PLVLSH (40, 35) ,PLT0T(35)

CO HHON/COaUNE/CLASIZ, TITLE, HOLT, PR0BLV,BAF,T3PANS, BEG ANS,A:RE

1,FEAHE, SAHHTD, LOGANS, IPILE
C0HH0N/REG/REGC0F,BEGINT, SPECIE,HINDBH,SEREG,NTR,RSCB,HEAHf

I.TOTNTR
COBHON/OOTVBASPAC, BASrOT,BDTOT,BDVLSa,CBTOT,CBVLSa,TREPAC,TR

1,SECnB,SEBBD,lJPHTS,PLVLSB,PLTOT,SEP0L,FINP0P, VARBRD,VABPIJL
DATA (SPECIE(I) , 1=1, 32) /• HP' , • BP* , ' PP* » ' Ha« , ' BF« , 'TA* , 'SP' , •

1,«SC«,«CE','0S« ,«SM«,«Ba«, 'HA', 'HI'.'AS'j'YP', 'BC, 'HB'

2,'lfB','SB', 'GB' ,'BE«, 'BH* ,'R0«, 'HO' ,'BO','EL', 'GO'

3,'HI','HH', 'OH'/
DATA ISTART(I) /1/,ISTART(2) /12/,IST0P(1) /1 1/,IST0P(2| /32/

ZERO ALL ACCOaULATORS

DO 10 ISPEC=1,35

BASTOT(ISPEC)=0:TBET0T(ISPEC) =0; CBTOT (ISPEC) =0; BDTOT (ISPEC)

PLTOT(ISPEC) =0:TOTHTR(ISPEC) =0

COBSQ(ISPEC) =0 ; BRDSQ (ISPEC) =0 ; POLSQ (ISPEC) =0

SECUB(ISPEC)=0:SEBED(ISPEC) =0 ; SEP OL (ISPEC) =0

DO 10 DBH=CLASIZ,40,CLASIZ

TREPAC(DBH, ISPEC) =0; BASP AC ( DBH, ISPEC) =0
CBVLSB(DBH, ISPEC) =0;BDVLSH(DBH, IS PEq=0
PLVLSB(DBH, ISPEC) =0

IF (IS PEC. LT- 33) TRECNT (DBH.ISPEC) =0
10 CONTINUE
NP»TS=0

O***** ******* *****«**************#******4>*
BEAD IN ALL TREE DATA AND ACCUBULATE SUBS
******************************************

INPUT POINT OR PLOT INFOBMATION

NS

20 READ(21,30,END=140) POINT, TREES
30 FORHAT (12,13)
CHECK TO SEE IF DATA IS IN SOHBARY FORB
IF (POINT- NE- -9) GOTO 40
NPNTS=TREES
TREES=1000000
INPUT INDIVIDUAL TREE DATA
40 LOOP=0
50 IF (LOOP. GE. TREES) GOTO 120

BEAD(21,60,END=120) S PEC, DEH, PROD, HEIGHT, IREP, PCTS AW
60 FORHAT(A2,X,I2,A1,F3. 0, I2,F3-2)
IF(IREP.EQ.O) IREP=1
LOOP=LOOP*IEEP
DECODE SPECIE CODES
DO 70 1=1,32
70 IP(SPEC.EQ.SPECIE(I)) GOTO 90

PRINT EBBOB HESSAGE

WRITE(5,80) SPEC, POINT
80 FORBAT('0***EBROB«** ILLEGAL SPECIE CODE: •,A2,' AT POINT:
V'ODATA FOR THIS THEE IGNORED! EXECUTION CONTINUES.')
GOTO 50

GUTS0010
SUTS0020
GUTS0030

***»*aOIS0040
G0TS0050
GOTS0060
GOTS0070
G0TS0080
GUTS0090
GUTS0100
GUIS0110
GnTS0120
G0TS0130

*****GUTS0140
GOrS0150
GUTSOI^O
GOTS0170
GOTSOISO
GUTS0 190
GUTS0200
GnTS0210
GUTS0220
G0TS0230
GaTS0240

S GOTS0250
GUTS0260
GUTS0270
GnTS0280

ETOI GUTS0290
GUIS0300
GOTS0310
GUTS0320
GOTS0330
GOTS0340
G0TS0350
GaTS03S0
GUTS0 37O
GUTS0330
GUTS0390
GOIS0400
GnTS0410
G0TS0420
SnTS0430
GOTS0440
GUTS0450
GUTS0460
GUTS0470
SnTS0480
G0TS0490
SUTS0500
GUTS0510
GUIS0520
GOrS0530
G0rS0540
GUTS0550
GnTS0560
30IS0570
G0TS0580
GOTS0590
GUIS0600
SUTSOeiO
GnTS0620
GUTS0630

soTsoeuo

GOrS0650
SOTS0660
GUTS0670
GnTS0680
GUTS069a
SUTS07aO
G0TS0710
GOTS0720
GOTS0730
GUrS07l»0
G0TS0750
GUTS0760
GOrS0770
GUTS0780

=0

.12,

-72-

C IF

C

c
c

c

CHECK DBH TO SEE IF HITHIH LEGAL LIHITS
90 IF(DBH-LT-CLASIZ) DBH=CLRSIZ

IF (DBH- GT. 40) DBB=aO
T&LLT TBEES FOB EACH SPECIES

ISPEC=I

T8ECNT (DBH,ISPEC)=TEECNT(DBH, ISPEC) ♦IBEP
CHECK TO SEE IF DIAHETER IS ABOVE aiHIMOM LIBIT

IF(DBH.LT.HIND6H(ISPEq . AND. PEGANS. EQ. • NO • ) GO
aCCOHOLATE HOHBEE OF HOH-COLL TBEES ABOVE BIN DBH S

TOTHrE(ISPEq=TOTNTB(ISPEC) +IBEP

IF (ISPEC. LE- 11) TOTNTE(33) =TOTHTE(33) -HEEP

IF (ISP EC. GE. 12) TOTNTE(3a)=TOTlITE(3a)+IBEP

T0TNTE(35)=T0rNTE(35) +IBEP
HEIGHT HAS NOT HEASUBED GENEBATE HEIGHT FEOH BEG

IF(HEIGHT.EQ.O. AND. REGANS. EQ.« NO' ) HEIGHr=BEGI

IF (HEIGHT. EQ.O) GOTO 50
COHPOTE VOLOHES
CHECK TO SEE IF TREE A COLL

IF(HEIGHT.EQ.999) GOTO 50

HEIGHT= HEIGHT *HULT
CHECK HEIGHT TO SEE IF BITHIN LEGAL LIHITS

IF(HEIGHT. LT. 0) HEIGHT=0

IF (BOLT. EO. 1. AND.HEIGHT.GT. 200) HEIGHT=200.

IF(flOLT.EQ. 16. AND. HEIGHT. GT. 128) HEIGHT=128

IF (TOPANS.EQ. TOTAL') CALL VOLTOT (DBH, HEIGHT, I

IF(TOPANS-EQ. 'BERCH') CALL VOLHER(DBH, HEIGHT, I

IF(TOPANS„EQ.'BIXED'. AND. ISPEC- GT. 11) CALL

1ISPEC,C0BVOL, BDVOL)

IF (TOPANS. EQ. ' flIXED • . AND. ISPEC. L£. 1 1 )
1 1S PEC, CO BVOL, BDVOL)

POLVOL=0
CHECK TO SEE IF % SAHLOG OPTION OSED

IF(PBOD. NE. 'X') GOTO 100
CHECK PCTSAW TO SEE IF WITHIN LEGAL LIHITS

IF(PCTSAH.LT-O) PCTSAH=0.0

IF (PCTSAW. GT. 1.0) PCTSAW=1-0

BDVOL =PCTSAW*BDVOL

P0LVOL= (1. -PCTSAW) ♦COBVOL

C0BV0L=PCTSAW*C0BVOL

GOTO 110
CHECK TO SEE IF TBEE IS POLPWOOD
100 IF(PBOD. NE. 'P') GOTO 110

P0LVOL=C0BVOL;C0BVOL=0:BDVOL=0
110 CBVLSHiDBH,ISPEC)=CBVLSfl(DBH,ISPEC)+C0BVOL*IEEP

BDVLSH(DBH, ISPEC) =BDVLSH(DBH, ISPEC) +BDVOL*IEEP

PLVLSa(DBH, ISPEC) =PLVLSB(DBH, ISPEC) +PDLVOL*IBEP

TO 50
ABPLED

BESSION

NI (I)>BEGCOF

SPEC, COB VOL,

SPEC, COBVOL,

VOLHER(DBH,

CALL VOLTOT (DBH,

ACCOBOLATE SOUS FOB VOLOBE STANDAED ERRORS

CHECK TO SEE IF DATA IS
IF (POINT. EQ. -9) GOTO 50

IN SOHHAEY FOEB

IF(FBABE-EQ. 'PLOT') FACTOB= 1.0/BAF

IF (FRABE. EQ. ' POINT' ) FACTOE=BAF/(. 005a5a*DBH**

COBPT (ISPEC) =C0BPT (ISPEC) +C0BVOL*IEEP*FACTOE

BRDPT(ISPEC)=BBDPT(ISPEC) ♦BDVOL ♦IEEP*FACTOB

POLPT(ISPEC)=P0LPT(ISPEC)+PULVOL»IEEP*FACTOB

IF(ISPEC„LE.11) ITYPE=33

IF(ISPEC.GT.II) IT?PE=Jl»

COBPT (ITrPE) = COBPT(ITTPE) +C0BVOL*IBEP*FACT0R

BBDPT(ITTPE)=BBDPT(ITYPE) +BDVOL«IEEP*FACTOE

PDLPT(ITJPE)=POLPT(ITYPE) ♦P0LVOL*IEEP*FACTOE

COBPT (35) =C0BPr (35) ♦C0BVOL*IEEP*FACTOR

BBDPT(35) =BRDPT(35) ♦BDVOL*IEEP»FACTOE

POLPT (35) =POLPr (35) ♦P0LVOL*IEEP*FACTOR

2)

GOTO 50
120 CONTINOE
IF IN SOHHART BODE DON'T ACCOHOLATE

IF (POINT. EQ. -9) GOTO 1U0
TALLY NOHBER OF POINTS SABPLED

NPNTS=HPNrS+1

VOLOHES OB INCREBENT NPNTS

G0TS0790

Gorsosoo

GOTS081Q
GDTS0820
GDTS0830
GOTSOBItO
G0TS0850
GOTSQBSa
GOTS0870
G0TS0880
60TS0890
60TS0900
GnTS0910
GOTS0920
(I)/DBHG0rS0930
SOTS09li0
G0rS0950
GOIS0960
G0rS0970
GOTS0980
GOrS0990
G0TS1000
GOTS1010
G0TS1020
BDVOL) aOTS1030
BDVOL) GOTS10H0
HEIGHT, GOTS1050
GOTS1060
HEIGHT, GOTS 1070
G0TS1080
GOTS1090
G0TS1100
GUTS1110
GOTS1120
G0TS1130

GOTsma

GOTS1150
GQTSIIBO
G0rS1170
GOTS1130
G0TS1190
S0TS1200
GOTS1210
G0TS1220
GOrS1230
GOrS12U0
G0TS1250
GDTS1260
G0rS1270
GaTS1280
G0TS1290

sorsi300

GOTS1310
G0TS1320
GOTS1330
G0TS1340
G0TS1350

GOTsneo

G0TS1370
GDTS1330
G0TS1390
5UTS1400
GOTSiaiO
G0TS1U20
GOTS1U30
GOIS1440
GOTS1450
50151460

Gorsia70

G0TS1430
G0TS1490
GOTS1500

-73-

ACCOHOLATE SOHS OP SQOAEES FOB VOLOHES
DO 130 ISPEC=1,35

COBSQ(ISPEC) =COBSQ(ISPEC) -fCOBPT (ISPEC) **2
BEDSQ(ISPEC) =BRDSQ{ISPEC) +BEDPT (ISPEC) **2
POLSQ(ISPEC) =POLSQ(ISPEC) +POLPT (ISPEC)**2
COBPT (ISPEC) =0
POI.PT(ISPEC) =0
130 BBDPT (ISPEC) =0

GOTO 20
1U0 CONTINOE

IF ONLY ONE POI HT SAHPLED, SET IN SOHHARY BODE
IP(HPNTS. EQ. 1) P0INT=-9

DETEBMINE SAHPLING METHOD
***««******«***«*«« ««>*• *******

IF(FRAHE-EQ.'PLOT«) GOTO 170

***************************************************
POINT SAHPLING

CALCULATE PER ACRE ESTIMATES FOR IHDIVIDOAL SPECIES
***************************************************

DO 160 ISPEC=1,32

DO 150 DBH=CLASIZ,l»0,CLASIZ
C CHECK TO SEE IF ANY TREES TALLIED

IF (TEECNICDBH, ISPEC) .EQ.O) GOTO 150
C BASAL AREA CALCULATIONS

BASPAC(DBH, ISPEC) =TRECNT (DBH, ISPEC) 'BAF/NPNTS
BASrOT (ISPEC) =BASTOT( ISPEC) +BASP AC (DB H, ISPEC)
C TREES PEB ACRE CALCULATIONS

TREPAC(DBH, ISPEC) =BASPAC (DBH, ISPEC) /(. 005l*5a*DBH**2)
TBETOT (I SPEC) =TRETOT (ISPEC) +TREPAC (DBH, ISPEC)
C CUBIC VOLUME CALCULATIONS

CBVLSH (DBH, ISPEC) =CBVLSM (DBH, ISPEC) *BAF/
1 (NPNTS*.005«5<»*DBH**2)

CBTOT (ISPEC) =CBTOT( ISPEC) +CBVLSM (DBH, ISPEC)
C BOARD FOOT VOLUME CALCULATIONS

BDVLSM(DBH, ISPEC) =BDVLSa (DBH, ISPEC) ♦BAF/
1 (NPNTS*.005«54*DBH**2)

BDTOT (ISPEC) =BDTOT (ISPEC) +BDVLSH (DBH, ISPEC)
C PULPWOOD CALCULATIONS

PLVLSM(DBH, ISPEC) =PLVLSfl(DBH, ISPEC) *BAF/
1 (NPNTS*.005a5tt*DBH**2)

PLTOT (I SPEC) =PLTOT (ISPEC) ♦PLVLSM (DBH, ISPEC)
150 CONTINOE
160 CONTINOE
GOTO 200

***************************************************
PLOT SAMPLING

CALCULATE PER ACRE ESTIMATES FOR INDIVIDUAL SPECIES
***************************************************

170 PLTSIZ=BAF

NPLOTS=NPNTS

DO 190 ISPEC=1,32

DO 180 DBH=CLASIZ,l»0,CLASIZ

CHECK TO SEE IF ANY TBEES TALLIED
IF(TRECNr (DBH, ISPEC). EQ-0)

TREES PEB ACRE CALCULATIONS

TBEPAC(DBH, ISPEC) =TRECNT(DB
TRETOT (ISPEC) =TRETOT(ISPEC)

BASAL AREA CALCULATIONS

BASPAC(DBH, ISPEQ =TEEPAC(DB
BASTOT (ISPEC) =BASTOT(ISPEC)

CUBIC VOLUME CALCULATIONS

CBVLSM(DBH, ISPEC) =CBVLSH(DB
CBTOT (I SPEC) =CBTOT (ISPEC) +C

BOARD FOOT VOLUME CALCULATIONS

BDVLSH(DBH, ISPEC) =BDVLSM(DB
BDTOT (ISPEC) =BDTOT( ISPEC) +B

PULPHOOD CALCULATIONS

PLVLSM (DBH, ISPEC) =PLVLSM (DB
PLTOT (I SPEC) =PLTOT (ISPEC) +P

183 CONTINOE

190 CONTINUE

GOTO 180

H,ISPEC) /PLTSIZ/NPLOTS
+TREPAC (DBH, ISPEC)

H, ISPEC) *.0054 5'»*DBH*»2
+BASPAC(DBH, ISPEC)

H, ISPEC) /PLTSIZ/NPLOTS
BVLSH(DBH, ISPEC)

H, ISPEC) /PLTSIZ/NPLOTS
DVLSM(DBH,ISPEC)

H, ISPEC) /PLTSIZ/NPLOTS
LVLSM(DBH, ISPEC)

GOTS1510
GUTS1520
GUTS1530

Gursisuo

GUTS1550
G0TS1560
G0TS1570
GUTSISBO
G0TS1590
GOTS1600

Gursieio

G0rS1620
GOTS1630
G0rS16l»0
GOrS1650

Guisieeo

G0TS1670
GUrsi680
G0TS1690
GOrS1700
G0rS1710
GOrS1720
G0TS1730
GUrS17l»0
GUTS1750
GUrS1760
GUrS1770
G0rS1780
GOrS1790
GOTS1800
G0TS1810
GOTS1820
GUrS1830
GOTS18U0
GOTS1850
G0rS1860
G0TS1870
G0TS1880
GUTS1890
G0TS1900
G0TS1910
GUTS1920
GUTS1930
GOrS19U0
GUTS1950
GUrS1960
G0rS1970
GUrS1980
GUTS1990
GUTS2000
GUTS2010
6UrS2020
GUTS2030
GniS20lt0
G0rS2050
G0TS2060
GUTS2070
G0IS2080
GUTS2090
GUrS2100
GUTS2110
GUTS2120
GOTS2130

Gors2no

G0TS2150
GUTS2160
GUrS2170
GUTS2180
GOTS2190
GOTS2200
GUTS2210
G0TS2220
GOrS2230
G0TS22a0
GUTS2250
G0rS2260
GUTS2270
GUrS2280
3UTS2290
GUIS2300

-74-

CALCULATE TOTALS FOB SOFTWOODS ARD HARDWOODS
^^^^^^*m*************************** **********

200 DO 230 ITIPE=33,34

DO 220 DBH=CLASIZ,UO,CLASIZ

DO 210 ISPEC=ISTABT(ITYPE-32| ,ISTOP

BASPAC(DBH,ITyPE) =BASPAC( DBH, ITYP

TBEPAC(DBB, ITYPE)=TEEPAC(DBH,ITIP

CBVLSH(DBH, ITSPE) =CBVLSH ( DBH, ITYP

BDVLSH(DBH,ITyPE)=BDVLSIl(DBH,ITYP

210 PLVLSH(DBH, ITtPE) =PLVLSH (DBH, ITYP

BASTOT(ITYPE)=BASTOT{ITYPE) ♦BASPACC

TEETOT(ITYPE)=TRETOT(ITYPE) +TEEPAC(

CBTOT (ITYPE) =CBTOT (ITYPE) -t-CBVLSH (DB

BDTOT(ITYPE)=BDTOT(ITYPE)+BDVLSH(DB

223 PLTOT (ITYPE) =PLTOT( ITYPE) ♦PLVLSa{DB

(ITYPE-32)

E) +BASPAC(DBH, ISPEC)

E) ♦TREPaC(DBH, ISPEC)

E) ♦CBVLSH(DBH, ISPEC)

E) +BDVLSa(DBH, ISPEC)

E)+PLVLSa(DBH,ISPEC)

DBH, ITYPE)

DBH, ITYPE)

H, ITYPE)

H, ITYPE)

H, ITYPE)

230 CONTINOE

*4t««4i* «*«*««*««*«««*»****«* ******

CALCULATE TOTALS FOE ALL SPECIES
*********************************

DO 2«0 DBH=CLASIZ,40,CLASIZ

BASPAC(DBH,35)=BASPAC (DBH, 33) ♦BASPiC (DBH, 34)
TBEPAC(DBH, 35)=TBEPAC(DBH, 33) +TEEPAC (DB H, 34)
CBVLSa(DBH,35)=CBVLSH(DBH,33) +CBVLSH (DBH, 34)
BDVLSH(DBH, 35)=BDVLSH(DBH, 33) +BDVLSH (DBH, 34)
PLVLSa(DBH,35)=PLVLSH(DBH,33)*PLVLSB(DBH,34)
240 CONTINUE

BASTOT(35) =BASTOT (33) ♦BASTOT(34 )

TBETOT (35) =TEErOT (33) +TBETOT (34)

CBTOT(35) =CBT0T(33) +CBTOT(34)

BDTOT(35)=BDT0T (33) ♦BDT0T(34)

PLT0T(35) =PLT0T{33) +PLT0T (34)

***************************

CALCULATE STANDARD ERRORS
***************************

FINPOP=1.0

IF(FRAHE.EQ.' PLOT') FINPOP= ( (ACRES/BAF) -NPNTS) /( ACEES/BAF)
IF(FINPOP. LT. 0.0) FINPOP=0.0
CHECK TO SEE IF DATA IS IN SUMHAEY FOEB
IF(POINT. EQ. -9) GOTO 270
DO 260 ISPEC=1,35

IF(TRETOT(ISPEC) .Eg.O) GOTO 260
IF(FBAHE. EQ. • PLOT') GOTO 250

PCTCEZ=100*NPNTS*BASTOT(ISPEC) /TRET OT (ISPEC) /ACEES/BAF
IF(PCTCBZ-GT. 100.0) PCTCRZ = 100.0
FINPOP= (100-PCTCRZ) /100
250 VAECnB= (CUBSQ (ISPEC) -NPNTS* (CBTOT (ISPEC) **2) ) /(NPNTS-I)

VAEBBD(ISPEC) = (BEDSQ(ISPEC)-NPNTS* (BDTOT (ISPEC) **2) ) /(HPNTS-1)
VAEPUL (ISPEC) = (PULSQ (IS PEC) -NPNTS *1PLT0T (ISPEC) •*2) )/(NPNTS-1)
SECOB(ISPEC) =SQBT((VAECUB/BPNTS)»FINPOP)
SEBED(ISPEC)=SQET((VAEBBD(ISPEC) /NPNTS) *FINPOP)
260 SEPUL(ISPEC) =SQRT ( (VARPOL(ISPEC) /NPNTS) *PINPOP)
CHECK TO SEE IF REGRESSION ROUTINE USED
270 IF(EEGANS.EQ. 'YES') RETURN
DO 280 ISPEC=1,35

IF(rOTNTE(ISPEC) . EQ.O) GOTO 280

EEGHGT = PLOAT(TOTNTR (ISPEC) -NTR(IS PEC) ) /FLOAT (TOTNTR (ISPEC) )
SECUB(ISPEC)=SQRT(SEC0B(ISPEC)»*2 *
1 REGHGT»(SEBEG (ISPEC) »CBTOT(ISPEC) /100.) **2)

SEBRD(ISPEC) =SQRT(SEBRD(ISPEC)**2 ♦
1 REGWGT* (SEEEG(ISPEC) *BDTOT( ISPEC) /I 00.) **2)

SEPOL (ISPEC) =SQRT(SEPUL(ISPEC) **2 *
1 EEGWGT*(SEREG(ISPEC)*PLTOT(ISPEC)/100.) ♦♦2)
280 CONTINUE
RETURN
END

GUTS2310
GBTS2320
GUrS2330
GUTS2340
G0rS2350
GaTS2360
GUr32370
GnTS2330
GUrS2390
GUTS2400
GOTS2410
GUrS2420
GUTS2430
GnTS2440
G0rS2450
GUTS24S0
GnTS2470
GUTS24B0
GUTS2490
GUTS2530
GUrS2510
GUTS2520
GUTS2530
G0TS2540
GOrS2553
GOTS2560
GUTS2570
G0rS2580
GnTS2590
GUTS2600
GUTS2610
GUTS2620
GOrS2630
GUTS2640
GUrS2650
GUrS2660
GUrS2670
G0TS2680
GUTS2690
Gnrs2700
GUTS2710
GUrS2720
GUTS2730
GnTS2740
G0TS2750
GUTS2760
GnTS2770
GUrS2780
GUTS2790
GUT528a0
GUTS2810
GOTS2820
G0rS2830
G0TS2840
GOTS2850
GUTS2860
SOrS2870
GnTS2880
GnTS2890
GUTS2900
GOTS2910
GUTS2920
G0TS2930
G0TS2940
G0TS2950
G0rS2960
GUrS2970
GOrS2980
GUTS2990
GOTS3030
GUTS3010
GUTS3020

-75-

C VLMROOIO

SOBBOOTINE VOLHEB (DBH, HEIGHT, ISPEC, COBVOL.BDVOL) VLHR0020

C VLHR0030

C VLBB9050

C THIS SOBPEOGBAM COHPOTES THE HERCHANTABLE COBIC AND BOARD F30r VLHR0060

C VOLOHES OF INDIVIDOAL TREES. THE INPOT PARAHETEES ABE DBH, VLME0070

C MERCHANTABLE HEIGHT, AND SPECIES. VLMBOOSO

C VLHB0090
C THE COBIC FOOT EQUATIONS ARE FROM BAENABD, JOSEPH E. ET AL. 1973. VLHE0100

C FOBEST SOEVET COBIC-FOOT VOLOHE EQUATIONS. OSDA FOBEST SEBV. VLHB0110

C NE. FOREST EXP. STA. RESEARCH NOTE NE-66. VLnR0120

C VL8R0130
C THE BOARD FOOT EQUATIONS ABE FBOH WIANT, H.V. AND CASTENADA, P. 1977 VLHEOIUO
C MESAVAGE AND GIRAED'S VOLUME TABLES FORMULATED. BESOURCE INVENTOBJ VLME0150

C NOTES, USDI, BLB, DENVEE, COLO. VLMR0160

C VLHB0170
C«*«* ♦♦*♦♦»*♦♦♦♦*****♦***♦****♦***♦*»*****♦♦****♦***♦*♦♦*♦♦****♦♦♦♦**♦♦» VLB EO 180

C VLHE0190

REAL HEIGHT, C0BIBT(17) ,C0BSLP(17) ,BDHT VLMEO2D0

INTEGEB DBH,GBP(32) ,D,FC(32) VLMR0210

COBBON/FOEBCL/FC VLaR0220

C VLBR0230

C DEFINE FORB CLASS FOE THE SPECIES VLBE02U3

C VLMR0250

DATA (FC(I) ,I=1,32)/3*76,70,7*76,21*7«»/ VLHE0260

C VLBB0270

C DEFINE VOLUME EQUATION GBOUPING FOB THE SPECIES VLaB0280

C VLBB0290

DATA (GEP(I) , 1= 1,32) /I, 1, 5, U, 3, 5, 2, 5, 5, 6, 1,7, 8, 9, 17,9,8,10 VLBR0 300

1,4*11,12,13,14, 15,14, 15,1 a, 16, 2*17/ VLMB0310

C VLaB0320

C DEFINE THE INTEBCEPTS FOE THE COBIC VOLOBE EQUATIONS VLSB0330

C VLBB0340

DATA (CUBINT(I) ,1=1 , 1 7) /3. 5142,2. 1998 , 1 . 4793, 2. 4784,-. 0496, 1.5817 VLSR0 350

1,1.6823,1. 1763,1.0067,1.1809, 1. 1339, 1-2 851 ,. 8976, 1 .2027,1 . 0009 VLBR0 35 3

2, 1.4438.. 8591/ VLHR0370

C VLBR0380

C DEFINE THE SLOPES FOR THE COBIC VOLOBE EQUATIONS VLMR0390

C VLBR34DD

DATA (COBSLP(I) ,1=1 ,1 7) /236 , 257 , 272, 242,3 03, 259, 310, 310,293 ,2 92 VLBR0410

1,281,326,349,3 01,300,316,309/ VLaR0420

C VLBB0433

C FONCTION TO COBPOTE BOARD FOOT VOLOBES FOR G. F. C 78 VLaR0440

C VLBR0450

HIANT(DBH,BDHT)=(-13.35212*9.58615»BDHT+1.52968*BDHT**2)+ VLaR0460

1 (1. 79620-2. 59995*BDHT-.27465*BDHT**2)*DBH+ VLHR0470

2 (.04482+.45997*BDHT-.00961*BDHT**2) *DBH**2 VLBB0480
C VLaE0490

C DETERBINE BOABD FOOT HEIGHT TO THE HALF-LOG VLMB0530

C VLBEOSIO

BDHT=HEIGHT/16-0 VLBE0520

CBHT=HEIGHT VLaE0530

C VLBB0540

C CALCULATE CUBIC AND BOARD FOOT VOLUBE VLBR0550

C VLMR0560
C0BVOL=COBINT(GEP(ISPEC)) +. 0000 1 ♦COBS LP (GBP (ISPEC) ) ♦CBHr*DBR**2 VLMR0570

BDVOL=HIANr(DBH, BDHT)*(((FC(ISPEC)-78)*.03)+1) VLaR0580

IF(BDHT.LT. 0.5) BDVOL=0 VLMR0590

IF{CUBVOL-LT.O) C0BVOL=0 VLaR0600

IF(BDVOL. LT.O) BDVOL=0 VLBE0610

HETOBN VLBB0620

EUO VLaB0630

-76-

C VLTT0013

SOBEOOTINE VOLTOT (DBH, HEIviHT.ISPEC.COBVOL, BDVOL) VLTT0020

C VLTT0030

C ******************************************************** ************«**VLTT0O(tO

C VLTT0050

C THIS SOBPBOGBAH COHPOTES THE HEECHANTftBLE COBIC AND BOABD FOOT VLTTOOeO

C VOLOHES OF INDIVIDUAL TE EES. THE INPDT PAEABETEES ARE DBH, VLTTOOTO

C TDTAL HEIGHT, AND SPECIES. VLTTOOaO

C VLTT0090

C ALL OF THE EQUATIONS OSED ABE FEOM HONEB, T. G. 1967. STANDAED VOLOBE VLTT0100

C TABLES AND BEECHANTABLE CONVEBSIOH FACTOES FOE THE COHBEECIAL TREE VLTTOIIO

C SPECIES OF CENTEAL AND EASTEBN CANADA. FHB-X-5, FBI, OTTAWA, CANADA VLTT0120

C VLrT0130

C**** ******************************************************* ************VLrTO lltO

C VLTT0150

INTEGEE GBP(32) ,DBH VLTTOieO

REAL C0BINr(16) ,C0BSLP(16) ,ASDBCF(3) ,ASDBBF(3) ,BEDT0P(2) .COBTOP {2 ) VLTTO 170

C VLTT0180

C DEFINE VOLOBE EQOATION GBOOPING FOB THE TBEE SPECIES VLTT0 190

C VLTT0200

DATA (GBP (I), 1=1, 32) /I, 2, 2, 6, a, 2, 3, 3, 2, 5, 2, 10, 10, 16, 15,7, 10, 1 3 VLTT0210

1,8,9,8,8, 12, 11, 3*16, H», 16,16,15,15/ VLTT0220

c vLrro230

C DEFINE THE INTEBCEPTS FOB THE TOTAL COBIC FT VOL EQOATIOH VLTT0240

C VLTT0250

DATA {CO BINT (I) , 1= 1, 16) /- 691 , . 710, 1.226,2. 139 ,4. 167, 1. 112,-. 3 12 VLTT0260

1,2.222,1.««9,1-0i»6, .9H8,.959,-033,.63i»,1.877,1.512/ VLrT0270

C VLrT0280

C DEFINE THE SLOPES FOB THE TOTAL COBIC FT VOL EQUATION VLTT0290

C VLTT0300

DATA (COBSLP(I) , 1= 1, 16) /363- 676,355.623 ,315. 832,30 1. 63U, 2au. 9 06 VLTT0310

1,350.09 2,iJ36.683,300-373, 3ai*.754,383-972,<»01.i»56,33i».829 VLTT0320

2, 393. 366, 4i»0.<»96, 332. 585, 336-509/ VLTT0330

C VLrT03i»0

C DEFINE STOMP HEIGHT AND BEECH TOPS VLTT0350

C VLTT0 350

DATA STOBP/0. 5/ VLTT0370

DATA (CaBTOP(I) ,I=1,2)/a.0,U.0/ VLTT0380

DATA (BEDTOP(I) ,1=1,2) /6.0, 8.0/ VLrT0390

C VLTTOUOO

C DEFINE ADJUSTED SQED DIA EATIO COBIC FT VOL CONVERSION COEFFICIENTS VLTTOIIO

C VLTT0420

DATA (ASDECF(I) ,1=1, 3)/.960«,-. 1660, -.7868/ VLTT0H30

C VLTTOUUO

C DEFINE ADJUSTED SQBD DIA EATIO BOAE D FT VOL CONVERSION COEFFICIENTS VLTTOttSO

C VLTT0460

DATA (ASDSBF(I) ,1=1, 3)/5.«332,- 1.6281,-4.4710/ VLTT0470

C VLTTOUSO

C DETERMINE SB OB HW TYPE VLTT0490

C VLTT0500

IF(ISPEC.LE.II) ITIPE=1 VLTTOSIO

IF(ISPEC. GE. 12) ITYPE=2 VLTr0520

c vLrro530

C COBPOTE TOTAL COBIC FOOT VOLOBE VLrT0540

C VLTTOSSO

TCBVOL=DBH**2/(C0BINT(GBP(ISPEC)) ♦COBSLP (GBP (ISPEC) ) /HEISHT) VLTT0560

C VLrr0570

C COBPOTE BEECH COBIC AND BOABD FT VOLUMES VLTT0580

C VLTr0590

XCUB=( (CUBTOP(ITYPE) /DBH) **2) * (1. O+STOMP/HEIGHT) VLrT0600

XBE&= ((BEDrOP(ITYPE)/DBH) **2) *( 1. 0 + STOMP/HEI3HT) VLTT0610

C0BVOL=TCBVOL*(ASDRCF(1)+ASDECF(2) *XCOB+ASDECF (3) *XC0B**2) VLTT0620

BDVOL =TCBV0L*(ASDBBF(1) ♦ASDEBF (2) •XBED+ASDRBF (3) »XBBD**2) VLrT0630

IF(COBVOL.LT. 0) C0BVOL=0 VLTT0640

IF(BDVOL.LT.O) BDVOL=0 VLTT0650

BETOBN VLTT0660

END VLrT0670

-77-

SOBROUTINE STRAT (ACRES, CLASIZ)
C

Q******1f^***********i***t ************* ******************m**^t.t.ttlr*m****^ti
C

C THIS SOBPROGEAH CALCULATES COMBINED STRATA STATISTICS
C
C**********************************************************************

c

INTEGER CLASIZ

REAL STSECB(35) ,STSEBD(35) .EFDFCB (35) ,EFDPBD(35)

REAL SrVHCB(35) ,STVHBD(35) , DFCB2(35) ,DFBD2(35)

REAL STV11PL(35) ,DFPL2(35) ,PLTST(35) ,STSEPL(35) , EFDFPL (35)

REAL PLVLST(l»0,35)

REAL CBTST(35) , BDTST(35) , CBVLST (1*0, 35) , BDVLST (40, 35)

COHMON/O0T1/BASPAC(«»0,35) ,BASTOT(35) ,BDTOT(35) , BDVLSH (40, 35) ,
1CBTOT(35) ,CBVLSn(40,35) , TREP AC (40, 35) , TRETOT(35)
2,SECOE(35) ,SEBRD(35) , NPNTS, PLVLSH (40, 35) ,PLT0T(35) ,SEP0L(35)
3,FISPOP,VARBED(35) ,VARPOL(35)

COHMON/STOOT/CBTST, BDTST,CBVLSr, BDVLST, STSECB ,STSEBD, EFDFCB
1 , EFDFBD, TA:RES, PLTST, PL VLST ,STS EPL, EFDFPL
C ZERO ACCOBOLATORS AND OTHER VARIABLES

DATA ((CBVLST(J,I) ,BDVLST(J,I) ,PLVLST(J,I) , J=1,40,1)
I.CBTST (I) ,BDTST (I) ,PLTST(I) ,STVHCB(I) ,STVHBD(I) ,STVnPL(I)
2,DFCB2 (I) ,DFBD2 (I) ,DFPL2(I) , EFDFCB (I) ,EFDFBD(I) , EFDFPL (I)
3,STSECB(I) ,STSEBD(I) ,STSEPL(I) , 1=1,35) / 4725»0.0 /

TACRES=TACRES+ ACRES
C

C ACCOaOLATE SOMS FOR VOLOHES
C

DO 20 ISPEC=1,35

00 10 DBH=CLASIZ, 40, CLASIZ

CBVLST (DBH, ISPEC) =CBVLST (DBH, ISPEC) +CBVLSfl (DBH,ISPEC) ♦ACRES
BDVLST(DBH,ISPEC)=BDVLST(DBH,ISPEC) +BDVLSa ( DBH, ISPEC) *ACRES
10 PLVLST{DBH, ISPEC) =PLVLST (DBH, ISPEC) +PLVLSB (DBH, ISPEC) ♦ACRES

CBTST (ISPEC) =CBTST(ISPEC) ♦CBTOT (ISPEC) ♦ACRES
BDTST (ISPEC) =BDTST(ISPEC) ♦BDTOT (ISPEC) ♦ACRES
PLTST (ISPEC) =PLTST(ISPEC) -fPLTOT (ISPEC) ♦ACRES
C

C ACCOnOLATE SOHS FOE VARIANCE AND D.F.
C

IF(NPNTS. LT.2) GOTO 20

ST VHCB (ISPEC) =STVBCB (ISPEC) ♦ (ACRES ^^2) ♦ (SEC OB (I SPEC) ^^2)
STVHBD ( ISPEC) =STVBBD (ISPEC) +(ACRES^^2)^(S EBRD (ISP EC) ^^2)
STVHPL (IS PEC) =STVBPL( ISPEC) ♦(ACRES^^2)^ (SEP OL (I SPEC) **2)
DFCB2(ISPEC) =DFCB2(ISPEC) + ACEES^^4^SEC0B(ISPEC) ♦♦4/(NPNTS-1 )
DFBD2(ISPEC) =DFBD2 (ISPEC) ♦ACRES^^4^SEBRD(ISPEC) ^^4/ (NPNTS- 1)
DFPL2 (ISPEC) =DFPL2(ISPEq ♦ACBES^^4^SEPaL (ISPEC) ♦♦4/(NPNTS-1 )
20 CONTINUE
RETURN
C

C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦•
C CALCULATE STANDAED EEROBS AND EFFECTIVE D.F.
C BASED OH PREVIOUSLY ACCOHOLATED SUflS
C*************** ***********************************
C

ENTRY ST EAT 2

DO 50 ISPEC=1,35

IF(PLTST(ISPEC) .EQ.O) GOTO 30
STSEPL(ISPEC)=SQET(SrVBPL(ISPEC) /TACSES^^2)
IF(DFPL2(ISPEC)-EQ.0) GOTO 30
EFDFPL(ISPEC) =STVBPL (ISPEC) ♦♦2/DFPL2 (ISPEC)
30 IF(CBTST(ISPEC) .EQ.O) GOTO 50

STSECB (ISPEC) =SQET(STVBCB (ISPEC) /TACRES^^2)
STSEBD(ISPEC) =SQBT (STVBBD (ISP EC) /TACRES^^2)
IF(DFCB2(ISPEC) .EQ.O) GOTO 40
EFDFCB(ISPEC) =STVBCB (ISPEC) ♦♦2/DFCB2 (ISPEC)
40 IF (DFBD2 (ISPEC). EQ.O) GOTO 50

EFDFBD(ISPEC) =STVaBD (ISPEC) ♦♦2/DFBD2 (IS PEC)
50 CONTINUE
RETURN
END

STET0010
STRr0020
STRT0030

♦STRT0040
STRT0050
STBT0050
STRT0070

♦STET0080
STRT0090
STRT0130
STRT0110
STRT0120
STRT0130
STRT0140
STRT0150
STRT0160
STBT0170
STBT0180
STRT0190
STRT0230
STRT0210
STRT0220
STRT0230
STRT0240
STET0250
STRT0260
STRT0270
STRT0280
STRT0290
STRT0300
STRT0310
STRT0320
STRT0330
STRT0340
STBT0350
STBr0360
STRT0370
STRT0380
STRT0390
STRT0400
STRT0410
STRT0420
STRT0430
STRT0440
STBT0450
STBT0460
STRT0470
STRT0480
STRT0490
STRT0500
STBr0510
STRT0520
STRT0530
STRT0540
STRT0550
STRT0560
STRT0570
STBT0580
STRT0590
STRT0600
STRT06 10
STRT0620
STRT0630
STRT0640
STBT0650
STET0660
STRT0670
STRT0680
STRT0690
STRT0700
STRT0710
STRT0720
STBT0730

-78-

C OTP10010

SOBEOOTIHE O0TP1 OTP 10020

C OTP10030
C***** ************************************************************** **$i**oip 10040

C OTP10050

C THIS SOBPROGRAH PRODOCES fiLL OF INVENTS TABLES AND SOHHAFIES OTP10060

C 0TP10070

C THE SOBPROGRAH HAS TWO ENTRY POINTS: OOTPI AND 00TP2 OTP10080

C OTP 10090

C OOTPI PRODOCES ALL OF THE INDIVIDOAL STRATA OOTPOT OTP10130

C 00TP2 PRODUCES THE STRATIFIED TOTAL OOTPOT OTP10110

C 0TP10120

INTEGER T0TACB,T0TABD,T0TAPL,IDATE(2) ,TOTNTR(35) .TOTPTS 0TP10U0

INTEGER DBH, TOPANS, HOLT, REGANS, FC(32J OTP10150

INTEGER SPETIT(35,3) ,TITLE(ia),CLASIZ,SPECIE(32) 0TP10160

INTEGER HINDBH(32) ,NTR(35) OTP10170

REAL VARBRD(35) ,VARP0L(35) ,REGC0F(33) ,REGINT(33) ,SEREG(35) OTP10180

REAL RSQR(33) ,HEANY(35) OTP10190

REAL TREPAC(<»0,35),TRETOT(35) ,CBT0T(35) ,BDT0T(35) ,HSD OTP10200

REAL SEC0B(35) ,SEBED(35) ,SEPUL(35) , PLVLSB(<»0 , 35) ,PLTOT(35) OTP10210

REAL CBVLSH(40,35) ,BDVLSH (40 ,35) ,BASPAC (40, 35) , BASTOT (35) OTP 10220

DOOBLE PRECISION IFILE,SAHHTD OTP10230

COMHON/REG/EEGCOF, REGIHT, SPECIE, MI NDBH,SEREG, NT R,RSQR,HEANr OTP 10 240

1,T0TNTR OTP10250

COHHON/FORflCL/FC OTP 10260

COHHON/STO0T/CBTST(35) ,BDTST(35) , CBVLST (40,35) , BD VLST(40,35( OTP10270

1,STSECB(35) ,STSEBD(35) ,EFDFCB(35) ,EFDFBD(35) ,TACRES OTP10280

2,PLTST(35) ,PLVLST(40, 35),STSEPL (35) ,EFDFPL(35) OTP 10290

COBHON/COHUNE/CL AS IZ, TITLE, nOLT,P80BLV,BAF,TOPA US, REGANS, ACRES OTP10300

1,FRAaE,SA[lMTD,L0GANS,IFILE OTP 10310

C0HH0N/0UT1/BASPAC, BASrOT,BDTOT,BDVLSH,CBT0T,CBVLSH,TREPAC,TRETor OTP10 320

1,SECUB,SEBRD,SPNTS,PLVLSn,PLTOT,SEP0L,FINPOP,VARBRD,VARP0L OTP 10330

DATA ((SPETIT (I,J) ,J=1,3) ,1=1,35) /'WHITE PINE RED PINE ',OTP10343

1'PITCH PINE HEMLOCK BALSAM FIR «, OTP10350

2«TAHARACK •, OTP10360

3'BD-BL-HH SPROCEHORHAY SPROCE SCOTCH PINE CEDAR ', OTP10370

4'OTHER SOFTWOODSSOGAR HAPLE RED MAPLE WHITE ASH «, OTP 10380

5'HILLOB ASPEN YELLOW POPLAR •, OTP10390

6«BLACK CHERRY WHITE BIRCH YELLOW BIRCH SWEEP BIRCH ', 0TP10400

7«GREY BIRCH BEECH BASSWOOD ', OTP10410

8«RED OAK WHITE OAK BLACK OAK •, OTP 10420

9«ELH GOH HICKORY HOPfiOENBEAM ', OTP10430

1 'OTHER HARDWOODSALL SOFTWOODS ', OTP 10440

2'ALL HARDWOODS ALL SPECIES '/ OTP10450

C ACCOMULATE NDMBER OF STRATA SAMPLED OTP10460

ISTRA=ISTRA*1 OTP 10470

C ACCOMOLATE TOTAL NUBBER OF POINTS SAHPLED IN ALL STRATA OTP10480

TOTPTS=TOTPTS + BPNTS OTP 10490

C OTPIOSOO

C CALCOLATE T- VALUE FOB THIS STRATA OTP 10510

C OTP10520

TVAL=0 OTP10530

DF=NPNTS-1 OTP 10540

IF{DF. GE. 1) CALL MDSTI (PROBLV, DF.TVAL) OTP10550

IPAGE=IPAGE+1 OTP 10560

C OTP10570

C******************lf ************** OTP 10580

C PRINT HEADING ON OOTPOT FILE OTP 10590

C********************************* OTP 106 00

C OTP10610

CALL DATE(IDATE) OTP10620

REGOSE='ALL' OTP10630

IF(EEGAMS.EQ.'NO') R EGOSE=' SOME' OTP 10640

WEITE(3,100) TITLE, IPAGE OTP10650

OFRAHE=FRAME OTP 10660

IF (FRAME. EQ. 'PLOT') OFRAHE=' PLOT' OTP10670

ALPHA=1.0-PROBLV OTP10680

WRITE (3, 10) IDATE,SAHBTD, BAF,0 FRAME, NPNTS,CLASIZ, rOPANS, LOS ANS, OTP 10690

1REG0SE,ISTRA,IFILE, ACRES, ALPHA OTP10700

10 FORflAT(//' INVENT VEH. 3 11/1/78 I. N. E. R. O.N.H-'/ 0TP10710

l/'ODATE OF RON: ' , 2 A5//' 0' , A9, ' =' , F6. 2 ,28X, ' NOBBER OF ',A5,'S ', OTP10720

2'SAHPLED =• ,I4/«0DIAMETER CLASS SIZE =',I3,21X OTP10730

3, 'HEIGHTS BEASORED BY ',A5,' HEIGHT IN ',A4/ 0TP10740

4«0HEIGHTS OF •,A4,' VOLOBE TREES BEASORED ' ,8X ,' STRATUM #',13/ OTP10750

5'OINPOT FILE NAME: • ,A1 0,18X, 'ACREAGE OF TRACT =',F7. 1/ OTP10760

6'OCONFIDENCE LEVEL = ' , F3. 2) OTP10770

-79-

C«** ********** **«**«****•***•****«•«» ****««*«**«««« OTP10780

C PRINT BEGBESSION ANRLYSIS SOMMAEY OTP 10790

Q:^^^^^**^^^t ************************* ************** OT PI 0800

C 0TP10810

IP(REGANS.EQ. 'lES') GOTO 90 OTP10820

IPiGE=IPAGE*1 OTP10830

HBITE(3,100) TITLE, IPAGE 0TP108U0

BBITE(3,20) OTP10850

20 FORaAT('0«,99(«*')/' *',97X,'*V' ♦' ,28X, • AHALJSIS OF THE HEIGHT' , OTP 10860

1* - DBH EELiriONSHIP',2eX,'*V' *',97X,«*V' *' ,35X, « MODEL: HEX • , OTP10870

2«GHT = BO ♦ B1/DBH',3l»X, '*•/' *',97X,'*V' ', 99 (•*•)/' 0« ,99 ('-« ) OTP108B0

3/« '.ISX,': •,5X,«B0« ,<«X,' :• ,5X,'B1' ,i»X,«:« ,5X,'R' ,5X,': STABDAE', OTP10890

IJ'D : AVEEAGE : VOLOHE : HEIGHT'/' SPECIES : I NTEBC • ,orP1 0930

5'EPT : SLOPE : SQOAEED : ERROR % : HEIGHT : TREES ',OTP10910

6' : TREES'/' ',99('-')/' ') OTP10920

IFLAG=0 OTP10930

DO 50 ISPEC=1,32 OTP1091J0

IF(TOTNTE(ISPEC)- EQ.O) GOTO 50 OTP10950

IF(STB(ISPEC) .LT. 3.AND. NTB(ISPEC)-LT. TOTNTR(ISPEC) ) IFLAG=1 OTP 10960

IF(NrB(ISPEC) .GE.3) HBITE(3,30) (SPETIT (ISPEC, J) , J= 1 , 3) OTP10970

1 ,EEGINT(ISPEC) ,B2GC0F(ISPEC) ,aSttE(ISPEC) ,SEBEG(ISPEC) , OTP10980

2 IJEAHY(ISPEC) ,TOTNTa(ISPEC),NTE(ISPEC) OTP10990
30 FOEMATC ' , 3 (A5) , 2( ' : ' , F9. 3, 2X) , 3 (' : ' ,F8. 1 ,3X) ,2 (• : ' ,18 ,3 X) > OTP11000

IF(NTR(ISPEC) .LT. 3.AND.BTE(ISPEC) . EQ. TOrNrB (ISPEC) ) OTP11010

1 HRITE(3,60) (SPETIT(ISPEC,J) ,J=1,3) ,SEBEG(ISPEC) , 0TP11020

2 flEABY (ISPEC) , TorNTR(ISPEC) , NTE (ISPEC) OTP11030
IF (BTR( ISPEC) .LT-3.ABD. HTE(ISPEC).LT. TorNrR(ISPEC) ) OTP 11 040

1 HRITE(3,«0) (SPETIT(ISPEC,J) ,J=1,3) ,SEREG (ISPEC) , OTP11050

2 MEAHY (ISPEC) ,TOTHTR(ISPEC) ,BTB(ISPEC) OTP11060
»0 FOEMATC ',3A5,':',10X,'*** HAENIBG ***',10X, 0TP11070

1 2{';',F8. 1,3X),2( •: ',I8,3X) ) OTP11030

50 COBTIBOE OTP 11090

WEITE(3,250) OTP11100

DO 70 ISPEC=33,35 OTP11110

WRITE (3, 60) (SPETIT(ISPEC,J),J=1,3) ,SEBEG (ISPEC) , MEAN J (IS PEC( OTP11120

1 ,TOTNTR(ISPEC) ,NTE(ISPEC) OTP11130

60 FOBHATC ' , 3A5 ,' : ' ,35X,2 (' : ' ,F8. 1 , 3X) ,2 (' : ' ,18, 3X) ) 0TP11140

70 IF(ISPEC. EQ. 34.0E.ISPEC. EQ. 35) BBITE(3,250) OTP11150

IF(IFLAG.EQ.I) WEITE(3,80) OTP11160

80 FOEHAT('0'/'0*** WAENING ♦*♦ INSOFFICIEBT HDMBER OF HEIGHT TB«,0TP11170

1'EES HEASOBED FOE THIS SPECIES!') OTP11180

C OTP11190

C ***************************************** OTP11200

C OOTPOT STAND AND STOCK TABLES OTP11210

C ***************************************** OTP 11 220

C OTP11230

90 DO 210 ISPEC=1,35 OTP11240

IF(TEEror (ISPEC). EQ.O. 0) GOTO 210 OTP11250

IPAGE=IPAGE+1 0TP11260

HBITE(3,100) TITLE, IPAGE OTP11270

100 FORMAT('1',14A5,26X,I3/'0') DTP11280

C OTP11290

C THE FOLLOWING SECTION COBTBOLS THE OOTPOT OF THE SPECIES FORM CLASS OTP11300

C PEIBTED IB THE STABD TABLE HEADING. THE FOEM CLASS IS NOT PRINTED OTP 11310

C ONDEE THE FOLLOWING CONDITIONS: OTP11320

C 1. STAND TABLES FOB SOFTWOOD, HARDWOOD, OB ALL SPECIES OTP 11330

C ( ISPEC >= 33 ) OTP113U0

C 2. TREES MEASORED BY TOTAL HEIGHT ( TOPANS = "TOTAL" ) 0TP11350

C 3. STAND TABLES FOE SOFTWOODS WHEN "MIXED" HEIGHTS OSED OTP11360

C ( TOPANS = "MIXED" AND ISPEC <= 11 ) OTP11370

C 0TP11380

AFcriT=' • ; AFCNOn=' • OTP11390

IF(ISPEC.GE.33) GOTO 120 OTP11400

IF(TOPANS.EQ. 'TOTAL') GOTO 120 OTP11410

IF(T0PAHS.EQ. 'MIXED'. AND.ISPECLE. 11) GOTO 120 OrP11420

AFCTIT='FC = ' OTP 11 430

C ENCODE WEITES THE INTEGEB VALOE OF FC INTO AFCBOP AS ASCII CHARACTERS OTP11440

EBCODE(2,110,AFCN01I) FC(ISPEC) 0TP11450

110 FORMAT (12) OTP11460

120 WEITE(3,130) (SPETIT (ISPEC, J) , J= 1 ,3) , AFCTIT, AFCNOM OTP11470

130 FOEHAT (•0',99(' *')/' ','*', 97X, '*' / OTP11480

1,' *',39X,3(A5) ,35X,A5,A2,X,'*'/' ♦',97X,'*'/ 0TP11490

2' ',99 (•*')/'0',99('-')/' ',10X,': BASAL ABBA :', 0TP11500

3' TEEES : SAWLOG : SAWLOG :', OTP11510

4' PULPWOOD'/' DIAMETER : '6X , 'PEE' , 5X , ' : ' ,6X, • PER* , 5X,' : • , 0TP11520

5' CUBIC FOOT VOLUME : BOARD FOOT VOLUME : CUBIC FOOT VOLUME'/ OTP11530

6' ',' CLASS', 3X, ': ',5X,'ACEE',5X, OTP11540

7':',5X,'ACRE' ,5X,3(':',5X,'PEE ACEE',6X)/' ',99('-')/' ') 0TP11550

-80-

DO 150 DBH=CLASIZ,aO,CLRSIZ

IF{TREPAC(DBH,ISPEC).EQ.O.O) GOTO 150

HRITE(3,140) DBH, BASPAC{DBH, ISPEC) ,TEEPAC (DBH.ISPEC)
1 ,CBVLSH(DBH,ISPEC) , BDVLSH( DBH, ISPEC) , PLVLSM (DBH, ISP

1IJ0 FOEMATC •,i»X,I2,aX,' :• ,F10.1,'JX,': ',?10. 1,«Z,' : SFIO

1 2(9X,':', F10. 1))

150 CONTINOE

IF (PLTOT (ISPEC). EQ.O) IEEPaL=0

IF(PLTOT(ISPEC) .NE.0)IEEPOL=(SEPUL (ISPEC) ♦TVAL*100/PLTOT (IS
1.5
IF (CBTOT (ISPEC) .EQ.O) IEECtrB=0

IF {CBTOT( ISPEC) -HE.O) IEECOB= (SECUB (ISPEC) *TVAL* 100/CBTOT ( IS
1.5
IF(BDTOT(ISPEC) -EQ.O) IERBBD=0

IF (BDTOT (ISPEC) -HE. 0) IERBED= (SEBED (ISPEC) ♦TV AL* 100/BDTOT ( IS
1.5

WRITE(3,160) BASTOT (ISPEC) ,TEETOr (ISPEC) .CBFOT (ISPEC) .lEECB
IBDTOT(ISPEC) ,IEEBBD,PLTOT(ISPEC) ,IEEPOL
160 PO RMAT('0', 99 ('-M/'O'rax, 'TOTAL', 2X,':',3X,F7. 1, <»X , • : • , 3X,
1,<»X,':',F10.1,' +',I3,'X',2(3X,':',F10.1,« +',13,'%')/
2' + ',52X,'_',2(19X,'_')/'0',99('-'))
C CALCOLATE AND OOTPOT HSD AND C. V.

BSD=SQET( BASTOT (ISPEC) /TEETOT (ISPEC) /-005U54)
IBDCV=0 : IPLCV=0

IF(BDTOr (ISPEC) -NE.O) IBDC?=SQET (VAE BR D (ISPEC) )*100/BDTOT(
1.5

IF (PLTOT (ISPEC)
1.5
HRITE(3,170) HSD.IBDCV
170 FOEHAT('0'/'0',5X,«nEAN STAND DIAHETEE =■.
1,9X,'B0AED FOOT COEFFICIENT OF VAEIATIOH
C CALCULATE AND OUTPOT MERCHANTABLE USD
RBSTOT=BASTOT (ISPEC)
ETRTOT=TRETOT(ISPEC)
DO 180 DBH=CLASIZ,7,CLASIZ

RBSTOT=RBSTOT-BASPAC (DBH, ISPEC)
180 ETflTOT=RrETOT-TREP AC (DBH, ISPEC)
IF(RTRT0T.LT„0-01) RHSD=0

IF(BTETOT.GE.0.01) EMSD=SQET (R BSTOT/BTRTOT/. 0051»5«)
WEITE(3,190) RHSD,IPLCV
190 FOEHAT('0'/'0' ,5X,' HERCHANTABLE H.S.D. = '.FH.I

1,9X,' POLP-WOOD COEFFICIENT OF VARIATION =',Ii»,'%'l
C
C CALCOLATE AND OOTPOT (ESTIHATED) % CROISE

IF(FRRHE.EQ.'PLOT') PCTCHZ=100 ♦( 1 .0-FINPOP)

IF( FRAME. EQ. 'POINT') PCTCRZ=- 5«5<»*NPNTS*MSD*HSD/ACRES/BAF
WRITE (3,200) PCTCRZ
200 FORMAT('0'/'0' ,5X, 'PERCENT CROISE =',F5-1,'X')
210 CONTINOE
C
C «*«**4i«i*«««**4i**4[** «******• **4i****4i ************ **4i**

C CALCOLATE AND OOTPOT SPECIES COMPOSITION PERCENTAGES

C

IPAGE=IPAGE+1
HRITE(3,100) riTLE,IPAGE
WRITE(3,220)
220 FORMAT('0',99(' *')/• *',97X,'*'/

1,' *',30X, 'SPECIES COMPOSITION B? PERCENT' , 37X, •*• /

2' *',97X,'*'/' ',99('*')/'0',99(«-')/

3* ',3(15X,':'), 2(5X,'SAWLOG',6X,':'),l»X,'P0LP»OOD'/

.NE.O) IPLCV=SQRT (VARPUL(ISPEC) ) *100/PLrOT(

' , F4 . 1
=',ia,'X')

1' •,3(15X,':') ,3X,«C0BIC FOOT', 4X , ' : ' , '
5,<»X,': ',3X,'CaBIC FOOT'/
6' ',I»X, 'SPECIES', «X,':',' BASAL AREA
73(': ',5X,'V0L0ME«,6X)/' ',99('-')/' ')
IF(CBTOr (35). EQ. 0) CBTOT (35)=-999

BOARD FOOT'

:',5X, 'TREES

IF(BDrOT (35)-EQ.O)
IF(PLTOT(35)-EQ.O)
DO 21*0 ISPEC=1,3a
IF (TEETOT (ISPEC)

BDTOT (35) =-999
PLTOT (35) =-999

EQ.0.0) GOTO 2U0
BASCOH=BASTOT (ISPEC) ♦lOO/BASTOT (35)
TRECOM=TRET0r (ISPEC) • 100/TRETOT (35)
C0BCOM=CBTOT(ISPEC) •100/CBTOT (35)
BDCOM=BDror (ISPEC) ♦100/BDTOT (35)
P0LCOM=PLTOT (ISPEC) ♦100/PLTOT (35)

OTP11560

0TP11570

OTP11580

EC) OTP11590

.1, OTP11600

OTP11610

OT?11620

OTP11530

PEC) )+ OTPIieUO

OTP 11650

OTP11660

PEQ ) + 0TP11670

OTP11630

OTP11690

PEC)) ♦ OTP 11700

OTP11710

B, DTP 11720

OTP11730

F7.1 0TP117U0

OTP11750

OTP11760

0TP11770

OTP11780

0TP11790

ISPEC) +OTP11800

OTP11810

ISPEC) +OTP 11820

OTP11830

0TP118«0

OTP11850

OTP11860

0TP11870

OTP11880

0TP11890

OTP11930

0TP11910

OTP11920

DTP11930

OTP119U0

OTP11950

0TP11960

OTP11970

0TP11980

OTP11990

OTP120D0

OTP12010

OTP12020

OTP 12030

OTP120IJ0

OTP12050

OTP 12060

OTP12070

OTP 12080

OTP12090

0TP12100

OTP12110

OTP12120

0TP12130

OTP121«0

OTP12150

0TP12160

OTP12170

0TP12180

OTP12190

OTP12200

0TP12210

OTP12220

0TP12230

OTP122U0

OTP12250

OTP12260

OTP12270

0TP12280

OTP12290

OTP12300

-81-

1

230
240
250

HBITE(3,230) (SPETIT (ISPEC,J) , J= 1 ,3) , BASCOM .TRECOH,

CO BCOM, BDCOH, POLCOM
FORHATC ',3 (AS) , ': ' ,2 (3X,F7. 2, 5X, • ; • ) . 2 (UX, F7. 2, 6X, • : ' )
,UX,F7. 2)

IF{ISPEC.EQ.32) HRITE(3,250)
WRITE(3,250)

FOBOATCO', 99 (•-•)/' ')
IF(CBTOT(35) .EQ.-999) CBTOT(35) =0
IF(BDTOT(35). EQ.-999) BDTOT(35)=0
IF(PLTOT{35).EQ--999) PI.TOT(35) =0

C

C ««***««***«'*• *****4'<)>«** ************************* **********

C OUTPUT VOLUHE TOTALS FOR ALL SPECIES

C*** ***************************************** **************

C

IPAGE=IPAGE+1
WRITE (3,100) TITLE.IPAGE
«RITE(3,260J
260 FORHAT(«0', 99 {•*•)/• *',97X,'*V

1' *«,30X, 'VOLUHE TOTALS FOR ALL S PECIES • , 38X, • *•/
2' *',97X,'*V' ' ,99('*«)/'0' ,99 ('-•)/
3« • ,15X,2(' :«,11X,'SAHL0G',11X) ,•: ', 10X, • PULP WOOD'/
Ui SPECIES : sex, 'CUBIC FOOT VOLUHE • ,5X, •:', 6X, • BOARD' ,

5' FOOT V0L[JBE',5X,':' ,6X, 'CUBIC FOOT VOLUHE'/
6' ',15X,2(': •,10X,'PER ACRE', 10X) ,'C',10X, 'PER ACRE'/
7' ',99('-')/' ')
DO 290 ISPEC=1,35

EQ. 0. AMD. ISPEC. LT. 33)

IF (CBTOT (ISPEC) .EQ.O.AND. PLTOT (ISPEC)
1 GOTO 280

HRITE(3,270) (SPETIT (ISPEC, J( , J= 1 , J) ,CBTOT{ISPEC) ,BDTOT (ISPEC) ,
1 PLTOT(ISPEC)

270 FORHATC ' ,3 (A5) , 3( ': ', F17. 1, 1 1X) )

280 IF(ISPEC.EQ.32.0R.ISPEC.EQ-31».OR.ISPEC. EQ.35) HBITE(3,250)
290 CONTINUE
C
C *****************************************************

C OUTPUT VOLUME TOTALS EXPANDED BY ACREAGE

C

IPAGE=IPAGE+1
WBITE(3,100) TITLE.IPAGE
HRITE(3,300)
300 FOBHAT('0',99('*')/' *',97X,'*'/

1' *' ,30X,'VOLOHE TOTALS EXPANDED BY ACREAGE' ,3«X, • *•/
2' *',97X,'*'/« ',99('*«)/'0', 99 ('-')/
3* ',15X,2(«:',11X,'SAHLOG',11X) ,': ' , 10X , ' PULPWOOD ' /
«• SPECIES :',6X, 'CUBIC FOOT VOLUME ', 5X ,':' ,6X, • BOARD' ,

5' FOOT VOLUME', 5X, •:' ,6X, 'CUBIC FOOT VOLUHE'/
6' ',99('-')/' ')
DO 330 ISPEC=1,35

IF (CBTOT (ISPEC) . EQ.O.AND. PLTOT (ISPEC) , EQ. 0. AND. ISPEC. LT. 33)
1 GOTO 320

TOTACB=CBTOT(ISPEC) *ACRES+.5
TOTABD=BDTOT (ISPEQ *ACRES*.5
TOTAPL=PLT0T (ISPEC) *ACRES+.5

WRITE(3,310) (SPETIT (ISPEC, J) ,J=1,3) , TOTACB ,T0TABD,T0TAPL
310 FORHATC ' ,3 ( A5) ,3 ( «: ' ,1 17, 1 1 X) )

320 IF(ISPEC. EQ.32.OR.ISPEC.EQ.3U.0R.ISPEC.EQ.35) BBITE(3,250|
330 CONTINUE
RETURN

OTP12310
OTP12320
OTP 12330
OTP123U0
OTP12350
OTP12360
OTP12370
OTP123B0
0TP12390
OTP12'»00
OTP12mO
OTP12II20
OTP12430
0TP12HH0
OTP12l*50
OTP 121160
OTP12U70
0TP12tt80
OTP12l»90
OTP12500
0TP12510
OTP12520
OTP 12530
OTP125a0
OTP12550
OTP 12560
OTP12570
OTP 12580
OTP12590
OTP12600
0TP12610
OTP12620
OTP 12630
OTP12640
OTP12650
OTP 12660
OTP12670
OTP 12680
OTP12690
0TP12700
OTP 127 10
OTP12720
OTP 12730
OTP127aO
3TP12750
OTP 12760
Orpi2770
OTP 12780
OTP12790
OTP12800
0TP12810
OTP12820
OTP12830
OTP12840
OTP 12850
OTP12860
OTP12870
OTP12880
OTP12890
OTP 12900

-82-

Q**^t* ************************ *******************************************

c

C OOTPOT TABLES FOE STRATIFIED TOTAL

C

C************ ******************* ****************************************

C

ENTRY O0rP2
WRITE(5,3«0)
3U0 FOEHAT(«0HHAT DO 100 WANT TO TITLE THE COHBINED STRATA OOTPOT?')

BEAD(5,350) TITLE
350 F0EHAT(1<JA5)
IPAGE=IPAGE+1
WRITE(3,100) TITLE, IPAGE
C PRINT HEADING ON OOTPOT FILE
OFRAME=FRAaE

IF(FRAHE. EQ. 'PLOT') OFBAHE=» PLOT*

WRITE (3,360) IDATE.CLASIZ, ISTEA,TACRES,OFRAHE ,TOTPTS, ALPS A
360 FORaAT(//' INVENT VER. 3 11/1/78 I. N. E. R. O.N.H. •//

VODATE OF RON: '.aAS//

2«0DIAMETER CLASS SIZE = • ,13 ,21X , • NOHBEB OF STRATA SAMPLED =•
3,I3/«0TOTAL ACEEAGE SAHPLED =• ,F9. 1 , 13X,« TOTAL •,A5
l»,'S SAMPLED =«,I5/«0CONFIDENCE LEVEL = ',F3.2)

C

C**** **********************************

C OOTPOT DIAMETER CLASS VOLOHE TABLES

c

DO U50 ISPEC=1,35

IF(CBTST{ISPEC) .EQ.O. AHD.PLTST(ISPEC) .EQ.O) GOTO 1*50
IPAGE=IPAGE+1
WRITE(3,100) TITLE, IPAGE
C

C THE FOLLOWING SECTION CONTROLS FOEB CLASS PRINTING AS PREVIOOSLY
C DESCRIBED. THE INPOTS OF THE LAST STRATOH ANALYSED CONTROLS WHAT
C IS PRINTED IN THE STRATIFIED TOTAL OOTPOT.
C

AFCTIT=' • : AFCNOa=« •
IF(ISPEC.GE. 33) GOTO 370
IF(TOPANS.EQ. "TOTAL') GOTO 370

IF(TOPANS.EQ."HIXED'. AND.ISPEC.LE. 11) GOTO 370
AFCTIT=«FC = '
C ENCODE WRITES THE INTEGER VALOE OF FC INTO AFCNOM AS ASCII CHARACTERS
ENCODE(2,110, AFCNOH) FC(ISPEC)
370 WRITE(3,380) (SPETIT(ISPEC, J) ,J=1 ,3) , AFCTIT, AFCNOH
380 FORMAT('0', 99 {•*•)/' *',97X,'*'/
*',39X,3(A5) ,35X,A5,A2,X,'*'/
• •,97X,'»'/' ', 99 C*')/'©', 99 {•-')/

DI AHETEH' , 1»X,2 ( • : • , 1 1X, 'SAWLDG • , 1 1X) , • : • , 1 OX , ' POLPWOOD' /
CLASS i',6X,'C0BIC FOOT VOLOME • ,5X, ':' ,6X, ' BOARD' ,

FOOT VOLOME', 5X, ':• ,6X,'C0BIC FOOT VOLOHE'/
',99('-')/' •)
DO UOO DBH=CLaSIZ,UO,CLASIZ

IF(CBVLST (DBH, ISPEC) . EQ- 0. AND- PLVLST(DBH, ISPEC) .EQ.O)

GOTO 400
ICBVST=CBVLST (DBH, ISPEC) +. 5
IBDVST=BDVLST(DBH,ISPEC)+.5
IPLVST=PLVLST (DBH, ISPEC) +.5
WRITE (3,390) DBH,ICBVST,IBDVST,IPLVST
390 FORHATC ',I8,7X,3(':',I17, 11X))

HOO CONTINOE

ICBrST=CBTST (ISPEC) ♦. 5

IPLTST=PLTST (ISPEC) +. 5

C CALCOLATE ERROR AS A PERCENT

IERCOB=0;IERBBD=0:IEBPOL=0;TVAL=0

IF(CBTST(ISPEC)-EQ,0) GOTO 410

IF(EFDFCB (ISPEC) .GE-1)CALL MDSTI (PEOBLV,EFDFCB(ISPEC) ,TVAL)

IERCOB= (STSECB(ISPEC) *TVAL*100*TACEES/CBTST (ISPEC) ) *.S

IF(BDTST(ISPEC) .EQ.O) GOTO 420

TVAL=0

IF(EFDFBD (ISPEC) . GE- 1 ) CALL MDSTI (PEOBLV, EFDFBD(ISPEC) , TVAL)

IERBRD= (STSEBD(ISPEC) *TVAL* 100*TACEES/BDTST (ISPEC) ) +.5

IF (PLTST (ISPEC). EQ.O) GOTO 430

TVAL=0

IF(EFDFPL(ISPEC).6E. 1) CALL MDSTI (PEOBLV, EFDFPL( ISPEC) ,TVAL)

IERPOL= (STSEPL (ISPEC) *TVAL*100*TACRES/PLrST (ISPEC) ) +. 5

WRITE (3,440) ICBTST, lERCOB ,IBDTSr,IEBBRD,IPLTSr,IERPOL

FORM AT ('0', 99 ('-' ) /'O ', 4X, 'TOTAL ',6X

1 ,3(':',I17,' +',I3,'X',5X)/'+*,15X,3(19X,'_',9X)/

2 '0',99('-'))

1*
2'
3'
U'
5'
6'

1

IBDTST=BDTST (ISPEC) +. 5

410

420

430

440

450 CONTINOE

OTP12910
DTP 12920
OTP12930
OTP12940
OTP 12950
OTP12960
OTP 12970
OTP12930
OTP 12990
OTP13000
OTP13010
OTP13020
OTP13030
OTP 13040
OTP13050
OTP13060
OTP 13070
OTP13080
OTP13090
OTP13100
OTP13110
0TP13120
OTP13130
0TP13140
0TP13 150
OTP13160
0TP13170
OTP13iaO
0TP13190
OTP13200
0TP13210
OTP13220
OTP 13230
OTP13240
OTP13250
OTP 13260
OTP13270
OTP 13280
OTP13290
OTP13300
OTP13310
OTP13320
0TP13330
OTP13340
OTP13350
OTP13360
OTP13370
OTP 13380
OTP13390
OTP13400
OTP 13410
OTP13420
OTP 13430
OTP13440
OTP 13450
OTP13460
OTP 13470
OTP13480
OTP 13490
OTP13500
OTP13510
OTP 13520
OTP13530
OTP 13540
OTP13550
OTP13560
0TP13570
OTP13580
OTP 13590
OTP13600
OTP13610
OTP13620
OTP13630
OTP 13640
OTP13650
OTP13660
OTP 13670
OTP13680
OTP 13690
OTP13700
OTP13710

-83-

C OTP 13720

C***!******* <■************* •****«*«** OTP13730

C OOTPOT &LL SPECIES VOLOME TABLE OTP 13740

C***«*«*4'***4<**** ********** ******** OT PI 3750

C OTP 13760

IPAGE=IPAGE + 1 OTP13770

WRITE (3, 100) TITLE, IPAGE OTP13780

(IBITE(3,<»60) OTP13790

460 FORHATCO', 99 (•*')/• •',97X,'*V OTP13830

1' *',30X,'VOL0HE TOTALS FOE ALL S PECIES • ,38X, • * V OTP13810

2« *',97X,«*V' •,99('*«)/'0',99('-')/ OTP13820

3« •,15X,2 Cx'.IIX.ISAHLOG'.IIX) ,•: •,10X,«POLPHOGDV OTP13830

«• SPECIES :«,6X,'C0BIC FOOT VOLO HE •, 5X ,«:•, 6X, • BOARD ', OTP138U0

5'FOOT VOLOHE'.SX, 'J'.ex, 'COBIC FOOT VOLOBE'/ OTP13850

6' ',99{'-')/' •) OTP13860

DO 480 ISPEC=1,35 OTP13870

IF(CBTST(ISPEC) . EQ. 0. AND. PLTST (ISPEC) . EQ- 0. AND. ISPEC. Lr.33) OTP 13880

1 GOTO 470 OTP13890

ICBTST=CBTST(ISPEC) +.5 ; IBDTST=BDTST (ISPEC) +.5 OTP139aO

IPLTST=PLTST(ISPEC) +. 5 OTP 13910

HRITE(3,310) (SPETIT(ISPEC,J) ,J=1,3) , ICBTST.IBDTST.IPLTST OTP13920

470 IF(ISPEC. EQ. 32. OR. ISPEC. EQ.34. OR. ISPEC. EQ.35) HRirE(3,250) OTP13930

480 CONTINUE OTP 13940

RETURN OTP13950

END OTP 13960

-8A-

C HDST0010

SOBPOOTINE HDSTI (Q,F,X) HDST0020

C f1DST0030

C »•*»•*»*••*****•*»•»*♦*»********»***************»♦********♦**»*»***♦♦♦♦ a DS TO OHO

C nDSTOOSO

C THIS SOBPROGBRM IS COMPOSED OF TWO SOBROOTINES ADAPTED FROH THE (1DST0060

C INTERNATIONAL HATHEWATICAL AND STATISTICAL LIBRARIES , INC. (IBSL) nDST0070

C HDSTI - INVERSE STODENT'S T DISTRIBOTION HDST0080

C BDNRIS - INVERSE NORHAL PROBABILITY DISTRIBOTION FUNCTION HDSr0090

C BOTH SUBROUTINES HAVE BEEN EDITED REHOVING CAPABILITIES NOT NEEDED nDSTOlOO

C IN THIS APPLICATION. THE PEHAINING CODE IS REPRODUCED HERE IN nDSTOIIO

C ACCORDANCE WITH THE IMSL POLICY ON RESEARCH WORK AS STATED 3N nDST0120

C PAGE INTRO-22 OF THE JULY, 1977 LIBRARY HANOAL. HDST0130

C HDSTOmO

Qm*m*********************************** **************************** *****fiQST:0 150

C SDSTOieO

C EXACT INTEGRAL FOP 2 D.F. HDST0170

IF (ABS(F-2.0) .GT. .000001) GO TO 10 HDST0180

X = SQET(2.0/(Q*(2.0-Q) )-2.0) (1DST0190

GO TO 50 nDST0200

10 HPI = 1.570796 MDST0210

C EXACT INTEGRAL FOE 1 D.F. nDST0220

IF (ABS(F-I.O) .GT. .000001) GO TO 20 HDST0230

A = Q*HPI HDST02U0

X = COS(A)/SIN(A) nDST0250

GO TO 50 HDST0260

C EXPANSION FOR N GREATER THAN 2 BDST0270

20 A = 1.0/(F-0.5) nDST0280

B = «b.O/(A*A) MDST0290

C = ( (20700. ♦A/B-98.) 'A-ie.) *A + 96. 36 BDST0300

D = ((9«.5/(B+C)-3.0) /B+1.0) •SQRT(A*HPI) ♦F MDST0310

XX = D*Q nDST0320

Y = XX«*(2.0/F) (1DST0330
IF (Y .GT. A*. 05) GO TO «0 BDST03H0

Y = { (1.0/( ((F+6.0)/(r*Y)-0.0«9*D-0.822)* (F+2.0)»3.0) ♦ BDST0350
1 0.5/(F+i».0) ) »Y-1.0) »(F+1.0)/(F + 2.0) -H.O/Y MDST0360

30 X = SQRT(F*Y) HDST0 370

GO TO 50 BDST0380

C ASYHPTOTIC INVERSE EXPANSION ABOUT BDSr0390

C HOBHAL BDSTOaOO

<»0 X = .5 * Q BDSTOHIO

CALL BDNRIS(X,XX) HDST0U2O

Y = XX * XX BDST0U30
IF (F .LT. 5.) C = C+0.3*{F-«.5) ♦(XX+0.6) HDSTOUHO
C= ({{.05*D*XX-5.0)»XX-7.0) *XX-2.0J*XX+B+C BDST0U50
f = (((((0.«»Y + 6. 3) *Y*36.) *Y + 9i«.5) /C-Y-3.0) /B + 1.0) •XX HDST0«60

Y = A»Y*Y HDST0l»70
D = Y BDST0U80

Y = .05*Y*Y+Y BDST0H90
IF (Y .GT. .002) Y = EXP{D) - 1.0 BDSTOSOO
GO TO 30 BDST0510

50 RETURN BDST0520

END MDST0530

-85-

c
c

SOBBOOTINE flDNR

DATA

DATA

DATA
1

DATA
1

DATA

DATA
1

DATA

DATA
1

DATA

GO TO 40
10 IF (X.LT. . 15)

IP {X.LT. 1.85)

SIGHA = -1.

A = 2.-X

B = X-1.

GO TO 50
20 Z = U-X

SIGHA = SIGN(1.

Z = ABS{Z)

GO TO 80
30 SIGHA = 1.

A = X

B = 1<.-X

GO TO 50

no INT = 3
X = 2- *P
GO TO 10

50 W = SQET(-ALOG (A+A*B) )
IF (W. LT. 2.5) GO TO 70
IF (W.LT- n. ) GO TO 60

IS(P,Y)

SQBT2/1. I»1«2136/

A 1,A2,A3/-. 575 17029, -1.89651 33, -.0 5496 2605/

BO, B1,B2,B3/-. 11 377303, -3. 2934740, -2. 3749959,

-1.1875145/
CO, C1,C2,C3/-. 11 466659, -.13147744, -.236820 10,

-050739749/
DO, D1,D2/-44. 279769, 21.985462,-7-5861027/
E0,E1,E2,E3/-. 0566842208,. 39370209, -.3 16650 10,

.062089629/
F0,F1,F2/-6. 2667859,4. 6662627,-2-96 28832/
GO, G1,G2,G3/. 000 1851 1591, -.0020281 520,

-. 14983844,. 010786386/
HO, HI, H2/. 099529751, .521 17329,-. 068883 009/

GO TO 30
GO TO 20

-Z)

WI = 1./H
SS = ((G3*iri+G2
SD = ((HI-t-H2)*H
F = » ■<■ H«(GO+S
GO TO 90

)*WI+G1)«BI
I+H1) ♦HI+HO
N/SD)

60 SB = ((E3«W+E2) *H*E1) *H
SD = ((H + F2)*H + F1)*ll*F0
F = H + H»(E0+SN/SD)
GO TO 90

70 SH = {(C3»8*C2) *B*C1) ♦«
SD = ((B+D2)*H+D1) *W+DO
F = H ♦ B»(CO+SH/SD)
GO TO 90

c

c

80

Z2 = Z*Z

F = Z*Z*{B0*A1«

c

c

90

Y = SIGHA*F
IF (INT.NE. 3)

Y = -Y*SQBT2

100

BETOEH
ESO

INVEBSE GAOSSIAN ENTBY

REDOCED ABGOHENT IS IN (.85,1.),
OBTAIN THE TBANSFOEHED VABIABLE

W GBEATER THAN 4., APPEOX. F BY A
BATIONAL FUNCTION IN 1./B

B BETBEEN 2.5 AND 4., APPEOX. F
BY A RATIONAL FONCTION IN B

B BETBEEN 1-13222 AND 2.5, APPEOX.
F BY A BATIONAL FaNCTION IN W

Z BETBEEN 0. AND -85, APPBOX. F
BY A EATION&L FONCTION IN Z

Z2/(B1*Z2+A2/(B2+Z2-fA3/(B3 + Z2) ) ) )
FOEH THE SOLUTION BY
THE PROPEE SIGN

GO TO 100

MOLT. F BY

HDNE0010
HDNB0020
HDNB0030
HDHB0040
HDHR0050
HDNE0060
HDNE0070
HDNEOOBO
HDNE0090
HDNR0100
HDNE0110
nDNE0120
HDHBOnO
HDNE0140
HDHB0150
HDNB0160
SDNR0170
HDNB01B0
HDNE0190
nDNB0200
HDNB0210
HDNB0220
HDNE0230
tlDNR0240
HDNE0250
HDNB0260
HDNE0270
HDNE02B0
HDNE0290
HDNB0300
aDNB0310
HDNB0320
MDNB0330
HDNB0340
HDNR0350
aDNB0360
HDNS0370
HDNB0330
HDNE0390
HDNB0400
HDNB0410
nDNB0420
HDNE0430
{1DNR0440
flDNR0450
HDNR0460
BDNE0470
aDNB0480
HDNE0490
HDNR0530
{1DNB0510
ODNE0520
HDNE0530
HDKB0543
HDNB0550
HDNB0560
HDNE0570
HDMB0580
HDNR0590
aDNB0600
nDNE0610
HDNE0620
flDNB0630
HDNE0640
HDNS0650
HDNE0660
HDNE0670
{IDNE0680

-86-

LITERATURE CITED

Abramowitz, Milton, and Irene Stegun. 1964. Handbook of Mathematical
Functions. U.S. Dept. Int. Nat. Bureau Stnds, 1044 pp.

Avery, T. Eugene. 1975. Natural Resource Measurements. New York: McGraw-
Hill, 339 pp.

Barnard, Joesph E. , et al. 1973. "Forest survey cubic-foot volume equations."
U.S. Dept. Agr, Forest Service, Research Note NE-66, 2 pp.

Barrett, James P., and Mary Nutt. 1975. Survey Sampling in the Environmental
Sciences: A Computer Approach. Hanover, N.H.: Project COMPUTe, 319 pp.

Bruce, Donald. 1961. "Prism cruising in the western United States and volume
tables for use therewith." Mason, Bruce and Girard. Portland, Oregon.

Cochran, William G. 1977. Sampling Techniques. New York: John Wiley &
Sons, 428 pp.

Cunia, Tiberius. 1959. "Notes on cruising intensity by the Bitterlich method."
Journal of Forestry 57:849-851.

Davis, Kenneth P. 1966. Forest Management: Regulation and Valuation. New
York: McGraw-Hill Book Co., 519 pp.

Dillworth, J.R., and J.F. Bell. 1976. Variable Probability Sampling.
Corvallis, Oregon: O.S.U. Bookstores, Inc., 130 pp.

Freese, Frank. 1962. Elementary Forest Sampling. U.S. Dept. Ag. Forest
Service, Ag. Hndbk, No. 232, 91 pp.

Freese, Frank. 1967. Elementary Statistical Methods for Foresters. U.S.
Dept. Agr. Forest Service, Ag. Hndbk. No. 317, 87 pp.

Honer, T.G. 1967. "Standard volume tables and merchantable conversion

factors for the commercial tree species of central and eastern Canada."
For. Mgt. Research and Services Inst., Ottawa, Ont . FMR-X-5, 154 pp.

Husch, Bertram, Charles I. Miller, and Thomas W. Beers. 1972. Forest Mensur-
ation. New York: John Wiley & Sons, 410 pp.

IMSL. 1977. Library 2 Reference Manual. Houston: International Mathematical
and Statistical Libraries, Inc.

Johnson, Floyd A. 1958. "Standard error of estimated average timber volume per
acre under point sampling when trees are measured for volume on a sub-
sample of all points." U.S. Dept. Ag. Forest Service, Research Note No. 201,
Portland, Oregon.

Mesavage, Clement, and James W. Girard. 1946. "Tables for estimating board-
foot volume of timber." U.S. Dept. Agr. Forest Service, 94 pp.

Neter, John, and William Wasserman. 1974. Applied Linear Statistical Models.
Homewood, 111.: Richard D. Irwin, Inc., 842 pp.

Satterthwaite, F.E. 1946. "An approximate distribution of estimates of variance
components." Biometrics, 2; 110-114.

Wiant, Harry V., and Froylan Castaneda. 1977. "Mesavage and Girard's volume
tables formulated." U.S. Dept. Int. B.L.M. Res. Inv. Nts. 4; 4 pp.

mG 1 3 2004

BioSci

"TJ JU . ) L.

■no. 501 516-

L.