Home > atmlab > mie > miecoated_S12.m

miecoated_S12

PURPOSE ^

Coated-sphere Mie-Scattering functions S1 and S2

SYNOPSIS ^

function result = miecoated_S12(m1, m2, x, y, u)

DESCRIPTION ^

 Coated-sphere Mie-Scattering functions S1 and S2
 for complex refractive index ratios m1,2=m1,2'+im1,2", 
 size parameters x=k0*a, y=k0*b, 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, July 2002

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

miecoated_S12.m

SOURCE CODE ^

0001 function result = miecoated_S12(m1, m2, x, y, u)
0002 
0003 % Coated-sphere Mie-Scattering functions S1 and S2
0004 % for complex refractive index ratios m1,2=m1,2'+im1,2",
0005 % size parameters x=k0*a, y=k0*b, 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, July 2002
0009 
0010 nmax=round(2+y+4*y^(1/3));
0011 abcd=miecoated_ab1(m1,m2,x,y);
0012 an=abcd(1,:);
0013 bn=abcd(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