Add a global nodal quadrature which should use J=1

Not sure if we need some other stuff to do this though, might not work
yet.
Rename argument -integration to -quadrature