// phase B — complete
A parameter nobody has measured
PFSS assumes the corona carries no currents, so the field is the
gradient of a potential satisfying Laplace’s equation. The model’s
real physical content is its upper boundary: a source surface
at radius Rss where the field is forced radial, standing in
for the solar wind dragging it open.
There is no such surface. Rss = 2.5 R☉ is a
convention from 1969, and recent work benchmarking against eclipse
images finds the best value moves with the solar cycle. So we do not
pick one — we run five, across all twelve realisations. Sixty solves,
124 seconds on a laptop.
Before trusting any of it: for a pure dipole boundary the potential
has a closed form, and the solver must reproduce it.
Br(r, θ) = b₀ cosθ · (Rss⁻³ + 2r⁻³) / (2 + Rss⁻³)
It does, to 5% at the source surface, with the polarity the right way
round, and it reproduces the input boundary map to 0.03% RMS. Sign
conventions, normalisation and the upper boundary condition are all
pinned by that single test.