vol eq 2

Volume_Norway_underbark.R

@@ -10,41 +10,37 @@ function (dbh.mm, height.dm, tree.sp)
         }
         x
     }
-    
+    if (any(is.na(height.dm))) stop("some heights are NA")
     if (any(dbh.mm < 50)) 
         stop("all dbh.mm should be at least 50mm")
     
-    i.tree.sp.1 <- num.species %in% c("Norway spruce")
-    i.tree.sp.2 <- num.species %in% c("Scots Pine")
-    i.tree.sp.3 <- num.species %in% c("Other Conif")
-    i.tree.sp.4 <- num.species %in% c("Hardwoods")
-    if ((sum(i.tree.sp.1) + sum(i.tree.sp.2) + sum(i.tree.sp.3) + 
-        sum(i.tree.sp.4)) != length(num.species)) {
-        print(num.species[which(!i.tree.sp.1 & !i.tree.sp.2 & 
+    i.tree.sp.1 <- tree.sp %in% c("Norway spruce")
+    i.tree.sp.3 <- tree.sp %in% c("Scots pine", "Other Conif")
+    i.tree.sp.4 <- tree.sp %in% c("Hardwoods")
+    if ((sum(i.tree.sp.1) +  sum(i.tree.sp.3) + 
+        sum(i.tree.sp.4)) != length(tree.sp)) {
+        print(tree.sp[which(!i.tree.sp.1 &
             !i.tree.sp.3 & !i.tree.sp.4)])
         stop("some species are not recognized")
     }
     vol.wo.bark.l <- rep(NA, length(dbh.mm))
 
     ## under bark
-       
-    vol.wo.bark.l[ i.tree.sp.1] <- picea.abies.vol(dbh = dbh.mm[i12], 
-            trh = height.dm[i12], bark = "ub", units = "l") * 
-            (1 - (nvl(vol.reduksjon[i12], 0)/100))
-    vol.wo.bark.l[i.tree.sp.2] <- sitka.vol(dbh = dbh.mm[i.tree.sp.2], 
-            trh = height.dm[i.tree.sp.2], bark = "ub", units = "l") * 
-            (1 - (nvl(vol.reduksjon[i.tree.sp.2], 0)/100))
-       
-    vol.wo.bark.l[i.tree.sp.3] <- pinus.sylvestris.vol(dbh = dbh.mm[i32], 
-            trh = height.dm[i32], bark = "ub", units = "l") * 
-            (1 - (nvl(vol.reduksjon[i32], 0)/100))
-    vol.wo.bark.l[i.tree.sp.4] <- harw.nor.vol(tsl = as.integer(tree.sp[i.tree.sp.4]), 
-            dbh = dbh.mm[i.tree.sp.4], trh = height.dm[i.tree.sp.4], 
-            bark = "ub", units = "l") * (1 - (nvl(vol.reduksjon[i.tree.sp.4], 
-            0)/100))
-    
-    volume.norway <- list()
+
+    vol.wo.bark.l[ i.tree.sp.1] <- picea.abies.vol(dbh = dbh.mm[i.tree.sp.1], 
+            trh = height.dm[i.tree.sp.1], bark = "ub", units = "l")     
+    vol.wo.bark.l[i.tree.sp.3] <-
+      pinus.sylvestris.vol(dbh = dbh.mm[i.tree.sp.3], 
+                           trh = height.dm[i.tree.sp.3],
+                           bark = "ub", units = "l")
+    ## assuming all hardwoods are birch
+    vol.wo.bark.l[i.tree.sp.4] <- harw.nor.vol(tsl = rep(30, sum(i.tree.sp.4)) ,
+                                               dbh = dbh.mm[i.tree.sp.4],
+                                               trh = height.dm[i.tree.sp.4], 
+                                               bark = "ub", units = "l") 
+
     
     vol.wo.tr.m3 <- vol.wo.bark.l/1000
+  
     return(vol.wo.tr.m3)
 }