anchor.diff:  real(8), allocatable :: rwfcr(:,:,:,:)
   ! radial charge density for core states
   real(8), allocatable :: rhocr(:,:)
  -
  +end module
  +module mod_energy
   !--------------------------!
   !     energy variables     !
   !--------------------------!
  
positive.code: ! (lmaxinr+1)^2
  integer lmmaxinr
  ! number of radial points to the inner part of the muffin-tin
  integer nrmtinr(maxspecies)
  ! index to (l,m) pairs
  integer, allocatable :: idxlm(:,:)
  
  !--------------------------------!
  !     spin related variables     !
  !--------------------------------!
  
negative.code: ! Lorentzian lineshape in optics
  logical :: optltz
  ! broadening for Lorentzian lineshape
  real(8) :: optswidth
  !</sag>
  logical lmirep
  ! spin-quantisation axis in Cartesian coordinates used when plotting the
  