anchor.diff:  
   contains
   
  +! Purpose: a concentration check routine. But it does not work with real(dp) concentrations
  +!          but rather with integers: For a given unit cell size, the range(2) is calculated
  +!          before the call to this function and contains the lower and upper limit of number 
  +!          of atoms of species el
   logical function check_labeling_numbers(labeling,el,range)
   integer, intent(in) :: labeling(:)
   integer, intent(in) :: el
  
positive.code:    !!!enddo
     do iHNF = 1, nHNF ! Write this labeling for each corresponding HNF
        jHNF = vsH(iHNF) ! Index of matching HNFs
        if (check_labeling_numbers(pplabeling,number_ElementN,number_Range)) then
          Tcnt = Tcnt + 1; Scnt = Scnt + 1
          write(14,struct_enum_out_formatstring) &
  
negative.code:      labIndx = labIndx - quot*multiplier(ilab) ! Take the remainder for the next step
     enddo
     if (postprocessLabeling) then
       ! see the comments at the beginning of the current routine
       call postprocess_labeling(n,nAllD,labeling,pplabeling,allD2LabelD) 
     else
       pplabeling = labeling ! nothing changes
  