S3L_reduce_axis uses the following arguments for output:
b - S3L array handle for the parallel array that will contain the result of the reduction.
ier (Fortran only) - When called from a Fortran program, S3L_reduce_axis returns error status in ier.