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:                      energyi,species(1)%n,species(1)%delta
           !
           nchars_ = 0
           !
           do kitem = 3,maxitems
             !
             l = len(trim(quantum_numbers(kitem,ilevelf)))
             !
             b_fmt = "(1x,a3)" ; if (l>3) b_fmt = "(1x,a8)"
             !
             write(b_fmt,"('(1x,a',i1,')')") nchars_quanta(kitem)
             !
             nchars_  = nchars_ + nchars_quanta(kitem)
             !
             if (nchars_>15) cycle
             !
             write(sunit,b_fmt,advance="no") trim(quantum_numbers(kitem,ilevelf))
             !
  
negative.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)
  