GeorgiaActon

> Brief experimentation with the input file for the implicitly treated MTM above suggests that the apparent bug with drifts_implicit =.true. is in fact caused by a bug with the...

### Notes on numerical instability: - Including parallel streaming in isolation seems stable, including mirror in isolation seems stable, but including them both without the presence of drifts can be...

Hi Yujia, I had a look at this and I think this problem is related to the initial condition. There are 2 options you have available to remove this problem:...

Here purple is your original input file, green is option 1. (using 'rh' for this initialisation) and blue is option 2 (removing the absolute value around the species charge in...

Sorry yes, in the new pull request of Stella (Stella-v1.0) in the file STELLA_CODE/arrays/initialise_distribution_function.f90, line 566

hmm - I have just tried re-running and I don't see any growth with all default unwinding values. Are there any differences that you can see in our input files...

If you run the initial condition 'default' with implicit_drifts = .true. does this also work (with fix 2. above)?

> I am happy with this, with a few minor exceptions. Here are my comments/questions: > > 1. How does include_electromagnetic work with the &electromagnetic namelist? E.g., if include_electromagnetic is...

Potential feature: - have maximal input file printed upon make with all default values in the code.

&extra_physics prp_shear_enabled = .false. hammett_flow_shear = .true. include_pressure_variation = .false. g_exb = 0.0 g_exbfac = 1.0 omprimfac = 1.0 changed to: & flow_shear g_exb = 0.0 --- EXB shearing rate...