[totork]

Right now in Fci, amax used for the dissipation in Div_par_mod / Div_par_fvv does not use the wave speed at the interface. This will lead to differences between cells that calculate the flux on the same interface, and hence spurious sources.

This change will require the parallel slices of the wave_speed_in to be set!

[dschwoerer]

Looks good 👍

Does this need changes in hermes-3 to work?

[totork]

Looks good 👍

Does this need changes in hermes-3 to work?

I think when sound_speed is on it communicates it, which means it should have parallel slices. Without it uses

Field3D fastest_wave; if (state.isSet("fastest_wave")) { fastest_wave = get<Field3D>(state["fastest_wave"]); } else { const Field3D T = get<Field3D>(species["temperature"]); const BoutReal AA = get<BoutReal>(species["AA"]); fastest_wave = sqrt(T / AA); }

Changing this from Field3D to Field3DParallel should do the trick