Description Usage Arguments Details Value Author(s) References See Also Examples
View source: R/compressive_properties3.R
It calculates the compressive parameters N, lambda and kappa using the pedo-transfer function from de Lima et al. (2018)
1 2 |
compressive_properties3(bulk.density, matric.suction, soil=c("SandyLoam","SandyClayLoam"))
|
bulk.density |
a numeric vector containing the values of bulk density, Mg m^{-3} |
matric.suction |
a numeric vector containing the values of matric suction, hPa |
soil |
the soil texture group 'SandyLoam' or 'SandyClayLoam'. See exemples |
Pedo-transfer function deveploed under no-till condition. See de Lima et al. (2018)
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>
de Lima, R. P., da Silva, A. P., Giarola, N. F., da Silva, A. R., Rolim, M. M., Keller, T., 2018. Impact of initial bulk density and matric suction on compressive properties of two Oxisols under no-till. Soil and Tillage Research, 175: 168-177.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 | # EXAMPLE 1
compressive_properties3(bulk.density=1.5, matric.suction=100, soil="SandyLoam")
compressive_properties3(bulk.density=1.5, matric.suction=100, soil="SandyClayLoam")
# EXAMPLE 2 for SandyLoam soil
matric.suction <- seq(from=30,to=1000,len=100)
out <- compressive_properties3(bulk.density=1.5,
matric.suction=matric.suction, soil="SandyLoam")
plot(x=matric.suction,y=out$N, ylab="N",
xlab="Matric suction (hPa)", log="x") # plot for N
# plot for lambda
plot(x=matric.suction,y=out$lambda, ylab="lambda",
xlab="Matric suction (hPa)", log="x")
# plot for kappa
plot(x=matric.suction,y=out$k, ylab="kappa",
xlab="Matric suction (hPa)", log="x")
# EXAMPLE 3 for SandyClayLoam soil
matric.suction <- seq(from=30,to=1000,len=100)
out <- compressive_properties3(bulk.density=1.5,
matric.suction=matric.suction,
soil="SandyClayLoam")
plot(x=matric.suction,y=out$N, ylab="N",
xlab="Matric suction (hPa)", log="x") # plot for N
# plot for lambda
plot(x=matric.suction,y=out$lambda,
ylab="lambda", xlab="Matric suction (hPa)", log="x")
# plot for kappa
plot(x=matric.suction,y=out$k,
ylab="kappa", xlab="Matric suction (hPa)", log="x")
# End (not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.