The presence of negative frequency and artificially small frequency cause overestimation of kappa.
Now, I add ad hoc treatment of cutoff of frequency in get_Sij
inv_omega=np.zeros(nmodes)
for i in range(nmodes):
#tentative
if omega[i] >0.1:
inv_omega[i]=1.0/np.sqrt(omega[i])
else:
inv_omega[i]=0.0
Consider better solution especially imposin the acoustic sum rule
The presence of negative frequency and artificially small frequency cause overestimation of kappa. Now, I add ad hoc treatment of cutoff of frequency in get_Sij
Consider better solution especially imposin the acoustic sum rule