docs: def sbootpars(Taus, Vs):
      """
       get bootstrap parameters for s data
      """
  #
      Tau1s, Tau2s, Tau3s = [], [], []
      V1s, V2s, V3s = [], [], []
      nb = len(Taus)
      bpars = {}
      for k in range(nb):
          Tau1s.append(Taus[k][0])
          Tau2s.append(Taus[k][1])
          Tau3s.append(Taus[k][2])
          V1s.append(Vs[k][0])
          V2s.append(Vs[k][1])
          V3s.append(Vs[k][2])
      x, sig = gausspars(Tau1s)
      bpars["t1_sigma"] = sig
      x, sig = gausspars(Tau2s)
      bpars["t2_sigma"] = sig
      x, sig = gausspars(Tau3s)
      bpars["t3_sigma"] = sig
      V1s=flip(V1s,combine=True)
      kpars = dokent(V1s, len(V1s))
      bpars["v1_dec"] = kpars["dec"]
      bpars["v1_inc"] = kpars["inc"]
      bpars["v1_zeta"] = (kpars["Zeta"] * np.sqrt(nb)) % 360.
      bpars["v1_eta"] = (kpars["Eta"] * np.sqrt(nb)) % 360.
      bpars["v1_zeta_dec"] = kpars["Zdec"]
      bpars["v1_zeta_inc"] = kpars["Zinc"]
      bpars["v1_eta_dec"] = kpars["Edec"]
      bpars["v1_eta_inc"] = kpars["Einc"]
      V2s=flip(V2s,combine=True)
      kpars = dokent(V2s, len(V2s))
      bpars["v2_dec"] = kpars["dec"]
      bpars["v2_inc"] = kpars["inc"]
      bpars["v2_zeta"] = (kpars["Zeta"] * np.sqrt(nb)) % 360.
      bpars["v2_eta"] = (kpars["Eta"] * np.sqrt(nb)) % 360.
      bpars["v2_zeta_dec"] = kpars["Zdec"]
      bpars["v2_zeta_inc"] = kpars["Zinc"]
      bpars["v2_eta_dec"] = kpars["Edec"]
      bpars["v2_eta_inc"] = kpars["Einc"]
      V3s=flip(V3s,combine=True)
      kpars = dokent(V3s, len(V3s))
      bpars["v3_dec"] = kpars["dec"]
      bpars["v3_inc"] = kpars["inc"]
      bpars["v3_zeta"] = (kpars["Zeta"] * np.sqrt(nb)) % 360.
      bpars["v3_eta"] = (kpars["Eta"] * np.sqrt(nb)) % 360.
      bpars["v3_zeta_dec"] = kpars["Zdec"]
      bpars["v3_zeta_inc"] = kpars["Zinc"]
      bpars["v3_eta_dec"] = kpars["Edec"]
      bpars["v3_eta_inc"] = kpars["Einc"]
      return bpars