Home > atmlab > mie > miecoated_prscan.m

miecoated_prscan

PURPOSE ^

Computation and plot of Mie Efficiencies of coated

SYNOPSIS ^

function result = miecoated_prscan(m1,m2,y,nsteps)

DESCRIPTION ^

 Computation and plot of Mie Efficiencies of coated
 spheres for given complex refractive-index ratios
 m1,2=m1',2'+im1",2" in kernel and coating, resp.
 and size parameters x=k0*a, y=k0*b, vs. x/y=1-wr (see below)
 with nsteps increments 
 a,b=inner,outer sphere radius, using complex Mie coefficients 
 an and bn, according to Bohren and Huffman (1983) BEWI:TDD122
 result: m1,m2,x,y,wr, efficiencies for extinction (qext), 
 scattering (qsca), absorption (qabs), backscattering (qb), 
 qratio=qb/qsca, asymmetry parameter (asy=<costeta>),
 volume fraction of coating w= 1-(a/b)^3 = 1-(x/y)^3
 C. M�zler, July 2002.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

miecoated_prscan.m

SOURCE CODE ^

0001 function result = miecoated_prscan(m1,m2,y,nsteps)
0002 
0003 % Computation and plot of Mie Efficiencies of coated
0004 % spheres for given complex refractive-index ratios
0005 % m1,2=m1',2'+im1",2" in kernel and coating, resp.
0006 % and size parameters x=k0*a, y=k0*b, vs. x/y=1-wr (see below)
0007 % with nsteps increments
0008 % a,b=inner,outer sphere radius, using complex Mie coefficients
0009 % an and bn, according to Bohren and Huffman (1983) BEWI:TDD122
0010 % result: m1,m2,x,y,wr, efficiencies for extinction (qext),
0011 % scattering (qsca), absorption (qabs), backscattering (qb),
0012 % qratio=qb/qsca, asymmetry parameter (asy=<costeta>),
0013 % volume fraction of coating w= 1-(a/b)^3 = 1-(x/y)^3
0014 % C. M�zler, July 2002.
0015 
0016 m1p=real(m1); m1pp=imag(m1);
0017 m2p=real(m2); m2pp=imag(m2);
0018 nx=(1:nsteps)';
0019 pr=(nx-0.5)/nsteps;               % note that pr=a/b=(1-w)^(1/3)
0020 x=y*pr;
0021 for j = 1:nsteps
0022     a(j,:)=miecoated(m1,m2,x(j),y,1);
0023 end;
0024 output_parameters='m1(real,imag), m2(real, imag), x, y, Qext, Qsca, Qabs, Qb, <costeta>, Qb/Qsca'
0025 
0026 % plotting the results
0027 plot(pr,a(:,1:5))
0028 legend('Qext','Qsca','Qabs','Qb','<costeta>')
0029 title(sprintf('Mie Efficiencies of coated sphere, y=%g, m1=%g+%gi, m2=%g+%gi ',y,m1p,m1pp,m2p,m2pp))
0030 xlabel('a/b')
0031 
0032 result=a;

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