### Janne Rämö ### # Bollandsås et al. (2008) growth model, age structured # param maxt := 250; #Time horizon param s := 3; #Harvesting interval param r := 0.03; #Interest rate param SI := 15; #Site index param LAT := 61.9; #Latitude param sp := 1; #Species, 1=Spruce, 2=Birch, 3=Pine param ac := 25; #No. of age classes ### For all matrices 1=Spruce, 2=Birch, 3=Pine, 4=Others ### param ps:= #Saw log price 1 55.463 2 48.551 3 55.897 4 0 ; param pp:= #Pulp price 1 23.708 2 15.58 3 16.105 4 0 ; ## 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 ; ## Initial stand ## param y0: 1 2 3 4 := #1=spruce, 2=birch, 3=pine, 4=others PBA0 addition: Spruce=10 Birch=8 Pine=0.3 Others=5 1 0 0 0 0 2 0 0 0 0 3 25 25 25 25 4 100 100 100 100 5 25 25 25 25 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 16 0 0 0 0 17 0 0 0 0 18 0 0 0 0 19 0 0 0 0 20 0 0 0 0 21 0 0 0 0 22 0 0 0 0 23 0 0 0 0 24 0 0 0 0 25 0 0 0 0 26 0 0 0 0 27 0 0 0 0 28 0 0 0 0 29 0 0 0 0 30 0 0 0 0 31 0 0 0 0 32 0 0 0 0 33 0 0 0 0 34 0 0 0 0 35 0 0 0 0 36 0 0 0 0 37 0 0 0 0 38 0 0 0 0 39 0 0 0 0 40 0 0 0 0 41 0 0 0 0 42 0 0 0 0 43 0 0 0 0 44 0 0 0 0 45 0 0 0 0 46 0 0 0 0 47 0 0 0 0 48 0 0 0 0 49 0 0 0 0 50 0 0 0 0 51 0 0 0 0 52 0 0 0 0 53 0 0 0 0 54 0 0 0 0 55 0 0 0 0 56 0 0 0 0 57 0 0 0 0 58 0 0 0 0 59 0 0 0 0 60 0 0 0 0 ; ## Initial diameters at breast height ## param DBH0 := 1 50 2 60 3 70 4 80 5 90 6 100 7 110 8 120 9 130 10 140 11 150 12 160 13 170 14 180 15 190 16 200 17 210 18 220 19 230 20 240 21 250 22 260 23 270 24 280 25 290 26 300 27 310 28 320 29 330 30 340 31 350 32 360 33 370 34 380 35 390 36 400 37 410 38 420 39 430 40 440 41 450 42 460 43 470 44 480 45 490 46 500 47 510 48 520 49 530 50 540 51 550 52 560 53 570 54 580 55 590 56 600 57 610 58 620 59 630 60 640 ;