function n=ref_ice(lambda)
% The complex refractive index of ice
%Reference: Warren, S. G., and R. E. Brandt (2008), Optical constants of ice from the ultraviolet to the microwave: A revised compilation.
%J. Geophys. Res., 113, D14220, doi:10.1029/2007JD009744. PDF link
%This research was supported by the National Science Foundation under grants ANT-00-03826 and OPP-06-36993.
%The 2008 compilation is a revision of the 1984 compilation:
%Warren, S.G., 1984: Optical constants of ice from the ultraviolet to the microwave. Applied Optics, 23, 1206-1225. PDF link
% input:
%      lambda in [um]
% output:
%      n - complex refractive index
x=load('IOP_2007_ASCIItable.dat');
m=interp1(x(:,1),x(:,2),lambda);
k=interp1(x(:,1),x(:,3),lambda);
n=complex(m,k);

