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