Test geocentric position call with scalar, array

Because this routine returns a 3-dimensional array, our vcall() routine
is hereby upgraded to handle the general case of n-dimensional result
sets; previously, it only worked with 2 dimensions.