User:Mwtoews/K

xlim <- c(3e-12,3) # reference to cm/s

drawblank <- function{ frame plot.window(xlim=xlim, ylim=c(1,2), log="x") # cm/s axis(1, 10^seq(-12,0), line=0)              # cm/s mtext("cm/s", 1, line=0, adj=1, las=1) plot.window(xlim=xlim*864, ylim=c(1,2), log="x") # m/day axis(1, 10^seq(-9,3), line=-3)                  # m/day mtext("m/day", 1, line=-3, adj=1, las=1) plot.window(xlim=xlim*2834.65, ylim=c(1,2), log="x") # ft/day axis(1, 10^seq(-9,3), line=-6)                      # ft/day mtext("ft/day", 1, line=-6, adj=1, las=1) plot.window(xlim=xlim*1164.958, ylim=c(1,2), log="x") # darcy axis(1, 10^seq(-9,3), line=-9)                       # darcy mtext("darcy", 1, line=-9, adj=1, las=1) }

fetter <- data.frame(material=c("clay", "silt", "silty sand", "well-sorted sand", "well-sorted gravel"), lower=c(1e-9, 1e-6, 1e-5, 1e-3, 1e-2), upper=c(1e-6, 1e-4, 1e-3, 1e-1, 1))
 * 1) "Applied Hydrogeology" Fetter, 3rd Ed., Table 4.6, Page 89

pos <- 1:5 drawblank plot.window(xlim=xlim, ylim=c(-3,6), log="x") # cm/s with(fetter, segments(lower, pos, upper, pos)) with(fetter, text(sqrt(upper*lower),pos+.3,material))

domenicoSeds <- data.frame(material=c("gravel","coarse sand","medium sand","fine sand","silt, loess","till","clay","marine clay"), lower=c(3e-4, 9e-7, 9e-7, 2e-7, 1e-9, 1e-12, 1e-11, 8e-13), upper=c(3e-2, 6e-3, 5e-4, 2e-4, 2e-5, 2e-6, 4.7e-9, 2e-9))
 * 1) "Physical and Chemical Hydrogeology" Domenico and Schwartz, 2nd Ed., Table 3.1, Page 38

pos <- rev(1:8) drawblank plot.window(xlim=xlim/100, ylim=c(-5,9), log="x") # m/s with(domenicoSeds, segments(lower, pos,upper, pos)) with(domenicoSeds, text(sqrt(upper*lower),pos+.3,material))