add area2radius.stack.R and implement in plots.r