ParametricPlot with bar legend

In this file, I have plotted region bounded by concentric circle of radius ‘rs’, which is a function of ‘rh’. Now I want to show the corresponding values of T in bar legend. The bar legend part is ok but the plot does not show the corresponding color for T values. for same rh, the color corresponding to T should show in the plot. Can anyone fix it? I have attached the notebook.

temp profile.nb (85.7 KB)

Use ColorFunction instead of PlotStyle.

I haven’t tried to reverse engineer your code to make sure this produces your desired result, and some of your code seems suspicious (like the FindMinimum/FindMaximum stuff), but ColorFunction is the option to use to color the plot according to the coordinates of points on the plot.