# NOT RUN {
top <- 2
bottom <- 6
wtr <- c(25.2,25.1,24.1,22.0,19.8,15.3,12.0,11.1)
depths <- c(0,1,2,3,4,5,6,7)
bthA <- c(10000,8900,5000,3500,2000,1000,300,10)
bthD <- c(0,1,2,3,4,5,6,7)
layer.density(top,bottom,wtr,depths,bthA,bthD)
# }
Run the code above in your browser using DataLab