Scripts underlying the publication: Predicting subsurface classification in 2D from cone penetration test data

This dataset contains the script underlying the manuscript titled 'Predicting subsurface classification in 2D from cone penetration test data'. The script can be used for predicting the most likely 1D classification of the subsurface at CPT locations. The uncertainties in the CPT measurements and transformation to soil units are accounted for by integrating the Robertson chart with the Bayesian approach. The classification is carried out in two steps: (a) for a given number of layers, identifying the most probable soil layer thicknesses, and (b) identifying the most likely number of layers which, together with the most probable thicknesses, maximises the likelihood of observing the CPT measurements. The script can be used with any CPT data in .gef format.