Home > atmlab > mie > mie_xscansmooth.m

mie_xscansmooth

PURPOSE ^

Computation and plot of Mie Efficiencies for given

SYNOPSIS ^

function result = mie_xscansmooth(m, nsteps, dx, Dx)

DESCRIPTION ^

 Computation and plot of Mie Efficiencies for given 
 complex refractive-index ratio m=m'+im" 
 and range of size parameters x=k0*a, 
 starting at x=0 with nsteps increments of dx
 a=sphere radius, using complex Mie coefficients an and bn 
 according to Bohren and Huffman (1983) BEWI:TDD122
 result: m', m", x, efficiencies for extinction (qext), 
 scattering (qsca), absorption (qabs), backscattering (qb), 
 qratio=qb/qsca and asymmetry parameter (asy=<costeta>).
 C. M�zler, May 2002.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

mie_xscansmooth.m

SOURCE CODE ^

0001 function result = mie_xscansmooth(m, nsteps, dx, Dx)
0002 
0003 % Computation and plot of Mie Efficiencies for given
0004 % complex refractive-index ratio m=m'+im"
0005 % and range of size parameters x=k0*a,
0006 % starting at x=0 with nsteps increments of dx
0007 % a=sphere radius, using complex Mie coefficients an and bn
0008 % according to Bohren and Huffman (1983) BEWI:TDD122
0009 % result: m', m", x, efficiencies for extinction (qext),
0010 % scattering (qsca), absorption (qabs), backscattering (qb),
0011 % qratio=qb/qsca and asymmetry parameter (asy=<costeta>).
0012 % C. M�zler, May 2002.
0013 nsmooth=Dx/dx;
0014 nx=(1:nsteps)';
0015 x=(nx-1)*dx;
0016 for j = 1:nsteps
0017     a(j,:)=mie(m,x(j));
0018 end;
0019 output_parameters='Qext, Qsca, Qabs, Qb, <costeta>, Qb/Qsca'
0020 for j=1:6,
0021     b(:,j)=smooth(a(:,j),nsmooth);
0022 end;
0023 b(:,3)=10000*b(:,3);
0024 % plotting the smoothed results
0025 m1=real(m);m2=imag(m);
0026 plot(x,b(:,1:6))
0027 legend('Qext','Qsca','10^4*Qabs','Qb','<costeta>','Qb/Qsca')
0028 title(sprintf('Smoothed (Dx=%g) Mie Efficiencies, m=%g+%gi',Dx,m1,m2))
0029 xlabel('x')
0030 
0031 result=b;

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