Add an option to include the Source field directly to the right hand side
for prognostic scalar fields. Adjust CG, DG, FV and CV to 
act on this as necessary. Include a very simple two field eigenvalue 
(using the time step loop as a power iteration loop) 1d test 
case for each discretisation checking order of converegence to test 
this new functionality.

This also includes a comment in the schema that variations in 
the Source field across subcycled time steps is not considered, 
a constant value across the whole time step is assumed (whether 
added directly or tested with shape_rhs).

In assemble/Advection_Diffusion_CG.F90 a redundant optional 
variable passed in called extra_discretised_source is removed, as 
well as the code that added it to the rhs.