phonons of Na2S
 &inputph
  recover = .true.,
   tr2_ph=1.0d-12,
   alpha_mix(1)=0.6,
    prefix='Na2S',
!  fildvscf='na2sdv',
  amass(1)=32.066,
  amass(2)=22.988,
  outdir='./OUT',
  fildyn='na2s.dyn',
!  elph=.true.,
!  trans=.true.,
   ldisp=.true.
   nq1=4, nq2=4, nq3=4
 /