m.e.ash   june 1969    subroutine prdhar
pick out derivatives w.r.t. gravitational potential harmonic coef.
also for target body initial conditions and harmonics.
this subroutine differs from prdbdy only in lhar being doubly
dimensioned and in allowing gaps in the mhar vector compared to
the lhar vector for target bodies.