"Calculate Flux Linkage" wrong result for non-unity current

[gijswl]

Hi all,

While trying to implement an inductance calculation in 3D magnetodynamics (WhitneyAVSolver), I found that the flux linkage calculated by the "Calculate Flux Linkage" option of MagnetoDynamicsCalcFields is incorrect when the coil current is not 1 A. Following the calculation in CalcFields.F90, it appears that the flux linkage per component is calculated as

$$\lambda = \iiint_{component} \textbf{J} \cdot \textbf{A} dV$$

This appears to be missing a division by the coil current, which is supported by the observation that the units have an extra ampere ($[H][A]^2$ versus the expected $[H][A]$). If I divide the result stored by SaveScalars by the actual coil current, the resulting number seems to be the correct flux linkage.

I have attached a MWE here, calculating the (incorrect) flux linkage for an E-shaped powder core inductor versus current: flux_linkage_example.zip