Leem is not the same as silt

[JobdePater]

For the soil workability A_SILT_MI is used. Is that correct? According to Locher and Backer leem = lutum + silt

https://github.com/AgroCares/Open-Bodem-Index-Calculator/blob/7b799211c348c3cd15e7544cec5cc36c627fc171/R/workability.R#L92C35-L92C35

[BrentWHR]

Do I understand correctly that you are wondering whether:

    # sandy soils with variable silt content    
    dt[soiltype.m == 'zand' & A_SILT_MI < 10, c(cols) := list(45,35)]
    dt[soiltype.m == 'zand' & A_SILT_MI >= 10 & A_SILT_MI < 20, c(cols) := list(55,30)]
    dt[soiltype.m == 'zand' & A_SILT_MI >= 20, c(cols) := list(60,30)]

should be changed to:

    # sandy soils with variable silt content    
    dt[soiltype.m == 'zand' & A_SILT_MI + A_CLAY_MI < 10, c(cols) := list(45,35)]
    dt[soiltype.m == 'zand' & A_SILT_MI + A_CLAY_MI >= 10 & A_SILT_MI+ A_CLAY_MI < 20, c(cols) := list(55,30)]
    dt[soiltype.m == 'zand' & A_SILT_MI + A_CLAY_MI >= 20, c(cols) := list(60,30)]

?

[JobdePater]

yes indeed. That will result in higher depths.

you forgot 1 in the second line:

dt[soiltype.m == 'zand' & A_SILT_MI + A_CLAY_MI >= 10 & A_SILT_MI + A_CLAY_MI < 20, c(cols) := list(55,30)]