function [ kh2o ] = continuum1000( v,T )
% Continuum absorption coefficient of water vapor (self-broadened) around
% 1000 cm-1
% From Pierrehumbert 2010, p260
% calculates capture cross section, kh2o in m2/kg for each wavenumber 
% in v between minv, maxv, at each temperature in T
% value of kh2o will be zero outside of minv - maxv

minv=500; maxv=1400;

[vv TT]=meshgrid(v,T);
kh2o=exp(12.167-0.050898.*vv+8.3207e-5.*vv.^2-7.0748e-8.*vv.^3+2.3261e-11.*vv.^4).*((296./TT).^4);

iv=find(v<minv | v>maxv);
kh2o(iv)=0;

end

