Function handle for the benchmark density profile \(\bar{\rho}(z)\).
Discussion
For user-specified gridded profiles this may be empty. For the
built-in analytical and benchmark cases it stores the
continuous background density used to initialize the concrete
solver.