-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 



