phonons of Si at Gamma
 &inputph
  tr2_ph=1.0d-14,
  prefix='si',
  epsil=.true.,
  amass(1)=28.08,
  outdir='./',
  fildyn='si.dynG',
 /
0.0 0.0 0.0
EOF