Assessing uncertainties in the numerical modelling of cosmogenic radionuclide production