clear
workdir=pwd;
dirs = {'plus2','plus1','neutral','minus1','minus2','minus3','minus4'};
param.n = [-1];
param.mu = [-163.887300];
param.efermi = 10.126067;
param.dielectric = 5.41*eye(3);
param.site = 1;
for ii = 1:numel(dirs)
   cd(workdir)
   ef(ii) = bareFormEnergy('bulk/results/scf',[dirs{ii},'/results/scf'],param);
   ec(ii) = corrKO('bulk/results/scf',[dirs{ii},'/results/scf'],param);
   cd(workdir)
end