anchor.diff:          call ArrayStart('ivib_state',info,size(ivib_state),kind(ivib_state))
           ivib_state = 1
           energies_vib = 0
  -        ivib1 = QN%vibcol(1)
  -        ivib2 = QN%vibcol(2)
  +        ivib1 = QN%vibcol(1)-3
  +        ivib2 = QN%vibcol(2)-3
         endif
         !
         if (trim(proftype)=='LIFETIME') THEN
  
positive.code:          energy = energies(i)
           !
           if (vibtemperature_do) then
             if (nint(jrot(i))==0) then 
                energies_vib(i) = energy
                ivib_state(i) = i
                ener_vib = energy
  
negative.code:             call readf(dtemp)
              call readf(ji)
              !
              if (info>0.and.int(ji)==0) then
                write(out,"('It is illegal to use not J-sorted states with vib-temperatures')")
                stop 'States file must be ordered by J to use with Tvib'
              endif
              !
              if (info>0) cycle
              !
              if (int(ji)==0) then 
                Nvib_states = Nvib_states + 1
              else 
                info = 1
  