/// <summary>
/// Returns the scattering coefficient for a given wavelength
/// </summary>
/// <param name="wavelength">Wavelength</param>
/// <returns>The scattering coefficient Mus</returns>
public double GetMus(double wavelength)
{
MieScattParams = new MieScatteringParameters();
BohrenHuffmanMie(wavelength);
double qSca = MieScattParams.Q[0];
//Let a = particle radius;
//fv = particle volume fraction;
//rho_s = fv/((4/3)*Pi/a^3) = particle number density;
//sigma_s = Qsca*A = cross-section
//A = Pi*a^2 = area
// Then mus = rho_s*sigma_s, which, upon simplifying, produces:
// mus = (fv*Qsca)/((4/3)*a)
double mus = VolumeFraction * qSca / ((4.0 / 3.0) * ParticleRadius * 1e-3); //radius in m
return(mus);
}
/// <summary>
/// Returns the anisotropy coefficient for a given wavelength
/// </summary>
/// <param name="wavelength">Wavelength</param>
/// <returns>The anisotropy coeffient g</returns>
public double GetG(double wavelength)
{
MieScattParams = new MieScatteringParameters();
BohrenHuffmanMie(wavelength);
return(MieScattParams.G);
}