dssatR: R Interface to DSSAT-CSM models

get.minBD <- function(sand,clay){
#    minBDtab = dget('~/R/source/rawlsbrakensiek1985fig1.r')
# The following table was taken from figure 1 of Rawls & Brakensiek (1985)
    minBDtab = 
structure(list(V3 = c(1.49, 1.52, 1.52, 1.52, 1.52, 1.52, 1.52, 
1.52, 1.5, 1.48, 1.47, 1.47, 1.45, 1.43, 1.42, 1.4, 1.33, 1.28, 
1.23, 1.2, 1.16), V4 = c(1.4, 1.45, 1.46, 1.46, 1.47, 1.47, 1.47, 
1.46, 1.45, 1.45, 1.44, 1.43, 1.42, 1.41, 1.4, 1.37, 1.33, 1.28, 
1.2, 1.17, -99), V5 = c(1.25, 1.35, 1.4, 1.4, 1.4, 1.4, 1.4, 
1.4, 1.39, 1.39, 1.38, 1.38, 1.38, 1.38, 1.37, 1.34, 1.32, 1.24, 
1.18, -99, -99), V6 = c(1.12, 1.2, 1.26, 1.3, 1.3, 1.3, 1.31, 
1.34, 1.35, 1.37, 1.37, 1.37, 1.37, 1.37, 1.33, 1.32, 1.25, 1.17, 
-99, -99, -99), V7 = c(1, 1.1, 1.18, 1.22, 1.25, 1.27, 1.29, 
1.32, 1.35, 1.37, 1.37, 1.37, 1.37, 1.35, 1.31, 1.26, 1.2, -99, 
-99, -99, -99), V8 = c(0.99, 1.12, 1.2, 1.23, 1.27, 1.3, 1.35, 
1.37, 1.39, 1.4, 1.41, 1.41, 1.4, 1.36, 1.3, 1.25, -99, -99, 
-99, -99, -99), V9 = c(1.08, 1.17, 1.25, 1.3, 1.35, 1.39, 1.42, 
1.43, 1.44, 1.44, 1.44, 1.43, 1.42, 1.41, 1.35, -99, -99, -99, 
-99, -99, -99), V10 = c(1.12, 1.21, 1.29, 1.35, 1.42, 1.45, 1.46, 
1.47, 1.48, 1.48, 1.48, 1.47, 1.45, 1.44, -99, -99, -99, -99, 
-99, -99, -99), V11 = c(1.14, 1.23, 1.31, 1.4, 1.45, 1.48, 1.5, 
1.51, 1.54, 1.54, 1.53, 1.52, 1.5, -99, -99, -99, -99, -99, -99, 
-99, -99), V12 = c(1.18, 1.25, 1.35, 1.44, 1.5, 1.52, 1.54, 1.56, 
1.57, 1.58, 1.58, 1.56, -99, -99, -99, -99, -99, -99, -99, -99, 
-99), V13 = c(1.15, 1.29, 1.4, 1.5, 1.53, 1.55, 1.58, 1.59, 1.61, 
1.61, 1.61, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99), 
    V14 = c(1.27, 1.33, 1.46, 1.52, 1.56, 1.58, 1.61, 1.62, 1.62, 
    1.62, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99
    ), V15 = c(1.3, 1.38, 1.51, 1.55, 1.59, 1.61, 1.62, 1.63, 
    1.63, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, 
    -99), V16 = c(1.35, 1.43, 1.53, 1.6, 1.62, 1.63, 1.65, 1.66, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, 
    -99), V17 = c(1.4, 1.47, 1.57, 1.62, 1.64, 1.66, 1.67, -99, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, 
    -99), V18 = c(1.44, 1.51, 1.6, 1.65, 1.68, 1.69, -99, -99, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, 
    -99), V19 = c(1.4, 1.55, 1.63, 1.7, 1.72, -99, -99, -99, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, 
    -99), V20 = c(1.5, 1.6, 1.64, 1.7, -99, -99, -99, -99, -99, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99
    ), V21 = c(1.52, 1.56, 1.64, -99, -99, -99, -99, -99, -99, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99
    ), V22 = c(1.52, 1.57, -99, -99, -99, -99, -99, -99, -99, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99
    ), V23 = c(1.53, -99, -99, -99, -99, -99, -99, -99, -99, 
    -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99, -99
    )), .Names = c("V3", "V4", "V5", "V6", "V7", "V8", "V9", 
"V10", "V11", "V12", "V13", "V14", "V15", "V16", "V17", "V18", 
"V19", "V20", "V21", "V22", "V23"), class = "data.frame", row.names = c(NA, 
-21L))
    minBD = vector(length=length(clay),mode='numeric')
    for (i in 1:length(minBD)){
        minBD[i]=minBDtab[ceiling(clay[i]/5),ceiling(sand[i]/5)]
    }
    return(minBD)
}