This subroutine behaves very like STRESS_MPI, except it calculates the current contribution from each pixel. A call to MPI_ALLREDUCE is needed to sum the individual contributions per processor onto the root (master) node. The 3 macroscopic current contributions (x,y,z) are available to each node, but root is the lone processor who uses these values for output purposes.