function [total, specular, volume] = gaussian_bscatter(x,diel,s,A)
%Input parameters:
%   diel = dielectric constant
%   rms = rms height (meters)
%   L = correlation length (meters)
%   theta = incidence angle (radians)

if x >= pi/2
    total = NaN;
    return
end

ref = ((1 - sqrt(diel))/(1 + sqrt(diel)))^2; %Fresnel reflection coefficient
specular = ((ref) ./ ((s.^2).* (cos(x).^4) )) .* exp( (-(tan(x).^2)) ./ ((s.^2) ));
volume = volume_scatter(A, x);
total = specular + volume;

end
