function beta=betcsh(r1,r2,rv1,rv2,gamma)
%  calculate coagulation kernel for shear use curvilinear kernel.   Need to use full blown approx
% to decrease the equal sizied particles as well use Hill 1992
%  input - r radius vector- cm  column vector
%          gamma   shear in 1/s
%  output- beta coagulation kernel
%                         11 jy 93
global ratrirg


con1=gamma*sqrt(8. * pi/15.);
rvec=[r1;r2];	% note: changed , to ; to match others
p=min(rvec)./max(rvec);
eff=1-(1+5*p+2.5*p.*p)./(1+p).^5;
beta=con1*eff.*(ratrirg*(r1+r2)).^3;

