I want to use `GLMakie`

for 3D plotting, but when I try to do the plotting I get a static image not in a Makie window:

https://i.imgur.com/R44xxJa.png

Earlier this morning it was working just fine, but I have since then upgraded to version 1.9. Why does the plotting seem to not use the GLMakie backend?

**MWE**

```
using GLMakie, DifferentialEquations, LinearAlgebra, SparseArrays
GLMakie.activate!()
Nā=31 # Number of waveguides / size of solution vector
Ī³=1 # Nonlinear term strength parameter
h=1 # Grid spacing
centerGrid = (Nā-1)/2;
x = -centerGrid:centerGrid;
# Coefficient matrix of second-order centered-difference operator (Ī“Ā²u)ā
M = spdiagm(-1 => fill(1,Nā-1), 0 => fill(-2,Nā), 1 => fill(1,Nā-1))
M[Nā,1] = 1; # Periodic boundary conditions
M[1,Nā] = 1;
# RHS of DNLS. The solution vector u is a Nāx1 complex vector
gā(u,p,t) = 1*im*(-1*M*u - @.(Ī³*((abs(u))^2).*u) )
# Julia is explicitly typed (e.g, cannot have Ixnt and Complex in same array) and so we must convert the object containing the initial data to be complex
u0 = Complex.(sech.(x))
tspan = (0.0,200)
prob = ODEProblem(gā,u0,tspan, [h])
sol = solve(prob, Tsit5(), reltol=1e-8, abstol=1e-8)
z= [abs(sol.u[i][j])^2 for j=1:Nā, i=1:size(sol)[2]] # |u|Ā²
y = sol.t
fig = Figure()
p1 = surface(fig[1,1], x, y, z,axis=(type=Axis3, xlabel="Space", ylabel="Time", zlabel="|u|Ā²"))
#xticks!(fig, -Nā:1:Nā)
p2 = contourf(fig[1,2], x, y, z, axis=(xlabel="Space", ylabel="Time"))
fig
#save("dnls_makie.png", fig)
#plt = plot(p1,p2,layout=(1,2), size=(1200,800))
```

julia> Makie.current_backend()

GLMakie