Mesh3d plot displays wrong element coordinate - Julia

I try to plot mesh3d to display heatmap data on polar coordinate. The plot skips the first row , e.g. it display from N30 to N90 , while it should display plot from N0 to N90.

r = 1.0:0.5:3.0
θ = 0.0:30.0:90
R = r'.*ones(size(θ))
T = ones(size(r))'.*θ  
 
FL = R.*sin.(deg2rad.(T));

fGrid(ai,aj,r,θ) = [length(θ)*(j+aj-1)+i+ai for i in 1:(length(θ)-2) for j in 1:(length(r)-1) ]
Ra = reshape(R, length(R))
Ta = reshape(T, length(T))
X = Ra.*sin.(deg2rad.(Ta))
Y = Ra.*cos.(deg2rad.(Ta))
Z = reshape(FL, length(FL))
I = [fGrid(0,0,r,θ); fGrid(1,1,r,θ)]   
J = [fGrid(1,0,r,θ); fGrid(0,1,r,θ)]    
K = [fGrid(0,1,r,θ); fGrid(1,0,r,θ)]   

NoAxes=attr(
        showgrid=false,
        zeroline=false,
        showline=false,
        showticklabels=false,
    )

lay = Layout(title="Hole Quality Plot", scene=attr(
        camera=attr(
            up=attr(x=1, y=0, z=1),
            eye=attr(x=0, y=0, z=1.5)
            ),
            aspectratio=attr(x=1, y=1, z=0.001),
            aspectmode="manual",
            xaxis_title="Easting",
            yaxis_title="Northing",
            xaxis =attr(showticklabels=false, showline=false, showgrid=false, zeroline=false),
            yaxis =attr(showticklabels=false, showline=false, showgrid=false, zeroline=false),
            zaxis =attr(showticklabels=false, showline=false, showgrid=false, zeroline=false),
       ),  

)

#t1=end
p1 = mesh3d(
    x=X,
    y=Y,
    z=Z,
    colorbar_title="FL",
    colorscale=[[0.0,"rgb(200,0,50)"],[0.5,"rgb(255,255,255)"],[1.0,"rgb(0,150,50)"]],
    #cmin = -100, cmax = 100, 
    intensity=reshape(FL, length(FL)), 
    i=I,
    j=J,
    k=K,
    name="y",showscale=true
)

tmp=0:0.01π:2π
Line = [scatter(x=minimum(r).*sin.(tmp), y=minimum(r).*cos.(tmp), z=0.0*(tmp), type="scatter3d", mode = "lines", name ="line", line_color="black", showlegend = false )]
tmp = [scatter(x=maximum(r).*sin.(tmp), y=maximum(r).*cos.(tmp), z=0.0*(tmp), type="scatter3d", mode = "lines", name ="line", line_color="black", showlegend = false )]
Line = [Line; tmp]

for t in 0:30:maximum(θ)
    tmp = scatter(x=[minimum(r), maximum(r)].*sin(deg2rad(t)), y=[minimum(r), maximum(r)].*cos(deg2rad(t)), z=[5000.0, 5000.0], type="scatter3d", mode = "lines", name ="line", line_color="#C6C7C1", showlegend = false )
    Line = [Line; tmp]
end


 
p = [p1; Line] 

plot(p,lay)

Capture

My work around is to add dummy elements by

replace

θ =  0.0:30.0:90

with

θ = [-2; 0.0:30.0:90]

Capture2