View source: R/compressive_properties5.R
compressive_properties5 | R Documentation |
It calculates the compressive parameteres N, lambda and kappa using the pedo-transfer function from O'Sullivan et al. (1999)
compressive_properties5(water.content, soil=c("SandyLoam","ClayLoam"))
water.content |
a numeric vector containing the values of gravimetric water content, \% |
soil |
the the soil texture group 'SandyLoam' or 'ClayLoam'. See exemples. |
See O'Sullivan et al. (1999).
N |
the specific volume at p = 1 kPa, N |
CI |
the compression index, lambda |
k |
the recompression index, kappa |
Renato Paiva de Lima <renato_agro_@hotmail.com> Anderson Rodrigo da Silva <anderson.agro@hotmail.com>
O'sullivan, M. F., Henshall, J. K., Dickson, J. W., 1999. A simplified method for estimating soil compaction. Soil and Tillage Research, 49: 325-335.
stressTraffic
# EXAMPLE 1 water.content <- 15 compressive_properties5(water.content=water.content, soil="SandyLoam") compressive_properties5(water.content=water.content, soil="ClayLoam") # EXAMPLE 2 - SandyLoam water.content <- seq(from=5,to=20,len=20) out <- compressive_properties5(water.content=water.content, soil="SandyLoam") plot(x=water.content ,y=out$N, ylab="N", xlab="Bulk density") # plot for N plot(x=water.content ,y=out$lambda, ylab="lambda", xlab="Bulk density") # plot for lambda plot(x=water.content ,y=out$kappa, ylab="kappa", xlab="Bulk density") # plot for kappa # EXAMPLE 3 - ClayLoam water.content <- seq(from=10,to=25,len=20) out <- compressive_properties5(water.content=water.content, soil="ClayLoam") plot(x=water.content ,y=out$N, ylab="N", xlab="Bulk density") # plot for N plot(x=water.content ,y=out$lambda, ylab="lambda", xlab="Bulk density") # plot for lambda plot(x=water.content ,y=out$kappa, ylab="kappa", xlab="Bulk density") # plot for rkappa # End (not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.