Added allocation, initialization, and basic stats for flux values.  Still need to write subroutine to actually tally, maybe give tally type options, call the subroutine, and normalize either in this routine or in the stats.