phonons of Carbon diamond on a grid
 &INPUTPH
  prefix = 'CD',
  asr = .true.
  ldisp = .true.
  nq1=4, nq2=4, nq3=4
 /