I know it is over two years later, but in case someone needs it, here is a solution I got to work with Julia 1.10
using Plots;
plotlyjs();
u1 = LinRange(0, 2 * pi, 100) .* ones(100)'
v1 = (LinRange(0, 2 * pi, 100) .* ones(100)')'
u2 = sin.(u1) .* sin.(v1)
nl, nc = size(u1)
mesh3d(vec(u1), vec(v1), vec(u2), intensity = vec(u2), color=:buda, showscale=true)
Where the colorscale is pulled from: https://docs.juliaplots.org/latest/generated/colorschemes/
I get the folowing display:
