Clear["Global`*"]
Kplus = wp/Sqrt[1 + wp/omega];
Kminus = wp/Sqrt[1 - wp/omega];
Plot3D[{Kplus, -Kplus, Kminus, -Kminus}, {wp, 0, 2}, {omega, -2, 2},
PlotRange -> {-4, 4}, BoxRatios -> {1, 1, 1}, BoundaryStyle -> None,
ClippingStyle -> None, PlotPoints -> 30, AxesLabel -> Automatic]

