Description Usage Arguments Details Value Author(s) References See Also Examples
View source: R/compressive_properties4.R
It calculates the compressive parameters N, lambda and kappa using the pedo-transfer function from de Lima et al. (2020)
1 2 |
compressive_properties4(matric.suction, soil=c("PloughLayer","PloughPan"))
|
matric.suction |
a numeric vector containing the values of matric suction, hPa. |
soil |
the soil compaction state 'PloughLayer' or 'PloughPan'. See the examples. |
Pedo-transfer function developed for a sandy loam soil texture. 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., Rolim, M. M., da C. Dantas, D., da Silva, A. R., Mendonca, E. A., 2020. Compressive properties and least limiting water range of plough layer and plough pan in sugarcane fields. Soil Use and Management, 00: 1-12.
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 37 | # EXAMPLE 1
compressive_properties4(matric.suction=100, soil="PloughLayer")
compressive_properties4(matric.suction=100, soil="PloughPan")
# EXAMPLE 2 for "PloughLayer"
matric.suction <- seq(from=10,to=10000,len=100)
out <- compressive_properties4(matric.suction=matric.suction, soil="PloughLayer")
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 "PloughPan"
matric.suction <- seq(from=10,to=10000,len=100)
out <- compressive_properties4(matric.suction=matric.suction,
soil="PloughPan")
# plot for N
plot(x=matric.suction,y=out$N,
ylab="N", xlab="Matric suction (hPa)", log="x")
# 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.