new_contents: #!/usr/bin/env python
  
  import yaml
  import numpy as np
  
  data = yaml.load(open("qpoints.yaml"))
  dynmat = []
  dynmat_data = data['phonon'][0]['dynamical_matrix']
  for row in dynmat_data:
      vals = np.reshape(row, (-1, 2))
      dynmat.append(vals[:, 0] + vals[:, 1] * 1j)
  dynmat = np.array(dynmat)
  
  eigvals, eigvecs, = np.linalg.eigh(dynmat)
  frequencies = np.sqrt(np.abs(eigvals.real)) * np.sign(eigvals.real)
  conversion_factor_to_THz = 15.633302
  print frequencies * conversion_factor_to_THz
  
  
  
subject: Add an example script in NaCle