Adds functionality to SPDP mode:
- all (field) summations done in double precision (end result still single precision)
- mesh face areas, face centres, cell volumes and cell centres calculated in double precision (but still stored in single precision)
- some mesh checks done in double precision