function y=edr(lidar,prof)

% lidar parameters calculation

c=3.0e-08; %speed of light
sigp=1.8e-7; % pulse width
dp=lidar.rangegatelength; %gate length of measurement
dr=sqrt(log(2))*c*sigp;  %full width at half max of lidar sensing volume

effective_resolution=dr+dp;

mu=(sqrt(2*log(2))*dp/dr;
chi=

g=inline('FF(x).*(UF(x,chi,m)-UF(x,chi,0))')

%select data subset matching profile extent

ind=find(lidar.zd(1,:)>prof.heights(1) & lidar.zd(1,:)<prof.heights(end)); %range index

diff_az=mean(diff(lidar.azimuth(:,1)));
