@@ -70,6 +70,7 @@ def analyze(data, dist, z_vec, wl_idx):
kappa1= np.sqrt(k_spp**2 - eps1*k_0**2)
kappa2= np.sqrt(k_spp**2 - eps2*k_0**2)
+ # TODO def check_field_match():
H1 = data_in_air[:,6,wl_idx]
H2 = data_in_gold[:,6,wl_idx]
E1 = data_in_air[:,4,wl_idx]
