diff --git a/moment_kinetics/src/load_data.jl b/moment_kinetics/src/load_data.jl index 9b5cdeb53..13f8ae0c3 100644 --- a/moment_kinetics/src/load_data.jl +++ b/moment_kinetics/src/load_data.jl @@ -629,8 +629,8 @@ function reload_evolving_fields!(pdf, moments, boundary_distributions, restart_p neutral_1V = (vzeta.n_global == 1 && vr.n_global == 1) restart_neutral_1V = (restart_vzeta.n_global == 1 && restart_vr.n_global == 1) - if geometry.bzeta != 0.0 && ((neutral1V && !restart_neutral_1V) || - (!neutral1V && restart_neutral_1V)) + if geometry.bzeta != 0.0 && ((neutral_1V && !restart_neutral_1V) || + (!neutral_1V && restart_neutral_1V)) # One but not the other of the run being restarted from and this run are # 1V, but the interpolation below does not allow for vz and vpa being in # different directions. Therefore interpolation between 1V and 3V cases