Consider a linear system $\lvert D \rvert \subset H^0(\mathbf{P}^2 \times \mathbf{P}^2, \mathcal{O}(a,b))$ with nonempty base locus $B$. This linear system defines a rational map $\mathbf{P}^2 \times \mathbf{P}^2 \dashrightarrow \mathbf{P}^n$, defined on the complement of $B$. (Here $n$ is of course one less than the dimension of the linear system.) By resolving the base locus of $D$, I mean resolving the indeterminacy locus of this rational map, i.e. finding a regular map $X \to \mathbf P^n$ from some smooth variety $X$ birational to $\mathbf{P}^2 \times \mathbf{P}^2$, extending the rational map. That such a resolution exists is a standard fact; I am seeking to describe it explicitly.

For instance, if the base locus of $B$ consists of a disjoint union of irreducible components, it seems clear to me that $X$ can be obtained as the blowup along these components; in fact, I believe the same should apply if the irreducible components intersect transversally. However, things seem more complicated if the irreducible components of the base locus do not intersect transversally.

For a simple example, consider the linear system $$\langle x_1y_2, x_2y_1 \rangle,$$ where the coordinates on $\mathbf{P}^2 \times \mathbf{P}^2$ are $((x_1,x_2,x_3),(y_1,y_2,y_3))$. This gives a rational map $\mathbf{P}^2 \times \mathbf{P}^2 \dashrightarrow \mathbf{P}^1$.

The base locus of this linear system has irreducible components $X_1 = V(x_1,x_2)$, $Y_1 = V(y_1,y_2)$, $Z_1 = V(x_1,y_1)$, and $Z_2 = V(x_2,y_2)$. Evidently we have $X_1 \cap Y_1 = Z_1 \cap Z_2$ is a point, and these intersections are transversal; however, the intersections $X_1 \cap Z_i$ and $Y_1 \cap Z_i$ are not transversal: each such intersection has dimension $1$ rather than the expected dimension $0$.

In this case, how can we explicitly resolve the base locus? Simply blowing up the irreducible components seems like it would cause problems, because of these non-transversal intersections, so my next guess would be to first blowup the intersections of the irreducible components (and, in general, the intersections of the intersections, and so on, in increasing order of dimension). Intuitively, each blowup will make the proper transforms of the next subvarieties disjoint, hence reducing the problem to the easy case.

I've attempted this approach for similar linear systems on $\mathbf{P}^n$ and it seems to be valid. However, I have reason to suspect that this does not work quite as well for $\mathbf{P}^2 \times \mathbf{P}^2$. Without getting into too much detail, computing the intersection theory for this sequence of blowups seems to give some incorrect intersection numbers for pullbacks of divisors from $\mathbf{P}^1$.

Is there anything obviously wrong with my approach to resolving this base locus, and if so, is there a better technique I should be looking at for performing such a resolution? Perhaps more to the point, after performing a sequence of blowups, how can I actually check that the resulting linear system is basepoint-free or not? (For this last question, it seems possible to explicitly compute coordinates for the blowup, but as you can imagine such computations would get unwieldy quite quickly.)