anchor.diff:                       energyi,species(1)%n,species(1)%delta
            !
            nchars_ = 0
  +         nchars_tot = 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)
  +           write(b_fmt,"('(a',i1,')')") nchars_quanta(kitem)
              !
              nchars_  = nchars_ + nchars_quanta(kitem)
              !
  -           if (nchars_>15) cycle
  +           if (nchars_>21) cycle
  +           nchars_tot = nchars_tot + nchars_quanta(kitem)
              !
              write(sunit,b_fmt,advance="no") trim(quantum_numbers(kitem,ilevelf))
              !
  
positive.code:            !
             !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,ileveli))
             !
           enddo
           !
           if ( mod(nint(2*Jf),2)==0 ) then
             !
             write(sunit,'(i7,1x,a1,1x,a1)',advance="no") nint(Jf),trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
             write(sunit,'(i7,1x,a1,1x,a1)',advance="no") nint(Ji),trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
             !
           else
             !
             write(sunit,'(f7.1,1x,a1,1x,a1)',advance="no") Jf,trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
             write(sunit,'(f7.1,1x,a1,1x,a1)',advance="no") Ji,trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
             !
           endif
           !
  
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)
  