### Janne Rämö ### # Bollandsås et al. (2008) kasvumalli # param maxt := 750; #Time horizon param s := 3; #Harvesting interval, 1=5y 2=10y etc. param r := 0.03; #Interest rate param SI := 15; #Site index param LAT := 61.9; #Latitude param w := 50; #Width of the size classes param sp := 1; #Species, 1=Spruce 2=Birch 3=Pine param sc := 12; #Number of size classes set SIs := 6 11 15 17; ### For all matrices 1=Spruce, 2=Birch, 3=Pine, 4=Others ### ### Puutavarahinnat, 1=kuusi, 2=koivu, 3=mänty, aikavälin 2000-2011 kantohinnat deflatoitu tukkuhintaindeksillä vuoteen 2011 ### param ps:= #Tukin hinta 1 55.463 2 48.551 3 55.897 4 0 ; param pp:= #Kuidun hinta 1 23.708 2 15.58 3 16.105 4 0 ; ### Price multiplier for sensitivity analysis. ### param pkerroin:= 1 1 #Kuusen kerroin 2 1 #Koivun kerroin 3 1 #Männyn kerroin 4 1 ; ### Growth multiplier for sensitivity analysis. ### param Ikerroin:= 1 1 #Kuusen kerroin 2 1 #Koivun kerroin 3 1 #Männyn kerroin 4 1 ; ### Ingrowth multiplier for sensitivity analysis. ### param dkerroin:= 1 1 #Kuusen kerroin 2 1 #Koivun kerroin 3 1 #Männyn kerroin 4 1 ; ## Recruit-function parameters ## # Probability # param rp: 1 2 3 4 := 1 -2.291 -0.904 -3.552 -3.438 #Intercept 2 -0.018 -0.037 -0.062 -0.029 #BA 3 0.066 0 0 0.123 #Site index 4 0.019 0.016 0.031 0.031 #PBA ; # Conditional # param rc: 1 2 3 4 := 1 43.142 64.943 67.152 31.438 # Intercept 2 -0.157 -0.161 -0.076 -0.1695 # Basal area 3 0.368 0.143 0 0.442 # Site Index 4 0.051 0.104 0 0.193 # Percentage of basal area of subject species (spruce) ; ## Diameter increment -function parameters ## param di: 1 2 3 4 := 1 17.839 11.808 25.543 2.204 #Intercept 2 0.0476 0 0.0251 0.063 #Diameter at breast height, DBH 3 -11.585 9.616 -5.660 -8.320 #DBH^2, *10^-5 4 0 -9.585 0 0 #DBH^3, *10^-8 5 -0.3412 0 -0.216 0 #Basal area or larger trees, BAL 6 0.906 0.519 0.698 0.359 #Site index 7 -0.024 -0.152 -0.123 -0.177 #Basal area 8 -0.268 -0.161 -0.336 0 #Latitude ; ## Mortality-function parameters ## param mp: 1 2 3 4:= 1 -2.492 -2.188 -1.808 -1.551 #Intercept 2 -0.020 -0.016 -0.027 -0.011 #DBH 3 3.200 2.700 3.300 1.400 #DBH^2 *10^-5 4 0.031 0.030 0.055 0.016 #Basal area ; #################################################################### ## Pulp volumes for different site indices, 1=Spruce, 2=Birch 3=Pine ## param Vp := [6,*,*]: 1 2 3 4:= 1 0.01214 0.01264 0.011693 0 2 0.05485 0.05384 0.055173 0 3 0.13494 0.12626 0.130175 0 4 0.06132 0.05713 0.042260 0 5 0.05286 0.04826 0.028843 0 6 0.04370 0.04047 0.025223 0 7 0.03921 0.03315 0.021293 0 8 0.03420 0.02912 0.019748 0 9 0.03224 0.02576 0.015598 0 10 0.03001 0.02644 0.012912 0 11 0.02723 0.02044 0.010455 0 12 0.02617 0.01808 0.008229 0 13 0 0 0 0 14 0 0 0 0 [11,*,*]: 1 2 3 4:= 1 0.01285 0.01445 0.033415 0 2 0.06061 0.06552 0.063695 0 3 0.15062 0.15522 0.096854 0 4 0.06857 0.07000 0.037383 0 5 0.06052 0.05743 0.029172 0 6 0.04872 0.04731 0.026897 0 7 0.04593 0.04769 0.025973 0 8 0.04370 0.04179 0.025512 0 9 0.03787 0.03290 0.025248 0 10 0.03573 0.03096 0.025085 0 11 0.03329 0.03058 0.024976 0 12 0.03035 0.02474 0.024900 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 [15,*,*]: 1 2 3 4:= 1 0.01374 0.01591 0.101943 0 2 0.06664 0.07464 0.128403 0 3 0.16690 0.18005 0.133533 0 4 0.08080 0.07854 0.111023 0 5 0.06482 0.06655 0.084151 0 6 0.05975 0.05827 0.064902 0 7 0.04978 0.04978 0.052703 0 8 0.05039 0.04865 0.044991 0 9 0.04324 0.04463 0.039959 0 10 0.03925 0.03891 0.036549 0 11 0.03317 0.03685 0.034153 0 12 0.03073 0.03268 0.032414 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 [17,*,*]: 1 2 3 4:= 1 0.01374 0.01591 0.101943 0 2 0.06664 0.07464 0.128403 0 3 0.16690 0.18005 0.133533 0 4 0.08080 0.07854 0.111023 0 5 0.06482 0.06655 0.084151 0 6 0.05975 0.05827 0.064902 0 7 0.04978 0.04978 0.052703 0 8 0.05039 0.04865 0.044991 0 9 0.04324 0.04463 0.039959 0 10 0.03925 0.03891 0.036549 0 11 0.03317 0.03685 0.034153 0 12 0.03073 0.03268 0.032414 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 ; ############################# ## Saw log volumes for different site indices, 1=Spruce, 2=Birch 3=Pine ## param Vs := [6,*,*]: 1 2 3 4:= 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 0.19151 0.17208 0.200336 0 5 0.35321 0.31225 0.354334 0 6 0.54620 0.47653 0.540876 0 7 0.76099 0.65933 0.759964 0 8 0.99508 0.85015 1.011596 0 9 1.24309 1.05669 1.295774 0 10 1.49430 1.26561 1.612496 0 11 1.75811 1.48208 1.961764 0 12 2.01921 1.69041 2.363576 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 [11,*,*]: 1 2 3 4:= 1 0 0 0 0 2 0 0 0 0 3 0 0 0.092443 0 4 0.21435 0.21483 0.256160 0 5 0.39553 0.39299 0.460942 0 6 0.61681 0.59908 0.709789 0 7 0.85638 0.82020 0.993701 0 8 1.11749 1.06492 1.321678 0 9 1.40218 1.32770 1.690720 0 10 1.68841 1.59244 2.100827 0 11 1.97974 1.85821 2.551200 0 12 2.28072 2.12589 3.044236 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 [15,*,*]: 1 2 3 4:= 1 0 0 0 0 2 0 0 0 0 3 0 0 0.097638 0 4 0.23419 0.25137 0.270343 0 5 0.44578 0.45137 0.485147 0 6 0.68392 0.69732 0.742052 0 7 0.96304 0.96304 1.041056 0 8 1.25313 1.24859 1.382161 0 9 1.57421 1.55035 1.765365 0 10 1.89981 1.86531 2.290670 0 11 2.21442 2.18117 2.658074 0 12 2.56544 2.49693 3.167579 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 [17,*,*]: 1 2 3 4:= 1 0 0 0 0 2 0 0 0 0 3 0 0 0.097638 0 4 0.23419 0.25137 0.270343 0 5 0.44578 0.45137 0.485147 0 6 0.68392 0.69732 0.742052 0 7 0.96304 0.96304 1.041056 0 8 1.25313 1.24859 1.382161 0 9 1.57421 1.55035 1.765365 0 10 1.89981 1.86531 2.290670 0 11 2.21442 2.18117 2.658074 0 12 2.56544 2.49693 3.167579 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 ; ## Initial stand ## param y0: 1 2 3 4 := #1=spruce, 2=birch, 3=pine, 4=others 1 2000 1750 1750 0 2 0 0 0 0 3 0 0 0 0 4 0 0 0 0 5 0 0 0 0 6 0 0 0 0 7 0 0 0 0 8 0 0 0 0 9 0 0 0 0 10 0 0 0 0 11 0 0 0 0 12 0 0 0 0 13 0 0 0 0 14 0 0 0 0 15 0 0 0 0 ; ## Diameters at breast height ## param DBH := 1 75 2 125 3 175 4 225 5 275 6 325 7 375 8 425 9 475 10 525 11 575 12 625 13 675 14 725 15 775 ;