Electron DoNothingBC's - Does the Energy Density follow?

[ghost]

I'm noticing that when I use the Electron DoNothingBC's the density profiles look as expected but I notice the electron temperature is quite high around the boundary at which the DoNothing is applied. I have a feeling that this is because mean_en is not leaving along with em, so it sees a no-flux boundary thus remaining constant, increasing the electron temperature.

Do I need to apply a separate energy BC at the boundaries where DoNothing is applied? My best guess is that I would need to write separate EnergyDoNothingBC to account for this, but the ElectronAdvectionDoNothingBC appears to only be a function of mean_en, so maybe this is the proper EnergyDoNothing?

Thanks!