How can I plot the geometry of my problem using pdegplot(g) on top of my three-dimensional solution using surf(qx,qy,qz) in the x-y plane?

In order to solve the Poisson equation in two-dimensions using finite element, I have created the geometry in the pdetoolbox. I can export the decomposed geometry into Matlab and plot it using the pdegplot(g) function. Separately, I have solved the problem using my own finite element code, and now want to plot my solution on top of the geometry layout. In order to plot the solution I will be using surf(qx,qy,qz), where qx & qy are defined as the meshpoint over my domain and qz includes the solution for each point respectively. I will be viewing this solution in the x-y plane (i.e. view(2))

There is no error in using

h=surf(qx,qy,qz)
shading interp
hold on
j=pdegplot(g)
hold off

But I am unable to send the surface plot to the back using:

uistack(h,'back')

Is there another alternative?

Thanks!