         coordinates. This matrix has the form

               If [1,1] = size(radius) then [3,3]   = size(jacobi)
               If [1,n] = size(radius) then [3,3,n] = size(jacobi)
                                             double = class(jacobi)

                   -                                -
                  |  dx/dr     dx/dlon     dx/dlat   |
                  |                                  |
                  |  dy/dr     dy/dlon     dy/dlat   |
                  |                       