Let $F:\mathbb{R} \to \mathbb{R}$ be locally Lipschitz, monotone and continuous. For the sake of concreteness only let us suppose it is of porous medium type (eg. $F(r) = r^{\frac 1m}$.)

Let $\Omega \subset \mathbb{R}^n$ be a bounded domain and let $Q=(0,T)\times\Omega$.

Given $u_0 \in L^\infty(\Omega)$ and $f \in L^\infty(Q)$, I have $F(u) \in L^\infty(Q) \cap L^2(0,T;H^{-1})$ with $u \in L^\infty(Q) \cap L^2(0,T;H^1)$ such that $$\int_0^T \langle (F(u))_t, \varphi \rangle + \int_0^T\int \nabla u \cdot \nabla \varphi = \int_0^T \int f\varphi$$ for all test functions $\varphi$. I also have the continuous dependence result for two solutions corresponding two two data: $$\lVert F(u_1) - F(u_2) \rVert_{L^1(0,T;L^1)} \leq C\left(\lVert u_{01} - u_{02}\rVert_{L^1} + \lVert f_1 - f_2 \rVert_{L^1(0,T;L^1)}\right).$$ Now I wish to extend my existence result to $L^1$ data satisfying a weaker formulation $$-\int_0^T \int F(u) \varphi_t - \int_0^T\int u \Delta \varphi = \int_0^T \int f\varphi$$ for smooth $\varphi$. Are there any standard tricks to do this using this continuous dependence result?

Of course we can approximate the $L^1$ data by $L^\infty$ data and using the above estimate (we obtain a Cauchy sequence and so) we find $F(u_n) \to F$ in $L^1(L^1)$ for some $F$ where $F(u_n)$ is the solution with data that approximates the $L^1$ data. From this I can obtain $u_n \to u$ for some $u$ pointwise a.e., but this is not enough pass to the limit.

**Edit**: I am aware of the book on PME by J Vazquez. If I recall correctly he handles $L^1$ data rather differently and I would like to know whether the above approach can work.