anchor.diff:  real(8), allocatable :: mommt(:,:)
   ! total muffin-tin moment
   real(8) mommttot(3)
  -
  +end module
  +module mod_APW_LO
   !-----------------------------------------!
   !     APW and local-orbital variables     !
   !-----------------------------------------!
   ! maximum allowable APW order
   integer, parameter :: maxapword=3
   ! APW order
  -integer apword(0:maxlapw,maxspecies)
  +integer apword(0:_MAXLAPW_,_MAXSPECIES_)
   ! maximum of apword over all angular momenta and species
   integer apwordmax
   ! APW initial linearisation energies
  -real(8) apwe0(maxapword,0:maxlapw,maxspecies)
  +real(8) apwe0(maxapword,0:_MAXLAPW_,_MAXSPECIES_)
   ! APW linearisation energies
   real(8), allocatable :: apwe(:,:,:)
   ! APW derivative order
  -integer apwdm(maxapword,0:maxlapw,maxspecies)
  +integer apwdm(maxapword,0:_MAXLAPW_,_MAXSPECIES_)
   ! apwve is .true. if the linearisation energies are allowed to vary
  -logical apwve(maxapword,0:maxlapw,maxspecies)
  +logical apwve(maxapword,0:_MAXLAPW_,_MAXSPECIES_)
   ! APW radial functions
   real(8), allocatable :: apwfr(:,:,:,:,:)
   ! derivate of radial functions at the muffin-tin surface
  
positive.code: real(8), allocatable :: tpgkc(:,:,:,:)
  ! structure factor for the G+k-vectors
  complex(8), allocatable :: sfacgk(:,:,:,:)
  
  !-------------------------------!
  !     q-point set variables     !
  !-------------------------------!
  
negative.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(:,:)
  
  end module
  module mod_spin
  !--------------------------------!
  !     spin related variables     !
  !--------------------------------!
  