code: def get_ptrm_dec_and_inc(self):
          
          PTRMS = self.PTRMS[1:]
          CART_pTRMS_orig = numpy.array([lib_direct.dir2cart(row[1:4]) for row in PTRMS])
          #B_lab_dir = [self.B_lab_dir[0], self.B_lab_dir[1], 1.] # dir
          tmin, tmax = self.t_Arai[0], self.t_Arai[-1]
          ptrms_dec_Free, ptrms_inc_Free, ptrm_best_fit_vector_Free, ptrm_tau_Free, ptrm_v_Free, ptrm_mass_center_Free, ptrm_PCA_sigma_Free = lib_direct.get_dec_and_inc(CART_pTRMS_orig, self.t_Arai, tmin, tmax, anchored=False)
          ptrms_angle = lib_direct.get_ptrms_angle(ptrm_best_fit_vector_Free, self.B_lab_cart)
          self.pars['ptrms_dec_Free'], self.pars['ptrms_inc_Free'] = ptrms_dec_Free, ptrms_inc_Free
          self.pars['ptrms_tau_Free'] = ptrm_tau_Free
          self.pars['ptrms_angle_Free'] = ptrms_angle
docs: not included in spd.