1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677 |
- from scattnlay import scattnlay
- import numpy as np
- dX = 3
- Xmax = 21.0
- m = np.array([1.46 + 0.0j], dtype = np.complex128)
- theta = np.linspace(0.0, np.pi, 500, dtype = np.float64)
- result = theta*180.0/np.pi
- for xl in np.arange(dX, Xmax, dX, dtype = np.float64):
- x = np.array([[xl]], dtype = np.float64)
- terms, Qext, Qsca, Qabs, Qbk, Qpr, g, Albedo, S1, S2 = scattnlay(x, m, theta)
- S11 = S1[0].real*S1[0].real + S1[0].imag*S1[0].imag + S2[0].real*S2[0].real + S2[0].imag*S2[0].imag
- result = np.vstack((result, S11/(2.0*xl*xl)))
- result = result.transpose()
- try:
- import matplotlib.pyplot as plt
- plt.plot(result[ : , 0], result[ : , 1:])
- ax = plt.gca()
- ax.set_xlim(0, 180)
- ax.set_yscale('log')
- plt.xlabel('Theta')
- plt.draw()
- plt.show()
- finally:
- np.savetxt("scattSiO2.txt", result, fmt = "%.5f")
- print result
|