anchor.diff:           !
            maxitems = max(nitems,maxitems)
            !
  -         !cycle
  +         if (vibtemperature_do) then 
  +            call readf(dtemp)
  +            call readf(ji)
  +            !
  +            if (info>0.and.nint(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 
  +              Nvib_states = Nvib_states + 1
  +            else 
  +              info = 1
  +            endif
  +            !
  +         endif
            !
            !10  exit
         end do
  
positive.code:           !
            if (trim(w)/="".and.trim(w)/="END") then
               !
               write (out,"('input: wrong last line in CONTRACTION =',a)") trim(w)
               stop 'input - illigal last line in CONTRACTION'
               !
            endif
            !
  
negative.code: 3module spectrum
    !
    use accuracy
    use timer
  