Yes, this works. I'll discuss the case $n=2$. I hope this can be generalized, but to be honest, I haven't thought that through properly.
Let $u(x,y)=\int_0^x f_1(s,y)\, ds$, $v(x,y)=\int_0^y f_2(x,t)\, dt$; note that these functions only change on a null set if a different representative of $f_j$ is chosen and in particular the distributions induced by $u,v$ are well defined. To avoid trouble later on, fix representatives once and for all.
We have $T_x=u_x$ in $\mathcal D'$. This implies that
$$
(T,\varphi) = (u,\varphi) + \int_{-\infty}^{\infty} S(\varphi(x,\cdot))\, dx \quad\quad\quad\quad (1)
$$
for some $S\in\mathcal D'(\mathbb R)$. See Hormander I, Theorem 3.1.4'.
Now consider test functions of product form $\varphi(x,y)=g(x)h(y)$. If we fix one of the factors, say $h$, then the composed map $g \mapsto (U,gh)$ gives us a distribution on $\mathbb R$ for any $U\in\mathcal D'(\mathbb R^2)$. We have
$$
-(T_y,gh)=(T,gh')=(u,gh')+S(h')\int g =-(v_y, gh)=(v,gh') .
$$
The third and fifth expressions, thought of as one-dimensional distributions acting on $g$, are functions. Hence
$$
S(h')=\int (v(x,y)-u(x,y))h'(y)\, dy
$$
for almost every $x$. Here the exceptional null set may depend on $h'$. However, since a countable set of $h$'s suffices to determine $S$ on the codimension $1$ subspace of derivatives $h'$, it is in fact true that $v(x,y)-u(x,y)=w(y)+c(x)$ for almost every $x$. See here.
Hence $S'=-w$ and thus also $S$ itself are functions. Now we can go back to (1) to confirm that $T$ is a function also.