 Solution

In Cartesian coordinates, $\nabla \Phi =({\Phi }_x,{\Phi }_y,{\Phi }_z)$. For orthogonal coordinates, let $\mathbf{r}(q_1,q_2,q_3)$ be the position, $h_i=|\partial \mathbf{r}/\partial q_i|$ the scale factors of orthogonal coordinates, and ${\mathbf{e}}_i=h_i^{-1}\partial \mathbf{r}/\partial q_i$ the corresponding unit vectors. Since $d\mathbf{r}={\sum }_i{h}_i{\mathbf{e}}_{i}d{q}_i$ and $d\Phi =\nabla \Phi \cdot d\mathbf{r}={\sum }_i{\Phi }_{q_i}d{q}_i$, comparison of independent increments gives the gradient in orthogonal coordinates:For the hyperbolic parametrization of oblate spheroidal coordinates, the three coordinate tangent vectors areTheir pairwise dot products are zero. Their squared Euclidean norms are ${\sinh }^{2}R+{\sin }^2\theta$, ${\sinh }^{2}R+{\sin }^2\theta$, and ${\cosh }^{2}R{\cos }^2\theta$, respectively. Thus the coordinate surfaces meet orthogonally wherever this is a regular coordinate chart, andThe last expression uses the usual latitude range $|\theta |<\pi /2$; without that range a scale factor is its absolute value. Axial or ring degeneracies are excluded from a regular chart.