--
&inputepw
  prefix      = 'MgB2',
  amass(1)    = 24.305,
  amass(2)    = 10.811
  outdir      = './'

  ep_coupling = .true.
  elph        = .true.
  kmaps       = .false.
  epbwrite    = .false.
  epbread     = .true.

  epwwrite    = .false.
  epwread     = .true.

  etf_mem     = .true.

  nbndsub     =  5,
  nbndskip    =  0

  wannierize  = .false.
  num_iter    = 500
  dis_froz_max= 8.8
  proj(1)     = 'B:pz'
  proj(2)     = 'f=0.5,1.0,0.5:s'
  proj(3)     = 'f=0.0,0.5,0.5:s'
  proj(4)     = 'f=0.5,0.5,0.5:s'

  iverbosity  = 0

  elinterp    = .true.
  phinterp    = .true.

  tshuffle2   = .true.

  parallel_k  = .true.
  parallel_q  = .false.

  eps_acustic = 5.0    ! Lowest boundary for the 
  ephwrite    = .true. ! Writes .ephmat files used when wliasberg = .true.

  fsthick     = 0.4  ! eV
  eptemp      = 300  ! K
  degaussw    = 0.10 ! eV
  degaussq    = 0.5 ! meV
  nqstep      = 500

  eliashberg  = .true.

  liso  = .true.
  lreal = .true.
  kerwrite = .true. 

  conv_thr_raxis = 1.0d-2

  wscut = 0.5   ! eV   Upper limit over frequency integration/summation in the Elisashberg eq
  wsfc  = 0.25  ! eV   Intermediate limit over frequency integration (only with lreal)

  nswc  = 100   ! frequency mesh from wsfc to wscut 
  nswfc = 200   ! frequency mesh up to wsfc 

  nstemp   = 1
  tempsmin = 25.00
  tempsmax = 30.00

  nsiter   = 500

  muc     = 0.16

  dvscf_dir   = './save'
 
  nk1         = 3
  nk2         = 3
  nk3         = 3

  nq1         = 3
  nq2         = 3
  nq3         = 3
 
  mp_mesh_k = .true.
  nkf1        = 6
  nkf2        = 6
  nkf3        = 6

  nqf1        = 6
  nqf2        = 6
  nqf3        = 6
 /
  6 cartesian
          0.000000000   0.000000000   0.000000000 1.0
          0.000000000   0.000000000   0.291867841 1.0
          0.000000000   0.384900179   0.000000000 1.0
          0.000000000   0.384900179   0.291867841 1.0
          0.333333333   0.577350269   0.000000000 1.0 
          0.333333333   0.577350269   0.291867841 1.0