Home > atmlab > mie > cloud_phasefunplot.m

cloud_phasefunplot

PURPOSE ^

Plot of Phase Function pm (unpolarised)

SYNOPSIS ^

function result = cloud_phasefunplot(lamda, rc, nsteps)

DESCRIPTION ^

 Plot of Phase Function pm (unpolarised)
 of clouds with normalisation to 'one' when integrated 
 over all directions/(4*pi), see Chandrasekhar 1960
 Eq. (28), 
 Input: wavelength lamda, mode radius rc, both in micron, 
 nsteps=number of scattering angles  
 Output: nstep pairs of: scattering angle, phase function  
 s. p. 111-114, Bohren and Huffman (1983) BEWI:TDD122
 C. M�zler, July 2003

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

cloud_phasefunplot.m

SOURCE CODE ^

0001 function result = cloud_phasefunplot(lamda, rc, nsteps)
0002 
0003 % Plot of Phase Function pm (unpolarised)
0004 % of clouds with normalisation to 'one' when integrated
0005 % over all directions/(4*pi), see Chandrasekhar 1960
0006 % Eq. (28),
0007 % Input: wavelength lamda, mode radius rc, both in micron,
0008 % nsteps=number of scattering angles
0009 % Output: nstep pairs of: scattering angle, phase function
0010 % s. p. 111-114, Bohren and Huffman (1983) BEWI:TDD122
0011 % C. M�zler, July 2003
0012 
0013 nsize=17;       % number of drop radii in size distribution
0014 dteta=pi/(nsteps-1);
0015 teta=(0:dteta:pi)';
0016 p=[];
0017 for j=1:nsteps,
0018     mu=cos(teta(j));
0019     c=cloud_phasefunction(lamda, rc, mu, nsize);
0020     p=[p; c(1)];
0021 end;
0022 semilogy(teta*180/pi,p,'r-')
0023 title(sprintf('Phase Function of Cloud with rc=%g micron, at lamda=%g micron ',rc,lamda))
0024 xlabel('Scattering Angle (deg)'),ylabel('Phase Function')
0025 result=[teta, p];

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