% Reproduces the solid line in Figure 3 of Westbrook (2013) 'Rayleigh scattering by hexagonal ice crystals and the interpretation of dual-polarisation radar measurements' QJRMS 

eps=3.15 % solid ice
V=1 % unit volume 

logw=-2:0.1:2; w=10.^logw; ZDR=logw*0; % create a list of aspect ratios w to loop through

for jj=1:length(w)
	[Xxx,Xzz]=hex_polarisabilities(V,eps,w(jj))
	if w(jj)>1.1
		ZDR(jj)=Xxx^2/Xzz^2;
	else
		ZDR(jj)=3/8*Xzz^2/Xxx^2+3/8+1/4*Xzz/Xxx;
	end
end

% Plot results
figure; semilogx(w,10*log10(ZDR),'k'); grid on; ylabel('ZDR [dB]'); xlabel('Aspect ratio w'); title('using shape factors from Westbrook (2013) QJRMS')


