s =   integration time variable (julian date+.5)
ktype =-1 determine rotation matrix for integrated body nplnt
ktype = 0 determine rotation matrix for central body ncentr
ktype = 1,2,3,4 determine rotation matrix for target body

if partials wrt target body rotation parameters should become
necessary, /rtpars/ can be expanded as follows:
common /rtpars/ pcnrot(4,3,3),ptrot(4,3,3,i_mxtrg),icrotp(4,2),
1 itrotp(4,2,i_mxtrg),icrot1,itrot1(i_mxtrg)