%Copyright WaVe 2012
f = 86; %frequency (Ghz)
lambda = (3*10^8)/(f*10^9);
e = 8.852 + j*15.41 ; %complex dielectric constant of water
rho = 1; %density of water (g/cm3)
cloud = [3 4 1; 2 3 0.6; 0.8 10 1; 0.6 10 0.4];

%Specific attenuation based on cloud type
for i = 1:4
    alphac(i) = 0.4343*((3*pi*cloud(i,3))/(32*lambda*rho))*abs(imag((1-e)/(2+e)));
    Ac(i) = alphac(i)*cloud(i,1); %Total zenith attenuation by cloud
end
Ac