Got some basics down for the jacobian construction.
The loops strucutre appears to give the correct index for within the jacobian.
However, since the memory is properly shared, only the root processor actually
saves the value to the root jacobian matrix.