diff --git a/R/FuruVol.R b/R/FuruVol.R index b6adcdc..9b29fab 100644 --- a/R/FuruVol.R +++ b/R/FuruVol.R @@ -16,17 +16,17 @@ function (dbh, trh, bark, enhet) f.dub <- round(dub, 1) } if (bark == "mb") { - i <- f.dmb < 11.1 + i <- f.dmb < 12.1 vol[i] <- 0.6716 + (0.075708 * (dmb[i]^2)) + (0.029679 * (dmb[i]^2) * h[i]) + (0.004341 * dmb[i] * (h[i]^2)) - i <- f.dmb > 11 + i <- f.dmb > 12 vol[i] <- -6.3954 + (0.178053 * (dmb[i]^2)) + (0.03317 * (dmb[i]^2) * h[i]) - (0.003008 * (dmb[i]^2) * b[i]) } else if (bark == "ub") { - i <- f.dub < 11.1 + i <- f.dub < 12.1 vol[i] <- 2.3393 + (0.010045 * (dub[i]^2)) + (0.038834 * (dub[i]^2) * h[i]) + (0.002732 * dub[i] * (h[i]^2)) - i <- f.dub > 11 + i <- f.dub > 12 vol[i] <- -3.5425 + (0.128182 * (dub[i]^2)) + (0.028268 * (dub[i]^2) * h[i]) + (0.008216 * dub[i] * (h[i]^2)) }