anchor.diff:                    ener_vib = energies_vib(ivib)
                     ener_rot = energyi-ener_vib
                   else
  -                  ivib = ivib_state(ileveli)
  +                  ivib = ivib_state(ilevelf)
                     ener_vib = energies_vib(ivib)
  -                  ener_rot = energyi-ener_vib
  +                  ener_rot = energyf-ener_vib
                   endif
                endif
                !
  
positive.code:         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)
        endif
        !
        if (trim(proftype)=='LIFETIME') THEN
  
negative.code:          energy = energies(i)
           !
           if (vibtemperature_do) then
             if (int(jrot(i))==0) then 
                energies_vib(i) = energy
                ivib_state(i) = i
                ener_vib = energy
  