T=6;
resol=250;
incid_angle=deg2rad(30);
%Vp=120;
Hp=2000;
Ha=500:15:1800;
R=(Hp -Ha)* tan(incid_angle);
delta_t_data= 5.625/360;
delta_t_plane=0.01;
delta_d_plane=delta_t*Vp;
N=T/delta_t;
az=(deg2rad(0:5.625:359))'*ones(1,T);
az=reshape(az,length(az)*T,1);
p_position_x=zeros(N,1);
p_position_y=((1:N)'.*delta_d_plane)-delta_d_plane;  

delta_x=250;
delta_y=250;
x_grid=-3000:delta_x:3000;
y_grid=-3000:delta_y:3000;
z_grid=0:30:1000;

[x_3dgrid,y_3dgrid,z_3dgrid]=meshgrid(x_grid,y_grid,z_grid);
[x_size,y_size,z_size]=size(x_3dgrid);
coverage=zeros(x_size,y_size,z_size);
es=[];

i=1;
Pc=1:length(t)
while d < resol 
    d=d+(Vp(i)/
    
    
    
    u=find(b>0);
temp=fix(b(u)/(pi));
temp2=b(u)-temp*pi;
bb=ones(100,1);
bb(u)=temp2;
u=find(b<0);
ess=-b(u);
temp=fix(-b(u)/(pi));
temp3=-(-b(u)-(temp+1)*pi);
bb(u)=temp3;



