use parallel_reduce to estimate min/max/average concentration

#51 (closed) The loop to estimate min/max/average concentration is updated.

update outputs for concentration status and energy