& n_dofs_process
INTEGER n_dofs_process (0:nPx*nPy-1)
local_dofs = n_dofs_process (mpiMyWid)
global_dofs = global_dofs + n_dofs_process (i)
& n_dofs_process(i-1)
local_dofs = n_dofs_process (0)
CALL MPI_ALLREDUCE(n_dofs_proc_loc,n_dofs_process,nPx*nPy,
n_dofs_process (0) = INT(DoFCount)+1