Question: how to plot z=f(r,phi) instead of z=f(x,y)?

so far i plotted functions where z values were functions z=f(x,y) in cartesian coordinates. this was easily done with the plot3d command. now i have a function described in polar coordinates, z=f(r,phi). so far i haven't succeded in plotting it. as far as i see it cylindrical coordinates would be required here, but maple expects r=f(z,angle) and i need z=f(r,angle). maybe a changecoords operation would help, but I'm not sure how to perfrom that (I've just started using maple). would be thankful for pointers on how to tackle my problem.

Here's the function I want to plot (in the preview of my post it looked kinda scrambled)

intensity := proc (rho, `ϕ`) options operator, arrow; Intensity*(rho/omega)^l*[L(m, l, 2*rho^2/omega^2)]^2*exp(-2*rho/omega^2)*cos(l*`ϕ`)^2 end proc