I don’t see color interpolation in this code although i am passing color as an array of vertices number size.
using Meshes, CoordRefSystems ,GLMakie
r = Float32[1.6, 1.754088, 1.9230156, 2.1082115, 2.311243, 2.5338273, 2.7778475, 3.0453684, 3.3386526, 3.6601818, 4.012676, 4.3991165, 4.8227735, 5.287231, 5.7964177, 6.354642, 6.9666257, 7.637547, 8.373081]
θ = Float32[0.0, 0.09817477, 0.19634955, 0.2945243, 0.3926991]
ϕ = Float32[0.0, 0.3926991, 0.7853982, 1.1780972, 1.5707964, 1.9634954, 2.3561945, 2.7488935, 3.1415927, 3.5342917, 3.9269907, 4.3196898, 4.712389, 5.105088, 5.497787, 5.8904862, 6.2831855]
fig = Figure()
ax::LScene = LScene(fig[1, 1], show_axis=false)
g = RectilinearGrid{𝔼,typeof(Spherical(0,0,0))}(r, θ, ϕ)
el = element(g,1)
clr = [-129.9932861328125, -123.21377563476562, -140.03717041015625, -145.09832763671875, -113.99639892578125, -110.56198120117188, -125.68373107910156, -126.87786865234375]
viz(el, color =clr, alpha=0.5)