anchor.diff:     !!!enddo
      do iHNF = 1, nHNF ! Write this labeling for each corresponding HNF
         jHNF = vsH(iHNF) ! Index of matching HNFs
  +      ! check if concentrations of this labeling match the user specification:
         if (check_labeling_numbers(pplabeling,number_ElementN,number_Range)) then
           Tcnt = Tcnt + 1; Scnt = Scnt + 1
           write(14,struct_enum_out_formatstring) &
  
positive.code: allocate(allD(nAllD)); allD = (/ (i,i=1,nAllD) /)
  allocate(allD2LabelD(nAllD));
  ! 1) Check whether we have to postprocess the labeling before writing it out
  postprocessLabeling = .not. (all(  abs( equivalencies-allD ) ==0))
  allocate(pplabeling(n*nAllD))  ! this is ok whether we do postprocessing (nAllD>nD) or not (nAllD==nD)
  ! 2) if yes:
  
negative.code: integer, intent(inout) :: Tcnt, Scnt, Hcnt ! counters for total number of labelings and number of this size and HNF running total 
  integer, intent(in) :: permIndx(:)
  character, intent(in) :: lm(:)
  integer, intent(in) :: number_ElementN     ! concentration check parameter: which element's concentration is checked?
  integer, intent(in) :: number_Range(2)     ! concentration check parameter: the lower and upper limit of the number of atoms
                                             !                                if spezies "number_ElementN" in the unit cell of 
                                             !                                size "n"
  
  integer nHNF ! Number of HNFs in the list that match the current block index, HNFi
  integer, allocatable :: vsH(:), vsL(:) ! Vector subscript for matching the HNF index to the HNF list
  