anchor.diff:            !
             call readf(temp)
             !
  +          if (Nitems>2) then 
  +             !
  +             call readu(w)
  +             !
  +             if (w(1:3)=='VIB') then
  +                 call readf(temp_vib)
  +                 vibtemperature_do  = .true.
  +               else
  +                 call report ("Illegal key, expected vib"//trim(w),.true.)
  +             endif 
  +             !
  +          endif
  +          !
           case ("PRESSURE")
             !
             call readf(pressure)
  
positive.code:          j = 2*int(jrot(i))+1
           energy = energies(i)
           !
           if (partfunc_do) then
             if (j/=j0) then
                !
  
negative.code: 3module spectrum
    !
    use accuracy
    use timer
  