issue #41: doubled rate of summing moments via array(nx,ny,nz,10)