Bulletin 574 November 1985 b~b~g g +Pautig erdaureu program to produce stand maps Alabama Agricultural Experiment Station, Auburn University David H. Teem, Acting Director, Auburn University. Alabama CONTENT S Page INTRODUCTION............................................................ PROGRAM APPLICATION ............................................... 3 4 A. Data File ...................................................... 5 B. Plotting Parameters .......................................... 5 C. Other Plotting Parameters................................. 6 EXAMPLE PLOTS ......................................................... SUMMARY..................................................................10 7 APPENDIX A. Computer Plotting: Basic Concepts............ 10 APPENDIX B. Program Listing.................................... 12 FIRST PRINTING 3M, NOVEMBER 1985 Information contained herein is available to all persons without regard to race, color, sex, or national origin. A GENERALIZED PLOTTING PROGRAM TO PRODUCE STAND MAPS J. S. Conway, B. L. Lanford, B. J. Stokes, and W. D. Greene 1 INTRODUCTION ALL THE TOOLS used by foresters, maps are among the most indispensable. One type of map which is becoming more popular is a stand map. Stand maps show the location of each tree on a plot of ground. Often tree diameter is drawn to scale and a tree identifying number may also be included. Stand maps have several uses: (1) documenting tree locations within a permanent research plot, (2)illustrating spatial relationships of trees within a stand, (3)illustrating thinning treatments, (4) moulding input for interactive simulation programs, or (5) adding input for harvesting or growth and yield simulation. Permanent field plots are often maintained by industry, government, or academic organizations for forest inventory and research measurements. Locating a tree which was previously measured is made easier with a stand map. Such a map should show tree location and any identifying number or species code. Forest stands generally occur in one of three spatial arrangements, clumped, random or uniform. Clumped stands are typically found where seed tree or shelterwood methods of regeneration are used. Random stands arise from seed or sprout regeneration while uniform arrangement results in planted stands such as pine plantations. Stand maps quickly illustrate the differences in spatial arrangement present in various stands. Stand maps can also illustrate the effect on a forest stand of different thinning methods (selective, third-row removal, 'The authors are former Graduate Research Assistant of Forestry, Associate Professor of Forestry, Forest Engineer, Southern Forest Experiment Station, USDA Forest Service, and Graduate Research Assistant of Forestry. OF 4 ALABAMA AGRICULTURAL EXPERIMENT STATION etc.). Maps can also be used with interactive simulation programs which are becoming increasingly popular for analyzing thinning treatments. In addition, most numerical harvesting simulation packages rely on stand information which is easily illustrated with a stand map. Many growth and yield simulation programs utilize similar stand information, which couldalso be represented with a stand map. While stand maps have many uses, they are time consuming to create by hand. To reduce the amount of time needed to produce a stand map, a general stand map plotting program was created. PROGRAM APPLICATION Most mainframe computer systems support some type of generalized plotting package. Often the package consists of a set of subroutines which perform a specific type of plotting function (e.g., draw a line, draw a symbol, etc.). These subroutines are called within a computer program in order to create a stand map or other complex graphic image. A more detailed discussion of plotting programs and their logic is found in Appendix A. The plotting program presented in this publication is named PLOTSTND (for "Plot-stand"). It was written in FORTRANG for use on the IBM 3033 mainframe computer available at Auburn University. The plotter supported by this system is the Versatec-07 electrostatic plotter. While the program presented here is specific to this computer and plotter combination, the logic and structure of the PLOTSTND program should be applicable to any plotting system. Appendix B contains a program flowchart and code. PLOTSTND is organized so that the user can quickly draw a plot of a study area which has been recorded in a data file. The format of data in the file will vary depending on the application, but the x- and y-coordinates and diameter at breast height (DBH) must be included. Tree number, species, and other variables also may be used. Before running the program, the user must specify the name of the data file to be used and set various plotting parameters within the program to create the desired plot. The program is then processed by the computer in a batch fashion. A discussion of each of the required inputs to the program follows. A PROGRAM TO PRODUCE STAND MAPS 5 Data File The data file is specified within the job control language (JCL) statements of the program. The data file is specified on the following JCL statement which is found near the end of the program: //GO.FT10F001 DD DSN=jobcode. liwords "jobThe brary.type(member),DISP =SHR. code.library.type(member)" are replaced by the name of the data file to be used. If the data are stored in a sequential file, the "(member)" portion is omitted. Plotting Parameters There are nine plotting parameters which should be specified by the user before using the program. The parameter values are changed within the program itself. Comment statements early in the program code explain the parameters, Appendix B. A brief description of each parameter follows. SCLMAP defines the scale of ground distances on the plot. This number must be in map inches per actual foot of ground distance. Using a value of SCLMAP which is a fraction of 10 (e.g., 0.10, 0.05, etc.) allows convenient use of an engineer's scale for measuring distance. DSCALE specifies the scale of tree diameters in map inches per actual DBH inch. It is often desirable to exaggerate tree diameter slightly in order to distinguish relative tree sizes. A DSCALE value of two times SCLMAP results in a six-fold exaggeration of DBH and is often useful. HTNUM specifies the character height in inches for printing tree numbers and species codes beside each tree location. Setting HTNUM = 0.05 provides clear but unobtrusive labelling of each tree symbol. TREEAC identifies the desired residual stand density for the plot. PLOTSTND will draw a grid overlay with TREEAC scaled blocks per acre on the plot. This feature provides a quick method for checking the number of trees remaining in a thinned stand. A stand which has been thinned to TREEAC trees per acre should have an average of one tree within each grid block. Specifying TREEAC = 0.0 will cancel the griddrawing procedure for a plot. THNDBH is used when thinning is performed within the plotting program. THNDBH specifies the minimum diameter of "leave" trees on the plot (i.e., all trees smaller than THNDBH will not be plotted). If thinning is not to occur within the 6 ALABAMA AGRICULTURAL EXPERIMENT STATION program, set THNDBH=0.0. This parameter must be used in conjunction with the ITHNOP parameter discussed below. XWIDTH is the actual length in feet of the study plot. In order to efficiently use the Versatec's plotting space, XWIDTH should equal the longest plot side. (If the Y dimension is the longest dimension on the ground, the axes will be transformed within the program to generate an accurate plot.) YWIDTH is the actual length in feet of the shorter side of the study plot. Because of limits on the paper sizes used with the Versatec plotter, the product of YWIDTH and SCLMAP should be less than 10 inches. ISYMB is an integer parameter which specifies the information to be plotted within the plot. If ISYMB = 1, only tree symbols will be drawn at each tree location. If ISYMB 2, the computer will plot tree symbols with tree number and species code numbers at each tree location. Tree number and species often clutter the appearance of a plot, therefore if they are not needed, ISYMB should be set equal to 1. ITHNOP defines the type of data file being input into the program. For files which have been thinned by another program and have "cut" and "leave" codes for each tree, set ITHNOP = 1. For "raw" data files which are to be thinned by this program, set ITHNOP = 2. When ITHNOP = 2 and THNDBH = 0.0, a plot of all trees before thinning is drawn. Other Plotting Parameters Two plotting parameters are found within the JCL statements rather than within the program code. Both of these parameters are found within the following JCL statement: &PLOT SCALE= 1.00 *YMAX=10.55 &END. SCALE determines the global scale of the plot. For example, when SCALE = 2.0, the map size will be doubled, and when SCALE = 0.5, the map size will be halved. The default value for SCALE is 1.0. The value of YWIDTH X SCLMAP X SCALE should be less than 10.0 inches for single strip maps. When SCLMAP is not an even factor of 10 (0.333 for example), the axes notation on the plot will take on non-integer values and be rounded before plotting. By specifying SCLMAP = 0.5, and SCALE = 0.667, a map with a scale of 1.0 inch equals 30 feet will be produced with complete axes. The user should note that when the global scale is changed, the scale label on the plot will be in error by a factor of scale times. A PROGRAM TO PRODUCE STAND MAPS 7 *YMAX defines the maximum width of the plot. As mentioned above, YWIDTH X SCLMAP must be less than or equal to 10.0 inches. This allows about half an inch along the map's edge for axis labels. Since the Versatec's paper is 10.55 inches wide, the default value for *YMAX is 10.55. If a map of greater width is desired, the parameter *YMAX must be increased in increments of 10.55 inches. Plotter output for larger map sizes will consist of two or more 10.55-inch wide strips which must be assembled to form a complete map. The value of YWIDTH X SCLMAP X SCALE should always be less than *YMAX- 0.50. EXAMPLE PLOTS To illustrate the types of plots available and some of their possible uses, four examples are discussed. Each example plot is shown along with the plotting parameters needed to produce the plot with PLOTSTND. This should aid beginning users in becoming familiar with the program and its possible uses. Figure 1 contains a stand map of a typical unthinned southern pine stand. Tree diameter is slightly exaggerated (3 times) so that a viewer can distinguish relative tree size. There are 594 trees per acre and row alignment is fairly well defined. Such a map would show the type of thinning to be performed or to illustrate how a stand appeared before thinning. The previous stand is shown after thinning in figure 2. There are now 350 trees per acre and row alignment is more difficult to discern due to the reduced number of trees. Comparison of the two stand maps provides an effective illustration of the changes in a stand due to thinning. Stand maps can also be used to represent a forest stand in an interactive simulation program, figure 3. The scale of tree diameters and ground distances are equivalent, allowing a scale model machine moved through the stand to accurately simulate spatial relationships present. This use of stand maps is relatively new, but has tremendous potential for use in analyzing logging problems. A more traditional use of stand maps is for documenting tree locations within a permanent (or long-term) study plot. A stand map for such a purpose would show trees in their respective locations with their identifying numbers alongside, figure 4. This reduces the difficulty of finding the tree on return trips to the plot for data collection. If space permits on the stand map itself, it can provide an effective tally sheet for field work. 8 ALABAMA AGRICULTURAL EXPERIMENT STATION X-DIRECTION 0.00 50.00 100.00 150.00 CD O 0 0 o O° o °° 0 0 o O° o °° o p ° 0 0. O° N " ° o 000 s® oop o o ° °° .0Os g o 0 0 0 o 0 0° o o , N 0 0 0° o o°° e 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 O o Oo op o o Oo 00 0 o Oo o"p 0 DSCLE O = HTNUM 000 .0 (inches per foot) inch) O o 0 ° 00 p o p ° ° 0 0 0 0 0 CD O (inches per O o ° o ° o °° (inches) 0 0 o O° o °° o O o O° o o° o p' o' z (trees per acre) TIONPH =.0 XWIDTH =l19. YWIDTH (inches) o 0 Oo , 0 0 ° O o 0 0 0 p f (feet) (feet) 1fet 1220. C o o O O O O 0 0 0 0 0 O 0 0 0 CD o Oo 00 o F- C-) CD O O o o Oo 00 0 0 0 0° O O o o O o 0o p o p o ITFINOP JCLL raetrs =2 BCD oo O o o 0 0 0 ° ° 0 0 0 ° ° O O o O G o 0 0 o O° o 0 0 O I 1. b Oo 000 OQ, ® o 000 Qp , 8 0 0 0 r SCALE = .0 1(1imes OCLIOPF' °'f o OP ° ° ° o ° 0 0 O 0 0 o f,°° °" 0 0 0 o *YO00 =1'05 ul o 0 o ° ff 0 0 O as p o Oo o"O ° ° o0 o 0 o ° Oo 0 0° O ° oo oo Oo 00 Oo O ° 0 u O ° o o0 , O0 ° 0 o ° 0° o ° ooO o e ° o O° o ° 0 0 ° o f O O O" 8 ° o ° oo O, 0 0 000 "o 0 , o 0 0 L0 0 40.00 50.00 100.00 150.00 X-DIRECTION FIG. 1. Stand map of a typical unthinned stand. X-DIPECTION 0°!00 O ° o 0 50.00 0° 0° 0 0 100.00 0 0 o 150.00 0 O O O N ° o° ° 0 0 0 0 ° ° o o ° o o O O C\i 0 0 0 ° p Plottinog Parameters~: " o O o " 0 o O o 0 0 o O p p o o 0 ° 0 ° v o ° O O o SCLMAP =002 DSCIOLE = 0.005 (inchs (iches per per inch) acre) z {fot) O O O ° o O O ° ° O o O O HTN011=00.4 TREEAC = 0.0 THNDBH = 0.0 XWIDTH =190. YWIDTH =220c. =1 000MB ITHNOP =2 JCL Parameters: SCALE *YOB (inches) (trees per (iche) (feet) Ifeet) O Ln0 - 0 0 ° ° 0 0 0 ° 0 O ° o O O o p O O o ° 0 p 0 o 0° op Z 0° ° ° 0 0 0 ° o ° 0 0 O° a O p ° 0 0 ° H O O O O O p p O O O O p p O H- U w O o O " O O o o O " O O " O CD U O LL! O Lr ° O v O 0 0 0 ° o O 0 ° o O 0 ° ° 0 0 =I1.00 = 10. 55 times (inches) OCLOOP) O O o 0 0 0 0 0 0 ° o ° " ° 0 0 0 0 0 0 0 0 0 0 0 o 0 O O o O "p v o0 O ov O o 0 ° o O 0 ° 0 O 0 ° ° 6 O 9 ° " O O ° p p O O ° " p p O ° o 0 0 0° p 0 0 O o 0 p o 0 0 0 0 O O O o° ° 0 00 0 0 ° 0 0 0° ° O °0.00 50.00 100.00 150.00 X-DIRECTION FIG. 2. Stand map of a typical thinned stand. A PROGRAM TO PRODUCE STAND MAPS X-DIRECTION 50.00 100130 150.00 9 0 00 C Plotting Prmeters: CD SCLMAF'= 0.02 (inches pea- foot) 0 LSCLLL= O. c:)6 a(nchesper (inachesm) HTNIUMI=a.04 ace =0.0 (tre pTF'EC T1-f'DBH =. (inches ) XWIDTH =191. (feet) inch) ja 0 in 0 O CD YWIDTH 220. oU (feet) o W orF .-a H ISYMB =1 ITHNOP_2 JCLParametersSCALE0=i.00aa(times 000550) *YMAO = 10.55(inches) ° 00 5000 100.00 150.00 X-D IRECT ION FIG. 3. Stand map used in interactive machine simulation. X-DIRECTION 0.00 50.00 100.00 100.00 a, 5 4212546 41 CDat 4410540;a1134it 6 41\ a J00st 55 Plotting POarametana: 05 O (anches (anches per pea a 05 m a 09 IyI 174 3 45 a15 4 1727 26 254 VS 6 SCLMAF = 0.0a2 0.a05 DSCOLE = foaot) inch) Om 7084 " 0 1 2 q s 1 HTNUM =a0.04 =a0a0 TREFOAC THNDBH =a0aa XWIDTH YWIDTH= 000MB=2 ITHN(IF =2 JC10CPaametera- =148. 22aa. =10.55 (inches) (trcam perracrae) (anches) (fect) (feet) 12e a 96 Cm 0 61 992Cm6267 02a062 Li 94 673 i'+ ass 610e26 619 64 C 3 5 60 7C4m 51 65 264 673CD L 348 73 a9e9 as em 455 e654 e803 em0 61 62 6564 083C3m' 1 SCALEt=1.00aa ims SCLMAF) *YMfO (inches QC m 811 a12e844 i65 ssa-86 082 Cm Z CmCZ) . a l 41 6 D6 a a4 a n 090732 ro o47am em a e m668 a2 em am830 52em m7 a67em1 66a ae W53 C LO 49 a4 a s Cm 45 479 m76aa4 3 53 a710 a5% a155aa assa5516a,283 am&5 a1a6a5C64e00a2a13 am55 Cm 537 a5 415 17 iaa a317 3565em 37 *6669 a °0.00 00.00 100.00 150.00 C X-DIRECTION FIG. 4. Stand map used for documentation of a study plot. 10 ALABAMA AGRICULTURAL EXPERIMENT STATION SUMMARY Maps of forest stands have a variety of uses. Computer generated stand maps can be used to keep records of how a plot's composition changes over time. Once a map has been constructed, it can be quickly updated. The stand map also can be used for data input in conjunction with a digitizer. More accurate distance measurements are also possible with computer-drawn maps. Maps can often save expensive trips to the field. Applications have also been developed recently for using stand maps as an integral part of machine modelling with interactive simulation. With the availability of sophisticated graphics software on today's microcomputers, the formation of stand or plot maps is a procedure that most managers or researchers can use to make their work more efficient. APPENDIX A Computer Plotting: Basic Concepts This discussion of plotting concepts is included to familiarize computer users with basic computer graphics. This is not a "how-to" plotting guide; rather it is an overview of the plotting process. Those users who may wish to modify the PLOTSTND program or write their own plotting program may want to review these general ideas. The Graphics Package Most computer graphics packages use a planar Cartesian coordinate system (perpendicular Xand Y-direction) as the operating scheme for plotting. The software will define the logical origin, axes convention, and plotting boundaries for ordinary plot generation. These plotting boundaries define what is referred to as the plotting "window." The user may specify plotting parameters and plotting window sizes other than the system-provided defaults when a special need arises. Extra JCL statements are required to change these default values. To draw diagrams using the computer graphics system, a program of plotter pen instructions is required. These instructions specify the X and Y coordinates to which the pen is to move from its present position. Whether the pen is to be up or down during the move must also be specified. The A PROGRAM TO PRODUCE STAND MAPS 11 A PROGRAM TO PRODUCE STAND MAPS 1 programmer must also instruct the pen to either return to the plotting origin or to remain at the present location after each move is made. The program of pen instructions must conform to the requirements of the plotting software. The programmer must make sure that the pen stays within the plotting window during all moves. Otherwise, lines may be unplottable or partially clipped if they fall outside the plotting window. When using the IBM-Versatec plotting system or any computer combination, system subroutines are available which greatly reduce the amount of programming required. These subroutines are sets of plotter pen directions which perform the most commonly used plotting functions, such as drawing axes, plotting graphs and bar charts, and printing labels and scale values. In the case of the IBM-Versatec plotting system, plotting subroutines may be accessed by entering the FORTRAN CALL statement: CALL subroutine name (argument). In the position labeled "subroutine name," the name of the subroutine being used is entered. Subroutines used in the PLOTSTND program include AXIS and SYMBOL. The "argument" of the subroutine is the set of parameters which dimension the plotting function. For example, in the argument for the CALL AXIS subroutine, the user must specify the axis starting coordinates, the axis label, the number of characters in the axis label, the length of the axis, the angle at which the axis line is to be drawn, the first value used in axis annotation, and the number of data units per map inch. A complete description of each plotting subroutine used in PLOTSTND may be found in the user guides for AUPLOT and VTEC PLOTTER available from the Division of University Computing at Auburn University. 12 ALABAMA AGRICULTURAL EXPERIMENT STATION APPENDIX B Listing of PLOTSTND.FORT //USERID# JOB (USERID,BIN),'USER NAME',MSGCLASS=P, NOTIFY= USERID# /*ROUTE PRINT LOCAL /*JOBPARM TIME= 05,LINES= 20K //STEP EXEC FTGVCLG //FORT.SYSIN DD * C C THIS IS A PROGRAM TO PLOT TREE LOCATIONS AND DIAMETERS ON THE AUCS VERSATEC PLOTTER. C C *PLOTTING PARAMETER DEFINITIONS: C C (SCLMAP) THE MAP SCALE IS ENTERED AS A FRACTION (INCHES PER FOOT). C (DSCALE) DIAMETER SCALE IS ENTERED AS A FRACTION (INCHES PER INCH). C (HTNUM) HEIGHT OF PRINTED TREE & SPECIES CODE NUMBERS (INCHES). C (TREEAC) TREES PER ACRE FOR THE GRID DENSITY. C **NOTE: 1. SET (TREEAC)= 0.0 TO DELETE GRID PLOTTING OPTION. C (THNDBH) MINIMUM DIAMETER OF "LEAVE" TREES. C (XWIDTH) THE X-DIRECTION DIMENSION OF ACTUAL STAND (FEET). C (YWIDTH) THE Y-DIRECTION DIMENSION OF ACTUAL STAND (FEET). C **NOTE: 1. (XWIDTH SHOULD BE THE LONGER OF THE TWO DIMENSIONS. C 2. (YWIDTH*SCLMAP) MUST BE .LE. 10.0. C (ISYMB) =1 FOR PLOTTING TREE SYMBOLS ONLY. C =2 FOR PLOTTING TREE SYMBOLS WITH NUMBERS. C (ITHNOP)= 1 FOR SOURCE FILES ALREADY THINNED. C =2 FOR SOURCE FILES NOT THINNED, THINNING IN THIS PROGRAM. C C C SCLMAP=0.10 DSCALE= 0.00833 HTNUM=0.05 TREEAC = 0.0 THNDBH=0.0 XWIDTH =66. YWIDTH = 66. ISYMB= 1 ITHNOP=2 C A PROGRAM TO PRODUCE STAND MAPS 13 A PROGRAM TO PRODUCE STAND MAPS C *COUNT # OF TREES IN DATA SET. C ITTLTR = 0 DO 10 1=1,100000 READ (10,100,END=20) ITNUM,DBHT,ISP,YLOC,XLOC,TLHT,LMDIAFLAG ITTLTR =ITTLTR + 1 10 CONTINUE 20 CONTINUE REWIND 10 WRITE (6,101) ITTLTR 101 FORMAT(1X,'NUMBER OF TREES IN FILE=',I6) C A SCALED AXES SYSTEM IN THE CORRECT DIRECTION ON THE VERSATEC'S PLOTTING PAPER. C C *THE FOLLOWING STATEMENTS ARE MARGINS, DO NOT CHANGE. C 1 *DRAW XSTART = 0.5 YSTART=0.5 C C C C C XLEN =XWIDTH*SCLMAP YLEN = YWIDTH*SCLMAP XAX = XLEN + XSTART YAX =YLEN + YSTART *INITIALIZE PLOTTING. CALL PLOTS (0,0,0) DVAL= 1/SCLMAP DBHV= 1/DSCALE XX=XLEN+0.6 C C C CALL PLOT (XSTART,YSTART, -3) C C C CALL AXIS (0.0, 0.0,'X-DIRECTION', 11,YLEN,90.0,0.0,DVAL) CALL AXIS (XLEN, 0.0,'X-DIRECTION',1 1,YLEN,90.0,0.0,DVAL) CALL AXIS (0.0, 0.0,'Y-DIRECTION',-1l1, * XLEN, 0.0,0.0,DVAL) CALL AXIS (0.0, YLEN,'Y-DIRECTLON', 11,XLEN, 0.0,0.0,DVAL) CALL SYMBOL (XX,0.0,0.10, * *THIS 'CALL PLOT' COMMAND PUTS THE LOGICAL CENTER ON THE MARGIN. *THE 'CALL AXIS' COMMAND DRAWS THE AXES TO SCALE. 'SCALES; DIST.: 1 INCH= FEET; DBH: 1 14 14 ALABAMA AGRICULTURAL EXPERIMENT STATION INCH= * INCHES', 90.0,55 ) CALL NUMBER ( XX,2.3,0.10,DVAL,90.0, 0) CALL NUMBER ( XX,4.5,0.10,DBHV,90.0, 0) C IF ( TREEAC .EQ. 0.0) GO TO 80 GSIZE = SQRT(43560/TREEAC) IXGRID=(XWIDTH/GSIZE)+ 1 IYGRID=(YWIDTH/GSIZE)+ 1 GRDSZE = GSIZE*SCLMAP YGCNTR=(GRDSZE/2(-1 ))+YSTART C C C CALL GRID (0.0,YGCNTR,IXGRID,GRDSZE,IYGRID,GRDSZE, -30584) C 80 CONTINUE C DO 200 I=1,ITTLTR READ (10,100) ITNUM,DBHT,ISP,YLOC,XLOC,TLHT,LMDIA,FLAG 100 FORMAT(13,F4.1,13,F4.0,F4.0,F4.0,F4.1,F4.0) C C *TRANSFORMATIONS FOR PLOTTING ARE WRITTEN HERE. C X = YLOC *THE 'CALL GRID' COMMAND DRAWS A TREE PER ACRE GRID OVERLAY. 90 95 110 111 C IF (XLOC .GT. 0) GO TO 90 Y= ABS(XLOC)+ 50.0 GO TO 95 Y=50.0 -XLOC CONTINUE TRNUM = ITNUM IF ( X.GT. XWIDTH) GO TO 200 IF ( Y.GT. YWIDTH) GO TO 200 IF ( ITHNOP .EQ. 1) GO TO 111 IF ( DBHT .GT. THNDBH) GO TO 110 FLAG=1 GO TO 111 CONTINUE FLAG=0 CONTINUE XDIST=X *SCLMAP YDIST =Y *SCLMAP DSIZE= DBHT*DSCALE C C C C *THE 'CALL NUMBER' COMMAND DRAWS THE TREE NUMBER AND SPECIES CODE AT THE TREE LOCATION. A PROGRAM TO PRODUCE STAND MAPS 15 A PROGRAM TO PRODUCE STAND MAPS IF ( ISYMB .EQ. 1) GO TO 195 CALL NUMBER (XDIST,YDIST,HTNUMTRNUM,0.0,-1) SPCODE= ISP C YDSTN =YDIST -(1.25*HTNUM) C CALL NUMBER (XDIST,YDSTN,HTNUMSPCODE,0.0,-1) C C *THE 'CALL SYMBOL' COMMAND PLOTS THE TREE LOCATIONS ON THE PLOTTING PAPER. THE SIZE OF EACH SYMBOL IS PROPORTIONAL TO THE TREE DBH. C **NOTE: 1. TREES WITH (FLAG=0) ARE 'LEAVE' TREES. C 195 CONTINUE IF (FLAG .EQ. 0) GO TO 199 IF (ITHNOP .EQ. 1) GO TO 199 198 CONTINUE CALL SYMBOL (XDIST,YDIST,DSIZE,3,O.,-1) GO TO 200 C 199 CONTINUE CALL SYMBOL (XDIST,YDIST,DSIZE,1,45.0,-i) C 200 CONTINUE C C END ALL PLOTTING C CALL PLOT (0,0.,999) C C STOP END //GO.PLOTPARM DD* &PLOT 1 SCALE= 1.00 *YMAX =10.55 &END //GO.FT1I OFOO 1 DD DSN =FYL56DG.STAND. MAPS 170PL9), DISP =SHR //GO.SYSIN DD DUMMY /7 Alabama's Agriculural xpcriment SLtmion SN stIIn AUBURN UNIVERSITY W\ ith an agricultural research unit in every major soil area, Auburn Unixersity serves the needs of field crop, livestock, forestry , and hor- 0 ticultural producers in each region in Alabama. Every citi zen of the State has a stake in this research progranm since anyx1a advantage from newx and more economical wxay s of produc- 0 o S O ing and handling farm procucts directly benefits the consuming pubiic. Research I 1h 1s 1 tiiiIdeItifi caitiin ® Main Agricultural Experiment Station, Auburn. ir E. V. Smith Research Center, Shorter. Tennessee Valley Substation, Belle Mina. Ssnd Mountain Substation, Crossville. North Alabama Horticulture Substation, Culiman. Upper Coastal Plain Substation, Winfield. Forestry Unit, Fayette County. Chilton Area Horticulture Substation, Clanton. Forestry Unit, Coosa County. Piedmont Substation, Camp Hill. Plant Breeding Unit. Tallassee Forestry Unit, Autauga County. Prattville Experiment Field, Prattville. Black Belt Substation, Marion Junction. The Turnipseed-Ikenberry Place, Union Springs. Lower Coastal Plain Substation, Camden. Forestry Unit, Barbour County. Monroeville Experiment Field, Monroeville. Wiregrass Substation, Headland. Brewton Experiment Field, Brewton. Solon Dixon Forestry Education Center, Covington and Escambia counties. 20. Ornamental Horticulture Substation, Spring Hill. 21 Gulf Coast Substation, Fairhope. 1. 2. 3. 4. 5. 6. 7. 8. 9 10 11 12 13 14. 15. 16. 17. 18. 19.