&input
  asr = 'crystal'
  flfrc = 'Si.fc'
  flfrq = 'Si.freq'
  q_in_band_form=.true.
/
3
0.5 0.5 0.5  20
0.0 0.0 0.0  20
1.0 0.0 0.0  1