Let $P$ be any partition of $[-1,1]$. Since $f$ contains positive values on every interval, the upper sum $U(P,f)$ is positive. Hence
$$\int f = \inf\limits_{P} U(P,f) \geq 0.$$
Similarly, $L(P,f)$ is negative, so that
$$\int f = \sup\limits_{P} L(P,f) \leq 0.$$
A more interesting question, which I am unsure of the answer for, is whether there is such a function that is Riemann integrable.
First of all, it seems like you might be thinking that all bounded functions are Riemann integrable, so I want to disavow you of that notion. Take $f(x) = 1$ for $x$ rational and $0$ for $x$ irrational. This is bounded on all compact intervals but is not Riemann integrable on any.
Strictly speaking, the definition of a Riemann integral is on a bounded function $f: [a, b] \longrightarrow \mathbb R$, so $f$ must be defined on all of $[a, b]$ no ifs ands or buts. If we had some $f$ defined on $[a, b] - A$ then it is technically incorrect to write $\int_a^b f(x) dx$. In fact, we cannot even ask if $f$ is Riemann integrable on $[a, b]$ as it is not a real function defined on $[a, b]$.
However, we are not always so strict in math. We can choose to abuse the notation of the Riemann integral and write $\int_a^b f(x) dx$ even when $f$ is not defined everywhere on $[a, b]$, but we must justify this abuse. Indeed, let $g_1, g_2: [a, b] \longrightarrow \mathbb R$ such that $g_1$ is Riemann integrable and $g_2$ differs from $g_1$ at only finitely many points. Then $g_2$ is Riemann integrable and $\int_a^b g_1(x) dx = \int_a^b g_2(x) dx$. To see why, we can work with the difference $h = g_1 - g_2$ and prove that $h$ is Riemann integrable with integral $0$. Roughly speaking, as $h$ is nonzero at only finitely many points, shrinking our partition around these points allows our Riemann sums to tend to $0$.
What this fact about $g_1, g_2$ says is that the Riemann integral does not care about the values of your function on any finite set. This is why we feel justified in talking about the Riemann integrability of a partially defined function $f: [a, b] - A \longrightarrow \mathbb R$. Whatever extension to the whole of $[a, b]$ we choose makes no difference - if $g_1, g_2: [a, b] \longrightarrow \mathbb R$ both extend $f$ then $g_1$ is Riemann integrable iff $g_2$ is, in which case their integrals coincide. We can therefore comfortably say that $f$ is Riemann integrable, when we mean that an extension is, and we can comfortably compute its integral by computing the integral of any extension.
To put this a bit more formally, let $R([a, b])$ be the set of all Riemann integrable functions on $[a, b]$. This is actually a real vector space and we have a linear map $\int_a^b: R([a, b]) \longrightarrow \mathbb R$ sending $f \mapsto \int_a^b f(x) dx$. Now, we define an equivalence relation on $R([a, b])$ via $f \sim g$ if $f(x) = g(x)$ at all but finitely many points. The quotient $R([a, b])/{\sim}$ is a real vector space, and the integral descends to a well defined map $R([a, b])/{\sim} \longrightarrow \mathbb R$ sending $[f] \mapsto \int_a^b f(x) dx$. This is just formalism for my above statement that the Riemann integral doesn't care about the values on a finite set. Additionally, take a partially defined function $f: [a, b] - A \longrightarrow \mathbb R$ where $A$ is finite. Suppose that $f$ can be extended to a Riemann integrable function on $[a, b]$. Then there is no canonical choice of $g \in R([a, b])$ to assign to $f$. However, there is a canonical equivalence class to assign to $f$ - that is, $[g]$ where $g$ is any extension of $f$. In fact, we can feel justified referring to this equivalence class as just $[f]$. Point being, this object $[f]$ is still specific enough to talk about Riemann integrability as much as we want. As an aside, this equivalence relation idea is extremely fruitful when dealing with the Lebesgue integral, where we replace "finite" with "measure 0."
Best Answer
An sketch for a proof, fill in the details.
Let $E_n:=\{x\in[a,b]:f(x)\geqslant 1/n\}$, then note that if $E_n$ is finite then $M:=\sup_{x\in[a,b]} f(x)<\infty $ and $M$ is attained, at most, at a finite number of points.
Now for any $x\in E_n$ you can set a closed interval containing it of length less or equal to $1/(n|E_n|)$, that is $[x-1/(2n|E_n|),x+1/(2n|E_n|)]\cap [a,b]$, where $|E_n|$ is the cardinality of $E_n$. Then for any partition $P$ with mesh smaller to $1/(2n|E_n|)$ we have that $$ 0\leqslant R(f,P)\leqslant \frac{b-a}n+\frac{M}n $$ where follows that the integral is zero.