Thank you very much for these explanations :)
Dsolve returns a rule. See this article for information on how to use rules in this situation:
http://support.wolfram.com/kb/12505
Plot[Evaluate[Abs[T[r]]] /. First[sol], {r, rmin, rmax}]
