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")