use the LU-decomposed response matrix and the contributions from the 'inhomogeneous' fields (phi, apar) to solve for (phi^{n+1}, apar^{n+1})
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| complex, | intent(inout), | dimension(:, :, -nzgrid:, :) | :: | phi | ||
| complex, | intent(inout), | dimension(:, :, -nzgrid:, :) | :: | apar | ||
| complex, | intent(inout), | dimension(:, :, -nzgrid:, :) | :: | bpar |