r=0:67:67*199;
z1=r*sin(deg2rad(0.8));
z2=r*sin(deg2rad(2.2));
z3=r*sin(deg2rad(3.8));
z4=r*sin(deg2rad(5.2));
z5=r*sin(deg2rad(6.8));
z=[z1 z2 z3 z4 z5];
[th,r,z]=meshgrid(deg2rad(azimuth),r,z);
% [x,y,z]=pol2cart(th,r,z);
% surf(x,y,z1,ppi1)
% hold on
% surf(x,y,z2,ppi2)
% surf(x,y,z3,ppi3)
% surf(x,y,z4,ppi4)
% surf(x,y,z5,ppi5)
% 
% shading flat