Context. In the extremely cold and dark environments of pre-stellar cores, methanol is formed on the surface of interstellar dust grains and released into the gas phase via non-thermal desorption mechanisms. Gaseous methanol constitutes the starting point for the formation of many massive complex organic molecules and is therefore of utmost importance for the build-up of chemical complexity. Aims. We aim to improve upon a previous model for the prediction of column densities and deuterium fractions of non-deuterated and singly deuterated methanol. Thereby, we try to identify crucial chemical and physical parameters for which the study of deuteration could provide valuable additional constraints. Methods. We employed a gas-grain chemical code to devise a model that is in agreement with the observed column density and deuterium fraction profiles of the innermost region of the pre-stellar core L1544. For that purpose, we developed a new treatment of reactive desorption, deriving an individual reactive desorption efficiency for every product species in a chemical reaction that depends on the reaction enthalpy and type of the underlying surface. Furthermore, we explored several options to promote the diffusion of hydrogen and deuterium atoms over the surface of interstellar dust grains in order to increase methanol formation. Results. Our fiducial model employs diffusion via the quantum tunnelling of hydrogen and deuterium atoms, resulting in CH3OH and CH2DOH column densities that are approximately an order of magnitude lower than the observed values, which is an improvement over the results of the previous model by a factor of 10. The N(CH2DOH)/N(CH3OH) ratio is reproduced within a factor of 1.2 for the centre and 1.8 for the position of the methanol peak. Given the large uncertainties that chemical models typically have, we consider our predictions to be in agreement with the observations. In general, we conclude that a diffusion process with a high diffusion rate needs to be employed to obtain methanol column densities that are in accordance with the observed values. Also, we find that the introduction of abstraction reactions into the methanol formation scheme suppresses deuteration when used in combination with a high diffusion rate.