# NOT RUN {
Area <- 1E4 ## 1 ha
HSIZE <- 30 ## 0--30 cm
ORCDRC <- 50 ## 5%
ORCDRC.sd <- 10 ## +/-1%
BLD <- 1500 ## 1.5 tonnes per cubic meter
BLD.sd <- 100 ## +/-0.1 tonnes per cubic meter
CRFVOL <- 10 ## 10%
CRFVOL.sd <- 5 ## +/-5%
x <- OCSKGM(ORCDRC, BLD, CRFVOL, HSIZE, ORCDRC.sd, BLD.sd, CRFVOL.sd)
x ## 20.25 +/-4.41 kg/m^2
## in tonnes per ha:
x[[1]] * Area / 1000
# }
Run the code above in your browser using DataLab