pi=Pi//N;
theta:=ArcTan[Tan[phi] - u/Cos[phi]];
bild=Plot3D[theta,{phi,0,2 pi},{u,-2,2},
PlotPoints->30,ImageSize->{400,300},
AxesLabel->{"phi","u","theta"}];