Tidying up Summer

[RyanOtago]

• Normalise initial conditions seperately

• Come up with a better name for 'grid_int' the volume per gridpoint for integral calculations

• Viscosity calculation. After determining ideal viscosity, change automatically with resolution of mesh. Need to check how it varies with other parameters

• Plot dissipated energy del_t E = avg(nukperp^(2n)(zeta^pm)^2). I.e the amount of energy dissipated at each time step. Should look kinda like a bell curve. 'Best' turbulence happens around peak of curve.

• Seems like there's a factor of 2pi missing in the time direction somewhere?

• Plot spectra/energy in terms of u and b (not zeta). Also useful to plot in terms of vorticity (emphasises small scale)

• Make sure there are no differences between the different RMHD_3D_* codes. Don't want previous mistakes to sneak back in. (Turbulence is most up to date!)

• Need to clean up spectrum and plotting code. Try and automate more of the process,

• In 'Initial Condition' section, add in 'SlowModes' condition instead of uncommenting slowmodes section

• Add switch in Spectrum code to choose between time averaged/instantaneous spectra