Plots.Surface rescaling gives incorrect color legend

Hello, I am trying to plot a surface via Julia Plots. The maximum of the series I would like to plot along the z axis is 10. But I want the z axis to rescale to (0,20). When I add in the argument: zlim = (0,20), I get what I want but the color scale on the right hand side is incorrect. It should go from 0 to 20 but it still does 0 to 10. Therefore, the colors don’t match with the numbers being graphed. Here is my output:

Is there a way to fix this? Thank you!

There is a related issue reported here for the default Plots.jl gr() backend.

You can use any of the other backends with 3D plotting capability: plotly(), plotlyjs(), pyplot(), gaston() or pgfplotsx().

The latter produces:

For such cliff-surfaces, plotlyjs() seems to handle color gradients quite well: