What's Changed

Bugfix for scenarios with non-zero Phi_b_dot (rate of change of toroidal flux at LCFS).

An erroneous implementation of the Phi_b_dot terms led to non-physical evolution for scenarios with large Phi_b_dot , e.g. rapid ramp-down including reduction of plasma volume. The current diffusion equation and Ohmic power calculations were particularly affected.

Previous tests did not have large enough Phi_b_dot to notice the effect. The present correct implementation was verified against RAPTOR for significantly larger Phi_b_dot values than previously used.

It is highly recommended to rerun any previous simulations which included time-varying geometry with Phi_b variations, and investigate the impact.