Create reusable routine for filenames
We have many calls of type
if (this%nspin == 1) then
write(filename, fmt='(a)') 'vhxc'
else
write(filename, fmt='(a,i1)') 'vhxc-sp', isp
end if
We should create a small function for this.