I am using the default of Plots, which is Gr. That command does not work there. I found this workaround:
my_cg = cgrad([:green,:green])
p_surf = surface(x,y,z,c=my_cg,colorbar=false)
Which produces a surface of uniform color in an annoying fashion by having to specify a custom colorbar.
Kind regards, thanks for your suggestion