==> bsdata.txt <==
       R[m]            Z[m]          Bp[T/A]           L[m]            Itor
   0.0000000e+00   0.0000000e+00   1.0204300e-06   0.0000000e+00   0.0000000e+00
   0.0000000e+00   2.9126403e-03   1.0180757e-06   2.9126403e-03   2.3624299e-03
   0.0000000e+00   5.8252807e-03   1.0107585e-06   5.8252807e-03   4.7136515e-03
   0.0000000e+00   8.7379210e-03   9.9847625e-07   8.7379210e-03   7.0421592e-03
   0.0000000e+00   1.1650561e-02   9.8056801e-07   1.1650561e-02   9.3356788e-03
   0.0000000e+00   1.4563201e-02   9.5682318e-07   1.4563201e-02   1.1580927e-02
   0.0000000e+00   1.7475842e-02   9.2577734e-07   1.7475842e-02   1.3762678e-02
   0.0000000e+00   2.0388482e-02   8.8698080e-07   2.0388482e-02   1.5863488e-02
   0.0000000e+00   2.3301123e-02   8.3728300e-07   2.3301123e-02   1.7861743e-02
   0.0000000e+00   2.6213763e-02   7.7544860e-07   2.6213763e-02   1.9730742e-02
   0.0000000e+00   2.9126402e-02   6.9483850e-07   2.9126402e-02   2.1434661e-02
   0.0000000e+00   3.2039020e-02   5.9082552e-07   3.2039020e-02   2.2924609e-02
   0.0000000e+00   3.4951635e-02   4.4646328e-07   3.4951635e-02   2.4126715e-02
   0.0000000e+00   3.7646666e-02   2.6479311e-07   3.7646666e-02   2.4889409e-02
   0.0000000e+00   4.0341698e-02   0.0000000e+00   4.0341698e-02   2.5171539e-02
   2.4937150e-03   4.0464410e-02   1.3225176e-07   4.2838430e-02   2.5301240e-02
   4.9515474e-03   4.0829007e-02   2.7838206e-07   4.5323158e-02   2.5707210e-02
   7.3876497e-03   4.1439790e-02   4.1098274e-07   4.7834662e-02   2.6396089e-02
   9.7179240e-03   4.2274233e-02   5.4281794e-07   5.0309833e-02   2.7335430e-02
   1.1999044e-02   4.3354582e-02   6.4794610e-07   5.2833850e-02   2.8531284e-02

==> floops.txt <==
probe_name     L[m]     Leff[m]    Bp[T/A]
FS05 & FS06  0.151897  0.137159  1.265287e-06
FS04 & FS07  0.202696  0.050581  1.298034e-06
FS03 & FS08  0.253496  0.050613  1.199657e-06
FS02 & FS09  0.304295  0.050641  9.952265e-07
FS01 & FS10  0.355109  0.102048  7.114703e-07
FL01 & FL10  0.596652  0.105056  4.530843e-07
FL02 & FL09  0.647234  0.050317  6.342801e-07
FL03 & FL08  0.697793  0.050306  7.379110e-07
FL04 & FL07  0.748375  0.050291  7.517815e-07
FL05 & FL06  0.798935  0.076356  6.672926e-07

==> surfpr.txt <==
probe_name     L[m]     Leff[m]    Bp[T/A]
S04P & S05P  0.177155  0.158882  1.299190e-06
S03P & S06P  0.227954  0.050665  1.264065e-06
S02P & S07P  0.278895  0.050639  1.109001e-06
S01P & S08P  0.329568  0.111561  8.627680e-07
L01P & L10P  0.621943  0.113361  5.530606e-07
L02P & L09P  0.672503  0.050321  6.963821e-07
L03P & L08P  0.723084  0.050305  7.567737e-07
L04P & L07P  0.773644  0.096471  7.218625e-07

The URL www.psicenter.org/marklin/hitsi_probe_info contains the following:
bb.txt              % Taylor state surface bfield data output from COMSOL
bsdata.txt          % see above - surface bfield data normalized to Itor=1
effective_lengths.m % MATLAB code - computes above data from COMSOL output
floops.txt          % see above - flux loop data fitting info
probe_info.txt      % this file
surfpr.txt          % see above - surface probe data fitting info

Peicewise Taylor fit:
Itor = (1/mu0) * sum_over_probes (Bp_measured * Leff)

Least square Taylor fit:
Error(Itor) = sum_over_probes (Bp_measured - Itor * Bp_Taylor)**2

Set dError/dItor=0 to get:
Itor = sum_over_probes (Bp_measured * Bp_Taylor) / sum_over_probes (Bp_Taylor**2)