anchor.diff:           j = 2*int(jrot(i))+1
            energy = energies(i)
            !
  +         if (vibtemperature_do) then
  +           if (nint(jrot(i))==0) then 
  +              energies_vib(i) = energy
  +              ivib_state(i) = i
  +              ener_vib = energy
  +              ener_rot = energy-ener_vib
  +           else
  +              loop_vib : do ivib  = 1,Nvib_states
  +                if (all(quantum_numbers(ivib1:ivib2,i)==quantum_numbers(ivib1:ivib2,ivib))) then
  +                  !
  +                  ivib_state(i) = ivib
  +                  ener_vib = energies_vib(ivib)
  +                  ener_rot = energy-ener_vib
  +                  exit loop_vib
  +                  !
  +                endif
  +              enddo loop_vib
  +           endif
  +         endif
  +         !
            if (partfunc_do) then
              if (j/=j0) then
                 !
  
positive.code:          !
           maxitems = max(nitems,maxitems)
           !
           !cycle
           !
           !10  exit
        end do
  
negative.code:       i = 0
        iline = 0
        maxitems = 5
        Nvib_states = 0
        info = 0
        !
        call input_options(echo_lines=.false.,error_flag=1)
        !
  