anchor.diff:              call readf(dtemp)
               call readf(ji)
               !
  -            if (info>0.and.nint(ji)==0) then
  +            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 (nint(ji)==0) then 
  +            if (int(ji)==0) then 
                 Nvib_states = Nvib_states + 1
               else 
                 info = 1
  
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:                   ener_vib = energies_vib(ivib)
                    ener_rot = energyi-ener_vib
                  else
                    ivib = ivib_state(ilevelf)
                    ener_vib = energies_vib(ivib)
                    ener_rot = energyf-ener_vib
                  endif
               endif
               !
  