function  d=refl_forward(m,xr)
% function  refl_forward
% Input: - Vector m0 with velocity and thickness
%        - Vector xr with station coordinates
%
% Output: Data vector d with nr arrival  times
% 
nr=length(xr);

for i=1:nr,
    d(i)= 2./(m(1))*sqrt( (xr(i)/2)^2 + m(2)^2);
end

% end