Same z-scale on plots3d with layout

Dear all,
I have the following code to plot several curves in a layout form.

using Plots
plotly()

pV = Vector{Any}(undef, 4);
n=[1,2,3,4];
p=1
for i=1:2, j=1:2
    z = n[p]*t
    pV[p] = scatter3d(x,y,z,lw=1,ratio=1,markersize = 0.5,xlabel="x",ylabel="y",zlabel="z",legend=false,zlims=(0,40))
    p=p+1
end
g = plot(pV...;legend=false,layout=(2,2),size=(600,600),zlims=(0,40)) 

and the associated output:

I would link of a same scale in all subplots, e.g, z \in [0,40]. However, when I plot the four curves, the z-scale is different.
How can I do this?