def getDamping(E, phase) :
damping_in = []
damping_lz = []
for ii in range(2**NE) :
damping_in.append(abs(dmp.IonNeutral_Damping(E[ii], phase, nu_n = 0, theta = 0).get("wi")))
damping_lz.append(abs(dmp.damping_lazarian_nopos(E[ii], phase)[0]))
return [damping_in, damping_lz]
wI_Alfven_o2 = np.zeros(len(E))
Ep = np.NaN
Em = np.NaN
Gamma_lz = np.zeros(len(E))
Gamma_nlld_inf = np.zeros(len(E))
Gamma_nlld_sup = np.zeros(len(E))
for e in range(len(E)) :
in_damping = dp.IonNeutral_Damping(E[e], phases[pi], nu_n = 0, theta = 0)
wR_Alfven[e] = in_damping.get("wr")
wI_Alfven[e] = in_damping.get("wi")
if (pi == 0) :
wI_Alfven[e] = np.NaN
in_damping = dp.IN_damping_approx_1(E[e], phases[pi], nu_n = 0, theta = 0)
wR_Alfven_o1[e] = in_damping.get("wr")
wI_Alfven_o1[e] = in_damping.get("wi")
in_damping = dp.IN_damping_approx_2(E[e], phases[pi], theta = 0)
wR_Alfven_o2[e] = in_damping.get("wr")
wI_Alfven_o2[e] = in_damping.get("wi")
Ep = in_damping.get("Ep")