The procedure for a more accurate calculation of the decay rate predictions within a specific GUT model (replacing the rough estimate in Eq. (REF )) is well-known, cf. e.g. [1]}, and has been implemented in the recently released Mathematica package ProtonDecay [2]} which (amongst its dimension five nucleon decay features) takes the dimension six nucleon decay operators as input and uses them to compute the nucleon decay rates of 13 different decay channels. With these predictions for the decay rates, a specific GUT model can be confronted with the current constraints [3]}, [4]}, [5]}, [6]}, [7]}, [8]}, [9]} and the future experimental sensitivities [10]}, [11]}, [12]}.
