The bhmie module needs significant improvement!

to check if your value makes sense compare to calculation with: http://omlc.org/calc/mie_calc.html



In [5]:

    
from atmPy.mie import bhmie



In [7]:

    
radius = .1 #micrometer
lambda_laser = .650 #micrometer
sizeParam = 2 * np.pi/lambda_laser * radius
n = complex(1.5,0.01) #refractive index
noOfAngles = 100 # this is the number of angles to be calculated between 0 and 90 degrees (i think)
mie = bhmie.bhmie_hagen(sizeParam, n, noOfAngles, 2*radius) # the last parameter is nesessary to calculate the scattering efficiency.
res = mie.return_Values_as_dict()



In [8]:

    
res.keys()









    Out[8]:





['scattering_crosssection',
 'phaseFct_S1',
 'phaseFct_S2',
 'extinction_efficiency',
 'asymmetry_parameter',
 'backscatter_efficiency',
 'extinction_crosssection',
 'scattering_efficiency']



In [9]:

    
res['extinction_crosssection']









    Out[9]:





0.0068098202800653228



In [ ]: