Init of nonlinear fluid variables
Bug description / feature request:
Similar to full-f and delta-f, we need an initialization strategy for nonlinear fluid variables. At the moment, we really only initialize perturbations.
Expected behavior:
Add perturbations to a background that is related to the MHD equilibrium.
Proposed solution: