Add new routine to read and summarize output

A fortran routine to read output is now available. This is way faster
than R and should be used if possible.
Change selectivity related variable names