|
@@ -3,6 +3,7 @@
|
|
|
import numpy as np
|
|
|
import matplotlib.pyplot as plt
|
|
|
import os
|
|
|
+import scipy.special.hankel2 as H2n
|
|
|
c = 299792458.0
|
|
|
eps_0 = 8.854187817e-12 # F/m
|
|
|
pi = np.pi
|
|
@@ -130,6 +131,8 @@ def analyze(data,wl):
|
|
|
if verbose > 5:
|
|
|
print("r =",r)
|
|
|
|
|
|
+ print(h_sp_m)
|
|
|
+
|
|
|
# print("S from full field",np.real(np.cross(E,np.conj(H))))
|
|
|
|
|
|
# print("H0 air (%5.4g %+5.4gj)"%(np.real(H1_0[wl_idx]), np.imag(H1_0[wl_idx])),
|