discretise_rootclasses: Discretise root area ratio over a number of discrete diameter...

In GJMeijer/FBMcw: Generic fibre bundle model for soil reinforcement by plant roots

Description

This function splits a root diameter range into a number of equal-width root diameter classes, and calculates the root area ratio that should be assigned to each class as well as the tensile strength. Power law distributions for root area ratio and tensile strength are assumed. The representative diameter for each class is chosen halfway the range of diameters in the class.

Usage

discretise_rootclasses(drmin, drmax, tru0, betat, phirt, betaphi, nc, dr0 = 1)

Arguments

drmin, drmax	minumum and maximum root diameter in bundle (numeric scalars)
tru0	Tensile strength of root with reference diameter (numeric array)
betat	Power coefficient for tensile strength (scalar)
phirt	Total root area ratio (numeric scalar)
betaphi	Power law coefficient for root diameter-root area ratio fit (numeric array)
nc	The number of discrete, equal-width root classes to use
dr0	Reference diameter (numeric scalar, default 1.0)

Value

a dataframe containing the fields 'dr', the diameter in the middle of each class, 'phir', the root area ratio that should be assigned to each class, and 'tru', the tensile strength assigned to each class. It also return 'drmini' and 'drmaxi', the lower and upper diameter in each class

Examples

discretise_rootclasses(1, 5, 10e3, -0.5, 0.01, -0.5, 10)

GJMeijer/FBMcw documentation built on Dec. 17, 2021, 9:23 p.m.