[FEAT] Embed the vertical ion drag in the viscostity calculation #164

The vertical velocity should be solved for in the viscosity solver, which is a semi-implicit solver. Because source terms can be included in the solver, they should be done here as opposed to spread out across the code.