Since $F_1,\dots,F_n$ are assumed to be functionally linearly independent at each point of $M_f$, this is also true near it, and on a open subset $A\subset M$ (containing $M_f$) one can find coordinate functions that complete the set given by $F_1,\dots,F_n$; let me call those functions $\varphi_1,\dots,\varphi_n$.
In my previous answer I mentioned that compact regular orbits are diffeomorphic to tori $\mathbb{T}^n$ and that
The foliation is a lagrangian fibration in a neighbourhood of each regular leaf; it defines a fibre bundle with lagrangian fibres.
The coordinate functions $F_1,\dots,F_n,\varphi_1,\dots,\varphi_n$ are exactly the ones used to show this claim (cf. problem of page 279 of Arnold's book).
As it is pointed out by Arnold (by the end of page 279), in these coordinates the symplectic form (over $A$) is not necessarily written as $\sum_{j=1}^{n}\mathrm{d}F_j\wedge\mathrm{d}\varphi_j$.
He proceeds by defining (page 282) functions $I_1,\dots,I_n$ and assumes that instead of $F_1,\dots,F_n,\varphi_1,\dots,\varphi_n$, one can use $I_1,\dots,I_n,\varphi_1,\dots,\varphi_n$ as coordinate functions (page 283, which reads: $\mathrm{det}(\partial\mathbf{I}/\partial\mathbf{f})|_{\mathbf{f}}\neq 0$).
Remark: Whilst Arnold proceeds by changing the coordinate functions $F_j$, in my first answer I change the coordinate functions $\varphi_j$; that is the main difference in our arguments, and this is why I can provide a semilocal proof directly.
Let me show here how he constructs these new coordinate functions. His constructions assume that the symplectic manifold is $\mathbb{R}^{2n}$ endowed with the canonical symplectic structure, here I am not going to assume this: as the Original Poster is interested in the global picture.
A Liouville torus is defined as the preimage of a point $\mathbf{f}\in\mathbb{R}^n$ by the momentum mapping $(F_1,\dots,F_n):M\to\mathbb{R}^n$,
$$M_\mathbf{f}:=\{x\in M \ ; \ (F_1(x),\dots,F_n(x))=\mathbf{f}\} \ , $$
and choosing $[\gamma_1],\dots,[\gamma_n]\in H_1^\infty(M_\mathbb{f};\mathbb{Z})$ generators of the first smooth singular homology group of the torus $M_\mathbb{f}$ (for details about this terminology, I refer the readers to Lee's Introduction to smooth manifolds, but this is what Arnold is doing with his cycles by the end of page 282) one can define functions $I_j:A\subset M\to\mathbb{R}$ by
$$I_j(x):=\int_{\gamma_j}\theta \ ,$$
where the $1$-form $\theta\in\Omega^1(A;\mathbb{R})$ is any differential form satisfying $\omega=\mathrm{d}\theta$ over the open subset $A\subset M$ (cf. page 198 of Arnold's book, and my comments on the Poincaré lemma with parameters; otherwise, one can simply assume that $M\cong\mathbb{R}^{2n}$, as Arnold does in page 282), and $\gamma_j$ are smooth singular cycles representing $[\gamma_j]$. For points $x\in M_\mathbf{f}\cap A$ this is well defined (problem of page 283 in Arnold's book); when $x\in A$ does not belong to a Liouville torus which is a preimage of $\mathbb{f}$, the fact that $A\cong\mathbb{T}^n\times\mathbb{R}^n$ implies that it belongs to another Liouville torus and the cycles are to be taken from that particular torus.
Now I can discuss the theorem of page 283 in Arnold's book whose proof is the reason why I am writing this. It states that the coordinate functions $I_1,\dots,I_n,\varphi_1,\dots,\varphi_n$ satisfy
$$\omega=\sum_{j=1}^{n}\mathrm{d}I_j\wedge\mathrm{d}\varphi_j \ , $$
over the whole open subset $A\subset M$.
I hope that the reader can appreciate the fact that $A\subset M$ is a neighbourhood of the torus $M_\mathbf{f}$ and that these coordinates are well defined on the whole of it (which is exactly what the Original Poster asked).
Now, for the proof Arnold observes that any $\theta$ satisfying $\omega=\mathrm{d}\theta$ over $A\subset M$ is closed when restricted to the submanifold $M_\mathbf{f}$, as the restriction of $\omega$ vanishes there (the torus is a lagrangian submanifold). This observation allows him to define a function $S:V\subset M_\mathbf{f}\to\mathbb{R}$ satisfying $\theta=\mathrm{d}S$ over a contractible neighbourhood $V\subset M_\mathbf{f}$ of a fixed point $x_0\in M_\mathbf{f}$ (cf. page 198 of Arnold's book). This is related to my use of the foliated version of the Poincaré lemma in my first answer.
By the definition of the coordinate functions $I_j$, and applying Stokes's Theorem, it "holds"
$$I_j=\int_{\gamma_j}\theta=\int_{\gamma_j}\mathrm{d}S=\int_{\partial\gamma_j}S=\Delta_jS \ .$$
The quotation marks on "holds" is due to $S$ not being well defined on the whole of $M_\mathbf{f}$ (a fact observed by Arnold at the beginning of page 284).
He, then, uses Darboux coordinate functions $p_1,\dots,p_n,q_1,\dots,q_n$ over a contractible open neighbourhood $W\subset M$ of the point $x_0\in M_{\mathbf{f}}$ (one can take it in such a way that its intersection with $M_{\mathbf{f}}$ is $V$), and applies the Hamilton-Jacobi method that he has developed in chapter 9 to conclude the proof of the theorem; with the generating function being
$$S(I_1,\dots,I_n,q_1,\dots,q_n)=\int_{\mathbf{q}_0}^{\mathbf{q}}\theta \ , $$
where the integral is computed from any curve in $M_{\mathbf{f}}$ joining $\mathbf{q}_0:=(q_1(x_0),\dots,q_n(x_0))\in M_{\mathbf{f}}$ and $\mathbf{q}:=(q_1,\dots,q_n)\in M_{\mathbf{f}}$.
The local coordinates $p_1,\dots,p_n,q_1,\dots,q_n$ are chosen by exploiting the fact that $M_{\mathbf{f}}\subset M$ is a lagrangian submanifold (a null manifold as Arnold calls it in page 274); thus, a small piece of the torus, $V=W\cap M_{\mathbf{f}}$, can be realised as a graph of a function $\mathbf{p}:\mathbb{R}^n\to\mathbb{R}^{2n}$. This graph has the property that the values of the coordinate functions $I_1,\dots,I_n$ are constant along it, and they vary as one moves from one small piece of Liouville torus $V\subset M_{\mathbf{f}}$ to a small piece of a neighbouring Liouville torus (not the same value as $\mathbf{f}$). This is why $I_1,\dots,I_n,q_1,\dots,q_n$ can be taken as independent coordinate functions.
It is important to remark that the use of the Hamilton-Jacobi method (at least the way it is developed in the book) implies that the proof is local, it does not automatically extends to a neighbourhood of a Liouville torus. However, if $\omega=\sum_{j=1}^{n}\mathrm{d}I_j\wedge\mathrm{d}\varphi_j$ holds near every point of $M_\mathbf{f}$, then one can conclude that this is true over a neighbourhood of $M_\mathbf{f}$, where the coordinate functions $I_1,\dots,I_n,\varphi_1,\dots,\varphi_n$ are well defined.
P.D: A user asked me to not delete a previous version of this answer, but I could not keep my promise to do so. The ones who have upvoted this answer are welcome to follow the proof of the linearisation theorem on my first answer to the original question.