PyMieScatt > example > ScatteringFunction() > SL, SR, SU

Environment

GeForce GTX 1070 (8GB)
ASRock Z170M Pro4S [Intel Z170chipset]
Ubuntu 16.04.4 LTS desktop amd64
TensorFlow v1.7.0
cuDNN v5.1 for Linux
CUDA v8.0
Python 3.5.2
IPython 6.0.0 -- An enhanced Interactive Python.
gcc (Ubuntu 5.4.0-6ubuntu1~16.04.4) 5.4.0 20160609
GNU bash, version 4.3.48(1)-release (x86_64-pc-linux-gnu)
scipy v0.19.1
geopandas v0.3.0
MATLAB R2017b (Home Edition)
ADDA v.1.3b6
gnustep-gui-runtime v0.24.0-3.1
PyMieScatt v1.7.0

Reference

• Bohren C.F. and D. R. Huffman, "Absorption and Scattering of Light by Small Particles", John Wiley, New York, NY (1983)
• PDF: MATLAB Functions for Mie Scattering and Absorption, Research Report No. 2002-08, June 2002 by Christian Mätzler

ScatteringFunction

code

On Jupyter Notebook

test_ScatteringFunction_180407.ipynb

%matplotlib inline
import PyMieScatt as PMS
import matplotlib.pyplot as plt
import numpy as np
from pylab import rcParams

# to enlarge figure
rcParams['figure.figsize'] = 10, 4
rcParams['figure.dpi'] = 110

m = 1.55 + 0.001j
wavelength = 10.
diameter = 5.

res = PMS.ScatteringFunction(m, wavelength, diameter)
thetas, SLs, SRs, SUs = res

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)

ax.plot(thetas, SLs, label='SL')
ax.plot(thetas, SRs, label='SR')
ax.plot(thetas, SUs, label='SU')

ax.set_title('SL, SR, SU (m=1.55 + 0.001j, wv=10nm, dia=5nm)')
ax.set_xlabel('Theta(rad)')
ax.set_xlim(0, np.pi)
ax.legend()