Xg = 100:.5:180;
Yg = -50:.5:0;
[Xg Yg]=meshgrid(Xg,Yg);
