anchor.diff:           !
            iloggf = nint(log(gf))
            !
  +         jf = jrot(ilevelf)
  +         ji = jrot(ileveli)
  +         !
  +         gamma1 = species(1)%gamma
  +         n1     = species(1)%n
  +         !
  +         Jpp = nint(Ji)
  +         Jp  = nint(Jf)
  +         !
  +         if ( trim(species(1)%filename)/="" ) then
  +           gamma1 = species(1)%gammaQN(Jpp,Jp-Jpp)
  +           n1     = species(1)%nQN(Jpp,Jp-Jpp)
  +         endif
  +         !
  +         if ( trim(species(2)%filename)/="" ) then
  +           gamma2 = species(2)%gammaQN(Jpp,Jp-Jpp)
  +         endif
  +         !
            write(out,"(i12,1x,f12.6,1x,i6,1x,f12.4,1x,f11.4,1x,f11.4,1x,f11.4,'||')") &
                        iso,lambda,iloggf,energyi,&
  -                     species(1)%gamma,species(2)%gamma,species(1)%n
  +                     gamma1,gamma2,n1
        enddo
        !
        !ielion, wl, gf-value, xi, igr, igs, igw
  
positive.code:    !
     if (trim(enrfilename)/="none".and..not.hitran_do) then
  
        if (verbose>=2) print('(a/)'),'Reading energies from .states ..'
        !
        write(ioname, '(a)') 'Energy file'
        call IOstart(trim(ioname),enunit)
  
negative.code:                write(out,"(4i12,2x,3es16.8)") ilevelf,ileveli,indexf,indexi,acoef,energyf,energyi
                 stop 'wrong order of indices'
                 cycle loop_swap
               elseif (energyf-energyi<small_) then
                 cycle loop_swap
               endif
               !
  