United States

Department of

Agriculture

Forest Service

Southern Forest

Experiment Station

New Orleans,

Louisiana

General Technical Report

so-52

December 1984

MERCHANDIZER:

Adapting USLYCOWG to

the Marketplace

Richard W. Guldin
SUMMARY

USLYCOWG is a computer model that estimates yields in cubic feet

per acre from unthinned slash and loblolly pine plantations established on cutover sites in the West Gulf states. However, most stump-

‘age in the West Gulf region is sold in other volumetric units, such as

cords, tons, and board feet, or by weight. A series of computer routines has been developed that converts USLYCOWG yields in cubic

feet to tons of fuel chips, cords of pulpwood and chip’n’saw logs, and

board feet of sawtimber, or weights of all products in tons.
MERCHANDIZER: Adapting

USLYCOWG to the Marketplace

Richard W. Guldin

INTRODUCTION

Deciding to invest in regeneration or other

forest management activities often depends

heavily on projected growth of a stand and on

yield of stumpage at harvest time. If stumpage

volumes are denominated in common marketing

units, it is easier to assign market values to the

volumes and calculate the economic consequences of management alternatives.

Dell et al. (1979) and Feduccia et al. (1979)

developed a computer model (USLYCOWG) that

projects the growth and yield of unthinned slash

and loblolly pine stands growing on cutover

sites in the West Gulf region. USLYCOWG requires the user to enter site index, stand age,

and seedling density. The specified density may

be either the number of trees surviving at the

given stand age or the number of trees planted.

If the latter, the user must also specify whether

Coile’s and Schumacher’s survival model or an

ideal survival model developed by USLYCOWG’s

authors is to be used to calculate the number of

stems surviving at the stand age desired. From

these parameters, the model can project stand

diameter and height distributions over time.

Diameters, heights, and ages are used, in turn,

by tree taper equations to calculate bole volumes in cubic feet for one-inch diameter classes.

This report assumes that the reader is familiar

with USLYCOWG. Copies of the research papers

describing USLYCOWG in detail are available

from the Southern Station.

Although the cubic foot is a commonly understood measure of volume, it is not a volume unit

of commerce in the marketplace for southern

timber. Wood is valued and sold in other volumetric units. For example, pulpwood and chip’n

‘saw logs are commonly sold by the cord. Sawlogs are most often sold by the board foot according to the Doyle, International Quarter Inch or

Scribner Decimal C Log Rules. Weight scaling

and whole tree logging have contributed to the

popularity of tons as a unit of commerce, for not

only fuelwood chips, but also for pulpwood and

logs. Market prices and stumpage bids are denominated in these merchantable units.

Simple multiplicative conversion of USLYCOWG cubic foot volumes to merchantable units

is not possible. USLYCOWG estimates cubic foot

volumes-say for the la-inch diameter classfrom stump height (0.5 feet) to a 2-, 3-, or 4-inch

top (both inside and outside bark). But USLYCOWG is unable to determine what portion of

the cubic foot volume of the trees in the 12-inch

class should be valued as sawtimber versus the

portion to be valued as pulpwood.

A set of computer routines, collectively called

MERCHANDIZER, has been developed to convert USLYCOWG yields into merchantable volume measures. These cover the various products

that can be obtained from trees in each diameter class. The routines are designed to be

embedded in USLYCOWG because they call up

existing USLYCOWG subroutines. The purpose

of this report is to describe how the MERCHANDIZER routines function and what additional input data are required by MERCHANDIZER. Appendix A is a listing of the USLYCOWG program written in PDP-11 FORTRAN

showing the additional program code (in italics)

needed to incorporate the MERCHANDIZER

r0utines.l Appendix B contains sample USLYCOWG output produced by the program after

MERCHANDIZER is incorporated.

’

‘The latest version of USLYCOWG incorporating MERCHANDIZER can be obtained in machine-readable form if

you send the author a blank g-track computer tape. The

program will be written on your tape in ASCII (1600 bpi, 80

characters per record, 1 record per block, unlabelled) and returned to you. If you require different tape characteristics,

please write to the author before sending your tape. The

mailing address is USDA Forest Service, T-10210 Postal Service Building, 701 Loyola Avenue, New Orleans, Louisiana

70113.

Richard W. Guldin is Research Forester, Southern Forest Experiment Station, Forest Service-USDA, New Orleans, La.
APPROACH

Data Required.-MERCHANDIZER requests a

second line of input data following the first line

called by USLYCOWG. The input specifications,

variable names, and default values are listed in

table 1.

Conversion Philosophy.-MERCHANDIZER

begins by calculating the diameter inside bark

(dib) at a height of 8,.75 feet from the ground,

given the diameter at breast height, outside bark

(dbh), total tree height, age, and associated taper

equation parameters. A height of 8.75 feet represents a stump height of 0.5 feet, a minimum

log length of 8.0 feet and a 0.25 foot trim allowance for bucking. The dib at 8.75 feet is first

compared to the minimum sawlog scaling diameter, SAWTBR (see table 1). If the dib at 8.7

feet exceeds SAWTBR, the program branches to

calculate sawlog volume. If dib does not exceed

SAWTBR, but is larger than the maximum pulpwood scaling diameter (PULPMX), the program

branches to calculate chip’n’saw log volume. If

dib is less than PULPMX, the program branches

to calculate pulpwood volume for the entire bole

up to the minimum pulpwood scaling diameter

(PULPMN).

When the program branches to the sawlog

portion of the routine, the height to SAWTBR is

calculated. If this height is less than 16.75 feet,

the height is rounded to the next lower evennumbered foot (BOLT length). The dib at (BOLT

+ 0.75 feet)2 becomes the scaling diameter for the

alternative sawlog volume equations. Volumes

are calculated for three log rules from:

Board Feet Doyle = (DIB-4.)2(BOLT) (1)

16

Board Feet International =

(0.796DIB2 + 1.375DIB - 1.23BOLT (2)

16

Board Feet Scribner = O.O494(BOLT)(DIB)2

- O.lB$(BOLT)(DIB)

- 0.269BOLT (3)

The volume from (3) is rounded to the nearest

10 feet to obtain board feet according to the

Scribner Decimal C Log Rule. Sawtimber volume

in cunits is calculated from the USLYCOWG

taper equations.

ZThe 0.75 feet includes a stump height of 0.5 feet and a 0.25

foot trim allowance.

If the height to SAWTBR exceeds 16.75 feet,

then progressive 16.25 foot bolts are cut and

scaled as 16 foot logs. When the height is less

than 16.25 feet, but more than 8.25 feet remains

until SAWTBR is reached, the length is reduced

to the next lowest even-numbered foot, the scaling dib calculated, and volume calculated. The

height to the top of this final bolt becomes the

base height for the next product-either chip’n

‘saw logs or pulpwood. When less than 8.25 feet

remains, the total height to the top of the last

18foot log becomes the base height for tGe next

product.

When PULPMX equals or is less than one

inch smaller than SAWTBR, no chip’n’saw volumes are calculated and the program branches

from the sawtimber volume section directly to

the pulpwood volume section. If SAWTBR and

PULPMX differ by more than one inch, the

height to PULPMX is calculated and the difference between that height and the height to

the top of the last sawtimber log becomes t‘he

length of the chip’n’saw log. The total chip’n

‘saw length is rounded to the next lower evennumbered foot and a cubic foot volume is calculated from the USLYCOWG taper equations.

This result is divided by CONV(3)-one of the

several conversion factors listed in table l-to

obtain chip’n’saw volume in cords. If the bole

length for a chip’n’saw log is less than 8.25 feet,

no chip’n’saw volumes are calculated. Instead,

the program branches to calculate pulpwood

volumes using the height to the top of the last

sawtimber log as the base of the pulpwood section of, the bole.

Once sawtimber and/or chip’n’saw volumes

are calculated and the height to the base of the

pulpwood section of the bole is determined,

height to PULPMN is calculated. The difference

between height to PULPMN and height to the

base of the pulpwood section of the bole is the

pulpwood bole length. The cubic foot volume of

the bole section is calculated from the USLYCOWG taper equations and divided by CONV(2)

to obtain the pulpwood volume in cords.

Fuel chip volume is calculated as the difference between the total cubic foot volume, ground

to tip outside bark (calculated by USLYCOWG)

and the sum of the sawtimber, chip’n’saw, and

pulpwood cubic foot volumes calculated by

MERCHANDIZER. Because the MERCHANDIZER volumes are all calculated based on dibs,

fuel chips include not only that portion of the

total stem above PULPMN, but also the volume

of the bark sheathing the stem below PULPMN.

If a clean chip volume (no bark) is desired, a

comment card in the program listing describes

2
Table L-MERCHANDIZER inputparameters

Column Variable Variable

number format name

Description and default values

l-4 F4.1

5-8 F4.1

9-12 F4.1

13-15 1 3

16-19 F4.0

20-23 F4.0

24-27 F4.0

28-31 F4.0

32-38 F7.2

39-45 F7.2

46-52 F7.2

53-59 F7.2

60-62 1 3

63-65 1 3

66-68 1 3

PULPMN

PULPMX

SAWTBR

IRULE

Minimum pulpwood diameter, inside bark, small end.

The default is 3.0 inches.

Maximum pulpwood diameter, inside bark, small end.

The default is 9.5 inches.

Minimum sawlog diameter, inside bark, small end.

The default is 9.5 inches.

Sawtimber log scaling rule. Enter -2 for weight

scaling; -1 for volume in cunits; 0 for board feet

Doyle, the default; +l for board feet International

Quarter Inch; or +2 for board feet Scribner Decimal

c;.

CONV(1)

CONV(2)

CONV(3)

CONV(4)

VALUE(l)

VALUE(S)

VALUE(B)

VALUE(4)

ICULLP

ICULLD

ICULLV

Fuelwood (dirty chips) conversion factor. Cubic feet

per ton. The default is 90.

Pulpwood conversion factor. If IRULE is -2, cubic feet

per ton with a default of 26. Otherwise, cubic feet

per cord with a default of 80.

Chip’n’saw conversion factor. If IRULE is -2, cubic

feet per ton with a default of 30. Otherwise, cubic

feet per cord with a default of 90.

Sawtimber conversion factor. If IRULE is -2, cubic

feet per ton with a default of 30. Not used otherwise, with a default of 0.

Fuelwood (dirty chip) value. Dollars per green ton. The

default is $0.00 for all VALUE(I).

Pulpwood stumpage value. Dollars per green ton if

IRULE = -2; dollars per cord otherwise.

Chip’n’saw stumpage value. Dollars per green ton if

IRULE = -2; dollars per cord otherwise.

Sawtimber stumpage value. If IRULE = -2, dollars per

green ton; if IRULE = -1, dollars per cunit; if

IRULE is 0, +l, or +2, dollars per thousand board

feet, applicable log rule.

Percentage of sawtimber volume in a diameter class

less than or equal to ICULLD that is unsuitable for

sawing lumber. Expressed as an integer between 5

and 100. The default is 0.

The maximum sawtimber diameter class where

ICULLP percent of the volume is unsuitable for

manufacturing sawn lumber. Default is SAWTBR

rounded down to the next whole inch.

The product class where the sawtimber cull volume is

added. Enter -1 for fuelwood chips; 0 for pulpwood,

the default; or +I for chip’n’saw logs.

how to save the USLYCOWG total volume

ground to tip inside bark instead of total volume

ground to tip outside bark for chip calculations.

The total chip volume in cubic feet is divided

by CONV(l) to obtain chip weight in tons.

If IRULE is -2, the program reports weights

instead of volumes. The total cubic foot volume

of each product in each diameter class is calculated, as just discussed. However, CONV(4) is

used to convert the volume into a green wood

tonnage. The result is an “inside bark” weight,

because inside bark taper equation coefficients

are used to calculate the volumes. The total

weight for each product is then increased by 10

percent to account for bark. Sample output is

included in Appendix B.

In some stands, a certain percentage of the

sawtimber-sized trees are unsuitable for manufacturing sawn lumber. For example, a loblolly

pine stand infested with fusiform rust at an

early age often has a certain percentage of butt

logs deformed by rust cankers. A cull deduction

option has been included in the routines to allow

the user to specify the percentage of the total

sawtimber volume in a particular diameter class

that is cull and the maximum sawtimber diameter class where this cull percentage applies.

When the user enters an integer between 5 and

100 and the ICULLP value, the routine subtracts

ICULLP percent of the sawtimber volume, whatever the scaling rule selected, from the diameter

class total. Then, ICULLP percent of the saw3
timber cubic foot volume is added to chip’n’saw

volume if ICULLV is set at 1, to pulpwood volume if ICULLV is set at 0, or to fuelwood chips

if ICULLV is set at -1.

DISCUSSION

MERCHANDIZER adds no precision to USLYCOWG yield estimates. It only converts the

estimated yields to merchantable units. The

flexibility inherent in being able to specify

diameter cutpoints and cubic foot conversion

factors for the various stumpage categories can

be detrimental if unrealistic diameter limits or

conversion factors are entered. The routine looks

for wildly inaccurate limits or factors (less than

50 percent of the default value), and, if found,

replaces them with the default values. If the

user has better estimates of the cord or tonnage

conversion factors than the defaults, they should

be used. The default values were obtained from

Koch’s (1972) review of the literature.

MERCHANDIZER assumes the midpoint of

the dbh class is the diameter of the average tree

in the class, calculates volume for that average

tree, and multiplies those volumes by the number of trees in the class to get diameter class

volume totals. If USLYCOWG were altered to

generate diameter distributions having smaller

dbh class widths, some marginal improvement

in the precision of volume estimates would

occur. Wider dbh class intervals would result in

poorer volume estimates, especially for sawtimber. However, one-inch classes should give

estimates that are sufficiently precise for most

users.

If one or more of the volume categories are

not of interest, the diameter limits for that class

should be set accordingly. For example, if a

market for chip’n’saw logs does not exist in the

user’s area, SAWTBR and PULPMX should be

equal. Merely setting chip’n’saw value equal to

$0 will not result in chip’n’saw volumes being

incorporated into either pulpwood or sawtimber

volumes. Fuel chip volumes are always calculated.

MERCHANDIZER provides no information on

log quality. Users desiring to know what proportion of chip’n’saw or sawtimber logs might

be suitable for plywood peelers must look elsewhere for that information. However, if one

knows, for example, that 60 percent ,of the logs

exceeding 7.5 inches dib small end are suitable

for veneer bolts and veneer bolts are bought on

an MBF Doyle basis, then one could set

SAWTBR at 7.5 inches and multiply the volume

result by 60 percent to estimate veneer bolt

volume. More precise estimates cannot be obtained from MERCHANDIZER.

The VALUES entered need not be stumpage

values. Prices for products delivered to the mill

(FOB the mill) could also be used. Landowners

who are vertically integrated wood growers,

harvesters, and transporters could use MERCHANDIZER to estimate wood value per acre

in terms of prices paid for wood delivered to the

mill. Total value per acre of products priced

F.O.B. the mill would allow cost analyses to

include harvesting and transportation cost components in addition to stand establishment and

management costs.

CONCLUSION

MERCHANDIZER breaks down USLYCOWG

cubic foot volumes by diameter class into four

merchantable product categories. When dollar

values are assigned to these volumes, land

managers have a better basis for evaluating

forest regeneration and stand management investments.

LITERATURE CITED

Dell, T.R.; Feduccia, D.P.; Campbell, T.E.; Mann,

W.F., Jr.; Polmer, B.H. Yields of unthinned

slash pine plantations on cutover sites in the

West Gulf region. Res. Pap. SO-147. New Orleans, LA: U.S. Department of Agriculture,

Forest Service, Southern Forest Experiment

Station: 1979.84 p.

Feduccia, D.P.; Dell, T.R.; Mann W.F., Jr.; Campbell, T.E.; Polmer, B.H. Yields of unthinned

loblolly pine plantations on cutover sites in

the West Gulf region. Res. Pap. SO-148. New

Orleans, LA: U.S. Department of .Agriculture,

Forest Service, Southern Forest Experiment

Station; 1979.88 p.

Koch, P. Utilization of the southern pines. Vol.

2 - Processing. Agric. Handb. 420. Washington, DC: U.S. Department of Agriculture; 1972.

2 vol.
Appendix A

USLYCOWG Program Listing Incorporating MERCHANDIZER

(MERCHANDIZER PDP-11FORTRAN statements added to USLYCOWG are in italics)
-~PROGRAM LISTING

C ---e--w

-------

; -- -- - --

C -e-e..--

C --se---

C -------

-------

c" -------

C -------

C -------

--em---

; - ----- -

C --..----

C -w----e

-------

; -------

-------

; ---- - - -

C -------

--se---

; -------

C -----mm

C -------

c---_---

c-- - - - - -

c---e--m

c-------

c -------

c--em---

c-------

c-----em

-------

c" -------

-------

c" ---- * --

----e--

c" ---- -- -

----e-e

; --- - - --

C -------

-------

c"-----

----e--

; ---- . - -

---wa--

c" - -- ----

C ----a--

YIELD TABLES AND STAND STRUCTURE FOR UNTHINNED SLASH

AND LOBLOLLY PINE PLANTATIONS ON CUT-OVER (UNPREPARED BUT

PROBLEM FREE) SITES IN THE WEST GULF COASTAL PLAIN.

REFERENCES:

SLASH--

DELL,T.R.,D.P.FEDUCCIA,T.E.CAMPBELL,W.F.MANN,JR.,

AND B.H.POLMER.

1979. YIELDS OF UNTHINNED SLASH PINE PLANTATIONS

ON CUTOVER SITES IN THE WEST GULF REGION,

U.S. DEP. AGRIC. FOR. SERV. RES. PAP. SO-147,

84 P. SOUTH. FOR. EXP. STN., NEW ORLEANS, LA.

LOBLOLLY--

FEDUCCIA,D.P.,T.R.DELL,W.F.MANN,JR.,T.E.CAMPBELL,

AND B.H.POLMER.

1979. YIELDS OF UNTHINNED LOBLOLLY PINE PLANTATIONS

ON CUTOVER SITES IN THE WEST GULF REGION,

U.S. DEP.,AGRIC. FOR. SERV. RES. PAP. SO-148,

88 P. SOUTH. FOR. EXP. STN., NEW ORLEANS, LA.

YIELD CONK!?RSIONS--

GULDIN, R.W.

1984. MERCHANDIZER: USLYCOWG-COMPATIBLE ROUTINES

FOR CONVERTING CUBIC FOOT YIELDS TO ST&'PAGE

VOLChlES, WEIGHTS, AND VALVES IN OTHER UNITS.

U.S. DEP. AGRIC. FOR. SERV. GEN. TECH. REPT.

SO-XX. XX P. SOUTH. FOR. EXP. STN. NEW ORLEANS, LA.

ACKNOWLEDGMENTS:

THIS COLLECTION OF ROUTINES IS ESSENTIALLY AN ADAPTATION

OF ROUTINES PROVIDED BY SMALLEY AND BAILEY IN CONJUNCTION

WITH USDA FS RESEARCH PAPERS SO-96 AND SO-97. WE DID

DIFFER FROM THEIR PROCEDURES IN THE APPLICATION OF TAPER

EQUATIONS TO DEFINE TREE VOLUME AND WERE ASSISTED IN THIS

ENDEAVOR BY F. THOMAS LLOYD AND EDWARD W. WHITEHORNE OF

THE SOUTHEASTERN FOREST EXPERIMENT STATION. THE YIELD

CONVERSION ROUTINES WERE PREPARED ESPECIALLY FOR VSLYCOWG.

THIS PROGRAM READS zwo CONTROL CARDSAND PRINTS TABLES. INPUT

UNIT IS INPT, OUTPUT UNIT IS IOUT. CONTROL CARD FORMAT :

200 FORMAT (914,2F4.0,213)

202 FORMAT (3F4.1,13,4F4.0,4F7.2,3I3)

C -w--L-e

c------- FIRST CARD

-------

~-~~~--c!!?!Y~~

VARIABLE MEANING

____-___ ____---_____-----------------------

C ------- ISB SITE INDEX: BEGINNING

-------

Em- - - - -

;: ; ISE ENDING

9-12 ISS STEP SIZE

C --e-s--

C ----..-- 13-16 IAB AGE: BEGINNING

C ----e-e 17-20 IAE (YEARS SINCE ENDING

-------

; -------

21-24 IAS PLANTING) STEP SIZE

(DO NOT USE PLANTATION

C e------ AGES LESS THAN 6 YEARS)

-------

; ------- 25-28 INB NUMBER STEMS: BEGINNING N

C ------- 29-32 INE (N) ENDING N

-------

E - -- ----

33-36 INS PER ACRE STEP SIZE FOR N

C ------- 37-40 BAGE BASE AGE FOR SITE INDEX(USUALLY 25)

6
-------

; -- -----

41-44

-------

; --- --- - 45-47

C -------

-------

;- -- -- --

C -------

-------

; --- ----

-------

; ------- 48-50

C -------

C -------

c------- SECOND CARD

c--------------

c-----em 1- 4

C-------

c-------

c-----me 5- 8

C-------

c------a

C------- 9-12

c-a-----

c-------

c -------

c -------

c -------

c -------

C -----mm

c -------

c -----me

c -------

c ------- 13-15

c ---- ---

c----e-e

c -------

c -------

c-- -- --

c ------- 16-19

c -------

C -------

c ------- 20-23

cm-e----

c---m-ec -------

c----e-e

c ------- 24-27

c-- -- -- -

c -------

c -------

c -------

c ------- 28-31

c -------

c -------

c-----me

c -------

ABAGE ALTERNATE BASE AGE FOR CORRESPONDING

SITE INDEX (OPTIONAL - USUALLY 50)

KEY =l USE N AS NUMBER PLANTED AND

CALCULATE SURVIVAL FROM IDEAL

SURVIVAL MODEL

=2 USE N AS NUMBER PLANTED AND

CALCULATE SURVIVAL FROM COILE +

SCHUMACHER 1964 MODEL

=3 USE N AS NUMBER SURVIVING

KSPCE =l SLASH

=2 LOBLOLLY

IRULE = -2 WEIGHT SCALE (ALL PRODUCTS IN TONS)

-1 SAWTIMBER VOLCME IN CVNITS

0 IN BOARD FEET, DOYLE RULE (THE DEFAULT')

+1 IN BOARD FEET, INTERNATIONAL l/4 INCH

+2 IN BOARD FEET, SCRIBNER DECIMAL C

CONVfl I = FUELWOOD (DIRTY CHIPS) CONVERSION FACTOR

CUBIC FEET PER TON. DEFAULT IS 90.

CONVf2) = PULPWOOD CONVERSION FACTOR fI.B.)

IF IRULE IS -2, CUBIC FEET PER TON WITH A

DEFAULT OF 26. OTHERWISE, CUBIC FEET PER

CORD WITH A DEFAULT OF 80.

CONVf3) = CHIP'N'SAW CONVERSION FACTOR (I.B.)

IF IRULE IS -2, CUBIC FEET PER TON WITH A

DEFAULT OF 30. OTHERWISE, CUBIC FEET PER

CORD WITH A DEFAULT OF 90.

CONV(4) = SAWTIMBER CONVERSION FACTOR

IF IRULE IS -2, CUBIC FEET PER TON WITH A

DEFAULT OF 30. NOT USED OTHERWISE WITH A

DEFAULT OF 0.

PULPMN =Xx.X MINIM&' PULPWOOD DIAMETER, I.B.,

,%ALL END. DEFAULT IS 3.0 INCHES

PULRllX =Xx.X MAXIMl&l PULWOOD DIAMETER, I.B.,

LARGE END. DEFAULT IS 9.5 INCHES

SAWTBR =Xx.X MINIM&l SAWLOG DIAMETER, I.B.,

SMALL END. DEFAULT IS 9.5 INCHES

IF PULRMX.NE.SAWLOG, AND THESE TWO

DIAMETERS ARE SEPARATED BY MORE THAN

8.25 FEET OF BOLE, CHIP'N'SAW LOG VOLCJ4ES

ARE ~A~~~LATED USING PUL~X AS MINIMUM

CHIP'N'SAW LOG DIAMETER, INSIDE BARK,

SMALL END. OTHERWISE, THE BOLE BELOW

SAWTBR DIAMETER IS CONVERTED TO PVLEWOOD.

c------- THE FOLLOWING DATA FIELDS FOR STMPAGE VALUES ALL HAVE 2 DECIMAL

c------- PLACES EMBEDDED TO READ COLLARS AND CENTS. WHEN THE WEIGHT SCALING

c------- OPTION IS SELECTED, ALL VALUES ARE IN DOLLARS PER GREEN TON.

cw-e----

c------- 32-38 VALUE(l) = FUEL (DIRTY) CHIP VALUE ($/GREEN TON)

C-------

c------- 39-45 VALUE(Z) = PULPWOOD STl&!PAGE VALUE ($/CORD)

C-------

c------- 46-52 VALUE(B) = CHIP'N'SAW STUMPAGE VALUE ($/CORD)

c-------

c------- 53-59 VALUE(4) = SAWTIMBER STUMPAGE VALUE

c------- IF IRULE = -I, $/CUNIT

C---a--- IF IRULE = 0, t2 OR t2, $/MBF
c---m---

C------- 60-62 ICULLP = PERCENTAGE OF SAWTlMBER-SIZED LOGS IN A

C---S- -- SPECIFIED DIAMETER CLASS LESS THAN OR

C-----em EQUAL TO DIAMETER CLASS 'ICULLD' THAT ARE

C-- -- -- UNSUITABLE FOR MANUFACTURING SAWN Ll&'BER.

C ---w--e EXPRESSED AS AN INTEGER BETWEEN 5 AND 100

C --m--s- WITH A DEFAULT OF 0.

C --w---m

C--m---m 63-65 ICULLD = THE MAXIM&f SAWTIMBER DIAMETER CLASS WHERE

C-----em 'ICULLP' PERCENT OF THE LOGS ARE UNSUITABLE

C --w---m FOR MANUFACTURING SAWN LUMBER. DEFAULT IS

C---e--s 'SAWTBR' ROUNDED DOWN TO THE NEXT WHOLE INCH.

C ---e--w

C--m---w 66-68 ICULLV = PRODUCT CLASS WHERE THE SAWTIMBER CULL

C------- VOLUME IS ADDED: -1 FOR FUELWOOD CHIPS;

C------- 0 FOR PULPWOOD, THE DEFAULT; +1 FOR

C------- CHIP'N'SAW LOGS.

C-mm----

INTEGER YIELD(2,9),CRATIO

DIMENSION BASAL(2),ITREES(41),VOLT(2),VOLM(6),COEFLM(2,4),CONVf4I.

lCOEFCR(2,4),BETA(4),STllMPVf2,91,VALUEf4~,ST~PG~40,6~,FYIELDf40~

-----em

; --- - - - - INITIALIZE LAMBDA COEFFICIENTS --

C --em--- COLUMN 1 OF COEFLM IS SLASH

C ------- COLUMN 2 OF COEFLM IS LOBLOLLY

C -----es

DATA COEFLMI

1 1.99830 , -5.6919 ,

: -;.g , 2.2974 ,

4 -0:24879 :

0.16988 ,

0.77763 /

E

---e--e

------- INITIALIZE MEAN CROWN RATIO COEFFICIENTS --

-------

; -- - - -- -

COLUMN 1 OF COEFCR IS SLASH

COLUMN 2 OF COEFCR IS LOBLOLLY

C -------

DATA COEFCR/

1 214.600 , 248.63 ,

2 -44.814 , -35.872 ,

3 -27.049 , -43.565

4 -23.453 , -33.424 i

DATA INPT,IOUT/5,1/

OPEN (uNIT=~,NAME=~USLY.~UT~,TYPE=~NE~I~,DISP=~SAYE~,

lCARRIAGECONTROL='FORTRAN')

READ (INPT,200) ISB,ISE,ISS,IAB,IAE,IAS,INB,INE,INS,BAGE,ABAGE,

1 KEY,KSPCE

READ (INPT,2021 PULI;MN,PULA41X,SAWTBR,IRULE,(CONVfI),I=lI,4),

1 (VALUE(I),I=l,41,ICULLP,ICULLD,ICULLV

C -------

C ------- CHECK FOR LIMITS ON AGE -- NOT LESS THAN 1

C --e-w--

IF (IAB.LT.~.~R.IAE.LT.~I Go To 220

GO TO (701,702),KSPCE

701 WRITE (IOUT,lOl)

GO TO 703

702 WRITE (IOUT,EOl)

703 WRITE (IOUT,301)

101 FORMAT (lH1///33X,66HYIELDS BY DIAMETER CLASSES FOR UNTHINNED SLAS

1H PINE IN PLANTATIONS)

201 FORMAT (lH1///32X,69HYIELDS BY DIAMETER CLASSES FOR UNTHINNED LOBL

1OLLY PINE IN PLANTATIONS)

301 FORMAT (1~ ,43x,47Ho~ CUT-OVER (UNPREPARED BUT PROBLEM FREE) SITES

1/63X,6HIN THE/55X,23HWEST GULF COASTAL PLAIN)

C----e-m

C------- CHECK INPUT PARAMETERS AND BUILD MERCHANDIZER HEADER

C-------

IFIPULRk'X.LT.4.95) PUL@lX=9.5

IFfPULMdN.LT.1.95) PULI141N=3.0

IFfSAWTBR.LE.5.95) SAWTBR=9.5
IF(SAWTBR-PULRb'X).LT.l.) PULEMX=SAWTBR

IFfIRULE.EQ.-2) GO TO 31

IF(CONV(l).LE.49.) CONV(l)=90.

IF(CONV(2J.LE.49.) CONVf2)=80.

IF(CONV(3J.LE.49.) CONV(3)=90.

IF(CONV(4j.NE.O.) CONV(I)=O.

GO TO 32

31 IF(CONV(l).LE.49.) CONV(l)=90.

IF(CONV(2).LE.l4.) CONV(Z)=ZS.

IF(CONV(3).LE.14.) CONV(3)=30.

IF(CONV(4j.LE.14.) CONV(4)=30.

32 IF(IABS(IRULEJ.GE.3) IRULE=O

LW 33 I=1,4

IF(VALUE(I).LT.O.) VALUE(I)=O.

33 CONTINUE

CULLP=FLOAT(ICULLP)/lOO.

IF(ICULLP.LT.5) ICULLD=INT(SAWTBR-1.)

IF(IABS(ICULLV).GT.l) ICULLV=O

WRITE(IOUT,119)

119 FORMAT(/////54X,26(lHm)/)

WRITE(IOUT,l20) SAWTBR

120 FORMAT(lH ,49X,34H*** MERCHANTABILITY PARAMETERS ***,/54X,26(lH-1,

l///lSX,67HSAWTIMBER VOLUMES ARE CALCULATED FROM STLhfP HEIGHT (0.5

BFEET) TO A ,F4.1,33H INCH TOP DIAMETER (INSIDE BARK)./lOX,lllHLOGS

3 ARE BUCKED INTO 16 FOOT LENGTHS, PLUS A 3 INCH ALLOWANCE. LOG S

4CALING DIAMETER IS THE DIAMETER INSIDE /lOX,ll2HBARX AT THE SMAL

5L END OF EACH LOG. IF THE LAST LOG BEFORE REACHING THE MINIMUM SC

6ALING DIAMETER EXCEEDS 8 FEET,/lOX,llOHBVT IS LESS THAN 16 FEET, I

7T IS BUCKED TO THE NEXT SMALLER EVEN-NUMBERED FOOT OF LENGTH, PLUS

8 TRIM ALLOWANCE. /lOX,lOSHTHE REMAINING LENGTH IS DROPPED INTO THE

9 NEXT LOWER STlbfPAGE VALUE CATEGORY. IF THE LAST LOG TO THE MINIM

Alk' /lOX,lOGHSCALING DIAMETER IS LESS THAN 8 FEET LONG, THE ENTIRE

BLENGTH IS DROPPED INTO THE NEXT LOWER STUMPAGE VALVE/lOX,9HCATEGOR

CY. /)

WRITE(IOVT,121) PVLA\IN,CONV(l~,VALUE(l~,PVLEMN,PVLA\IX,CONV(2~,

1 VALUE(2),PULI%lX,SAWTBR,CONV(3),VALVE(3I

121 FOEZMAT(25X,8HSTOMPAGE,lSX,8HDIAMETER,20X,6HVOL~E,20X,8HST~PACE/

127X,5HVALUE,l7X,6HLIMITS,22X,4HVNIT,22X,5HVALVE/27X,SHCLASS,l8X,4H

Z(IB),22X,7HFACTORS,19X,BHPER VNIT/2SX,88ilH )//ZSX,l4HFUELWOOD CHI

3PS,3X,lOHLESS THAN ,F4.1,7H INCHES,8X,F4.O,ilH CUBIC FEET,9X,2H $,

4F7.2,8H PER TON/46Y,l3H NO BRANCHES ,17X,7HPER TON//25X,lSHPULPIJOO

SD BOLTS ,6X,F4.1,2H -,F4.1,15X,F4.O,llH CUBIC FEET,7X,4H $,F7.2,

69H PER CORD/76X,BHPER CORD//ZSX,lSHCHIP'N'SAW LOGS,6X,F4.1,2H -,F4

7.1,15X,F4.O,llH CUBIC FEET, 7X,4H $.F7.2,9H PER CORD/76X,BHPER CO

8RD//I

IF(IRVLE) 342,343,344

342 WRITE(IOVT,125) VALUE(B), SAWTBR

125 FORMAT(25X,lSHSAWTlMBER LOGS ,6X,l2HGREATER THAN,13X,l6H CUNIT

1s SX,3H $,F7.2,1OH PER CVNIT/46X,F4.1, 7H INCHES/25X,88(lHe))

GO TO i49

343 WRITE(IOVT,l26/ VALVE(II,SAWTBR

126 FORMAT(25X,lSHSAWTIMBER LOGS ,6X,lZHGREATER THAN,13X,

116HBOARD FEET DOYLE, 7X,3H $,F7.2,8H PER MBF/46X,F4.1, 7H INCHES/25

2X,88(lH )I

GO TO 379

344 IF(IRVLE.EQ.2) GO TO 345

WRITE(IOUT,127) VALUE(4),SAWTBR

127 FORMAT(25X,lSHSAWTlMBER LOGS ,6X,lBHGREATER THAN,l3X,

116HBOARD FEET INT'L,7X,3H $,F7.2.8H PER MBF/46X,F4.1,7H INCHES/25

2X,88(lH ))

GO TO 379

345 WRITE(IOVT,128) VALVE(4),SAWlBR

128 FORMAT(25X,lSHSAWTlMBER LOGS ,6X,l2HGREATER THAN,6X,

129HBOARD FEET SCRIBNER DECIMAL C,X,3H $,F7.2,8H PER MBF/46X,

2F4.1,7H INCHES/25X,88(lH_I)

GO TO 349

372 WRITE(IOUT,373) PVLEMN,CONV(l),VALVE(lI,PUL~N,PVLpMX,CONV(2) 1 VALUE(2),PVLI;MX,SAWTBR,CONV(3).VALUE(31,CONV(4) 2 VALVE(IJ,SAWlBR

9
373 FORMAT(25X,8HST~PAGE,l5X,8HDIAMETER,20X,6HWEIGHT,20X,8HST~PACE/

127X,5HVALUE,l7X,6HL~ITS,22X,4HUNIT,22X,SHVALUE/27X,5HCLASS,l8X,4H

2tIB),22X,7HFACTORS,19X,BHPER UNIT/25X,88(lH )//25X,lCHFUELWOOD CHI

BPS,3X,lOHLESS THAN ,F4.1, 7H INCHES,8X,F4.0,ilH CUBIC FEET,9X,2H $,

4F7.2,8H PER TON/46X,l3H NO BRANCHES ,17X, 7HPER TON//25X,l5HPVLFWOO

SD BOLTS ,6X,F4.1,2H -,F4.1,15X,F4.O,llH CUBIC FEET,9X,2H $,F7.2,9H

6 PER TON /76X, 7HPER TON//25X,lSHCHIP'N'SAW LOGS,6X,F4.1,2H -,F4.1,

715X,F4.O,llH CUBIC FEET,9X,2H $,F7.2,8H PER TON/76X,7HPER TON//ZSX

8,lSHSAWTIMBER LOGS ,6X,lBHGREATER THAN,l3X,F4.O,llH CUBIC FEET,9X,

92H $,F7.2,8H PER TON/46X,F4.1,7H INCHES,lgX, 7HPER TON/25X,88flHe)/

A/25X,SOHWEIGHT UNIT FACTORS FOR PRODUCTS OTHER THAN FUELWOOD CHIPS

B ARE IN CUBIC FEET OF GREEN WOOD/25X,85HPER TON. TOTAL WEIGHT OF

CTHE GREEN WOOD OF THESE PRODUCTS IS INCREASED BY 10 PERCENT/BSX,27

DHTO ACCOUNT FOR BARK WEIGHT.)

349 IF(ICULLP.NE.0) WRITEfIOUT374) IClJLLP,ICULLD

374 FORMAT~//lOX,81HSPECSAL NOTE: THIS RUN HAS ACTIVATED THE SAWTIMBER

1 CULL FEATURE OF MERCHANDIZER. ,I3,25H PERCENT OF THE SAWTIMBER/

2X,53HVOLlhME IN EACH DIAMETER CLASS BELOW AND INCLUDING THE,I3, 36H

3 INCH CLASS HAS BEEN ASSIGNED TO THEI

IF(ICULLV) 380,381,382

380 WRITE(IOUT,383)

383 FORMAT(24X,79HFUELWOOD CHIP VOLUME CATEGORY BECAUSE, AS CULL VOLLfM

lE, IT WAS DEEMED lJNSlJITABLE/24X,23HTO PRODUCE SAWN LUMBER.)

GO TO 386

381 WRITE(IOUT,384)

384 FORMAT(24X,74HPULPWOOD VOLUME CATEGORY BECAUSE, AS CULL VOLUME, IT

1 WAS DEEMED lJNSUITABLE/24X,23HTO PRODUCE SAWN LVMBER.)

GO TO 386

382 WRITE(IOVT,385)

385 FORMAT(24X,76HCHIP1N'SAW VOLUME CATEGORY BECAUSE, AS CULL VOLUME,

1IT WAS DEEMED lJNSVITABLE/24X,23HTO PRODUCE SAWN LUMBER.)

C----e-m

C------- END OF MERCHANDIZER HEADER

C---we--

386 IBAGE=BAGE

E

---m--e

------- ANY STEP SIZES LESS THAN 1 ARE SET EQUAL TO 1 --

C ------- ALLOW BLANK STEP SIZES AS INPUT WHEN OUTPUT IS

C ------- ONE SITE INDEX, AGE, OR NUMBER VALUE

C ---_---

IF (.ISB.EQ.ISE) ISS=l

IF (INB.EQ.INE) INS=1

IF (IAB.EQ.IAE) IAS=~

i

---_--_

------- SET UP DO LOOPS TO

C -----em AGE, AND NUMBER

C ---m--m

DO 100 IS=ISB,ISE,ISS

SI=IS

DO 100 IN=INB,INE,INS

ANO=IN

ANNOW=IN

DO 100 IA=IAB,IAE,IAS

AGE=IA

AS=AGE

C ---_---

GENERATE TABLES FOR EACH COMBINATION OF SITE,

THIS IS DONE TO

GENERATED FOR ONLY

C - - - - - - - CALL ROUTINE TO COMPUTE MEAN HEIGHT OF DOMS AND CODOMS PLUS THE

-------

c" ---- ---

SITE INDEX FOR THE ALTERNATE BASE AGE (ABAGE)

CALL HEIGHT (H,AS,BAGE,SI,KSPCE,ABAGE,IASI)

C ----e-m

C ------- CALL ROUTINE TO CALCULATE SURVIVAL IF INPUT IS NUMBER PLANTED

C -------

IF (KEY.EQ.l.OR.KEY.EQ.2) ANNOW=SURVIV(H,AGE,ANO,KSPCE,KEY)

IF (KSPCE.EQ.l) WRITE (IOUT,114)

IF (KSPCE.EQ.2) WRITE (IOUT,llS)

WRITE (IOUT,lOZ) IS,IBAGE

102 FORMAT (lHO,SPX,lOHSITE INDEX,I4,11H (BASE AGE,I3,1H))

114 FORMAT (lHl//GlX,lOHSLASH PINE)

115 FORMAT (lH1//6OX,13HLOBLOLLY PINE)

10
IF (KEY.EQ.31 GO TO 10

5 WRITE (IOUT,103) IN

103 FORMAT (lH0,52X,I4,23H STEMS PER ACRE PLANTED)

GO TO 11

10 WRITE (IOUT,112) IN

112 FORMAT (lH0,5X,I4,25H STEMS PER ACRE SURVIVING)

11 IF (KEY.EQ.~) WRITE (1ouT,116)

IF (KEY.EQ.2) WRITE (IOUT,117)

116 FORMAT (lH0,56X,20HIDEAL SURVIVAL MODEL)

117 FORMAT (lH0,57X,18HC-S SURVIVAL MODEL)

IH=Ht.5

WRITE (IOUT,104)

104 FORMAT (lHO/58X,24(1H-),17HCUBIC FOOT VOLUME,24(1H-)/6X,7HGROWING,

145X,9HALL TREES,9X,47H5-INCH CLASS AND GREATER, STUMP HEIGHT .5 FE

2ET,/6X,12HSEASONS AV.,8X,5HSTEMS,28X,6HGROUND,llX,19HFOR O.B. TOP

3s OF --/7X,23HSINCE DtC DBH PER,5X,17HBASAL CR AV.,7X,6

4HT0 TIP,llX,8H2 INCHES,llX,8H3 INCHES,llX,8H4 INCHES/8X,lOHEST.

5HT.,9X,4HACRE,5X,4HAREA,9X,3HHT.,6X,llHO.B. I.B.,6X,llHO.B. I.

6B.,8X,llHO.B. I.B. ,8X,ll:iO.B. I.B.)

WRITE (IOUT,105)

105 FORMAT (lX,3H***,122X,3H***)

DO 12 J1=1,2

DO 12 J2=1,9

12 YIELD(Jl,J2)=0

C -----we

C ------- CALL ROUTINE TO CALCULATE DIAMETER DISTRIBUTION USING WEIBULL

C ------- PARAMETER ESTIMATES

C -------

CALL DIADST (AGE,H,ITREES,DMAX,ANNOW,A,B,C,KSPCE)

nL - - - - - - -

C -m-e--.- EXTRACT COEFFICIENTS FROM COEFLM (LAMBDA) AND COEFCR (MEAN CROWN

C --m---e RATIO) ACCORDING TO SPECIES BEING SUMMARIZED AND CALCULATE

-------

; -- - - - - -

LAMBDA AND MEAN CROWN RATIO

DO 18 K1=1,2

00 15 K2=1,4

GO TO (13,14),Kl

13 BETA(K2)=COEFLM(KSPCE,K2)

GO TO 15

14 BETA(K2)=COEFCR(KSPCE,K2)

15 CONTINUE

DUM=BETA(1)tBETA(2)*ALOG1O(H)tBETA(3)*ALOGlO(AGE)

1 tBETA(4)*ALOGlO(ANNOW)

IF (DuM.LT.~.~) DuM=o.~

GO TO (16,17),~1

16 ALAMDA=DUM

GO TO 18

17 ACROWN=DUM

18 CONTINUE

C ----e-e

C ------- IF CROWN RATIO IS GREATER THAN 100, THEN SET IT EQUAL TO 100

C --e---w

IF (ACROWN-100.) 20,20,19

19 ACROWN=99.9995

C -------

------- CALCULATE ARITHMETIC MEAN DBH (DA) AND QUADRATIC MEAN DBH (DQ)

;- - -- - - -

20 DA=AtB*GAMMA(C**(-l.)tl.)

DQ=SQRT(A**2t2.*A*B*GAMMA(C**(-l.)tl.)tB**2*GAMMA(2.*C**(-l.)tl.))

BASAL(2)=0.

KDFST=O

LINCK=O

ITREES(41)=0

C -------

------- FOR EACH DBH CLASS (l-40) MAKE VARIOUS CALCULATIONS. THE ARRAY

; ----- -- YIELD HAS TWO ROWS AND 9 COLUMNS. ROW 1 CONTAINS INFORMATION FOR

C ------- THE DBH CLASS, ROW 2 HOLDS TOTALS OVER ALL CLASSES. COLUMN 1

C ------- CONTAINS AVERAGE TREE HEIGHT, COLUMNS 2-3 HOLDS VOLUME FROM

GROUND TO TIP, AND COLUMNS 4-9 HOLDS MERCHANTABLE VOLUME IN

CUBIC FEET.

1 1
DO 90 IDBH=1,40

DBH=IDBH

-------

;---+- -- CALL FUNCTION ROUTINE TO COMPUTE AVERAGE HEIGHT OF TREE

C -------

yIELD(l,l)=MEANHT(DBH,DMAX,AGE,H,ANNOW,KSPCE)+*5

FREQ=FLOAT(ITREES(IDBH))

IF (FREQ) 22,22,23

22 IF (ITREES(41)) 90,90,95

E -------

------- SET AVERAGE HEIGHT TO 5 IF IT IS LESS THAN 5

C ------- 23 IF (YIELD(l,l)-5) 24,25,25

24 YIELD(l.1)=5 ------..

------- ACCUMULATE TOTAL NUMBER STEMS PER ACRE IN ITREES(41) AND TOTAL

C ------- BASAL AREA IN BASAL(2). ALSO CALCULATE BASAL AREA IN BASAL(l)

C -----a- AND CROWN RATIO.

25 ITREES(41)=ITREES(41)+ITREES(IDBH)

BASAL(1)=.00545415*DBH*DBH*FREQ

BASAL(z)=BASAL(2)+BASAL(l)

CRATIO=1OO.*(l.-EXP((DBH/DA)**ALAMDA*ALOG(1.-ACROWN/100.)))~.5

C -e----s

C ------- SELECT CROWN RATIO CLASS

C -------

IF (CRATIO-35) 26,26,27

26 ICRG=l

GO TO 30

27 IF (CRATIO-50) 28,28,29

28 ICRG=E

GO TO 30

29 ICRG=J

30 HT=FLOAT(YIELD(l,l))

LINCK=LINCK+l

LsLIN=(FLOAT(LINCK)/5.-FLOAT(LINCK/S))*lO.

C -------

------- CALL ROUTINE TO COMPUTE TOTAL VOLUME FOR THIS DBH CLASS

;- --- --

CALL TVOL (DBH,ICRG,HT,VOLT,KSPCE)

:

----_-_

------- STORE TOTAL VOLUME FOR DBH CLASS AND ACCUMULATE OVERALL TOTALS

C ---e-e- I

DO 35 J1=2,3

YIELD(l,Jl)=VOLT(Jl-l)*FREQ+.5

35 YIELD(2,Jl)=YIELD(2,Jl)+YIELD(l,Jl)

C -----em

C ------- IF DBH CLASS IS 5 OR GREATER, CALCULATE MERCHANTABLE VOLUMES AND

e-s----

: --- -- --

STORE, OTHERWISE SKIP CALCULATIONS

IF (IDBH-5) 40,45,45

40 JT1=3

GO TO 70

i

-------

------- CALL ROUTINE TO COMPUTE MERCHANTABLE VOLUMES FOR THIS DBH CLASS

C .e------

45 CALL MVOL (DBH,ICRG,HT,VOLM,KSPCE)

DO 65 J1=4,9

YIELD(l,Jl)=VOLM(Jl-3)*FREQ+.5

65 YIELD(2,J1)=YIELD(2,J1)+YIELD(1,Jl)

JT1=9

70 IF (KDFST) 75,75,80

75 WRITE (IOUT,106) IA,IH,IDBH,ITREES(IDBH),BASAL(1),CRATIO,

1 (~IELD(~,J~),J~=~,JT~)

KDFST-1

GO TO 85

80 WRITE (IOUT,107) IDBH,ITREES(IDBH),BASAL(l),CRATIO,

(YIELD(~ ~1) ~l=l JT~)

85lIF ((L5LIN.EQ.O).AND.(IT~EESiIDBH;2).NE.O).AND.(IDBH.LT.38))

108lFORMAT (1H )

WRITE (IOUT,108)

12
c-- -- -- -

c-------

C-------

c-------

c-------

CALL 'MERCHAWDIZER' ROUTINE TO CONVERT CUBIC FOOT MERCHANTABLE

VOLlRdES TO MERCHANTABLE VOLlJi'ES FOR FUEL CHIPS (TOWS)

PULPWOOD AND/OR CHIP'W'SAW LOGS (CORDS OR TOWS), AND SAWTIMBER

(TOWS, CUWITS, OR MBF DOYLE, INTERNATIONAL QUARTER INCH, OR

c------- SCRIBWER DEClMAL C LQG RULES)

c-------

CALL MRCHDZ(DBH,HT,ICRG,KS~E;ST~PG,COWV,IRULE,PULEMW,PUL~X,

1SAWTBR)

c-------

c------- FUEL CHIPS ARE THE DIFFERENCE BETWEEW THE TOTAL VOLME OB AND

c------- THE PULEWOOD, CHIP'W'SAW, AND SAWTlMBER VOLlJdES (ALL 3 IB AND

c-----mm ALL MEASURED IN CUBIC FEET). TO OMIT BARK VOLlhi'E FROY THE FUEL

c------- CHIP CALCULATIONS, CHANGE YIELD(l,2) TO YIELD(l,3/ BELOW.

C-------

FYIELDfIDBH)=FWAT(YIELD(l,2))/FREQ

c---w-sC-------END OF MERCHAWDIZER CALCULATIONS IN LOOP "Do 90"

c-------

90 CONTINUE

C -------

C ------- WRITE OUT OVERALL TOTALS FOR STEMS, BA, AND VOLUMES

C -----em

96 WRITE(IOUT,109) ITREES(41),BASAL(2),(YIELD(2,Jl),Jl=2,9)

TS=ITREES(41)

C -------

-------

; -- -- -- -

IF INPUT IS NUMBER PLANTED, CALCULATE PERCENT SURVIVAL,

OTHERWISE SKIP THIS CALCULATION

C ---e-eIF (KEY.EQ.3) GO TO 98

PCTSUR=TS/ANO*lOO.

WRITE (IOUT,llO) DA,DQ,A,B,C,PCTSUR,ACROWN,ALAMDA

GO TO 99

98 WRITE (IOUT,113) DA,DQ,A,B,C,ACROWN,ALAMDA

PCTSUR=O.

C -------

------- IF THERE IS AN ALTERNATE BASE AGE GIVEN AS INPUT, PRINT IT,

; - -- - -- - OTHERWISE SKIP THIS WRITE STATEMENT

C -------

99 IF (ABAGE-0.1 97,97,96

96 IF (AS.LT.6) GO TO 97

IABAGE=ABAGE

WRITE (IOUT,lll) IABAGE,IASI

97 WRITE (IOUT,105)

c-------

c------- WRITE 'MERCHAWDIZER' OUTPUT

c--es---

IF(KSECE.EQ.1) WRITE(IOUT,114)

IF(KSPCE.EQ.2) WRITE(IOUT,llS)

WRITE(IOUT,lOZ) IS, IBAGE

IF(KEY.EQ.3) GO TO 340

WRITE(IOUT,103) IN

GO TO 341

340 WRITE(IOUT,l12) IN

341 IFfKEY.EQ.1) WRITE(IOUT,116)

IF(KEY.EQ.2) WRITEfIOUT,117)

IFF(IRULE.EQ.-2) GO TO 370

c-I-----

c------- BUILD A TABLE HEADER FOR WG SCALE OPTION

c-------

WRITEIIOUT,l29)

129 FORMAT(//8X,7HGROIJING,l5X,38H-----------VOLUMES PER ACRE----------

ls7X,49H-------------STLrMPAGE VALUE PER ACRE-------------/8X,7HSEA

ZSOWS,SX, SHSTEMS,2X,4HFUEL,4X,8HPUPLWOOD,X,lOHCHIP'W'SAW,2X,9HSAWT

31MBER,7X,4HFIIEL,3X,8HPUL~OOD,2X,lOHCHIP'W'SAW,2X,9HSAWTIMBER,2X,9

4HDBH CLASS/SX,SHSIWCE,4X,3HDBH,3X,3HPER,3X, SHCUIPS,SX,SHBOLTS,SX,4

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

6TAL

IFCITRULEI 346,347,347

346 WRITE(IOUT,132)

1 3
132 FOHMAT(IOX,IHEST.,3X,SHCLASS,ZX,4HACRG;ZX,6HITONS),3X,7H(CORDS),3X

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

2(4),2X//)

GO TO 348

347 WRITE(IOUT,l34)

134 FORMAT(lOX,4HE'ST.,3X,SHCLASS,2X,4HAcR~,2X,6H(ToNs),3X, 7H(coRDS),3X

1, ?H(CORDS),3X,9H (MBF)

2H($),X/)

,SX,3H($),6X,3H($).8X,3H($),8X,3H($),8X,3

GO TO 348

C----___

C------ BUILD A TABLE HEADER FOR WEIGHT SCALE OPTION

C-----__

370 WRITE(IOUT,l40)

140 FORMAT(//8X,7HGROWING,lSX,38H-----------WEIGHTS PER ACRE----------

l-,7X,49H -------------STUMPAGE VALUE PER ACRE-------------/8X,7HsEA

2SONS,8X, 5HST~MS,2X,4HFUEL,4X,8HPULWOOD,X,lOHCHIP'N'SAW,2X,9HSAWT

3~BER,7X,4HFVEL,3X,8HPULpWOOD,2X,lOHCHIP'N'SAW,2X,9HSAWTIMBER,2X,9

4HDBH cLASS/9X,SHSINC$4X,3HDBH,3X,3HP~R,3X,SHCHIPS,5X,SHBOLTS,5X,4

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

~TAL/~OX,~HEST..~X,~HCLASS,~X,~HACRE,~X,~H(TONS),~X,~HITONS),~X,~H(

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

348 WRITE(IOUT,l35) IA

135 FORMAT(8X,I4)

C- -- - - - -

c------- CALCULATE MERCHANDILED STUMPAGE VALUES

C-- - - - - -

Do 357 J=l,9

STUMPV(2,J)=O.O

357 CONTINUE

Do 360 IDBH=l,lMAX

IF(ITREES(IDBH).EQ.O) GO TO 360

Do 358 J=1,9

STUMPV(l,J)=O.O

358 CONTINUE

TREEz=FLOAT(ITREES(IDBH))

STUMPV(~,~)=TREEZ*(FYIELD(IDBH)-STUMPG(IDBH,~)

1 -STUMPG(IDBH,2)-STUMPG(IDBH,3))/CONV(l)

IF(STUMPV(l,l).LT.O.0) STUMPV(l,l)=0.0

IF(IRULE.EQ.-2) GO TO 371

~TUMP~(~,~)=TREEZ~STUMPG(IDBH,~)/~~N~(~)

STUMPV(l,3)=TREEZ*STUMPGfIDBH,2)/CONV(3)

IF(IRULE) 351,352,353

351 STUMPV(l,4)=TREEZ*STUMPG(IDBH,3)/100.

GO TO 354

352 STUMPV(1,4)=TREEZxST&'PG(IDBH,4)/1000.

GO TO 354

353 STUMPV(l,4)=TREEZ"ST~PG(IDBH,5)/1000.

IF(IRULE.E&.~) STUMPV(l,4)=TREEZ*STUMpG(IDBH,6)/1000.

GO TO 354

c------

C ------- CALCULATE WEIGHTS AND VALVES FOR WEIGHT SCALE OPTION

C------- WEIGHT OF THE WOOD INT HE DIAMETER CLASS Is CALCULATED BY

C------- PRODUCT AND THE TOTAL WEIGHT INCREASED 10 PERCENT To ACCOUNT

C-- -- -- - FOR BARK WEIGHT.

Cm------

371 STUMP~(~,~)=~.~*TREEZ*STL~~~~IDBH,~)/~~N~~~)

STUMPV(1,3)=l.l*TREEZ*ST~~(IDBH,2)/CONv(3)

STUMPV(l,4)=l.l*TREEZ*STU'@GfIDBH,3)/cONv(4)

C-------

c------- DEDUCT CULL FRC&l SAWTIMBER, CALCULATED AS A PERCENTAGE OF

C ------- SAWTlMBER CUBIC FOOT VOLUME, AND ASSIGN THE CULL VOLUME TO

c------- THE NEXT LOWER-VALUED PRODUCT MERCHANDIZED

C -------

354 IF(ST~PV(l,4).LE.0.000) GO To 361

IF(IDBH.GT.ICVLLD) GO TO 361

STUMPV(l,4)=ST~PV(l,4)"(1.0-CULLP)

CULVOL=TREEZ%S'Tlk'PG(IDBH,3)xCULLp

IF(ICVLLV) 375,376,377

375 IF(IRULE.EQ.-2)STMPV(l,l)=ST~PV(l,l)+(l.l*cULvoL/coNv(l))

IF(IRULE.NE.-2) STUMPV(l,l)=STl&lPV(l,l)+(cULvOL/coNv(l))

Go TO 361

1 4
376 IF(IRVLE.EQ.-2) STUMPV(l,2)=ST&‘PV(l, 2/+(1 .I *CVLVOL/CONVfl))

IF(IRVLE.NE.-2) STt!kfPV(l,2I=ST~PV(l,2~+(CVLVOL/CONV/l)~

GO TO 361

377 IF(IRVLE.EQ.-2) STUMPV(l,3)=STUk’PV(l,3)+(1 .I “CVLVOL/CONVfB))

IF(IRVLE.NE.-2) STUMPV(l,3)=STUMPVil,3~+fCVLVOL/CONV(2))

C m-v----

C m-w---- CALCULATE STUMFAGE VALVES AND AGGREGATE VOLUMES AND VALVES

C ---w-e- BY DBH CLASS

C-------

361 DO 355 J=5,8

STUMPV(1, J)=STUMPVIl, J-4) “VALVE (J-4)

IF(STUMPV(l,J) sLT.O.0) STlhlPV(l,JJ=O.

355 STUMPV(l,9~=STuMPV(l, 9)+STUMPV(l, J)

Do 356 J=1,9

356 STi%lPV(B, J)=STUMPV’(2, J)+STUMPVIl, J)

C-----

C -----me END OF MERCHANDIZER OUTPUT CALCULATIONS

C -m-w---

WRITE(IOVT, 136) IDBH, ITREESfIDBHI, (STUMPVfl, JI, J=l, 9)

136 FORMAT(l9X,I2,2X,I4,4X,F4.1,4X,F5.1,5X,F5.1,4X,F7.3,9X,F6.2,2X,F7.

12, 4X, F7.2,4X, F7.2,4X, F7.2/)

360 CONTINUE

WRITE (IOUT, 13 7)

137 FORMAT(/23X. 4H 3X, ----, 4H 5X, ----, 5H 5X, 5H-7y-;6X, _____, 5H--- -, 9X, 6

IH

W~~&~ :87----- ITRE-i?4;;TSTUifP$>, ----- _----- J=l, 97’ ’ H

-------71)

138 FO~AT(23X,I4,4X,F4.1,4X,F~.l,5X,F5.1,4X,F7.3,9X,F6.2,2X,F7.2,4X,F

17.2, 4X, F7.2,4X, F7.21

C-----

C ------- END OF MERCHANDIZER OUTPUT

C ----me100 CONTINUE

GO TO 230

220 WRITE (IOUT,~~~)

106 FORMAT (8X,I3,4X,I4,2X,I3,3X,I4,3X,F6.1,3X,I3,3X,I3,2(5X,I5,2X,I5~

107lFORMAT ('

2(7X,15,2X,15))

21X,I3,3X,I4,3X,F6.1,3X,I3,3X,I3,2(5X,I5,2X,I5~

1 ,2(7X,15,2X,15))

109 FORMAT (1H ,26X,4H----,3X,sH------,17X,2(5H-----,2X,5H-----,5X),

1 2X,2(5H-----,2X,SH-----,7X)/27X,I4,3X,F6.1,17X,

2(15 2X 15 5X) 2X 2(15 2X 15 7X))

1lO'FORMAT (1HO ;2X:16iARI:H.‘MEAN’DBti= ;5 1 25H

lBH=,F5.1,27; WEIBULL PARAMETER: A-I f4.1,;H

QUADRATIC MEAN D

2=,F5.2/30X,17;PERCENT SURVIVAL=,F5.1,25; ,

B=,F6.2,6H C

MEAN CROWN RAT;O=,

3F5.1,14H , LAMBDA=,F6.3)

111 FORMAT (1H ,43X,37HCORRESPONDING SITE INDEX FOR BASE AGE,IS,lH=,I3

113l:ORMAT (1HO lPX,16HARITH.

lBH=,F5.1,27;

MEAN DBH= F5 1 25H

WEIBULL PARAMETERS A-I ;4.1,;H

QUADRATIC MEAN D

2=,F5.2/46X,17bMEAN CROWN RATIO=,F5.1,14;

B=,F6.2,6H , C

LAMBDA= F6.3)

118 FORMAT (lH1,43HAGE LIMITS ARE TOO SMALL -1 PROGRAM ABOiTED)

230 STOP

END

C -------

\ BLOCK DATA

COMMON /TAPER/COEFSS(3,14),COEFSL(3,14)

C -------

C******* TAPER COEFFICIENTS

C -------

------- INITIALIZE SLASH TAPER COEFFICIENTS -- BY CROWN RATIO GROUP

BETAlOB - B(1)

C ------- BETAPOB - B(2)

BETA30B - B(3)

BETA4OB - B(4)

ETAOB - ETAOB

GAMMAlOB - G(1)

C ----e-w GAMMA20B - G(2)

C ------- BETAlIB

15
----e-m

; - - -- - --

BETAZIB

BETA3IB

-e----s

; - ----- -

BETA4IB

ETAIB

C ------- GAMMAlIB

C -----mm GAMMA2 IB

c"

------- COLUMN 1

------- COLUMN 2

C ------- COLUMN 3

C -------

DATA COEFSS/

1 19.685000,

5 oo:Ei?

4 -13.065000:

:

-0.135320,

0.000000,

7 1.000000,

8 10.991000,

9 0.524640,

A 0.022500,

F

-6.630700,

-0.109430,

0 -0.481480,

E 0.892630,

C -------

UF COEFS IS CROWN RATIO LESS THAN 36

OF COEFS IS CROWN RATIO 36-50

OF COEFS IS CROWN RATIO GREATER THAN 50

9.142700, 6.062700,

1.634400, 1.919300,

-0.041715, -0.068244,

-8.588700, -7.181500,

-0.135320, -0.135320,

0.000000, 0.000000,

1.000000, 1.000000,

2.739100, 0.774720,

1.495400, 1.610900,

-0.036476, -0.053570,

-2.873900, -1.801600,

-0.109430, -0.109430,

-0.481480, -0.481480,

0.892630, 0.8926301

C ------- INITIALIZE LOBLOLLY TAPER COEFFICIENTS -- BY CROWN RATIO GROUP

C -------

DATA COEFSL/

1 21.901200, 14.118000, 9.610800,

5 0.059713, 0.130381, -0.015343, 1.164817, -0.556774, 0.062742,

4 -15.313400, -12.350600, -4.239100,

5 -0.101050, -0.110650, -0.122970;

6 0.000000, 0.000000, 0.000000,

; lEE?

1.000000, 1.000000,

9.573100, 6.278900,

9 -0.150774: 1.153531, -0.567733,

A 0.089512, -0.012421, 0.077350,

c"

-9.852540, -7.888530, -1.695430,

-0.102870, -0.103780, -0.103700,

0 -0.462060, -0.462060, -0.462060,

E 0.931570, 0.931570, 0.931570/

END

C -se----

C******* NUMBER TREES SURVIVING, IDEAL AND C-S

-------

; --- - -- - ARGUMENTS:

-------

; --- - - - -

(1) H = AVERAGE HEIGHT OF DOMS AND CODOMS

(2) AGE = AGE OF TREES

C ------- (3) AN0 = NUMBER STEMS PLANTED

C ---w---

-----mm

c"- --- - - -

["S', ;;pyC = SPECIES

= SELECTED SURVIVAL MODEL

FUNCTION SURVIV(H,AGE,ANO,KSPC,KEY)

DIMENSION BETA(3,2),BETAG(3,2)

C -------

------- IDEAL MODEL COEFFICIENTS

; -- - -- -- ROW 1 IS SLASH

C ------- ROW 2 IS LOBLOLLY

C -------

DATA BETA/

l-0.015038,-0.00063965, +0.0095738,

2-0.013480,-0.00060783, +0.0084124/

C --w--m-

------- C-S MODEL COEFFICIENTS

c" ---- - - - ROW 1 IS SLASH

C ---e--e ROW 2 IS LOBLOLLY

16
C -------

DATA BETAG/ :

2.3949, -1.2505,

2.1346, -1.1103,

GO TO (5,10),KEY

C -----mm C ------- COMPUTE SURVIVAL FOR IDEAL MODEL C -------

5 SURVIV= ANO*lO.**( BETA(l,KSPC)*AGE*ALOGlO(ANO) 1 +BETA(2,KSPC)*AGE*H t BETA(3,KSPC)*AGE*SQRT(H))

RETURN

C -------

------- COMPUTE SURVIVAL FOR C-S MODEL ; - --- ---

10 SURVIV=10.**(ALOG10(ANO)t(BETAG(1,KSPC)tBETAG(2,KSPC) 1 *ALOGlO(ANO)tBETAG(3,KSPC))*AGE/lOO.)

RETURN

END

C -------

C******* MEAN HEIGHT OF DOMS AND CODOMS AND SITE INDEX

------- ; ------ 

- ARGUMENTS:

-------

;- -- ---

(1) H = AVERAGE HEIGHT OF DOMS AND CODOMS

(2) AS = AGE OF TREES

------- ; - - - - - -

-

I.3; ;;GE = BASE AGE

= SITE INDEX

C ------- (5) KSPC = SPECIES C ------- (6) ABAGE = ALTERNATE BASE AGE

------- ; - - - - -- -

(7) IAST = ALTERNATE SITE INDEX

SUBROUTINE HEIGHT(H,AS,BAGE,SI,KSPC,ABAGE,IASI)

DIMENSION BETA(4,2),COL(4)

C -------

------- BETA COEFFICIENTS

c" -- - -- -- ROW 1 IS SLASH

C ------- ROW 2 IS LOBLOLLY C ----..--

DATA BETA/ :

-8.80405 , 22.7952 , 0.0000 , 0.0000

-21.0977 , 316.282 ,-2443.85 , 6318.86 ; C --e---e

------- IF AGE LT 6 THEN USE INTERPOLATION ; - - -- - --

IF (AS.GE.6) GO TO 4

DO 2 J=1,4

2 COL(J)=BETA(J,KSPC)

CALL HTSI(COL,BAGE,AS,SI,1,6.O,H)

IF (ABAGE-0.) 10,10,5

C ------- C ------- COMPUTE AVERAGE HEIGHT C -e---e4 H=SI*lO.**(BETA(1,KSPC)*(1./AS-1./BAGE)tBETA(2,KSPC)*((l./AS)**2

l-(l./BAGE)**2)tBETA(3,KSPC)*((l./AS)**3-(l./BAGE)**3)

2tBETA(4,KSPC)*((l./AS)**4-(l./BAGE)**4))

C ------- C ------- COMPUTE ALTERNATE SITE INDEX IF ALTERNATE BASE AGE IS GIVEN C -------

IF (ABAGE-0.) 10,10,5

5 IASI=H*lO.**(BETA(1,KSPC)*(1./ABAGE-1./AS)tBETA(2,KSPC)*((l./ABAGE

1)**2-(l./AS)**2)tBETA(3,KSPC)*((l./ABAGE)**3-(l./AS)**3)

2tBETA(4,KSPC)*((l./ABAGE)**4-(1./AS)**4))t.5

10 RETURN

END

SUBROUTINE HTSI(CO,AI,X,SH,I,AL,H)

REAL C0(4),AI,X,SH,I,AL

Vl=CO(l)/AItCO(2)/AI**2tCO(3)/AI**3tCO(4)/AI**4

IF (X.GE.AL) GO TO 10

IF (I.LT.0) GO TO 20

V2=CO(l)/ALtCO(2)/AL**2tCO(3)/AL**3tCO(4)/AL**4

TSH=(X*(SH*lO.**(V2-Vl)))/AL

H=TSH

RETURN

1 7
10 v2=CO(1)/X +CO(2)/X**2 + co(3)/x**3 + co(4)/x**4

TSH=(SH*lO.**(I*(V2-Vl)))

H=TSH

RETURN

2 0 H=O.

RETURN

END

C -----em

C******* MERCHANTABLE VOLUME 06 AND IB FROM .5 STUMP TO 2, 3, AND

--e---w

; - - - - - - -

4 INCH TOP OB

__-----

; - - -- - --

ARGUMENTS:

( 1) D = DBH OB

C ---..-a- = CROWN RATIO GROUP

-------

; - - - - - - -

I i; ERG = TOTAL HEIGHT OF TREE

( 4) VOLM(1) = OB VOLUME TO 2 INCH TOP OB

C -----em ( 5) VOLM(2) = IB VOLUME TO 2 INCH TOP OB

-------

c"- - - -- - -

( 6) VOLM(3) = OB VOLUME TO 3 INCH TOP OB

( 7) VOLM(4) = IB VOLUME TO 3 INCH TOP OB

C----e-m ( 8) VOLM(5) = OB VOLUME TO 4 INCH TOP OB

C --..---- ( 9) VOLM(6) = IB VOLUME TO 4 INCH TOP 08

-------

; -- ..--- -

(la) KSPCE = SPECIES

SUBROUTINE MVOL(D,ICRG,HT,VOLM,KSPCE)

DIMENSION VOLM(6),B(4),G(2)

COMMON /TAPER/COEFSS(3,14),COEFSL(3,14)

HL=.5

C -----em

; ---- - --

------- HL (STUMP HEIGHT) COULD BE DEFINED AS A FUNCTION, SAY OF DBH

DO 30 L=1,3

GO TO (5,6,7),L

5 GD=2

GO TO 10

6 GD=3

GO TO 10

7 GD=4

C --e---e

C ------- ARRAY B CONTAINS BETA COEFFICIENTS, FIRST FOR OB, THEN FOR IB.

C --w---e ETA HOLDS THE OB ETA, THEN THE IB ETA.

-------

c"- - - ---

ARRAY G CONTAINS GAMMAS FOR OB, THEN FOR IB.

10 DO 30 1=1,2

K=(I-1)*7

IF(KSPCE.EQ.2) GO TO 21

DO 15 J=1,4

15 B(J)=COEFSS(ICRG,K+J)

ETA=COEFSS(ICRG,K+5)

DO 20 J=6,7

20 G(J-5)=COEFSS(ICRG,K+J)

GO TO 24

21 DO 22 J=1,4

22 B(J)=COEFSL(ICRG,K+J)

ETA=COEFSL(ICRG,K+5)

DO 23 J=6,7

23 G(J-5)=COEFSL(ICRG,K+J)

24 GO TO (25,26),1

C -------

C ------- CALL ROUTINE TO CALCULATE HEIGHT OF TREE, GIVEN THE DIAMETER

C -------

25 HU=HGD(D,HT,GD,ETA,G,B)

26 M=2*(L-1)

cc

-------

------- CALL ROUTINE TO CALCULATE VOLUME BETWEEN TWO HEIGHTS. LOWER

---m---

;--L - - -

HEIGHT IN THIS CASE IS ALWAY

f

.5 (STUMP HEIGHT).

VOLM(M+I)=VIOB(D,HT,HL,HU,ETA,G,B)

30 CONTINUE

RETURN

END

1 8

, /
------- cc******* TOTAL VOLUME FROM GROUND TO TIP OB AND ALSO IB

-------

;------ ARGUMENTS:

-------

~----

(1) D =

C -----mm

------- ; - - - - - -

-

-------

E --- ----

DBH OB

CROWN RATIO GROUP

TOTAL HEIGHT OF TREE

VOLUME OB

VOLUME IB

SPECIES

SUBROUTINE TVOL(D,ICRG,HT,VOLT,KSPCE)

DIMENSION VOLT(2),B(4),G(2)

COMMON /TAPER/COEFSS(3,14),COEFSL(3,14)

HL=O.

HU=HT

C ------- C ------- ARRAY B CONTAINS BETA COEFFICIENTS, FIRST FOR OB, THEN FOR IB.

------- ; ---- - - -

ETA HOLDS THE OB ETA, THEN THE IB ETA.

ARRAY G CONTAINS GAMMAS FOR OB, THEN FOR IB.

C -------

DO 30 1=1,2

K=( I-1)*7

IF(KSPCE.EQ.2) GO TO 20

DO 10 J=1,4

10 B(J)=COEFSS(ICRG,KtJ)

ETA=COEFSS(ICRG,Kt5)

DO 15 J=6,7

15 G(J-5)=COEFSS(ICRG,KtJ)

GO TO 28

20 DO 22 J=1,4

22 B(J)=COEFSL(ICRG,KtJ)

ETA=COEFSL(ICRG,Kt5)

DO 24 J=6,7

24 G(J-5)=COEFSL(ICRG,KtJ)

c"

-------

------- CALL ROUTINE TO CALCULATE VOLUME BETWEEN TWO HEIGHTS

C -------

28 VOLT(I)=VIOB(D,HT,HL,HU,ETA,G,B)

30 CONTINUE

RETURN

END

----m-e

;******* MEAN HEIGHT OF TREE

-------

E -- --- -- ARGUMENTS:

C ------- (1) DBH = DIAMETER CLASS

-------

;---- --

(2) DMAX = LARGEST DIAMETER CLASS WITH COMPUTED VALUE

(3) AGE = AGE OF TREES

-------

cc--- --- -

(4) H = AVERAGE HEIGHT OF DOMS AND CODOMS

C -------

1;; ;li;;I li;M3f3;SSTEMS,SURVIVING

C ---..---

FUNCTION MEANHT(DBH,DMAX,AGE,H,ANOW,KSPC)

DIMENSION BETA(6,2)

-------

------- BETA COEFFICIENTS

C ------- ROW 1 IS SLASH C ------- ROW 2 IS LOBLOLLY

DATA BETA/

1+0.050341,3.1868,1.5708E-5,0.0114942,-2.0981,

2 1.4034

3 0.032876:1.9930,2.5047E-5,0.0043249,-1.0360,

/ 4-0.026038/

MEANHT=10.**(ALDG10(H)tBETA(1,KSPC)-(1./DBH-l./DMAX)*

1 (BETA(2,KSPC)tBETA(3,KSPC)*AGEXANOWtBETA(4,KSPC)*ANOW/AGE 5 tBETA(5,KSPC)*ALOGlO(ANOW/AGE)tBETA(6,KSPC)*ALOGlO(H/AGE)))t .5

RETURN

END

19
C******* WEIBULL PARAMETER ESTIMATES

C --e-"-m

-----me

;- - - -- - -

ARGUMENTS:

-----es

;-- - - - - -

1;; ;GE = AGE OF TREES

= AVERAGE HEIGHT OF OOMS AND CODOMS

-----me

c"- - - -- -_

I;; ;NT = NUMBER STEMS SURVIVING

= WEIBULL LOCATION PARAMETER

(5) B = WEIBULL SCALE PARAMETER

C -----em (6) C = WEIBULL SHAPE PARAMETER

---m--m

; - - - - - - -

(7) KSPC = SPECIES

SUBROUTINE PARMS(AGE,H,ANT,A,B,C,KSPC)

DIMENSION BETA(4,3,2),PRE(3)

----w-m

~--- - BETA COEFFICIENTS

C ----e-e 1ST DIMENSION OF ARRAY IS 1 FOR SLASH, 2 FOR LOBLOLLY

----s-m

; ---- - - -

2ND DIMENSION OF ARRAY IS 1 FOR PARAMETER A, 2 FOR AtB,

3 FOR C

DATA BETA/

1 1.3986 , 2.9217 , -1.1126 , -1.8477 ,

2 2.5800 , 10.1380 , -3.6275 , -2.5005 ,

3 9.1471 , 6.5959 , -2.4479 , -7.6706 ,

4 0.13391 , -0.076245, -0.70442 , 2.3569 ,

5 2.14570 , 8.1415 , -4.13600 , 2.1264

6 2.69350 , 5.3594 , -0.53123 , -5.1274 ;

DO 10 K=1,3

C -------

C ------- PARAMETERS A, AtB, AND C ARE COMPUTED BELOW

C L------

10 PRE(K)=BETA(l,K,KSPC)tBETA(2,K,KSPC)*ALOGlO(H)tBETA(3,K,KSPC)*ALOG

llO(ANT)+BETA(4,K,KSPC)*ALOGlO(AGE)

IF (PRE(2).LT.0.02) PRE(2)=0.02

IF (PRE(3).LT.l.l) PRE(3)=1.1

A=PRE(l)

IF (A.LT.o.) A=.OI

B=PRE(2)-PRE(l)

C=PRE(3)

RETURN

END

C ..------

C******* EVALUATES THE GAMMA FUNCTION (GENERALIZED FACTORIAL) OF DN

C ------_

FUNCTION GAMMA(DN)

DLGGM=O.

IF (DN) 90,90,91

91 DTERM=O.

DNX=DN

7 IF (DNX-20.) ll,lO,lO

11 DTERM=DTERMtALOG(DNX)

DNX=DNXtl.

GO TO 7

10 DLGGM=(DNX-.5)*ALOG(DNX)-DNXt1./(12.*DNX)

:

-1./(36.El*DNX**3)+1./(1.26E3*DNX**5)-1./(1.68E3*DNX**7)

+1./(1188.E1*DNX**9)-6.91E2/(36036.El*DNX**ll)

3 +1./(1.56E2*DNX**13)+0.918938533204673-DTERM

GAMMA=EXP(DLGGM)

99 RETURN

90 WRITE (6,100) DN

100 FORMAT (lHO,lOX,57HARGUMENT GIVEN TO GAMMA FUNCTION .LE. ZERO, ARG

1UMENT = ,F15.8)

GO TO 99

END

C -------

C******* DIAMETER AT GIVEN HEIGHT ABOVE GROUND EITHER OUTSIDE OR

C ----e-m INSIDE BARK DEPENDING ON ARGUMENTS 4, 5, AND 6

-------

c" - -- - - - - ARGUMENTS:

C -----mm 1 1) D = DBH OB

C ------- ( 2) HT = TOTAL HEIGHT OF TREE

20
-------

~---~

( 3) H = HEIGHT ABOVE GROUND

( 4) ETA = PARAMETER OF LOWER STEM CURVE

-----me

;- ---- - -

( 5) G = PARAMETERS OF BARK RELATIONS AT 4.5 FEET

( G(l)=O, G(2)=1 FOR OB DIAMETER)

-------

; - - -- -- -

( 6) B = PARAMETERS OF UPPER STEM CURVE

FUNCTION DIO(D,HT,H,ETA,G,B)

DIMENSION G(2),8(4)

IF(H.GT.4.5) GO TO 10

C -------

------- CALCULATE DIAMETER IF TREE HEIGHT IS'BELOW BREAST HEIGHT

;- - -- - --

DIO=(G(l)tG(2)*D)*(H/4.5)**ETA

RETURN

---a---

; --- - - - - CALCULATE DIAMETER IF TREE HEIGHT IS ABOVE BREAST HEIGHT

C ----m-s

10 FD=G(l)tG(2)*D

FHl=(HT-H)/(HT-4.5)

FH2=(HT-H)*(H-4.5)/(HT*HT)

FH3=FD*FH2

FH4=FD*FH3

FH5=(FH2*(HTtHT-H-4.5))/HT

DIO=FD*FH1tB(1)*FH2tB(2)*FH3tB(3)*FH4tB(4)*FH5

RETURN

END

C -------

C******* HEIGHT TO GIVEN DIAMETER EITHER OUTSIDE OR INSIDE BARK

C -m-s--- DEPENDING ON ARGUMENTS 3-6

-------

; - - -- - - -

HEIGHT IS CALCULATED BY AN ITERATIVE PROCESS, COMPUTING

DIAMETER MIDWAY BETWEEN GRADUALLY DECREASING UPPER AND

C ----e-e GRADUALLY INCREASING LOWER HEIGHT VALUES UNTIL THE COMPUTED

--w-m--

; -- - -- - -

DIAMETER IS FINALLY VERY CLOSE TO THE GIVEN DIAMETER

C -----em ARGUMENTS:

-------

; --- - - - -

(11 D = DBH OB

(2) HT = TOTAL HEIGHT OF TREE

C ---e-v- (3) GD = GIVEN DIAMETER

-------

c" - - - - - - -

I;; ;TA = PARAMETER OF LOWER STEM CURVE

= PARAMETERS OF BARK RELATIONS AT 4.5 FEET

-------

; - - - - - --

( G(l)=O, G(2)=1 FOR OB DIAMETER)

(6) B = PARAMETERS OF UPPER STEM CURVE

C -------

FUNCTION HGD(D,HT,GD,ETA,G,B)

DIMENSION G(2),8(4)

C -------

C ------- DETERMINE WHETHER GIVEN DIAMETER IS ABOVE OR BELOW BREAST HEIGHT

C -e--.--e

DD=G(l)tG(2)*D

IF(GD-DD) 20,80,10

C -------

C ------- GIVEN DIAMETER IS BELOW BREAST HEIGHT HERE

C ---I---

10 HU=4.5

HL=O.

GO TO 30

C -------

-------

; - -- - - - -

GIVEN DIAMETER IS ABOVE BREAST HEIGHT HERE

20 HU=HT

HL=4.5

30 H=(HU-HL)/2.+HL

C e-e----

-------

; - - -- - --

CALL ROUTINE TO CALCULATE DIAMETER AT A GIVEN HEIGHT

DB=DIO(D,HT,H,ETA,G,B)

IF(ABS(DB-GD)-.05) 70,70,40

40 IF(DB-GD) 50,50,60

50 HU=H

GO TO 30

-1

21
60 HL=H

GO TO 30

70 HGD=H

RETURN

80 HGD=4.5

RETURN

END

C ---^---

C******* VOLUME BETWEEN TWO HEIGHTS ABOVE GROUND EITHER OUTSIDE

----_--

; -- - - - - -

OR INSIDE BARK DEPENDING ON ARGUMENTS 5, 6, AND 7

C ----m-w ARGUMENTS:

e------

;-- - -- -

( 1) D = DBH OB

( 2) HT = TOTAL HEIGHT OF TREE (

C ------- ( 3) HL = LOWER HEIGHT ABOVE GROUND

w------

; -- -- - - -

( 4) HU = UPPER HEIGHT ABOVE GROUND

-------

;- - -- - --

; ;; ;TA = PARAMETER OF LOWER STEM CURVE

= PARAMETERS OF BARK RELATIONS AT 4.5 FEET

( G(l)=O, G(2)=1 FOR OB VOLUME)

s------

c"- - - - - - -

( 7) B = PARAMETERS OF UPPER STEM CURVE

FUNCTION VIOB(D,HT,HL,HU,ETA,G,B)

DIMENSION G(2),8(4)

VIOB=O.

HNP=HU

HNV=HL

C -------

c" ..- - - -_ -

------- IF UPPER HEIGHT IS BELOW BREAST HEIGHT, THEN VOLUME IS

CALCULATED FOR LOWER PORTION OF TREE

C ------- IF LOWER HEIGHT IS BELOW BREAST HEIGHT AND UPPER HEIGHT IS

C ------- ABOVE BREAST HEIGHT, VOLUME IS CALCULATED FOR LOWER PORTION,

-------

;-----L -

THEN ADDED TO UPPER PORTION VOLUME

IF(HU.LT.4.5) GO.TO 10

IF(HL.LT.4.5) GO TO 20

C -------

C ------- CALCULATE VOLUME ABOVE BREAST HEIGHT

C -----L30 HTSQ=HT*HT

HT45=HT-4.5

HT3=HTSQ*HT

DSQ=D*D

DHT=D/HT

D2HT=DSQ/HT

DHT2=D/HTSQ

D2HT2=DSQ/HTSQ

RO=G(l)*HT/HT45+G(2)*HT*D/HT45-4.5*B(l)/HT-4.5*G(l)*B(2)/HT-4.5*G(

12)*B(2)*DHT-4.5*G(l)*G(l)*B(3)/HT-g.*G(l)*G(2)*B(3)*DHT-4.5*G(2)*G

2(2)*9(3)*D2HT-9.*B(4)/HTt20.25*B(4)/HTSQ

Rl=-G(l)/HT45-G(2)*D/HT45tB(l)/HTt4.5*B(l)/HTSQtG(l)*B(2)/HTt4.5*G

l(l)*B(2)/HTSQtG(2)*B(2)*DHTt4.5*G(2)*B(2)*DHT2tG(l)*G(l)*B(3)/HTt4

2.5*G(l)*G(l)*B(3)/HTSQt2.*G(l)*G(2)*B(3)*DHTt9.*G(l)*G(2)*B(3)*DHT

32tG(2)*G(2)*B(3)*D2HTt4.5*6(2)*6(2)*8(3)*D2HT2t2.*B(4)/HT+g.*B(4)/

4HTSQ-20.25*B(4)/HT3

R2=-(B(l)/HTSQtG(l)*B(2)/HTSQtG(2)*B(2)*DHT2tG(l)*G(l)*B(3)/HTSQt2

l.*G(l)*G(2)*B(3)*DHT2+G(2)*G(2)*B(3)*D2HT2t3.*B(4)/HTSQ)

R3=B(4)/HT3

HNP2=HNP*HNP

HNP3=HNP*HNP2

HNP4=HNP*HNP3

HNP5=HNP*HNP4

HNP6=HNP*HNP5

HNP7=HNP*HNP6

HNV2=HNV*HNV

HNV3=HNV*HNV2

HNV4=HNV*HNV3

HNV5=HNV*HNV4

HNV6=HNV*HNV5

HNV7=HNV*HNV6

VIOB=5.454l538E-3*(RQ*RO*(HNP-HNV)tRO*R1*(HNP2-HNV2)t(Rl*Rlt2.*RO*

2 2
lR2)/3.*(HNP3-HNV3)t(RO*R3tRl*R2)/2.*(HNP4-HNV4)t(R2*R2t2.*Rl*R3)/5

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

RETURN

C -------

C ------- CALCULATE VOLUME OF TREE BELOW BREAST HEIGHT ONLY

C -------

10 VIOB=.00545415*(G(1)tG(2)*D)**2/(4.5**(2.*ETA)*(2.*ETAtl.))*(HNP

1**(2.*ETAtl.)-HNV**(2.*ETAtl.))

RETURN

C -------

C ------- CALCULATE VOLUME OF TREE BELOW BREAST HEIGHT AND TRANSFER TO

C ------- UPPER VOLUME COMPUTATION

C -------

20 HNP=4.5

VIOB=.00545415*(G(1)tG(2)*D)**2/(4.5**(2.*ETA)*(2.*ETAtl.))*(HNP

1**(2.*ETAtl.)-HNV**(2.*ETAtl.))

HNV=4.5

HNP=HU

GO TO 30

END

C --e-e..-

C******* DIAMETER DISTRIBUTION

C ----emARGUMENTS:

(1) AGE = AGE OF TREES

C ----se- (2) H = AVERAGE HEIGHT OF DOMS AND CODOMS

(3) IDST = DIAMETER DISTRIBUTION

(4) DMAX = LARGEST DIAMETER CLASS COMPUTED VALUE

;6"; YNOW = NUMBER STEMS SURVIVING

= WEIBULL LOCATION PARAMETER

C ------- (7) B = WEIBULL SCALE PARAMETER

C ------- (8) C = WEIBULL SHAPE PARAMETER

C ------- (9) KSPC = SPECIES

C -------

SUBROUTINE DIADST(AGE,H,IDST,DMAX,ANNOW,A,B,C,KSPC)

DIMENSION AIDST(41),IDST(41)

C -------

------- CALL ROUTINE TO GENERATE WEIBULL ESTIMATES

; ---- ---

CALL PARMS(AGE,H,ANNOW,A,B,C,KSPC)

cuM=o.o

NNOW=ANNOW+.5

ICUM=O

MODE=0

C -----em

------- COMPUTE DIAMETER DISTRIBUTION

;

1

4

11

DO 11 1=1,41

X=FLOAT(I)t.5

AIDST(I)=O

IF(ABS(CUM-ANNOW)-.001)1,3,3

IF(X-A)l,l,E

AIDST(I)=(l.-EXP(-((X-A)/B)**C))*ANNoW-CUM

CUM=CUMtAIDST(I)

IDST(I)=AIDST(I)t.S

ICUM=ICUMtIDST(I)

IF(IDST(I)-MODE)11,11,4

K=I

MoDE=IDST(I)

CONTINUE

IDST(K)=IDST(K)tNNOW-ICUM

DMAX=O.

C -------

C ------- DETERMINE LARGEST DIAMETER CLASS HAVING A COMPUTED VALUE

C -------

DO 34 J=1,40

K=41-J

ITWEN=AIDST(K)+.5

IF(ITWEN.EQ.0) GO TO 34

IF(DMAX)35,35,34
35 DMAX=K

34 CONTINUE

IDST(41)=NNOW

RETURN

END

c-------

C"""**"" MERCHANDIiiER: A ROUTINE TO CONVERT USLYCOWG YIELDS IN CUBIC

c------- FEET TO STUMPAGE VOLVMES IN OTHER UNITS

c-----

c----- ARGUMENTS:

c----- 11) D = DBH OB

c----- (2) HT = TOTAL HEIGHT OF TREE

c---- 131 ICRG = CROEJN RATIO GROUP

c- -- - - - - (4) KSPCE = SPECIES, l=SLASH, 2=LOBLOLLY

c----- (5) STVMPG(D,l) = PVLFWOOD STVMPAGE VOLVME,'TREE (CUEIC FT IBI

c--e--e_ STVMPGID,2) = CHIP'N'SAW STUMPAGE VOLmE/TREE (CUBIC FT IB)

c----- STVMPGfD,3) = SAWTIMBER VOLVME/TREE (CVNITS IB)

c----- STVMPG(D,4) = SAWTIMBER VOLUME/TREE (BOARD FEET DOYLE)

c------ STU'%?GfD,S/ = SAWTIMBER VOLUME/TREE (BOARD FEET INT'L l/4")

c-m - - - - - sTvMPG(D,6) = SAWTIMBER V~LVME~TREE (BOARD FEET SCRI~NER DEC.C.

c- - - -- - - (6) CONV(l) = FVE'L CHIP CONVERSION FACTOR (CVBIC FEET TO TONS)

c------_ CONVl2) = PVLpkrOOD CONVERSION FACTOR (CVBIC FEET TO CORDS OR

c---I - -- OR TONS, DEPENDING ON IRULE)

c- - - - - - - CONV13) = CHIP'N'SAW CONVERSION FACTOR (CVBIC FEET TO CORDS OR

c------ OR TONS, DEPENDING ON IRVLE)

c----- CONVf41 = SAWTIMBER CONVERSION FACTOR (CUBIC FEET TO TONS,

c-----mm USED ONLY WHEN IRVLE = -2)

c- - - - - - - (71 IRVLE = SAWTIMBER VOLUME RULE (-2, WEIGHT SCALE; -1, CVNITS;

C-em..--- (0, DOYLE; +I, INT'L l/4 INCH; +2, SCRIBNER DEC. Cl

cm-- - - - - (8) PULPMN = MINIMUM PULPWOOD DIAMETER, IB

c--m--e- (9) PULEMX = MAXIMUM PULFWOOD (MINIMUM CKIP'N'SAWJ DIAMETER, IB

c------- (10) SAWTBR = MINIMUM SAWTIMBER (MAXIMM CHIt"N'SAW) DIAMETER, IB

c-------

SUBROUTINE MRCHDZ(D,HT,.TCRG,KSPCE,ST~PG,CONV,IRVLE,PvL~N,

lPVLEMX,SAWTBR)

c-------

c---- ARRAY B CONTAINS BETA COEFFICIENTS FOR INSIDE BARK EQNS

c------- ETA HOLDS THE INSIDE BARK ETA

c-----we ARRAY G CONTAINS GAMMAS FOR INSIDE BARK

c------

DIMENSION Bi41,G(2),ST~PC(40,6/,CONV(II

CC?&'ON/TAPER/COEFFS~3,14),COEFSL(3,14)

.lD=D

DO 5 J=l,6

5 STVMPGIID,J)=O.OOO

IFfKS~E.EQ.2) GO TO 20

Do 10 J=1,4

10 BfJ)=COEFSSIICRG,7+5)

ETA-+?OEFSSIICRG,12)

Do 15 J=6,7

15 G1J-S)<OEFSSIICRG,7-kJ)

GO TO 28

20 Do 22 J=1,4

22 BfJI=COEFSLfICRG, 7+J)

ETA-&'OEFSL(ICRG,12)

DO 24 J=6,7

24 G(J-~~~oEFSL~I~RG,~+J)

C---_---

C ------- DIB AT THE TOP OF THE a-FOOT BUTT LOG DETERMINES INITIAL. PRODUCT

c------- IS THE BUTT LOG BIG ENOUGH FOR SAWTIMBER? IF SO, GO TO 310.

c------- IF NOT, IS IT BIG ENOUGH FOR A CHIP'N'SAW LOG? IF SO, W TO 320.

c ---w--v IF NOT, THE ENTIRE STEM IS PVLRJOOD AND FUEL CHIPS.

c ---m--m

28 H=8.75

DIB=DIOiD,HT,H,ETA,G,B)

IFIDIB.GE.SAWTBRl GO TO 310

IFIF(DIB.GT.PVLpMXI TO TO 320

PBAZ=O.S

GO TO 331

330 PBAZ=CBAZ

24
c

-----me

c

-we---- CALCULATE PULPWOOD CUBIC FOOT VOLUME AND CONVERT IT TO CORDS

C_------

331 HITE=HGD(D,HT,PULmN,ETA,G,B)

STUMpG(ID,l )=VIOB(D,HT,PBAZ,HITE,ETA,G,B)

RETURN

C---_-mm

C----m-m CALCULATE CHIP'N'SAW CUBIC FOOT VOLU%'E AND CONVERT IT TO CORDS

C---m-s320 CBAZ=O.S

321 HITE=HGDfD,HT,PULmX,ETA,G,B)

IF(HITE-CBAZ-8.25)330,322,322

322 BOLT=2"(AINT((HITEdBAZ)/2))

HU=CBAZ+BOLT

STUMPG(ID,2)=VIOB(D,HT,CBAZ,HU,ETA,G,B)

PBAZ=CBAZ+BOLTtO.25

GO TO 331

C_------

C ------- CALCULATE SAWTIMBER VOLL!blES

C-me----

310

311

312

HITE=HGD(D,HT,SAWTBR,ETA,G,B)

IF(HITE.LT.8.75) GO TO 320

BOLT=2"(AINTIHITE/B))

SBAZ=O.S

ST&%(ID,3)=VIOB(D,HT,SBAZ,BOLT,ETA,G,B)

CBAZ=BOLTtO.75

IF(IRULE) 321,312,312

IF(BOLT.GT.16.0) GO TO 313

C-------

C ------- BOARD FOOT VOLUME OF THE LAST SAWTIMBER BOLE SECTION (8 TO 16 FT)

C-----

H=BOLT+O.75

315 DIB=DIO(D,HT,H,ETA,G,B)

STUMPG(ID,4)=STUk'PG(ID,4)+((DIB-4.) "*21"BOLT/l6

ST~~(ID,5)=ST~~(ID,5)tf(BOLT/l6.)"f(0.796ADIB"*2)t

1(1.375*DIBI-1.23))

ST~C(ID,6)=ST~PC~ID,6/+(0.0494"BOLT~DIBA*2~-f0.124*DIB*BOLT~-

l(0.269"BOLT)

STUMPG(ID,6)=INT((STUMPG~ID,6)/10)+0.5)*10

CBAZ=H

W TO 321

C ----m-e

C------- BOARD FOOT VOL@JES OF 16-FOOT LOGS PROGRESSING UP THE BOLE

C-------

313 NLOGS=IFIX((HITE-0.5)/16.25)

Do 314 N=l,NLOGS

SLOGS=FLOAT(N)

H=O.St(SLOGS*l6.25)

DIB=DIO(D, HT, H, ETA,G, BI

STIYMPCIID,~)=STIIMPG(ID,~~~IDIB-~.I""~

ST~~(ID,5)=ST~~(ID,5)+((0.796~DIB~*"2)+(1.375~DIB~-l.23

314 ST~~(ID,6)=ST~~(ID,6)+(0.7904"DIB"~2)-(1.984"DIB)-4.304

BOLT=2*(AINT((HITE-H)/2))

H=HtBOLT

GO TO 315

END

2 5

Appendix B

Sample USLYCOWG Output With MERCHANDIZER Activated

Output from three runs is included. The first two batches of output

show sawtimber volumes in board feet International Quarter Inch Log

Rule. The initial run did not activate the cull option feature, the second

did. The third batch of output is from a run where the weight scaling

option was activated.

27
YIELDS BY DIANETER CLASSES FOR ONTHINNED LOBLOLLY PINE IN PLANTATIONS

ON CUT-OVER (UNPREPARED BUT PROBLEN FREE) SITES

IN THE

tiEST GULF COASTAL PLAIN

SAWTIMBER VOLUMES ARE CALCULATED FROM STUNP HELGHT (0.5 FEET) TO A 9.5 INCH TOP DIAMETER (fNS:DE BARK).

LOGS ARE BUCKED INTO 16 FOOT LENGTHS, PLUS A 3 INCH ALLOWANCE. LOG SCALING DIAMETER IS THE DIAMETER INSIDE

BARK AT THE SMALL END OF EACH LOG. IF THE LAST LOG BEFORE REACHING THE NINIAUH SCALING DIAWETEH EXCEEDS 8 FEET,

BUT IS LESS THAN i6 FEET, IT IS BUCKED TO THE NEXT SMALLER EVEN-NUMBERED FOOT OF LENGTH, PLUS TRIM ALLONANCE.

THE REMAINING LENGTH 1s DROPPED L&TO THE NEXT LOWER STUMPAGE VALUE CATEGORY. If THE LAST LOG TO THE MINIMUIll

SCALING DIAMETER IS LESS THAN 8 FEET LONG, THE ENTIRE LENGTH IS DROPPED INTO THE NEXT LOWER STUHPAGE VALUE

CATEGORY.

PULPWOOD BOLTS 3.0 - 7.5 80. CUBIC FEET

PER CORD

s 17.00 PER CORD

CHIP'N'SAW LOGS 7.5 - 9.5 90. CUBIC FEET

PER CORD

s 45.00 PER CORD

SAWTIMBER LOGS GREATER THAN

9.5 INCHES

BOARD FEET INT'L S 150.00 PER MBF
GROhUZNG

SEASONS AV.

SINCE Dtc

EST. HT.

***

58

ARITH.

STEMS

DBH PER

ACRE

3 3

4 1 1

5 24

6 42

7 59

8 68

9 65

10 49

1 1 29

12 13

13 4

14 1

"--*

368

BASAL

AREA

CR AV.

HT.

0.1 24 33

1.0 27 41

3.3 31 46

8.2 33 50

15.8 36 53

23.7 38 55

20.7 40 57

26.7 42 59

19.1 44 60

10.2 46 61

3.7

1.1

LIIIIl

141.7

40 62

49 63

LQBLOLLY PIKE

SITE INDEX 85 (F5bSE AGE 50)

726 STEMS PER bCi(E PLANTED

C-S SURVIVAL UUDEL

.-cII---.---------~-.~~-CUBIC FOQT VQLU~E--~-----r-rr-r----------

ALL TREES S-INCH CLASS AND GREATER, STUMP HEIGHT .5 FEET,

GROUlD FUR O.H. TOPS OF --

TO TIP 2 INCHES 3 INCRES 4 INCHES

0.8. 1.B. U.B. 1.B. O.B. 1.8. 0.B. 1.B.

***

3

20 1;

73 57 68 54 47 34 25

199 159 188 151 1:;: 143 145 115

408 329 398 314 376 305 349 282

643 526 614 505 604 497 579 475

810 670 777 646 769 639 747 621

783 653 753 631 747 626 734 615

573 481 552 465 548 462 541 456

311 263 300 254 299 253 296 251

114 97 110 94 110 94 109 93

34 29 33 28 32 28 32 28

-Wm.- -w-c- -c-". m-a-3 -w-w" -"m-w --I)-- w--w3971 3281 3783 3142 3722 3094 3566 2961

UEAN DBH= 8.1 , QUADRATIC MEAN DBH= 8.4 , WEIBULL PARAMETEHS A= 1.5 , 0= .7.39 , C = 3.62

PERCENT SURVIVAL= 50.7 , MEAN CROYM RATIO= 38.6 , LAMBDA= 0.597

CORRESPONDING SITE IIQDFX FOH BASE AGE 25= 58

***
,

GROLiXNG

SEASONS STE@!S

SINCK DBH

EST. CLASS

25

3

4

5

6

7

8

9

10

1 3

12

1 3

14

PCR

ACdE

FUEL PULP@000 CHIY'ti'SAW SAMTIMBER FUEL PULP#OOD CHIP'W'SAW

CHIPS bc)LTS LUGS LUGS CHIPS BOLTS LUGS

(TONS) CCOHDS) (cuRosl (#HF) (S) (Sl CS)

3 0.0 0.0

1 1 0.2 0.1

24 0.3 0.S

42 0.7 1.7

59 1.2 3.8

68 1.7 6.1

65 2.0 5.5

45 1.8 3.4

29 1.3 1.5

13 0.7 0.6

4 0.2 0.2

1 0.1 0.0

w--I) --m-w

3bi3 10.2 23.5

LOBLOLLY PINE

STTF. INDEX RS (BASE AGE 50)

726 STEMS PLR ACRE PLAPJTED

C-S SURVIVAL HODEL

-9a--------VOLUVES PER ACRE q---w---m-m -------------STUtiPAGE VALUE PER ACRE ~-~~-~~-~c1II

0.0

0.0

0.0

0.0

0 . 0

0.0

2.1

3.9

3.7

0.8

0.2

0 . 0

-m--w

10.8

SAWTIRlr3ER DBH CLASS

LOGS TOTAL

($1 (S)

0.000 0.17 0.00 0.00 0.00 0.17

0.000 0.83 1.06 0.00 0.00 1.89

0.000 1.69 9.03 0.00 0.00 10,73

0.000 3.40 29.28 0.00 0.00 32.68

0.000 5.93 64.02 0.00 0.00 69.95

0.000 8.40 104.51 0.00 0.00 112.91

0.000 10.06 93.48 94.54 0.00 198.07

0.000 9.06 57.16 175.46 0.00 241.68

0.000 6.36 25.98 168.14 0.00 200.48

1.084 3.56 10.65 37.30 162.61 214.12

0,453 1.24 2.90 10.08 67.90 82.12

0.158 0.37 0.60 2.06 23.65 26.68

.e-w---

1.694 51.08 398.67

-m---e487.58------- -w--v--

254.16 1191.48
YIELDS BY D1ArrlETEA CLASSES FUR UNTHINt’JED LOBLOLLY PINE IM PLANTATIONS

Or.4 CUT-CIVEH (UNPREPARED BUT PRO8LER FREE) SITES

XN ‘THE

*EST GULF COASTAL PLAIN

SAWTIMUER VOLUMES ARE CAkCULATED FROM STUMP tfEIGBT (0.5 FEET) TO A 9.5 INCH TOP DIAWETER (INSIDE BARK).

COGS AtiE BI!CI(ED IHTO 16 FCIDT LEKGTHS, PLUS A 3 INCti ALLOWANCE. LOG SCALING DIAfiETEA IS THE DIAMETER INSIDE

dAili( AT THE SMALL EPlU OF EACH LOG. IF THE LAST kDG BEFORE REACHLNG THE WINIKUM SCALING DIAMETER EXCEEDS 8 FEET,

RliT IS LFSS THAiV 14 FEET, IT X.5 BUCKED TO THE KEXT Sn4LLER EVEN-NUMBERED FUDT OF LENGTH, PLUS TRIW ALLOlwANCE.

Tnc: HEMA~~TNG LEKGTH Is DRDPPEIJ I&TO THF NEXT kDwCR STUKPAGE VALU1E CATECOitY. IF THE LAST LOG TO THE MIb11MiJti

SCALIi’dG DIAf4ETb.H IS LESS TtiAN H FEET LONG, THE EkTIRE LENGTH 1;s DROPPED INTO THE EvEXT LOWER STUMPAGE VALUE

CATEGUKX.

I; iit-itWOOD CHIPS LESS THAK 3.0 IWCHES 90. CUBIC FEET s 5.00 PER TON

NQ BRAkCHES PER TON

PULPWOOD BOLTS 3,O - 7.5 80. CUBPC FEET s 17.00 PER CORD

PER CORD

CrfIP’N’S43 LOtiS 7.5 - 9.5 90. CUBIC FEET 5 45.00 PER CORD

PER CURD

SPECIAL NOTE: THIS HUN HAS ACTIVATED THE SAdTIMBER CULL FEATURE OF HERCHANDIZER. 25 PERCERT OF THE SAWTIMBER

VOLlJhE IFI *:ACli DXAMfiTER CLASS i\F:LOr AND INCLIJDING THE 13 INCH CLASS HAS MEN ASSIGNEO TO THE

PUkPirODD VOLUME CATEGOHY BECAUSE, AS CULL VilLUfiE’, XT WAS DEEMED UNSUITABLE

‘t’fl PHUDUCE SAdf4 LUrtcER.
LOBLOLLY PINE

SITE INDEX 85 (BASE AGE 59)

726 STEMS PEK ACRE PLANTED

GROWING

SEASONS AV.

SINCE rJ+c

EST. HT.

***

25 58

STEMS

DRH PER

ACRE

BASAL

AREA

CR AV.

HT.

-,,.-,-----,-c-----.....CUB~C T‘OI-JT V()Ll.jME ---.-------“---.-.-.....

ALL TREUS S-INCH CLASS AEID GREATER, STUYP dElGtiT .5 FEET,

GROUND FOR 0.B. TOPS OF -y

TO TIP 2 S.NCHES 3 INCHES 4 INCHES

0.0. I.-e. 0.B. X.0. O.B, 1.B. O.B. 1.0,

3 3 0.1 24 33 3 2

4 11 I.0 27 41 20 15

5 24 3.3 31 46 73 57 68 54 59 47 34 2s

6 42 8.2 33 50 199 159 188 151 178 143 145 115

7 59 15.8 36 53 408 329 388 314 376 305 349 282

t 68 65

10 49

11 29

12 13

23.7 38 5 5 643 526 614 505 604 497 579 475

28.7 40 57 RIO 670 777 646 769 639 747 621

26.7 42 59 783 653 753 631 747 626 734 615

19.1 44 60 573 481 552 465 548 462 541 456

10.2 46 61 311 263 300 254 299 253 296 251

13 4

14 1

--"-

368

3.7

I.1

.--0-0

141.7

48

49

62 114 97 110 94 110 94 109 93

63 34 2\9 33 28 32 28 32 28

-m--- ---(I- -11-.- ----* r-*-o --I)-- ----- o-"-w

3971 3281 3783 3142 3722 3094 3566 2961

***

C-S SdRVXUAL NClDEL

***

ARITH. LlEAb DBH= Cl.1 QUADRATIC MEAN DBH= 8.4 ,

PERCEkT SURVIVAL= 50.7

WEIBULL PARAMETERS A= 1.5 , f3= 7.39 , C= 3.62

MEAN CROWN RATIO= 38.6

CORRESPONhING SITE INDEX FOR BASE AGE 2;=

LAMBDA= 0.597

5 8

***
GROWING

SEASONS s1 EHS

SINCE DBH

EST. CLASS

25

3

4

5

6

7

(3

9

10

11

12

13

14

PER

ACRE

3

11

24

42

s9

68

65

49

29

13

4

1

I---

368

LOBLOLLY PINE

SITE INDEX 85 (BASE AGE 50)

726 STEMS PER ACRE PLANTEU

C-S SURVIVAL MODEL

---------.,-VOLUMES PER ACHE-----------

FUEL PIJLP~OOD CYIP*W'SAW SAvlTIMBER

CHIPS BOLTS LOGS LOGS

(TONS) (CORDS) (CORDS) CMBF)

0.0

0.2

0.3

0.7

1.2

1.7

2.0

1.8

1.3

0.7

0.2

0.1

e-w"

10.2

0.0 0.0 0.000 0.17 0.00 0.00 0.00 0.17

0.i 0.0 0.000 0.83 1.06 0.00 0.00 1.89

0.5 0.0 0.000 1.69 9.03 0.00 0.00 10.73

1.7 0.0 0.000 3.40 29.28 0.00 0.00 32.68

3.8 0.0 0.000 5.93 64.02 0.00 0.00 69.95

6,l 0.0 0.000 8.40 104.51 0.00 0.00 112.91

5.5 2.1 0.000 10.06 93.48 94.54 0.00 198.07

3.4 3.9 0.000 9.06 57.16 175.46 0.00 241.68

1.5 3.7 0.000 6.36 25.98 168.14 0.00 200.40

4.0 0.8 0.813 3.56 68.24 37,30 121.96 231.06

1.6 0.2 0.339 1.24 26.95 10.08 so.92 89.19

0.0 0.0 0.158 0.37 0.60 2.06 23.65 26.68

28.3

w-w-1

1.310

0.-0-c-r..a--STURPAGE VALUE

FUEL PULIwOOD CHIP'N'SAW

CHIPS

($1

BOLTS

(5)

--

LOGS

(5)

-ww----

51.08 480.31 487.58 196.53 1215.50

SAWTIWBER DEH CLASS

LOGS TOTAL

(S) ($1

m-111111
YIELDS YY OIAMETER CLASSES FUR UNTHINNED LOALOLLY PINE IN PLANTATIONS

ON CUT-OVFH (UN"HEPAWED BUT PROBLEM FREE) SITES

IN THE

WEST GULF COASTAL PLAIN

--------------------------

*** .MERCHANTABILITY PARAMETERS ***

--------------------------

SAkJTI~fiFH VDLrJMFS ARE CALClJLATED FHDM STUMP HEIGHT CO.5 FEET) TO A 15.0 INCH TOP DIAMETER (INSIDE BARK).

LUGS ARE BUCKED INTO 16 FOOT LENGTHS, PLUS A 3 INCH ALLO*ANCE. LOG SCALING DIAMETER TS THE DIAMETER INSIDE

HARK AT THE SMALL END OF EACH LOG. IF THE LAST LGG BFFCJHE REACHXNG THE MINZMUM SCALING DIAMETER EXCEEDS 8 FEET,

BUT IS LESS THPQ 16 FEET, IT IS BUCKED TO THE NEXT SMALLER EVEQ-NUMBERED FOOT OF LENGTH, PLUS TRIM ALLOYANCE.

THE PEHATr\‘LQC LENGTH IS OHUPPED TkTD TFiF NEXT LOWFH STIIMPACE VALUE CATEGORY. IF THE LAST LOG TO THE MINIMUM

SCALlNG DIAMETER ZS LFSS THAN 8 FEET LO&G, THE ENTIRE LENGTH IS DROPPED INTO THE NEXT LOWER STUMPAGE VALUE

CATEGOkY.

STUWPAGE DIAHETER

VALUE LIMITS

CLASS CiH)

#EIGHT

UNIT

FACTORS

STUMPAGE

VALUE

PER UNIT

FUELrlOUD CHIPS LESS THAN 3.0 INCHES 90. CUBIC FEET S 5.00 PER TON

MO BKBNCHES PER TON

PULPWOOD BOLTS 3.0 - 7.5 26. CUBIC FEET

PER TON

8 6.25 PER TON

CHIY'N'SAW LOGS 7.5 -15.0 30. CURLC FEET

PER TON

s 16.50 PER TON

SAwTJMHUR LOGS GREATER THAN

15.0 INCHES

30. CUBIC FEET

PER TON

S 31.00 PER TON

*EIGHT UNIT FACTORS FCR PRODUCTS OTHER THAN FUELNDDD CHIPS ARE IN CUBIC FEET OF GREEN YOOD

PER TON. TOTAL WEIGHT OF THE GREEN 1000 OF THESE PQODUCTS IS INCREASED RY 10 PERCENT

TO ACCOUNT F‘f)R HAWK HEIGHT.
LORLOLLY PlNE

SITE INDEX 85 (RASE ACE 50)

726 STEMS PER ACRE PLANTED

C-S SlJRVIVAlr MODEL

GROWING

SEASONS AV.

SIWE D+C

EST. HT.

***

20 52

DBH

2

3

4

5

6

7

R

9

10

11

12

ARITH. MEAN

***

STEMS

PER

ACRE

BASAJr

AREA

CR

1 0.0 25

6 0.3 30

19 1.7 34

41 5.6 37

hR 13.4 40

8EI 23.5

86 30.0

64 ad.3

34 18.5

12 7.9

3 2.4

--(I- -w---w

422 131.5

42

44

4 6

48

50

51

-,,~,,,,,,~,,~,-,,-,,,--coBIcF’JOT VOLUNE------------------------

ALL TREtS S-INCH CLASS AND GREATER, STUMP HEIGHT .5 FEET,

GRIJUkD FOR O.R. TOPS OF --

AV. TT) TIP 2 INCHES 3 INCHES 4 INCHES

HT. O.B. I.H. 0.B. I.B. 0.8. I.B. O.B. 1.6.

23 0 0

33 5 4

39 33 25

44 119 91 110 85 96 74 59 43

47 308 242 289 229 274 218 228 178

49 571 459 542 437 526 425 480 394

41 764 623 728 597 716 587 687 562

5 3 751 619 718 595 711 589 691 572

54 so4 419 483 403 480 400 471 393

55 220 184 211 177 210 176 207 174

56 63 52 61 50 60 50 6 0 4 9

----I) B---L -1111-e 1-w-- --m-. -*me- --*-m --w-w

3338 2718 3142 2573 3073 2519 2891 2365

***

DBH= 7.4 , QUADRATIC YE4N DBH= 7.6 , lJEXtjULL PARAMETERS A= 1.2 , B= 6.78 , C- 3.81

PFiRCEnST SURVIVAL= 58.1 , MEAN CROWN RATIO= 42.8 , I>AUBDA= 0.506

CORRRSQONDXNG SITE INDEX FOR BASE AGE 25= 58

***
GROWING

SEASONS STEMS

PER

ACHE

SINCE DISH

EST. CLASS

20

2

3

4

5

6

7

8

9

1 (I

1 1

12

F !j E L PlJLf’dWD CHIP’N’SAW

CHIPS Wlbl’S LOGS

(TUNS) (TONS) (TONS)

1 0.0

6 0.1

19 0.3

41 9.6

68 1.1

88 1.7

86 2.0

64 1.3

34 1.2

12 9.5

3 9.1

-w-e

422

v-m9.5

0.0

0.0

9.4

2.8

8.9

17.7

24.6

16.6

7.5

1.9

0.4

w-w--

89.9

9.0

9.9

0.0

0 . 0

0.0

0.0

0.0

6.8

8.1

4.7

1.5

-w--e

21.1

SAuiTIi4BER

LOGS

( T 0 N s 1

0.090

0.000

rl.000

0.000

0.009

0.000

0.000

0.000

0.000

0.000

0.000

0.000

FUEL

CHIPS

CS)

0.00 0.00 0.00

0.28 0.00 0.00

1.37 2.23 0.00

2.R9 17.70 0.00

5.48 55.37 0.00

8.42 110.93 0.00

10.15 153.71 0.09

9.53 104.05 112.45

5.98 46.67 132.96

2.50 12.12 78.15

0.74 2.56 24.15

47.34

PULPvlOOD

BOLTS

(S)

--m----

505.35

CHIP’N’SAW

LOGS

($1

w-a-1-w

347.72

SAYTIWRER

LOGS

($1

0.00

0.00

0.09

0.00

0.09

0.00

0.00

0.00

0.00

0.00

0.00

-I)-----

0.00

DBH CLASS

TOTAL

($1

0.00

0.28

3.59

20.59

60.85

119.35

163.86

226.04

185.62

92.77

27.46

--e--w900.41
Guldin, Richard W. MERCHANDIZER: Adapting

USLYCOWG to the marketplace. U.S. Department of

Agriculture Forest Service, General Technical Report

SO-52. South For. Exp. Stn. New Orleans, La. 36 p.

1984.

USLYCOWG is a computer model that estimates

yields in cubic feet per acre from unthinned slash and

loblolly pine plantations established on cutover sites

in the West Gulf states (see Research Papers SO-147

and SO-148). A series of computer routines called

MERCHANDIZER has been developed that convert

USLYCOWG yields in cubic feet to tons of fuel chips,

cords of pulpwood and chip’n’saw logs, and board feet

of sawtimber, or weights of all products in tons.

Keywords: Economics, stumpage, forest products,

volume predictions, volume yield, Pinus taeda, Pinus

elliottii var. elliottii, unthinned plantation yields, Weibull distribution, taper curves, crown ratio, survival

prediction, problem-free sites.