anchor.diff:       integer(ik),intent(in) :: nchars_quanta(maxitems),gtot(nlines)
        character(len=20),intent(in) :: quantum_numbers(0:maxitems,nlines)
        integer(ik) :: nchars_,kitem,ierror_nu,iE,ierror_i,ierror,ierror_S,l,qni,qnf,ierror_f,&
  -                    ilevelf,ileveli,iswap
  +                    ilevelf,ileveli,iswap,nchars_tot
        character(9) b_fmt
        !
        !
  
positive.code:    !
     !Prepare VoigtKampff
     if(trim(proftype(1:5))=='VOI-F') then
          !
       if(trim(proftype(1:6))=='VOI-FN') then
        call fast_voigt%construct(dpwcoef,offset,dfreq,.true.)
        !call initalize_voigt_kampff(dpwcoef,dfreq,offset,voigt_index,1)
  
negative.code:     endif
      !
      if (if_species_defined.and.trim(proftype)=="VOI-FAST") then
        write(out,"('Warning: VOI-FAST is now used together with SPECIES')")
        !stop 'Input Error: VOI-FAST cannot be used together with Species'
      endif
      !
      if (vibtemperature_do.and..not.if_QN_defined) then 
  