[FR] zcolor specifictation surfaces

[rrhiemstra]

zcolor specification seems to not be implemented, although I have the recollection that it used to be in Plots around 2016.

using Plots
plotly() # or pyplot() - gr() does not work

X(r,theta,phi) = r * sin(theta) * sin(phi)
Y(r,theta,phi) = r * sin(theta) * cos(phi)
Z(r,theta,phi) = r * cos(theta)

thetas = range(0, stop=pi,   length=50)
phis   = range(0, stop=pi/2, length=50)

xs = [X(1, theta, phi) for theta in thetas, phi in phis] 
ys = [Y(1, theta, phi) for theta in thetas, phi in phis]
zs = [Z(1, theta, phi) for theta in thetas, phi in phis]
zc = sqrt.(xs.^2 .+ ys.^2 .+ zs.^2) 

surface(xs, ys, zs, zcolor=zc,
    seriestype=:surface, 
    aspect_ratio=:equal)

[BoundaryValueProblems]

How about using fill_z instead of zcolor?

[BeastyBlacksmith]

zcolor is an alias for marker_z and thus only works for points