# load test data
data(sp2)
# promote to SPC
depths(sp2) <- id ~ top + bottom
# assign abrupt textural change to a profile
newdata <- data.frame(id = c("hon-21"),
restrkind = "abrupt textural change",
restrdep = 46)
# do left join
restrictions(sp2) <- newdata
# inspect site table: newvalue TRUE only for horizons
# with top depth equal to zero
restrictions(sp2)
Run the code above in your browser using DataLab