Home > atmlab > mie > miecoated_wrscan.m

miecoated_wrscan

PURPOSE ^

Computation and plot of Mie Efficiencies of coated

SYNOPSIS ^

function result = miecoated_wrscan(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. 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_wrscan.m

SOURCE CODE ^

0001 function result = miecoated_wrscan(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. 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 opt=1;
0017 m1p=real(m1); m1pp=imag(m1);
0018 m2p=real(m2); m2pp=imag(m2);
0019 nx=(1:nsteps)';
0020 wr=0.1*(nx-0.5)/nsteps;            % note that wr=(b-a)/b=1-(1-w)^1/3)
0021 x=y*(1-wr);
0022 for j = 1:nsteps
0023     a(j,:)=miecoated(m1,m2,x(j),y,opt);
0024 end;
0025 output_parameters='m1(real,imag), m2(real, imag), x, y, Qext, Qsca, Qabs, Qb, <costeta>, Qb/Qsca'
0026 
0027 % plotting the results
0028 loglog(wr,a(:,1:5))
0029 legend('Qext','Qsca','Qabs','Qb','<costeta>')
0030 title(sprintf('Mie Efficiencies of coated sphere, y=%g, m1=%g+%gi, m2=%g+%gi ',y,m1p,m1pp,m2p,m2pp))
0031 xlabel('(b-a)/b')
0032 
0033 result=a;

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