anchor.diff:             !
            case ('JJKK-X')
              !
  -           read(bunit,*,end=15) gamma_,n_,Jpp,Jp,Kp,Kpp
  +           read(bunit,*,end=15) ch_broad(1:3),gamma_,n_,Jpp,Jp,Kp,Kpp
              !
              if (abs(Jp-Jpp)>1) stop 'Jp-Jpp in broadening file'
              if (Kp>Jp.or.Kpp>Jpp) stop 'K > J in broadening file'
  
positive.code:            !
           case ('J','A0')
             !
             read(bunit,*,end=15) gamma_,n_,J
             species(i)%gammaQN(J,:) = gamma_
             species(i)%nQN(J,:) = n_
             !
           case ('JJ','A1')
             !
             read(bunit,*,end=15) gamma_,n_,Jpp,Jp
             !
             if (abs(Jp-Jpp)>1) stop 'Jp-Jpp in broadening file'
             !
  
negative.code:      real(rk)  :: Jf,Ji,abscoef,acoef,tranfreq,energyf,energyi
       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
       character(9) b_fmt
       !
       !
  