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