&INPUTPP
  outdir = './tmp'
  prefix = 'C.scf'
  filplot = 'wave'
  plot_num = 8
/

#KS_write yes

bphase yes

ks_block 1 4

bphase_dir all # "all" for computation of maximmally joint WCs (MJWC) (step 6b)
# bphase_dir x (or y or z) for hermaprodite WCs (HWC) along x (or y or z)  (step 6a)  

scf_output nscf.out

#only_gamma yes # if berry phase matrices to be calculated only for the gamma point

# single_point yes # if only single point berry phase matrices to be calculated