--
&inputepw
  prefix      = 'si'
  amass(1)    = 28.0855
  outdir      = './'

  iverbosity  = 0

  elph        = .true.
  kmaps       = .true.
  epbwrite    = .false.
  epbread     = .false.
  epwwrite    = .false.
  epwread     = .true.
  use_ws      = .false.

  scattering  = .true.
  scattering_serta = .true.
  int_mob     = .false.
  carrier     = .true.
  ncarrier    =  -1E13
  mp_mesh_k   = .true.
  iterative_bte = .true.
  epmatkqread = .false.
  mob_maxiter = 20
  broyden_beta= 0.7

  nstemp      =   1
  temps       = 100

  restart      = .true.
  restart_step = 5
  selecqread   = .false.

  etf_mem     = 1
  lpolar      = .false.
  vme         = .true.
  lphase      = .true.

  nbndsub     =  8

  wannierize  = .false.
  num_iter    = 1500
  iprint      = 2
  dis_win_max = 18
  dis_froz_max= 8.5
  proj(1)     = 'Si : sp3'   
  wdata(1) = 'bands_plot = .true.'
  wdata(2) = 'begin kpoint_path'
  wdata(3) = 'L 0.50 0.00 0.00 G 0.00 0.00 0.00'
  wdata(4) = 'G 0.00 0.00 0.00 X 0.50 0.50 0.00'
  wdata(5) = 'end kpoint_path'
  wdata(6) = 'bands_plot_format = gnuplot'
  wdata(7) = 'guiding_centres = .true.'
  wdata(8) = 'dis_num_iter      = 500'
  wdata(9) = 'num_print_cycles  = 10'
  wdata(10) = 'dis_mix_ratio     = 1.0'
  wdata(11) = 'conv_tol = 1E-9'
  wdata(12) = 'conv_window = 4'

  elecselfen  = .false.
  phonselfen  = .false.
  a2f         = .false.

  fsthick     = 0.4
  degaussw    = 0.0 ! eV

  dvscf_dir   = './save/'

  efermi_read  = .true
  fermi_energy = 6.355486 ! 6.255486
  
  nkf1        = 12
  nkf2        = 12
  nkf3        = 12

  nqf1        = 12
  nqf2        = 12
  nqf3        = 12

  nk1         = 6
  nk2         = 6
  nk3         = 6

  nq1         = 2
  nq2         = 2
  nq3         = 2
 /
