add populus as the 8th species

Author: BayesianLearner

R/parameters.r

@@ -1,6 +1,6 @@
 ###default parameters
-speciesNam <- c("pisy","piab","beal","fasy","pipi","eugl","rops")
-pCROB <- matrix(NA,43,7,dimnames = list(NULL,speciesNam))
+speciesNam <- c("pisy","piab","beal","fasy","pipi","eugl","rops",'popu')
+pCROB <- matrix(NA,43,8,dimnames = list(NULL,speciesNam))
 rownames(pCROB) <-  c("cR","rhow","sla","k","vf","vr","c","mf","mr","mw",
                       "z","beta0","betab","betas","rhof2","s1","kRein","s0scale","x","aETS",
                       "alfar1","alfar2","alfar3","alfar4","alfar5","sarShp","S_branchMod","conifers=1 or decidous=2","p0_ref","ETS_ref",
@@ -53,6 +53,15 @@ pCROB[,7] <- c(0.24661449,257.6439127,45.75561945,0.287913095,0.902917082,
                0.331516085,0.566954087,0.694243495,0.768673015,0.843102534,
                1,1,2,1.5,2855,0.047822508,21.98403388,0.061856656,0.000105076,
                -3.324,0.135,-0.03,0.011098704,45.75561945,0,0.958581963,0,1.1)
+pCROB[,8] <- c(0.0552574805543716,146.1815811,59.47661807,0.293587959,0.989521945,
+                  0.766409722,0.207791519,0.228413058,0.281946456,0.021923979,
+                  1.971003606,0.299262023,0.348366895,0.445408413,157.6778994,
+                  0.016181779,1315.101353,0.318597647,0.944661382,
+                  0.003504272,0.297756481,0.380447231,0.428789134,
+                  0.468428346,0.597870642,1,1,2,1.5,2855,0.048016582,
+                  24.97203048,0.038961807,0.000412673,-3.324,0.135,-0.03,
+                  0.031652341,59.47661807,0,0.818727858,0,1.1) #Populus is based on a small dataset
+
 pPREL = c(413.000000, 0.450000, 0.118000, 3.000000, 0.745700, 10.930000, -3.063000, 17.720000,
           -0.102700, 0.036730, 0.777900, 0.500000, -0.364000, 0.271500, 0.835100, 0.073480,
           0.999600, 0.442800, 1.200000, 0.330000, 4.970496, 0.000000, 0.000000, 160.000000,
@@ -67,7 +76,7 @@ pYAS = c(4.897147e-01, 4.913873e+00, 2.419735e-01, 9.487642e-02, 4.362893e-01, 2
          -7.913602e-05, 3.518549e-02, -2.089906e-04, -1.808920e+00, -1.172547e+00, -1.253595e+01,
          4.596472e-03, 1.302583e-03, -4.389227e-01, 1.267467e+00, 2.569142e-01)
 
-parsAWEN <- matrix(NA,12,7,dimnames = list(NULL,speciesNam))
+parsAWEN <- matrix(NA,12,8,dimnames = list(NULL,speciesNam))
 parsAWEN[,1] <- c(0.518000,0.177300,0.088700,0.216000,0.474660,0.019012,0.078308,
                   0.430248,0.670000,0.022500,0.007500,0.285000)
 parsAWEN[,2] <- c(0.482600,0.131700,0.065800,0.319900,0.474660,0.019012,
@@ -85,8 +94,10 @@ parsAWEN[,6] <- c(0.360000,0.281200,0.098800,0.260000,
                   0.750000,0.030000,0.000000,0.220000)
 parsAWEN[,7] <- c(0.407900,0.198000,0.099000,0.295100,0.474660,0.019012,
                   0.078308,0.430248,0.715000,0.015000,0.000000,0.275000)
+parsAWEN[,8] <- c(0.407900,0.198000,0.099000,0.295100,0.474660,0.019012,
+                  0.078308,0.430248,0.715000,0.015000,0.000000,0.275000)
 
-pHcM <- matrix(NA,7,7,dimnames = list(NULL,speciesNam))
+pHcM <- matrix(NA,7,8,dimnames = list(NULL,speciesNam))
 ###parameters Hcmodel pisy
 pHcM[,1] <- c(1.4628,-0.256007,0.017293,-0.089891,0.213387,-0.079094,-0.062191)
 ###parameters Hcmodel piab
@@ -101,9 +112,11 @@ pHcM[,5] <- c(1.4628,-0.256007,0.017293,-0.089891,0.213387,-0.079094,-0.062191)
 pHcM[1:5,6] <- c(-1.067271, -0.017684, 0.07708,	-0.619978,0.045046)
 ###parameters Hcmodel rops
 pHcM[1:3,7] <- c(0.04237,-0.13308,0.31382)
+###parameters Hcmodel poplus
+pHcM[1:3,8] <- c(0.04237,-0.13308,0.31382)
 
-litterSizeDef <- matrix(0.,3,7,dimnames = list(NULL,speciesNam))
-litterSizeDef[1,] <- c(30,30,10,30,30,20,20)
+litterSizeDef <- matrix(0.,3,8,dimnames = list(NULL,speciesNam))
+litterSizeDef[1,] <- c(30,30,10,30,30,20,20,20)
 litterSizeDef[2,] <- 2
 
 ClCut_birch <- matrix(NA,2,4)