Nothing
# get soil depth based on morphology
estimateSoilDepth <- function(f, name='hzname', top='hzdept', bottom='hzdepb', p='Cr|R|Cd', no.contact.depth=NULL, no.contact.assigned=NULL) {
# extract horizons
h <- horizons(f)
# extract possible contact
contact.idx <- grep(p, h[[name]], ignore.case=TRUE)
# everything else
no.contact.idx <- grep(p, h[[name]], ignore.case=TRUE, invert=TRUE)
# no contact defined, use deepest hz bottom depth
if(length(contact.idx) < 1) {
d <- max(h[[bottom]][no.contact.idx], na.rm=TRUE)
# is there a user-specified depth at which we assume a standard depth?
if(!missing(no.contact.depth) & !missing(no.contact.assigned)) {
if(d >= no.contact.depth & !is.null(no.contact.assigned))
res <- no.contact.assigned
else
res <- d
}
# otherwise use depth of deepest horizon
else
res <- d
}
# contact defined
else
res <- min(h[[top]][contact.idx], na.rm=TRUE)
return(res)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.