Made solver dimension independant.
The fluxes can now be avaluated dimensionwise