anchor.diff:        enddo 
      endif
      !
  -   if (hitran.and.partfunc<0.0) then
  +   if (hitran_do.and.partfunc<0.0) then
        write(out,"('For HITRAN partition function must be defined using PF or QSTAT keywords')")
        stop 'Undefined PF'
      endif
  
positive.code:           !
          case ("HITRAN")
            !
            hitran = .true.
            !
          case ("SELECT","FILTER")
            !
            filter  = .true.
  
negative.code:               case ("N")
                  !
                  call readf(species(i)%N)
                  !
                case("DELTA")
                  !
                  call readf(species(i)%delta)
                  !
                case ("RATIO")
                  !
  