 &inputph
  tr2_ph=1.0d-14,
  prefix='Si',
  epsil=.false.,
  ldisp = .true.,
  nq1=4,
  nq2=4,
  nq3=4,
  amass(1)=28.0855,
  fildyn='Si.dynG',
 /
