Home > atmlab > mie > mie_S12.m

mie_S12

PURPOSE ^

Computation of Mie Scattering functions S1 and S2

SYNOPSIS ^

function result = mie_S12(m, x, u)

DESCRIPTION ^

 Computation of Mie Scattering functions S1 and S2
 for complex refractive index m=m'+im", 
 size parameter x=k0*a, and u=cos(scattering angle),
 where k0=vacuum wave number, a=sphere radius;
 s. p. 111-114, Bohren and Huffman (1983) BEWI:TDD122
 C. M�zler, May 2002

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

mie_S12.m

SOURCE CODE ^

0001 function result = mie_S12(m, x, u)
0002 
0003 % Computation of Mie Scattering functions S1 and S2
0004 % for complex refractive index m=m'+im",
0005 % size parameter x=k0*a, and u=cos(scattering angle),
0006 % where k0=vacuum wave number, a=sphere radius;
0007 % s. p. 111-114, Bohren and Huffman (1983) BEWI:TDD122
0008 % C. M�zler, May 2002
0009 
0010 nmax=round(2+x+4*x^(1/3));
0011 ab=mie_ab(m,x);
0012 an=ab(1,:);
0013 bn=ab(2,:);
0014 
0015 pt=mie_pt(u,nmax);
0016 pin =pt(1,:);
0017 tin =pt(2,:);
0018 
0019 n=(1:nmax);
0020 n2=(2*n+1)./(n.*(n+1));
0021 pin=n2.*pin;
0022 tin=n2.*tin;
0023 S1=(an*pin'+bn*tin');
0024 S2=(an*tin'+bn*pin');
0025     
0026 result=[S1;S2];

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