routine plpndx - j.f.chandler - 1982 feb
cross-reference partials requested for integration with list of
partials available on central-body or target-body tape
quantities of interest ('possible partials' + motion) are (in
the order of the output arrays):
1:   motion
2-7: partials w.r.t. initial conditions
8... any other partials

nctl = number of partials to be integrated
ictrl= list of descriptors for integrated partials
nlist= number of possible partials besides central/target-body i.c.'s
klist= list of descriptors to be compared with 'ictrl'
nplp=  planet number of central or target body
nq =   (on return) number of partials (+ motion) to be interpolated
kq =   array of 'nq' ptrs to list of possible partials (maps the
interpolation y-vectors to the output arrays)
kr =   array of 'nq' ptrs into partials on integration tape
kp =   array of 'nctl' ptrs from integrated partials to output arrays
numki= number of integration partials controls on integration tape
ki =   array of integration controls