there are a couple useless commits in there when I tried to merge what you had previously done with the postprocessing. landmark test cases are in landmark folder and the results of 3 simulations to 1 millisecond with each case are saved.
What I did manually to make method uncoupled:
change source term in source definition.
multiply pe by 0 in compute fluxes
commented out smooth_if in HLLE flux
added diffusion to both ion continuity and momentum, with 1mm diffusive length (should only be momentum but that breaks)
I should wrap this in some function and make it an option in the simulation struct.
there are a couple useless commits in there when I tried to merge what you had previously done with the postprocessing. landmark test cases are in landmark folder and the results of 3 simulations to 1 millisecond with each case are saved.
What I did manually to make method uncoupled: change source term in source definition. multiply pe by 0 in compute fluxes commented out smooth_if in HLLE flux added diffusion to both ion continuity and momentum, with 1mm diffusive length (should only be momentum but that breaks)
I should wrap this in some function and make it an option in the simulation struct.