Reconcile the code with the method from [Vejmelka et al 2016](https://doi.org/10.1071/WF14085) and assess its error