% #this funtion returns refraction times
% # input parameters: t_intercept is intercept time; v_refractor; x source receiver distance
% #	  x - source-receiver distance (can be an scalar or array, 
function [t_refract]=trefr(t_intercept,v_refractor,x)
	
	t_refract=t_intercept+abs(x)./v_refractor;
end