Home > atmlab > mie > mie_abs.m

mie_abs

PURPOSE ^

Computation of the Absorption Efficiency Qabs

SYNOPSIS ^

function result = mie_abs(m, x)

DESCRIPTION ^

 Computation of the Absorption Efficiency Qabs
 of a sphere of size parameter x, 
 complex refractive index m=m'+im", 
 based on nj internal radial electric field values
 to be computed with Mie_Esquare(nj,m,x)
 Ref. Bohren and Huffman (1983) BEWI:TDD122,
 and my own notes on this topic;
 k0=2*pi./wavelength;
 x=k0.*radius;
 C. M�zler, May 2002, revised July 2002.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

mie_abs.m

SOURCE CODE ^

0001 function result = mie_abs(m, x)
0002 
0003 % Computation of the Absorption Efficiency Qabs
0004 % of a sphere of size parameter x,
0005 % complex refractive index m=m'+im",
0006 % based on nj internal radial electric field values
0007 % to be computed with Mie_Esquare(nj,m,x)
0008 % Ref. Bohren and Huffman (1983) BEWI:TDD122,
0009 % and my own notes on this topic;
0010 % k0=2*pi./wavelength;
0011 % x=k0.*radius;
0012 % C. M�zler, May 2002, revised July 2002.
0013 
0014 nj=100*round(2+x+4*x.^(1/3))+300;
0015 e2=imag(m.*m);
0016 dx=x/nj;
0017 x2=x.*x;
0018 nj1=nj+1;
0019 xj=(0:dx:x);
0020 en=Mie_Esquare(m,x,nj);
0021 en1=0.5*en(nj1).*x2;     % End-Term correction in integral
0022 enx=en*(xj.*xj)'-en1;    % Trapezoidal radial integration
0023 inte=dx.*enx;
0024 Qabs=4.*e2.*inte./x2;
0025 result=Qabs;

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005