function del=refraction(el)

%refraction in Earth atmosphere
% corection in elevation angle
del=3.51561*(0.1594+0.0196*el+0.00002*el.^2)./(1.+0.505*el+0.0845*el.^2);
a=find(el<-0.56);
del(a)=0.56;

