In our paper, we propose a robust and operational algorithm to determine the height of the melting level that can be applied to either quasi-vertical profiles (QVPs) or vertical profiles (VPs) of polarimetric radar variables. The algorithm is applied to 1 year of rainfall events that occurred over southeast England and validated using radiosonde data. The algorithm proves to be accurate as the errors (mean absolute error and root mean square error) are close to 200 m.